Differences in whole-body fat oxidation kinetics between cycling and running.

This study aimed to quantitatively describe and compare whole-body fat oxidation kinetics in cycling and running using a sinusoidal mathematical model (SIN). Thirteen moderately trained individuals (7 men and 6 women) performed two graded exercise tests, with 3-min stages and 1 km h(-1) (or 20 W) increment, on a treadmill and on a cycle ergometer. Fat oxidation rates were determined using indirect calorimetry and plotted as a function of exercise intensity. The SIN model, which includes three independent variables (dilatation, symmetry and translation) that account for main quantitative characteristics of kinetics, provided a mathematical description of fat oxidation kinetics and allowed for determination of the intensity (Fat(max)) that elicits maximal fat oxidation (MFO). While the mean fat oxidation kinetics in cycling formed a symmetric parabolic curve, the mean kinetics during running was characterized by a greater dilatation (i.e., widening of the curve, P < 0.001) and a rightward asymmetry (i.e., shift of the peak of the curve to higher intensities, P = 0.01). Fat(max) was significantly higher in running compared with cycling (P < 0.001), whereas MFO was not significantly different between modes of exercise (P = 0.36). This study showed that the whole-body fat oxidation kinetics during running was characterized by a greater dilatation and a rightward asymmetry compared with cycling. The greater dilatation may be mainly related to the larger muscle mass involved in running while the rightward asymmetry may be induced by the specific type of muscle contraction.

Uncertainties regarding dengue modeling in Rio de Janeiro, Brazil

Dengue fever is currently the most important arthropod-borne viral disease in Brazil. Mathematical modeling of disease dynamics is a very useful tool for the evaluation of control measures. To be used in decision-making, however, a mathematical model must be carefully parameterized and validated with epidemiological and entomological data. In this work, we developed a simple dengue model to answer three questions: (i) which parameters are worth pursuing in the field in order to develop a dengue transmission model for Brazilian cities; (ii) how vector density spatial heterogeneity influences control efforts; (iii) with a degree of uncertainty, what is the invasion potential of dengue virus type 4 (DEN-4) in Rio de Janeiro city. Our model consists of an expression for the basic reproductive number (R0) that incorporates vector density spatial heterogeneity. To deal with the uncertainty regarding parameter values, we parameterized the model using a priori probability density functions covering a range of plausible values for each parameter. Using the Latin Hypercube Sampling procedure, values for the parameters were generated. We conclude that, even in the presence of vector spatial heterogeneity, the two most important entomological parameters to be estimated in the field are the mortality rate and the extrinsic incubation period. The spatial heterogeneity of the vector population increases the risk of epidemics and makes the control strategies more complex. At last, we conclude that Rio de Janeiro is at risk of a DEN-4 invasion. Finally, we stress the point that epidemiologists, mathematicians, and entomologists need to interact more to find better approaches to the measuring and interpretation of the transmission dynamics of arthropod-borne diseases.

A Novel Pulse-Chase SILAC Strategy Measures Changes in Protein Decay and Synthesis Rates Induced by Perturbation of Proteostasis with an Hsp90 Inhibitor.

Standard proteomics methods allow the relative quantitation of levels of thousands of proteins in two or more samples. While such methods are invaluable for defining the variations in protein concentrations which follow the perturbation of a biological system, they do not offer information on the mechanisms underlying such changes. Expanding on previous work [1], we developed a pulse-chase (pc) variant of SILAC (stable isotope labeling by amino acids in cell culture). pcSILAC can quantitate in one experiment and for two conditions the relative levels of proteins newly synthesized in a given time as well as the relative levels of remaining preexisting proteins. We validated the method studying the drug-mediated inhibition of the Hsp90 molecular chaperone, which is known to lead to increased synthesis of stress response proteins as well as the increased decay of Hsp90 "clients". We showed that pcSILAC can give information on changes in global cellular proteostasis induced by treatment with the inhibitor, which are normally not captured by standard relative quantitation techniques. Furthermore, we have developed a mathematical model and computational framework that uses pcSILAC data to determine degradation constants kd and synthesis rates Vs for proteins in both control and drug-treated cells. The results show that Hsp90 inhibition induced a generalized slowdown of protein synthesis and an increase in protein decay. Treatment with the inhibitor also resulted in widespread protein-specific changes in relative synthesis rates, together with variations in protein decay rates. The latter were more restricted to individual proteins or protein families than the variations in synthesis. Our results establish pcSILAC as a viable workflow for the mechanistic dissection of changes in the proteome which follow perturbations. Data are available via ProteomeXchange with identifier PXD000538.

