Bayesian inference of agent-based models: a tool for studying kidney branching morphogenesis

The adult mammalian kidney has a complex, highly-branched collecting duct epithelium that arises as a ureteric bud sidebranch from an epithelial tube known as the nephric duct. Subsequent branching of the ureteric bud to form the collecting duct tree is regulated by subcellular interactions between the epithelium and a population of mesenchymal cells that surround the tips of outgrowing branches. The mesenchymal cells produce glial cell-line derived neurotrophic factor (GDNF), that binds with RET receptors on the surface of the epithelial cells to stimulate several subcellular pathways in the epithelium. Such interactions are known to be a prerequisite for normal branching development, although competing theories exist for their role in morphogenesis. Here we introduce the first agent-based model of ex vivo kidney uretic branching. Through comparison with experimental data, we show that growth factor-regulated growth mechanisms can explain early epithelial cell branching, but only if epithelial cell division depends in a switch-like way on the local growth factor concentration; cell division occurring only if the driving growth factor level exceeds a threshold. We also show how a recently-developed method, "Approximate Approximate Bayesian Computation", can be used to infer key model parameters, and reveal the dependency between the parameters controlling a growth factor-dependent growth switch. These results are consistent with a requirement for signals controlling proliferation and chemotaxis, both of which are previously identified roles for GDNF

Topology and Computational Performance of Attractor Neural Networks

To explore the relation between network structure and function, we studied the computational performance of Hopfield-type attractor neural nets with regular lattice, random, small-world and scale-free topologies. The random net is the most efficient for storage and retrieval of patterns by the entire network. However, in the scale-free case retrieval errors are not distributed uniformly: the portion of a pattern encoded by the subset of highly connected nodes is more robust and efficiently recognized than the rest of the pattern. The scale-free network thus achieves a very strong partial recognition. Implications for brain function and social dynamics are suggestive.Comment: 2 figures included. Submitted to Phys. Rev. Letter

Exclusive Queueing Process with Discrete Time

In a recent study [C Arita, Phys. Rev. E 80, 051119 (2009)], an extension of the M/M/1 queueing process with the excluded-volume effect as in the totally asymmetric simple exclusion process (TASEP) was introduced. In this paper, we consider its discrete-time version. The update scheme we take is the parallel one. A stationary-state solution is obtained in a slightly arranged matrix product form of the discrete-time open TASEP with the parallel update. We find the phase diagram for the existence of the stationary state. The critical line which separates the parameter space into the regions with and without the stationary state can be written in terms of the stationary current of the open TASEP. We calculate the average length of the system and the average number of particles

A dual point description of mesoscopic superconductors

We present an analysis of the magnetic response of a mesoscopic superconductor, i.e. a system of sizes comparable to the coherence length and to the London penetration depth. Our approach is based on special properties of the two dimensional Ginzburg-Landau equations, satisfied at the dual point $(\kappa = \frac{1}{\sqrt{2}}).$ Closed expressions for the free energy and the magnetization of the superconductor are derived. A perturbative analysis in the vicinity of the dual point allows us to take into account vortex interactions, using a new scaling result for the free energy. In order to characterize the vortex/current interactions, we study vortex configurations that are out of thermodynamical equilibrium. Our predictions agree with the results of recent experiments performed on mesoscopic aluminium disks.Comment: revtex, 20 pages, 9 figure

Structure comparison of binary and weighted niche-overlap graphs

In ecological networks, niche-overlap graphs are considered as complex systems. They represent the competition between two predators that share common resources. The purpose of this paper is to investigate the structural properties of these graphs considered as weighted networks and compare their measures with the ones calculated for the binary networks. To conduct this study, we select four classical network measures : the degree of nodes, the clustering coefficient, the assortativity, and the betweenness centrality. These measures were used to analyse different type of networks such as social networks, biological networks, world wide web, etc. Interestingly, we identify significant differences between the structure of the binary and the weighted niche-overlap graphs. This study indicates that weight information reveals different features that may provide other implications on the dynamics of these networks

Electronic Collective Modes and Superconductivity in Layered Conductors

