When the optimal is not the best: parameter estimation in complex biological models

Background: The vast computational resources that became available during the past decade enabled the development and simulation of increasingly complex mathematical models of cancer growth. These models typically involve many free parameters whose determination is a substantial obstacle to model development. Direct measurement of biochemical parameters in vivo is often difficult and sometimes impracticable, while fitting them under data-poor conditions may result in biologically implausible values. Results: We discuss different methodological approaches to estimate parameters in complex biological models. We make use of the high computational power of the Blue Gene technology to perform an extensive study of the parameter space in a model of avascular tumor growth. We explicitly show that the landscape of the cost function used to optimize the model to the data has a very rugged surface in parameter space. This cost function has many local minima with unrealistic solutions, including the global minimum corresponding to the best fit. Conclusions: The case studied in this paper shows one example in which model parameters that optimally fit the data are not necessarily the best ones from a biological point of view. To avoid force-fitting a model to a dataset, we propose that the best model parameters should be found by choosing, among suboptimal parameters, those that match criteria other than the ones used to fit the model. We also conclude that the model, data and optimization approach form a new complex system, and point to the need of a theory that addresses this problem more generally

The K2K SciBar Detector

A new near detector, SciBar, for the K2K long-baseline neutrino oscillation expe riment was installed to improve the measurement of neutrino energy spectrum and to study neutrino interactions in the energy region around 1 GeV. SciBar is a 'fully active' tracking detector with fine segmentation consisting of plastic scintillator bars. The detector was constructed in summer 2003 and is taking data since October 2003. The basic design and initial performance is presented.Comment: 7 pages, 4figures, Contributed to Proceedings of the 10th Vienna Conference on Instrumentation, Vienna, February 16-21, 200

Study of Bc+B_c^+ decays to the K+K−π+K^+K^-\pi^+ final state and evidence for the decay Bc+→χc0π+B_c^+\to\chi_{c0}\pi^+

A study of $B_c^+\to K^+K^-\pi^+$ decays is performed for the first time using data corresponding to an integrated luminosity of 3.0 $\mathrm{fb}^{-1}$ collected by the LHCb experiment in $pp$ collisions at centre-of-mass energies of $7$ and $8$ TeV. Evidence for the decay $B_c^+\to\chi_{c0}(\to K^+K^-)\pi^+$ is reported with a significance of 4.0 standard deviations, resulting in the measurement of $\frac{\sigma(B_c^+)}{\sigma(B^+)}\times\mathcal{B}(B_c^+\to\chi_{c0}\pi^+)$ to be $(9.8^{+3.4}_{-3.0}(\mathrm{stat})\pm 0.8(\mathrm{syst}))\times 10^{-6}$. Here $\mathcal{B}$ denotes a branching fraction while $\sigma(B_c^+)$ and $\sigma(B^+)$ are the production cross-sections for $B_c^+$ and $B^+$ mesons. An indication of $\bar b c$ weak annihilation is found for the region $m(K^-\pi^+)<1.834\mathrm{\,Ge\kern -0.1em V\!/}c^2$, with a significance of 2.4 standard deviations.Comment: All figures and tables, along with any supplementary material and additional information, are available at https://lhcbproject.web.cern.ch/lhcbproject/Publications/LHCbProjectPublic/LHCb-PAPER-2016-022.html, link to supplemental material inserted in the reference

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

Search for direct pair production of the top squark in all-hadronic final states in proton-proton collisions at s√=8 TeV with the ATLAS detector

The results of a search for direct pair production of the scalar partner to the top quark using an integrated luminosity of 20.1fb−1 of proton–proton collision data at √s = 8 TeV recorded with the ATLAS detector at the LHC are reported. The top squark is assumed to decay via t˜→tχ˜01 or t˜→ bχ˜±1 →bW(∗)χ˜01 , where χ˜01 (χ˜±1 ) denotes the lightest neutralino (chargino) in supersymmetric models. The search targets a fully-hadronic final state in events with four or more jets and large missing transverse momentum. No significant excess over the Standard Model background prediction is observed, and exclusion limits are reported in terms of the top squark and neutralino masses and as a function of the branching fraction of t˜ → tχ˜01 . For a branching fraction of 100%, top squark masses in the range 270–645 GeV are excluded for χ˜01 masses below 30 GeV. For a branching fraction of 50% to either t˜ → tχ˜01 or t˜ → bχ˜±1 , and assuming the χ˜±1 mass to be twice the χ˜01 mass, top squark masses in the range 250–550 GeV are excluded for χ˜01 masses below 60 GeV

Evidence for the h_b(1P) meson in the decay Upsilon(3S) --> pi0 h_b(1P)

