Automated multi-objective calibration of biological agent-based simulations

Computational agent-based simulation (ABS) is increasingly used to complement laboratory techniques in advancing our understanding of biological systems. Calibration, the identification of parameter values that align simulation with biological behaviours, becomes challenging as increasingly complex biological domains are simulated. Complex domains cannot be characterized by single metrics alone, rendering simulation calibration a fundamentally multi-metric optimization problem that typical calibration techniques cannot handle. Yet calibration is an essential activity in simulation-based science; the baseline calibration forms a control for subsequent experimentation and hence is fundamental in the interpretation of results. Here, we develop and showcase a method, built around multi-objective optimization, for calibrating ABSs against complex target behaviours requiring several metrics (termed objectives) to characterize. Multi-objective calibration (MOC) delivers those sets of parameter values representing optimal trade-offs in simulation performance against each metric, in the form of a Pareto front. We use MOC to calibrate a well-understood immunological simulation against both established a priori and previously unestablished target behaviours. Furthermore, we show that simulation-borne conclusions are broadly, but not entirely, robust to adopting baseline parameter values from different extremes of the Pareto front, highlighting the importance of MOC's identification of numerous calibration solutions. We devise a method for detecting overfitting in a multi-objective context, not previously possible, used to save computational effort by terminating MOC when no improved solutions will be found. MOC can significantly impact biological simulation, adding rigour to and speeding up an otherwise time-consuming calibration process and highlighting inappropriate biological capture by simulations that cannot be well calibrated. As such, it produces more accurate simulations that generate more informative biological predictions

Monocytes regulate the mechanism of T-cell death by inducing Fas-mediated apoptosis during bacterial infection.

Monocytes and T-cells are critical to the host response to acute bacterial infection but monocytes are primarily viewed as amplifying the inflammatory signal. The mechanisms of cell death regulating T-cell numbers at sites of infection are incompletely characterized. T-cell death in cultures of peripheral blood mononuclear cells (PBMC) showed 'classic' features of apoptosis following exposure to pneumococci. Conversely, purified CD3(+) T-cells cultured with pneumococci demonstrated necrosis with membrane permeabilization. The death of purified CD3(+) T-cells was not inhibited by necrostatin, but required the bacterial toxin pneumolysin. Apoptosis of CD3(+) T-cells in PBMC cultures required 'classical' CD14(+) monocytes, which enhanced T-cell activation. CD3(+) T-cell death was enhanced in HIV-seropositive individuals. Monocyte-mediated CD3(+) T-cell apoptotic death was Fas-dependent both in vitro and in vivo. In the early stages of the T-cell dependent host response to pneumococci reduced Fas ligand mediated T-cell apoptosis was associated with decreased bacterial clearance in the lung and increased bacteremia. In summary monocytes converted pathogen-associated necrosis into Fas-dependent apoptosis and regulated levels of activated T-cells at sites of acute bacterial infection. These changes were associated with enhanced bacterial clearance in the lung and reduced levels of invasive pneumococcal disease

Protecting eyewitness evidence: Examining the efficacy of a self-administered interview tool

Given the crucial role of eyewitness evidence, statements should be obtained as soon as possible after an incident. This is not always achieved due to demands on police resources. Two studies trace the development of a new tool, the Self-Administered Interview (SAI), designed to elicit a comprehensive initial statement. In Study 1, SAI participants reported more correct details than participants who provided a free recall account, and performed at the same level as participants given a Cognitive Interview. In Study 2, participants viewed a simulated crime and half recorded their statement using the SAI. After a delay of 1 week, all participants completed a free recall test. SAI participants recalled more correct details in the delayed recall task than control participants

Chedoke Arm and Hand Activity Inventory-9 (CAHAI-9): Perceived clinical utility within 14 days of stroke

Purpose: The Chedoke Arm and Hand Activity Inventory-9 (CAHAI-9) is an activity-based assessment developed to include relevant functional tasks and to be sensitive to clinically important changes in upper limb function. The aim of this study was to explore both therapists' and clients' views on the clinical utility of CAHAI-9 within 14 days of stroke. Method: Twenty-one occupational therapists actively working in stroke settings were recruited by convenience sampling from 8 hospitals and participated in semistructured focus groups. Five clients within 14 days of stroke were recruited by consecutive sampling from 1 metropolitan hospital and participated in structured individual interviews. The transcripts were analyzed thematically. Results: Six themes emerged from the focus groups and interviews: collecting information, decisions regarding client suitability, administration and scoring, organizational demands, raising awareness, and clients' perceptions of CAHAI-9 utility. All therapists agreed CAHAI-9 was suited for the stroke population and assisted identification of client abilities or difficulties within functional contexts. Opinions varied as to whether CAHAI-9 should be routinely administered with clients who had mild and severe upper limb deficits, but therapists agreed it was appropriate for clients with moderate deficits. Therapists made suggestions regarding refinement of the scoring and training to increase utility. All clients with stroke felt that the assessment provided reassurance regarding their recovery. Conclusion: The findings indicate that CAHAI-9 shows promise as an upper limb ability assessment for clients within 14 days of stroke

Deuteron and antideuteron production in Au+Au collisions at sqrt(s_NN)=200 GeV