A distinctive feature of layered conductors is the presence of low-energy electronic collective modes of the conduction electrons. This affects the dynamic screening properties of the Coulomb interaction in a layered material. We study the consequences of the existence of these collective modes for superconductivity. General equations for the superconducting order parameter are derived within the strong-coupling phonon-plasmon scheme that account for the screened Coulomb interaction. Specifically, we calculate the superconducting critical temperature Tc taking into account the full temperature, frequency and wave-vector dependence of the dielectric function. We show that low-energy plasmons may contribute constructively to superconductivity. Three classes of layered superconductors are discussed within our model: metal-intercalated halide nitrides, layered organic materials and high-Tc oxides. In particular, we demonstrate that the plasmon contribution (electronic mechanism) is dominant in the first class of layered materials. The theory shows that the description of so-called ``quasi-two-dimensional superconductors'' cannot be reduced to a purely 2D model, as commonly assumed. While the transport properties are strongly anisotropic, it remains essential to take into account the screened interlayer Coulomb interaction to describe the superconducting state of layered materials.Comment: Final version (minor changes) 14 pages, 6 figure

Tomato: a crop species amenable to improvement by cellular and molecular methods

Tomato is a crop plant with a relatively small DNA content per haploid genome and a well developed genetics. Plant regeneration from explants and protoplasts is feasable which led to the development of efficient transformation procedures. In view of the current data, the isolation of useful mutants at the cellular level probably will be of limited value in the genetic improvement of tomato. Protoplast fusion may lead to novel combinations of organelle and nuclear DNA (cybrids), whereas this technique also provides a means of introducing genetic information from alien species into tomato. Important developments have come from molecular approaches. Following the construction of an RFLP map, these RFLP markers can be used in tomato to tag quantitative traits bred in from related species. Both RFLP's and transposons are in the process of being used to clone desired genes for which no gene products are known. Cloned genes can be introduced and potentially improve specific properties of tomato especially those controlled by single genes. Recent results suggest that, in principle, phenotypic mutants can be created for cloned and characterized genes and will prove their value in further improving the cultivated tomato.

Modelling thirty-day mortality in the acute respiratory distress syndrome (ARDS) in an adult ICU

Publisher's copy made available with the permission of the publisher © Australian Society of AnaesthetistsVariables predicting thirty-day outcome from Acute Respiratory Distress Syndrome (ARDS) were analysed using Cox regression structured for time-varying covariates. Over a three-year period, 1996-1998, consecutive patients with ARDS (bilateral chest X-ray opacities, PaO₂/FiO₂ ratio of <200 and an acute precipitating event) were identified using a prospective computerized data base in a university teaching hospital ICU. The cohort, 106 mechanically ventilated patients, was of mean (SD) age 63.5 (15.5) years and 37% were female. Primary lung injury occurred in 45% and 24% were postoperative. ICU-admission day APACHE II score was 25 (8); ARDS onset time from ICU admission was 1 day (median: range 0-16) and 30 day mortality was 41% (95% CI: 33%-51%). At ARDS onset, PaO₂/FiO₂ ratio was 92 (31), 81% had four-quadrant chest X-ray opacification and lung injury score was 2.75 (0.45). Average mechanical ventilator tidal volume was 10.3 ml/ predicted kg weight. Cox model mortality predictors (hazard ratio, 95% CI) were: APACHE II score, 1.15 (1.09-1.21); ARDS lag time (days), 0.72 (0.58-0.89); direct versus indirect injury, 2.89 (1.45-5.76); PaO₂/FiO₂ ratio, 0.98 (0.97-0.99); operative versus non-operative category, 0.24 (0.09-0.63). Time-varying effects were evident for PaO₂/FiO₂ ratio, operative versus non-operative category and ventilator tidal volume assessed as a categorical predictor with a cut-point of 8 ml/kg predicted weight (mean tidal volumes, 7.1 (1.9) vs 10.7 (1.6) ml/kg predicted weight). Thirty-day survival was improved for patients ventilated with lower tidal volumes. Survival predictors in ARDS were multifactorial and related to patient-injury-time interaction and level of mechanical ventilator tidal volume.J. L. Moran, P. J. Solomon, V. Fox, M. Salagaras, P. J. Williams, K. Quinlan, A. D. Berstenhttp://www.aaic.net.au/Article.asp?D=200332