Evaluation of rate law approximations in bottom-up kinetic models of metabolism.

BackgroundThe mechanistic description of enzyme kinetics in a dynamic model of metabolism requires specifying the numerical values of a large number of kinetic parameters. The parameterization challenge is often addressed through the use of simplifying approximations to form reaction rate laws with reduced numbers of parameters. Whether such simplified models can reproduce dynamic characteristics of the full system is an important question.ResultsIn this work, we compared the local transient response properties of dynamic models constructed using rate laws with varying levels of approximation. These approximate rate laws were: 1) a Michaelis-Menten rate law with measured enzyme parameters, 2) a Michaelis-Menten rate law with approximated parameters, using the convenience kinetics convention, 3) a thermodynamic rate law resulting from a metabolite saturation assumption, and 4) a pure chemical reaction mass action rate law that removes the role of the enzyme from the reaction kinetics. We utilized in vivo data for the human red blood cell to compare the effect of rate law choices against the backdrop of physiological flux and concentration differences. We found that the Michaelis-Menten rate law with measured enzyme parameters yields an excellent approximation of the full system dynamics, while other assumptions cause greater discrepancies in system dynamic behavior. However, iteratively replacing mechanistic rate laws with approximations resulted in a model that retains a high correlation with the true model behavior. Investigating this consistency, we determined that the order of magnitude differences among fluxes and concentrations in the network were greatly influential on the network dynamics. We further identified reaction features such as thermodynamic reversibility, high substrate concentration, and lack of allosteric regulation, which make certain reactions more suitable for rate law approximations.ConclusionsOverall, our work generally supports the use of approximate rate laws when building large scale kinetic models, due to the key role that physiologically meaningful flux and concentration ranges play in determining network dynamics. However, we also showed that detailed mechanistic models show a clear benefit in prediction accuracy when data is available. The work here should help to provide guidance to future kinetic modeling efforts on the choice of rate law and parameterization approaches

Lithic technological responses to Late Pleistocene glacial cycling at Pinnacle Point Site 5-6, South Africa

There are multiple hypotheses for human responses to glacial cycling in the Late Pleistocene, including changes in population size, interconnectedness, and mobility. Lithic technological analysis informs us of human responses to environmental change because lithic assemblage characteristics are a reflection of raw material transport, reduction, and discard behaviors that depend on hunter-gatherer social and economic decisions. Pinnacle Point Site 5-6 (PP5-6), Western Cape, South Africa is an ideal locality for examining the influence of glacial cycling on early modern human behaviors because it preserves a long sequence spanning marine isotope stages (MIS) 5, 4, and 3 and is associated with robust records of paleoenvironmental change. The analysis presented here addresses the question, what, if any, lithic assemblage traits at PP5-6 represent changing behavioral responses to the MIS 5-4-3 interglacial-glacial cycle? It statistically evaluates changes in 93 traits with no a priori assumptions about which traits may significantly associate with MIS. In contrast to other studies that claim that there is little relationship between broad-scale patterns of climate change and lithic technology, we identified the following characteristics that are associated with MIS 4: increased use of quartz, increased evidence for outcrop sources of quartzite and silcrete, increased evidence for earlier stages of reduction in silcrete, evidence for increased flaking efficiency in all raw material types, and changes in tool types and function for silcrete. Based on these results, we suggest that foragers responded to MIS 4 glacial environmental conditions at PP5-6 with increased population or group sizes, 'place provisioning', longer and/or more intense site occupations, and decreased residential mobility. Several other traits, including silcrete frequency, do not exhibit an association with MIS. Backed pieces, once they appear in the PP5-6 record during MIS 4, persist through MIS 3. Changing paleoenvironments explain some, but not all temporal technological variability at PP5-6.Social Science and Humanities Research Council of Canada; NORAM; American-Scandinavian Foundation; Fundacao para a Ciencia e Tecnologia [SFRH/BPD/73598/2010]; IGERT [DGE 0801634]; Hyde Family Foundations; Institute of Human Origins; National Science Foundation [BCS-9912465, BCS-0130713, BCS-0524087, BCS-1138073]; John Templeton Foundation to the Institute of Human Origins at Arizona State Universit

Bioinformatics challenges and potentialities in studying extreme environments