Radiation induced apoptosis and initial DNA damage are inversely related in locally advanced breast cancer patients.

Background. DNA-damage assays, quantifying the initial number of DNA double-strand breaks induced by radiation, have been proposed as a predictive test for radiation-induced toxicity. Determination of radiation-induced apoptosis in peripheral blood lymphocytes by flow cytometry analysis has also been proposed as an approach for predicting normal tissue responses following radiotherapy. The aim of the present study was to explore the association between initial DNA damage, estimated by the number of double-strand breaks induced by a given radiation dose, and the radio-induced apoptosis rates observed. Methods. Peripheral blood lymphocytes were taken from 26 consecutive patients with locally advanced breast carcinoma. Radiosensitivity of lymphocytes was quantified as the initial number of DNA double-strand breaks induced per Gy and per DNA unit (200 Mbp). Radio-induced apoptosis at 1, 2 and 8 Gy was measured by flow cytometry using annexin V/propidium iodide. Results. Radiation-induced apoptosis increased in order to radiation dose and data fitted to a semi logarithmic mathematical model. A positive correlation was found among radio-induced apoptosis values at different radiation doses: 1, 2 and 8 Gy (p < 0.0001 in all cases). Mean DSB/Gy/DNA unit obtained was 1.70 ± 0.83 (range 0.63-4.08; median, 1.46). A statistically significant inverse correlation was found between initial damage to DNA and radio-induced apoptosis at 1 Gy (p = 0.034). A trend toward 2 Gy (p = 0.057) and 8 Gy (p = 0.067) was observed after 24 hours of incubation. Conclusions. An inverse association was observed for the first time between these variables, both considered as predictive factors to radiation toxicity.

The demographic benefits of belligerence and bravery: defeated group repopulation or victorious group size expansion?

Intraspecific coalitional aggression between groups of individuals is a widespread trait in the animal world. It occurs in invertebrates and vertebrates, and is prevalent in humans. What are the conditions under which coalitional aggression evolves in natural populations? In this article, I develop a mathematical model delineating conditions where natural selection can favor the coevolution of belligerence and bravery between small-scale societies. Belligerence increases an actor's group probability of trying to conquer another group and bravery increase the actors's group probability of defeating an attacked group. The model takes into account two different types of demographic scenarios that may lead to the coevolution of belligerence and bravery. Under the first, the fitness benefits driving the coevolution of belligerence and bravery come through the repopulation of defeated groups by fission of victorious ones. Under the second demographic scenario, the fitness benefits come through a temporary increase in the local carrying capacity of victorious groups, after transfer of resources from defeated groups to victorious ones. The analysis of the model suggests that the selective pressures on belligerence and bravery are stronger when defeated groups can be repopulated by victorious ones. The analysis also suggests that, depending on the shape of the contest success function, costly bravery can evolve in groups of any size.

The risk of acquiring the new influenza A(H1N1) for Brazilian travelers to Chile, Argentina and the USA

We estimate the risk of acquiring the new influenza A(H1N1) for Brazilian travelers to Chile, Argentina and the USA. This is done by a mathematical model that quantifies the intensity of transmission of the new virus in those countries and the probability that one individual has of acquiring the influenza depending on the date of arrival and time spent in the area. The maximum estimated risk reached 7.5 cases per 10,000 visitors to Chile, 17 cases per 10,000 travelers to Argentina and 23 cases per 10,000 travelers to the USA. The estimated number of imported cases until 27 July is 57 ± 9 from Chile, 136 ± 27 from the USA and 301 ± 21 from Argentina, which are in accord with the official figures. Estimating the number of imported cases was particularly important for the moment of the disease introduction into this country, but it will certainly be important again as a tool to calculate the number of future imported cases from northern countries in our next inter-epidemic season, were imported cases can constitute again the majority of the new influenza burden to the Brazilian health services.