Using a sample of 122 million Upsilon(3S) events recorded with the BaBar detector at the PEP-II asymmetric-energy e+e- collider at SLAC, we search for the $h_b(1P)$ spin-singlet partner of the P-wave chi_{bJ}(1P) states in the sequential decay Upsilon(3S) --> pi0 h_b(1P), h_b(1P) --> gamma eta_b(1S). We observe an excess of events above background in the distribution of the recoil mass against the pi0 at mass 9902 +/- 4(stat.) +/- 2(syst.) MeV/c^2. The width of the observed signal is consistent with experimental resolution, and its significance is 3.1sigma, including systematic uncertainties. We obtain the value (4.3 +/- 1.1(stat.) +/- 0.9(syst.)) x 10^{-4} for the product branching fraction BF(Upsilon(3S)-->pi0 h_b) x BF(h_b-->gamma eta_b).Comment: 8 pages, 4 postscript figures, submitted to Phys. Rev. D (Rapid Communications

Search for new phenomena in final states with an energetic jet and large missing transverse momentum in pp collisions at √ s = 8 TeV with the ATLAS detector

Results of a search for new phenomena in final states with an energetic jet and large missing transverse momentum are reported. The search uses 20.3 fb−1 of √ s = 8 TeV data collected in 2012 with the ATLAS detector at the LHC. Events are required to have at least one jet with pT > 120 GeV and no leptons. Nine signal regions are considered with increasing missing transverse momentum requirements between Emiss T > 150 GeV and Emiss T > 700 GeV. Good agreement is observed between the number of events in data and Standard Model expectations. The results are translated into exclusion limits on models with either large extra spatial dimensions, pair production of weakly interacting dark matter candidates, or production of very light gravitinos in a gauge-mediated supersymmetric model. In addition, limits on the production of an invisibly decaying Higgs-like boson leading to similar topologies in the final state are presente

Human oral viruses are personal, persistent and gender-consistent.

Viruses are the most abundant members of the human oral microbiome, yet relatively little is known about their biodiversity in humans. To improve our understanding of the DNA viruses that inhabit the human oral cavity, we examined saliva from a cohort of eight unrelated subjects over a 60-day period. Each subject was examined at 11 time points to characterize longitudinal differences in human oral viruses. Our primary goals were to determine whether oral viruses were specific to individuals and whether viral genotypes persisted over time. We found a subset of homologous viral genotypes across all subjects and time points studied, suggesting that certain genotypes may be ubiquitous among healthy human subjects. We also found significant associations between viral genotypes and individual subjects, indicating that viruses are a highly personalized feature of the healthy human oral microbiome. Many of these oral viruses were not transient members of the oral ecosystem, as demonstrated by the persistence of certain viruses throughout the entire 60-day study period. As has previously been demonstrated for bacteria and fungi, membership in the oral viral community was significantly associated with the sex of each subject. Similar characteristics of personalized, sex-specific microflora could not be identified for oral bacterial communities based on 16S rRNA. Our findings that many viruses are stable and individual-specific members of the oral ecosystem suggest that viruses have an important role in the human oral ecosystem

Indole is an essential herbivore-induced volatile priming signal in maize

Herbivore-induced volatile organic compounds prime non-attacked plant tissues to respond more strongly to subsequent attacks. However, the key volatiles that trigger this primed state remain largely unidentified. In maize, the release of the aromatic compound indole is herbivore-specific and occurs earlier than other induced responses. We therefore hypothesized that indole may be involved in airborne priming. Using indole-deficient mutants and synthetic indole dispensers, we show that herbivore-induced indole enhances the induction of defensive volatiles in neighbouring maize plants in a species-specific manner. Furthermore, the release of indole is essential for priming of mono- and homoterpenes in systemic leaves of attacked plants. Indole exposure markedly increases the herbivore-induced production of the stress hormones jasmonate-isoleucine conjugate and abscisic acid, which represents a likely mechanism for indole-dependent priming. These results demonstrate that indole functions as a rapid and potent aerial priming agent that prepares systemic tissues and neighbouring plants for incoming attacks

Observation of Bc+ →j /ψD (∗)K (∗) decays

A search for the decays B+c→J/ψD(*)0K+ and B+c→J/ψD(*)+K*0 is performed with data collected at the LHCb experiment corresponding to an integrated luminosity of 3 fb−1. The decays B+c→J/ψ0K+ and B+c→J/ψD*0K+ are observed for the first time, while first evidence is reported for the B+c→JψD*+K*0 and B+c→J/ψD+K*0 decays. The branching fractions of these decays are determined relative to the B+c→J/ψπ+ decay. The B+c mass is measured, using the J/ψD0K+ final state, to be 6274.28±1.40(stat)±0.32(syst) MeV/c2. This is the most precise single measurement of the B+c mass to date