Cold environments are populated by organisms able to contravene deleterious effects of low temperature by diverse adaptive strategies, including the production of ice binding proteins (IBPs) that inhibit the growth of ice crystals inside and outside cells. We describe the properties of such a protein (EfcIBP) identified in the metagenome of an Antarctic biological consortium composed of the ciliate Euplotes focardii and psychrophilic non-cultured bacteria. Recombinant EfcIBP can resist freezing without any conformational damage and is moderately heat stable, with a midpoint temperature of 66.4 degrees C. Tested for its effects on ice, EfcIBP shows an unusual combination of properties not reported in other bacterial IBPs. First, it is one of the best-performing IBPs described to date in the inhibition of ice recrystallization, with effective concentrations in the nanomolar range. Moreover, EfcIBP has thermal hysteresis activity (0.53 degrees C at 50 mu M) and it can stop a crystal from growing when held at a constant temperature within the thermal hysteresis gap. EfcIBP protects purified proteins and bacterial cells from freezing damage when exposed to challenging temperatures. EfcIBP also possesses a potential N-terminal signal sequence for protein transport and a DUF3494 domain that is common to secreted IBPs. These features lead us to hypothesize that the protein is either anchored at the outer cell surface or concentrated around cells to provide survival advantage to the whole cell consortium

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

Jet energy measurement with the ATLAS detector in proton-proton collisions at root s=7 TeV

The jet energy scale and its systematic uncertainty are determined for jets measured with the ATLAS detector at the LHC in proton-proton collision data at a centre-of-mass energy of √s = 7TeV corresponding to an integrated luminosity of 38 pb-1. Jets are reconstructed with the anti-kt algorithm with distance parameters R=0. 4 or R=0. 6. Jet energy and angle corrections are determined from Monte Carlo simulations to calibrate jets with transverse momenta pT≥20 GeV and pseudorapidities {pipe}η{pipe}<4. 5. The jet energy systematic uncertainty is estimated using the single isolated hadron response measured in situ and in test-beams, exploiting the transverse momentum balance between central and forward jets in events with dijet topologies and studying systematic variations in Monte Carlo simulations. The jet energy uncertainty is less than 2. 5 % in the central calorimeter region ({pipe}η{pipe}<0. 8) for jets with 60≤pT<800 GeV, and is maximally 14 % for pT<30 GeV in the most forward region 3. 2≤{pipe}η{pipe}<4. 5. The jet energy is validated for jet transverse momenta up to 1 TeV to the level of a few percent using several in situ techniques by comparing a well-known reference such as the recoiling photon pT, the sum of the transverse momenta of tracks associated to the jet, or a system of low-pT jets recoiling against a high-pT jet. More sophisticated jet calibration schemes are presented based on calorimeter cell energy density weighting or hadronic properties of jets, aiming for an improved jet energy resolution and a reduced flavour dependence of the jet response. The systematic uncertainty of the jet energy determined from a combination of in situ techniques is consistent with the one derived from single hadron response measurements over a wide kinematic range. The nominal corrections and uncertainties are derived for isolated jets in an inclusive sample of high-pT jets. Special cases such as event topologies with close-by jets, or selections of samples with an enhanced content of jets originating from light quarks, heavy quarks or gluons are also discussed and the corresponding uncertainties are determined. © 2013 CERN for the benefit of the ATLAS collaboration

Stage-Specific Expression Profiling of Drosophila Spermatogenesis Suggests that Meiotic Sex Chromosome Inactivation Drives Genomic Relocation of Testis-Expressed Genes

In Drosophila, genes expressed in males tend to accumulate on autosomes and are underrepresented on the X chromosome. In particular, genes expressed in testis have been observed to frequently relocate from the X chromosome to the autosomes. The inactivation of X-linked genes during male meiosis (i.e., meiotic sex chromosome inactivation—MSCI) was first proposed to explain male sterility caused by X-autosomal translocation in Drosophila, and more recently it was suggested that MSCI might provide the conditions under which selection would favor the accumulation of testis-expressed genes on autosomes. In order to investigate the impact of MSCI on Drosophila testis-expressed genes, we performed a global gene expression analysis of the three major phases of D. melanogaster spermatogenesis: mitosis, meiosis, and post-meiosis. First, we found evidence supporting the existence of MSCI by comparing the expression levels of X- and autosome-linked genes, finding the former to be significantly reduced in meiosis. Second, we observed that the paucity of X-linked testis-expressed genes was restricted to those genes highly expressed in meiosis. Third, we found that autosomal genes relocated through retroposition from the X chromosome were more often highly expressed in meiosis in contrast to their X-linked parents. These results suggest MSCI as a general mechanism affecting the evolution of some testis-expressed genes