Projecting the effectiveness of RotaTeq® against rotavirus-related hospitalisations in Brazil

RotaTeq® (Merck & Company, Inc, Whitehouse Station, NJ, USA) is an oral pentavalent rotavirus vaccine (RV5) that has shown high and consistent efficacy in preventing rotavirus gastroenteritis (RGE) in randomised clinical trials previously conducted in industrialised countries with high medical care resources. To date, the efficacy and effectiveness data for RV5 are available in some Latin American countries, but not Brazil. In this analysis, we projected the effectiveness of RV5 in terms of the percentage reduction in RGE-related hospitalisations among children less than five years of age in four regions of Brazil, using a previously validated mathematical model. The model inputs included hospital-based rotavirus surveillance data from Goiânia, Porto Alegre, Salvador and São Paulo from 2005-2006, which provided the proportions of rotavirus attributable to serotypes G1, G2, G3, G4 and G9, and published rotavirus serotype-specific efficacy from the Rotavirus Efficacy and Safety Trial. The model projected an overall percentage reduction of 93% in RGE-related hospitalisations, with an estimated annual reduction in RGE-related hospitalisations between 42,991-77,383 in the four combined regions of Brazil. These results suggest that RV5 could substantially prevent RGE-related hospitalisations in Brazil.

Potential impact of a presumed increase in the biting activity of dengue-virus-infected Aedes aegypti (Diptera: Culicidae) females on virus transmission dynamics

Recently, we showed that infection with dengue virus increases the locomotor activity of Aedes aegypti females. We speculate that the observed increased locomotor activity could potentially increase the chances of finding a suitable host and, as a consequence, the relative biting rate of infected mosquitoes. We used a mathematical model to investigate the impact of the increased locomotor activity by assuming that this activity translated into an increased biting rate for infected mosquitoes. The results show that the increased biting rate resulted in dengue outbreaks with greater numbers of primary and secondary infections and with more severe biennial epidemics.

Stochastic optimal bid to electricity markets with environmental risk constraints

There are many factors that influence the day-ahead market bidding strategies of a generation company (GenCo) in the current energy market framework. Environmental policy issues have become more and more important for fossil-fuelled power plants and they have to be considered in their management, giving rise to emission limitations. This work allows to investigate the influence of both the allowances and emission reduction plan, and the incorporation of the derivatives medium-term commitments in the optimal generation bidding strategy to the day-ahead electricity market. Two different technologies have been considered: the coal thermal units, high-emission technology, and the combined cycle gas turbine units, low-emission technology. The Iberian Electricity Market and the Spanish National Emissions and Allocation Plans are the framework to deal with the environmental issues in the day-ahead market bidding strategies. To address emission limitations, some of the standard risk management methodologies developed for financial markets, such as Value-at-Risk (VaR) and Conditional Value-at-Risk (CVaR), have been extended. This study offers to electricity generation utilities a mathematical model to determinate the individual optimal generation bid to the wholesale electricity market, for each one of their generation units that maximizes the long-run profits of the utility abiding by the Iberian Electricity Market rules, the environmental restrictions set by the EU Emission Trading Scheme, as well as the restrictions set by the Spanish National Emissions Reduction Plan. The economic implications for a GenCo of including the environmental restrictions of these National Plans are analyzed and the most remarkable results will be presented.

Identificació i caracterització funcional dels complexes proteics dels receptors BRL1 i BRL3 en el teixit vascular d’Arabidopsis Thaliana