The production of deuterons and antideuterons in the transverse momentum range 1.1 < p_T < 4.3 GeV/c at mid-rapidity in Au + Au collisions at sqrt(s_NN)=200 GeV has been studied by the PHENIX experiment at RHIC. A coalescence analysis comparing the deuteron and antideuteron spectra with those of protons and antiprotons, has been performed. The coalescence probability is equal for both deuterons and antideuterons and increases as a function of p_T, which is consistent with an expanding collision zone. Comparing (anti)proton yields p_bar/p = 0.73 +/- 0.01, with (anti)deuteron yields: d_bar/d = 0.47 +/- 0.03, we estimate that n_bar/n = 0.64 +/- 0.04.Comment: 326 authors, 6 pages text, 5 figures, 1 Table. Submitted to PRL. Plain text data tables for the points plotted in figures for this and previous PHENIX publications are (or will be) publicly available at http://www.phenix.bnl.gov/papers.htm

An XMM-Newton view of M101 - III. Diffuse X-ray emission

We present a study of the X-ray properties of the nearby face-on Scd spiral galaxy M101 based on recent XMM–Newton observations. In this third and final paper in the present series, we focus on the spatial and spectral properties of the residual emission, after excluding bright X-ray sources with LX > 1037 erg s−1. Within a central region of radius 10 arcmin (21 kpc), the X-ray emission broadly traces the pattern of the spiral arms, establishing a strong link with recent star formation, but it also exhibits a radial scalelength of ≈2.6 arcmin (5.4 kpc) consistent with optical data. We estimate the soft X-ray luminosity within the central 5 arcmin (10.5 kpc) region to be LX ≈ 2.1 × 1039 erg s−1 (0.5–2 keV), the bulk of which appears to originate as diffuse emission. We find a two-temperature thermal model best fits the spectral data with derived temperatures of keV which are very typical of the diffuse components seen in other normal and starburst galaxies. More detailed investigation of the X-ray morphology reveals a strong correlation with images recorded in the far-ultraviolet through to V band, with the best match being with the U band. We interpret these results in terms of a clumpy thin-disc component which traces the spiral arms of M101 plus an extended lower halo component with large filling factor

Effective Rheology of Bubbles Moving in a Capillary Tube

We calculate the average volumetric flux versus pressure drop of bubbles moving in a single capillary tube with varying diameter, finding a square-root relation from mapping the flow equations onto that of a driven overdamped pendulum. The calculation is based on a derivation of the equation of motion of a bubble train from considering the capillary forces and the entropy production associated with the viscous flow. We also calculate the configurational probability of the positions of the bubbles.Comment: 4 pages, 1 figur

Single Electrons from Heavy Flavor Decays in p+p Collisions at sqrt(s) = 200 GeV

The invariant differential cross section for inclusive electron production in p+p collisions at sqrt(s) = 200 GeV has been measured by the PHENIX experiment at the Relativistic Heavy Ion Collider over the transverse momentum range $0.4 <= p_T <= 5.0 GeV/c at midrapidity (eta <= 0.35). The contribution to the inclusive electron spectrum from semileptonic decays of hadrons carrying heavy flavor, i.e. charm quarks or, at high p_T, bottom quarks, is determined via three independent methods. The resulting electron spectrum from heavy flavor decays is compared to recent leading and next-to-leading order perturbative QCD calculations. The total cross section of charm quark-antiquark pair production is determined as sigma_(c c^bar) = 0.92 +/- 0.15 (stat.) +- 0.54 (sys.) mb.Comment: 329 authors, 6 pages text, 3 figures. Submitted to Phys. Rev. Lett. Plain text data tables for the points plotted in figures for this and previous PHENIX publications are (or will be) publicly available at http://www.phenix.bnl.gov/papers.htm

Measurement of the inclusive and dijet cross-sections of b-jets in pp collisions at sqrt(s) = 7 TeV with the ATLAS detector

The inclusive and dijet production cross-sections have been measured for jets containing b-hadrons (b-jets) in proton-proton collisions at a centre-of-mass energy of sqrt(s) = 7 TeV, using the ATLAS detector at the LHC. The measurements use data corresponding to an integrated luminosity of 34 pb^-1. The b-jets are identified using either a lifetime-based method, where secondary decay vertices of b-hadrons in jets are reconstructed using information from the tracking detectors, or a muon-based method where the presence of a muon is used to identify semileptonic decays of b-hadrons inside jets. The inclusive b-jet cross-section is measured as a function of transverse momentum in the range 20 < pT < 400 GeV and rapidity in the range |y| < 2.1. The bbbar-dijet cross-section is measured as a function of the dijet invariant mass in the range 110 < m_jj < 760 GeV, the azimuthal angle difference between the two jets and the angular variable chi in two dijet mass regions. The results are compared with next-to-leading-order QCD predictions. Good agreement is observed between the measured cross-sections and the predictions obtained using POWHEG + Pythia. MC@NLO + Herwig shows good agreement with the measured bbbar-dijet cross-section. However, it does not reproduce the measured inclusive cross-section well, particularly for central b-jets with large transverse momenta.Comment: 10 pages plus author list (21 pages total), 8 figures, 1 table, final version published in European Physical Journal