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

Can Genetic Programming Do Manifold Learning Too?

Exploratory data analysis is a fundamental aspect of knowledge discovery that aims to find the main characteristics of a dataset. Dimensionality reduction, such as manifold learning, is often used to reduce the number of features in a dataset to a manageable level for human interpretation. Despite this, most manifold learning techniques do not explain anything about the original features nor the true characteristics of a dataset. In this paper, we propose a genetic programming approach to manifold learning called GP-MaL which evolves functional mappings from a high-dimensional space to a lower dimensional space through the use of interpretable trees. We show that GP-MaL is competitive with existing manifold learning algorithms, while producing models that can be interpreted and re-used on unseen data. A number of promising future directions of research are found in the process.Comment: 16 pages, accepted in EuroGP '1

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

Epidemiology of astrovirus infection in children

Human astrovirus (HAstV) is a major cause of acute diarrhea among children, resulting in outbreaks of diarrhea and occasionally hospitalization. Improved surveillance and application of sensitive molecular diagnostics have further defined the impact of HAstV infections in children. These studies have shown that HAstV infections are clinically milder (diarrhea, vomiting, fever) than infections with other enteric agents. Among the 8 serotypes of HAstV identified, serotype 1 is the predominant strain worldwide. In addition to serotype 1, the detection rate of HAstV types 2 to 8 has increased by using newly developed assays. HAstV is less common compared with other major gastroenteritis viruses, including norovirus and rotavirus; however, it is a potentially important viral etiological agent with a significant role in acute gastroenteritis. A better understanding of the molecular epidemiology and characteristics of HAstV strains may be valuable to develop specific prevention strategies

Chronic Melatonin Administration Reduced Oxidative Damage and Cellular Senescence in the Hippocampus of a Mouse Model of Down Syndrome

Previous studies have demonstrated that melatonin administration improves spatial learning and memory and hippocampal long-term potentiation in the adult Ts65Dn (TS) mouse, a model of Down syndrome (DS). This functional benefit of melatonin was accompanied by protection from cholinergic neurodegeneration and the attenuation of several hippocampal neuromorphological alterations in TS mice. Because oxidative stress contributes to the progression of cognitive deficits and neurodegeneration in DS, this study evaluates the antioxidant effects of melatonin in the brains of TS mice. Melatonin was administered to TS and control mice from 6 to 12 months of age and its effects on the oxidative state and levels of cellular senescence were evaluated. Melatonin treatment induced antioxidant and antiaging effects in the hippocampus of adult TS mice. Although melatonin administration did not regulate the activities of the main antioxidant enzymes (superoxide dismutase, catalase, glutathione peroxidase, glutathione reductase, and glutathione S-transferase) in the cortex or hippocampus, melatonin decreased protein and lipid oxidative damage by reducing the thiobarbituric acid reactive substances (TBARS) and protein carbonyls (PC) levels in the TS hippocampus due to its ability to act as a free radical scavenger. Consistent with this reduction in oxidative stress, melatonin also decreased hippocampal senescence in TS animals by normalizing the density of senescence-associated â-galactosidase positive cells in the hippocampus. These results showed that this treatment attenuated the oxidative damage and cellular senescence in the brain of TS mice and support the use of melatonin as a potential therapeutic agent for age-related cognitive deficits and neurodegeneration in adults with DS

G-Quadruplex Visualization in Cells via Antibody and Fluorescence Probe

G-quadruplexes (G4s) are noncanonical nucleic acids structures involved in key regulatory and pathological roles in eukaryotes, prokaryotes, and viruses: the development of specific antibodies and fluorescent probes represent an invaluable tool to understand their biological relevance. We here present three protocols for the visualization of G4s in cells, both uninfected and HSV-1 infected, using a specific antibody and a fluorescent G4 ligand, and the effect of the fluorescent ligand on a G4 binding protein, nucleolin, upon binding of the molecule to the nucleic acids structure

High-resolution genomics identifies pneumococcal diversity and persistence of vaccine types in children with community-acquired pneumonia in the UK and Ireland.

BACKGROUND: Streptococcus pneumoniae is a global cause of community-acquired pneumonia (CAP) and invasive disease in children. The CAP-IT trial (grant No. 13/88/11; https://www.capitstudy.org.uk/ ) collected nasopharyngeal swabs from children discharged from hospitals with clinically diagnosed CAP, and found no differences in pneumococci susceptibility between higher and lower antibiotic doses and shorter and longer durations of oral amoxicillin treatment. Here, we studied in-depth the genomic epidemiology of pneumococcal (vaccine) serotypes and their antibiotic resistance profiles. METHODS: Three-hundred and ninety pneumococci cultured from 1132 nasopharyngeal swabs from 718 children were whole-genome sequenced (Illumina) and tested for susceptibility to penicillin and amoxicillin. Genome heterogeneity analysis was performed using long-read sequenced isolates (PacBio, n = 10) and publicly available sequences. RESULTS: Among 390 unique pneumococcal isolates, serotypes 15B/C, 11 A, 15 A and 23B1 were most prevalent (n = 145, 37.2%). PCV13 serotypes 3, 19A, and 19F were also identified (n = 25, 6.4%). STs associated with 19A and 19F demonstrated high genome variability, in contrast to serotype 3 (n = 13, 3.3%) that remained highly stable over a 20-year period. Non-susceptibility to penicillin (n = 61, 15.6%) and amoxicillin (n = 10, 2.6%) was low among the pneumococci analysed here and was independent of treatment dosage and duration. However, all 23B1 isolates (n = 27, 6.9%) were penicillin non-susceptible. This serotype was also identified in ST177, which is historically associated with the PCV13 serotype 19F and penicillin susceptibility, indicating a potential capsule-switch event. CONCLUSIONS: Our data suggest that amoxicillin use does not drive pneumococcal serotype prevalence among children in the UK, and prompts consideration of PCVs with additional serotype coverage that are likely to further decrease CAP in this target population. Genotype 23B1 represents the convergence of a non-vaccine genotype with penicillin non-susceptibility and might provide a persistence strategy for ST types historically associated with vaccine serotypes. This highlights the need for continued genomic surveillance

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

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