Els esteroids juguen papers clau en el creixement I el desenvolupament d’eucariotes multicel•lulars. En plantes, aquestes hormones, anomenades Brassinosteroides (BRs), estan involucrades en una gran varietat de processos biològics essencials per a les plantes. S’han descrit anteriorment dos receptors de BRs del tipus Leucine Rich Repeat Receptor Like Kinase LRR-RLK, BRASSINOSTEROID RECEPTOR LIKE 1 i 3 (BRL1 i BRL3 respectivalemt) que són homòlegs al receptor principal BRI1 i són necessaris pel desenvolupament vascular. Tot i que els principals components de la senyal ja han estat identificats pel seu homòleg més pròxim, el receptor BRI1, els complexes de BRL1 i BRL3 juntament amb els candidats co-receptors així com els components de la ruta de senyalització encara no han sigut identificats. Per tal d’entendre millor la funció molecular d’aquests receptors de BRs en la planta aquesta tesis doctoral planteja dues aproximacions: com a primera aproximació, vaig realitzar un estudi fenotípic del desenvolupament del teixit vascular a la planta model Arabidopsis thaliana (Arabidopsis). Disposant d'una amplia bateria de mutants de síntesis de la hormona i senyalització del receptor BRI1, vam analitzar quantitativament el seu patró vascular a la tija d'Arabidopsis. Vam establir els paràmetres en les plantes silvestres [Col-0 wild type, (WT)] i els vam analitzar a tots i cadascun dels mutants. Això conjuntament amb una col•laboració amb la Dr. Marta Ibañes, física de la Universitat de Barcelona que va construir un model matemàtic per simular la formació del patró vascular ens va permetre el•laborar una hipòtesis que vam demostrar experimentalment i va ser publicada a la revista PNAS. Posteriorment vam observar que les plantes knock-out d'aquests dos receptors BRL1 y BRL3 a diferència de BRI1, no tenien cap fenotip obvi en el teixit vascular de la planta adulta. Així, a continuació, per entendre quina necessitat té la planta de disposar de tres receptors tant altament homòlegs que poden percebre la mateixa hormona, vam utilitzar una aproximació bioquímica en col•laboració amb el Prof. de Vries de la Universitat de Wageningen (Holanda) per tal de purificar els complexes dels receptors in vivo i els seus interactors. Això ens ha permès entendre millor el paper funcional d'aquests receptors en la planta. Els resultats d’aquests experiments estan resumits en un article en preparació que aviat estarà en revisió.

Chronic potassium depletion increases adrenal progesterone production that is necessary for efficient renal retention of potassium.

Modern dietary habits are characterized by high-sodium and low-potassium intakes, each of which was correlated with a higher risk for hypertension. In this study, we examined whether long-term variations in the intake of sodium and potassium induce lasting changes in the plasma concentration of circulating steroids by developing a mathematical model of steroidogenesis in mice. One finding of this model was that mice increase their plasma progesterone levels specifically in response to potassium depletion. This prediction was confirmed by measurements in both male mice and men. Further investigation showed that progesterone regulates renal potassium handling both in males and females under potassium restriction, independent of its role in reproduction. The increase in progesterone production by male mice was time dependent and correlated with decreased urinary potassium content. The progesterone-dependent ability to efficiently retain potassium was because of an RU486 (a progesterone receptor antagonist)-sensitive stimulation of the colonic hydrogen, potassium-ATPase (known as the non-gastric or hydrogen, potassium-ATPase type 2) in the kidney. Thus, in males, a specific progesterone concentration profile induced by chronic potassium restriction regulates potassium balance.

Metabolic Flux and Compartmentation Analysis in the Brain In vivo.

Through significant developments and progresses in the last two decades, in vivo localized nuclear magnetic resonance spectroscopy (MRS) became a method of choice to probe brain metabolic pathways in a non-invasive way. Beside the measurement of the total concentration of more than 20 metabolites, (1)H MRS can be used to quantify the dynamics of substrate transport across the blood-brain barrier by varying the plasma substrate level. On the other hand, (13)C MRS with the infusion of (13)C-enriched substrates enables the characterization of brain oxidative metabolism and neurotransmission by incorporation of (13)C in the different carbon positions of amino acid neurotransmitters. The quantitative determination of the biochemical reactions involved in these processes requires the use of appropriate metabolic models, whose level of details is strongly related to the amount of data accessible with in vivo MRS. In the present work, we present the different steps involved in the elaboration of a mathematical model of a given brain metabolic process and its application to the experimental data in order to extract quantitative brain metabolic rates. We review the recent advances in the localized measurement of brain glucose transport and compartmentalized brain energy metabolism, and how these reveal mechanistic details on glial support to glutamatergic and GABAergic neurons.