Ecological Adaptation of Diverse Honey Bee (Apis mellifera) Populations

BACKGROUND: Honey bees are complex eusocial insects that provide a critical contribution to human agricultural food production. Their natural migration has selected for traits that increase fitness within geographical areas, but in parallel their domestication has selected for traits that enhance productivity and survival under local conditions. Elucidating the biochemical mechanisms of these local adaptive processes is a key goal of evolutionary biology. Proteomics provides tools unique among the major 'omics disciplines for identifying the mechanisms employed by an organism in adapting to environmental challenges. RESULTS: Through proteome profiling of adult honey bee midgut from geographically dispersed, domesticated populations combined with multiple parallel statistical treatments, the data presented here suggest some of the major cellular processes involved in adapting to different climates. These findings provide insight into the molecular underpinnings that may confer an advantage to honey bee populations. Significantly, the major energy-producing pathways of the mitochondria, the organelle most closely involved in heat production, were consistently higher in bees that had adapted to colder climates. In opposition, up-regulation of protein metabolism capacity, from biosynthesis to degradation, had been selected for in bees from warmer climates. CONCLUSIONS: Overall, our results present a proteomic interpretation of expression polymorphisms between honey bee ecotypes and provide insight into molecular aspects of local adaptation or selection with consequences for honey bee management and breeding. The implications of our findings extend beyond apiculture as they underscore the need to consider the interdependence of animal populations and their agro-ecological context

Cognition based bTBI mechanistic criteria; a tool for preventive and therapeutic innovations

Blast-induced traumatic brain injury has been associated with neurodegenerative and neuropsychiatric disorders. To date, although damage due to oxidative stress appears to be important, the specific mechanistic causes of such disorders remain elusive. Here, to determine the mechanical variables governing the tissue damage eventually cascading into cognitive deficits, we performed a study on the mechanics of rat brain under blast conditions. To this end, experiments were carried out to analyse and correlate post-injury oxidative stress distribution with cognitive deficits on a live rat exposed to blast. A computational model of the rat head was developed from imaging data and validated against in vivo brain displacement measurements. The blast event was reconstructed in silico to provide mechanistic thresholds that best correlate with cognitive damage at the regional neuronal tissue level, irrespectively of the shape or size of the brain tissue types. This approach was leveraged on a human head model where the prediction of cognitive deficits was shown to correlate with literature findings. The mechanistic insights from this work were finally used to propose a novel helmet design roadmap and potential avenues for therapeutic innovations against blast traumatic brain injury

Analysis of ecological thresholds in a temperate forest undergoing dieback.

Positive feedbacks in drivers of degradation can cause threshold responses in natural ecosystems. Though threshold responses have received much attention in studies of aquatic ecosystems, they have been neglected in terrestrial systems, such as forests, where the long time-scales required for monitoring have impeded research. In this study we explored the role of positive feedbacks in a temperate forest that has been monitored for 50 years and is undergoing dieback, largely as a result of death of the canopy dominant species (Fagus sylvatica, beech). Statistical analyses showed strong non-linear losses in basal area for some plots, while others showed relatively gradual change. Beech seedling density was positively related to canopy openness, but a similar relationship was not observed for saplings, suggesting a feedback whereby mortality in areas with high canopy openness was elevated. We combined this observation with empirical data on size- and growth-mediated mortality of trees to produce an individual-based model of forest dynamics. We used this model to simulate changes in the structure of the forest over 100 years under scenarios with different juvenile and mature mortality probabilities, as well as a positive feedback between seedling and mature tree mortality. This model produced declines in forest basal area when critical juvenile and mature mortality probabilities were exceeded. Feedbacks in juvenile mortality caused a greater reduction in basal area relative to scenarios with no feedback. Non-linear, concave declines of basal area occurred only when mature tree mortality was 3-5 times higher than rates observed in the field. Our results indicate that the longevity of trees may help to buffer forests against environmental change and that the maintenance of old, large trees may aid the resilience of forest stands. In addition, our work suggests that dieback of forests may be avoidable providing pressures on mature and juvenile trees do not pass critical thresholds