N-dimensional dynamical systems exploiting instabilities in full

We report experimental and numerical results showing how certain N-dimensional dynamical systems are able to exhibit complex time evolutions based on the nonlinear combination of N-1 oscillation modes. The experiments have been done with a family of thermo-optical systems of effective dynamical dimension varying from 1 to 6. The corresponding mathematical model is an N-dimensional vector field based on a scalar-valued nonlinear function of a single variable that is a linear combination of all the dynamic variables. We show how the complex evolutions appear associated with the occurrence of successive Hopf bifurcations in a saddle-node pair of fixed points up to exhaust their instability capabilities in N dimensions. For this reason the observed phenomenon is denoted as the full instability behavior of the dynamical system. The process through which the attractor responsible for the observed time evolution is formed may be rather complex and difficult to characterize. Nevertheless, the well-organized structure of the time signals suggests some generic mechanism of nonlinear mode mixing that we associate with the cluster of invariant sets emerging from the pair of fixed points and with the influence of the neighboring saddle sets on the flow nearby the attractor. The generation of invariant tori is likely during the full instability development and the global process may be considered as a generalized Landau scenario for the emergence of irregular and complex behavior through the nonlinear superposition of oscillatory motions

Beyond GDP: Methodological and measurement issues in redefining “wealth”

In the last decades; a growing stock of literature has been devoted to the criticism of GDP as an indicator of societal wealth. A relevant question is: what are the perspectives to build, on the existing knowledge and consensus, alternative measures of prosperity? A starting point may be to connect well-being research agenda with the sustainability one. However, there is no doubt that there is a lot of complexity and fuzziness inherent in multidimensional concepts such as sustainability and well-being. This article analyses the theoretical foundations and the empirical validity of some multidimensional technical tools that can be used for well-being evaluation and assessment. Of course one should not forget that policy conclusions derived through any mathematical model depend also on the conceptual framework used, i.e. which representation of reality (and thus which societal values and interests) has been considered.

Minimizing recombinations in consensus networks for phylogeographic studies

Background: We address the problem of studying recombinational variations in (human) populations. In this paper, our focus is on one computational aspect of the general task: Given two networks G1 and G2, with both mutation and recombination events, defined on overlapping sets of extant units the objective is to compute a consensus network G3 with minimum number of additional recombinations. We describe a polynomial time algorithm with a guarantee that the number of computed new recombination events is within ϵ = sz(G1, G2) (function sz is a well-behaved function of the sizes and topologies of G1 and G2) of the optimal number of recombinations. To date, this is the best known result for a network consensus problem.Results: Although the network consensus problem can be applied to a variety of domains, here we focus on structure of human populations. With our preliminary analysis on a segment of the human Chromosome X data we are able to infer ancient recombinations, population-specific recombinations and more, which also support the widely accepted 'Out of Africa' model. These results have been verified independently using traditional manual procedures. To the best of our knowledge, this is the first recombinations-based characterization of human populations. Conclusion: We show that our mathematical model identifies recombination spots in the individual haplotypes; the aggregate of these spots over a set of haplotypes defines a recombinational landscape that has enough signal to detect continental as well as population divide based on a short segment of Chromosome X. In particular, we are able to infer ancient recombinations, population-specific recombinations and more, which also support the widely accepted 'Out of Africa' model. The agreement with mutation-based analysis can be viewed as an indirect validation of our results and the model. Since the model in principle gives us more information embedded in the networks, in our future work, we plan to investigate more non-traditional questions via these structures computed by our methodology.