Dynamical modeling and analysis of large cellular regulatory networks.

The dynamical analysis of large biological regulatory networks requires the development of scalable methods for mathematical modeling. Following the approach initially introduced by Thomas, we formalize the interactions between the components of a network in terms of discrete variables, functions, and parameters. Model simulations result in directed graphs, called state transition graphs. We are particularly interested in reachability properties and asymptotic behaviors, which correspond to terminal strongly connected components (or "attractors") in the state transition graph. A well-known problem is the exponential increase of the size of state transition graphs with the number of network components, in particular when using the biologically realistic asynchronous updating assumption. To address this problem, we have developed several complementary methods enabling the analysis of the behavior of large and complex logical models: (i) the definition of transition priority classes to simplify the dynamics; (ii) a model reduction method preserving essential dynamical properties, (iii) a novel algorithm to compact state transition graphs and directly generate compressed representations, emphasizing relevant transient and asymptotic dynamical properties. The power of an approach combining these different methods is demonstrated by applying them to a recent multilevel logical model for the network controlling CD4+ T helper cell response to antigen presentation and to a dozen cytokines. This model accounts for the differentiation of canonical Th1 and Th2 lymphocytes, as well as of inflammatory Th17 and regulatory T cells, along with many hybrid subtypes. All these methods have been implemented into the software GINsim, which enables the definition, the analysis, and the simulation of logical regulatory graphs.

The value of selected in vitro and in silico methods to predict acute oral toxicity in a regulatory context: results from the European Project ACuteTox.

ACuteTox is a project within the 6th European Framework Programme which had as one of its goals to develop, optimise and prevalidate a non-animal testing strategy for predicting human acute oral toxicity. In its last 6 months, a challenging exercise was conducted to assess the predictive capacity of the developed testing strategies and final identification of the most promising ones. Thirty-two chemicals were tested blind in the battery of in vitro and in silico methods selected during the first phase of the project. This paper describes the classification approaches studied: single step procedures and two step tiered testing strategies. In summary, four in vitro testing strategies were proposed as best performing in terms of predictive capacity with respect to the European acute oral toxicity classification. In addition, a heuristic testing strategy is suggested that combines the prediction results gained from the neutral red uptake assay performed in 3T3 cells, with information on neurotoxicity alerts identified by the primary rat brain aggregates test method. Octanol-water partition coefficients and in silico prediction of intestinal absorption and blood-brain barrier passage are also considered. This approach allows to reduce the number of chemicals wrongly predicted as not classified (LD50>2000 mg/kg b.w.).

Time-Course In Vivo Trafficking Pattern and Effector Function of Donor-Specific Regulatory T Cells in Response to an Allograft.

Background: Experimental data have suggested that adoptive transfer of CD4+CD25+Foxp3+ regulatory T cells (Tregs), capable of controlling immune responses to specifi c auto- or alloantigens, could be used as a therapeutic strategy to promote specifi c tolerance in T-cell mediated diseases and in organ transplantation (Tx). However, before advocating the application of immunotherapy with Tregs in Tx, we need to improve our understanding of their in vivo homeostasis, traffi cking pattern and effector function in response to alloantigens. Methods : Donor-antigen specifi c murine Tregs were generated and characterized in vitro following our described protocols. Using an adoptive transfer and skin allotransplantation model, we have analyzed the in vivo expansion and homing of fl uorescent-labeled effector T cells (Teff) and Tregs, at different time-points after Tx, using fl ow-cytometry as well as fl uorescence microscopy techniques. Results: Tregs expressed CD62L, CCR7 and CD103 allowing their homing into lymphoid and non-lymphoid tissues (gut, skin) after intravenous injection. While hyporesponsive to TCR stimulation in vitro, transferred Tregs survived, migrated to secondary lymphoid organs and preferentially expanded within the allograft draining lymph nodes. Furthermore, Foxp3+ cells could be detected inside the allograft as early as day 3-5 after Tx. At a much later time-point (day 60 after Tx), graft-infi ltrating Foxp3+ cells were also detectable in tolerant recipients. When transferred alone, CD4+CD25- Teff cells expanded within secondary lymphoid organs and infi ltrated the allograft by day 3-5 after Tx. The co-transfer of Tregs limited the expansion of alloreactive Teff cells as well as their recruitment into the allograft. The promotion of graft survival observed in the presence of Tregs was in part mediated by the inhibition of the production of effector cytokines by CD4+CD25- T cells. Conclusion: Taken together, our results suggest that the suppression of allograft rejection and the induction of Tx tolerance are in part dependant on the alloantigendriven homing and expansion of Tregs. Thus, the appropriate localization of Tregs may be critical for their suppressive function in vivo.

Analysis of a finite volume element method for the Stokes problem

In this paper we propose a stabilized conforming finite volume element method for the Stokes equations. On stating the convergence of the method, optimal a priori error estimates in different norms are obtained by establishing the adequate connection between the finite volume and stabilized finite element formulations. A superconvergence result is also derived by using a postprocessing projection method. In particular, the stabilization of the continuous lowest equal order pair finite volume element discretization is achieved by enriching the velocity space with local functions that do not necessarily vanish on the element boundaries. Finally, some numerical experiments that confirm the predicted behavior of the method are provided.

A hormonal regulatory module that provides flexibility to tropic responses.

Plants orient their growth depending on directional stimuli such as light and gravity, in a process known as tropic response. Tropisms result from asymmetrical accumulation of auxin across the responding organ relative to the direction of the stimulus, which causes differential growth rates on both sides of the organ. Here, we show that gibberellins (GAs) attenuate the gravitropic reorientation of stimulated hypocotyls of dark-grown Arabidopsis (Arabidopsis thaliana) seedlings. We show that the modulation occurs through induction of the expression of the negative regulator of auxin signaling INDOLE-3-ACETIC ACID INDUCIBLE19/MASSUGU2. The biological significance of this regulatory mechanism involving GAs and auxin seems to be the maintenance of a high degree of flexibility in tropic responses. This notion is further supported by observations that GA-deficient seedlings showed a much lower variance in the response to gravity compared to wild-type seedlings and that the attenuation of gravitropism by GAs resulted in an increased phototropic response. This suggests that the interplay between auxin and GAs may be particularly important for plant orientation under competing tropic stimuli.

Sedimentology, geochemistry and mineralogy of the Paleocene-Eocene thermal maximum (PETM) : sediment records from Egypt, India and Spain

A gradual increase in Earth's surface temperatures marking the transition from the late Paleocene to early Eocene (55.8±0.2Ma), represents an extraordinary warming event known as Paleocene-Eocene Thermal Maximum (PETM). Both marine and continental sedimentary records during this period reveal evidences for the massive injection of isotopically light carbon. The carbon dioxide injection from multiple potential sources may have triggered the global warming. The importance of the PETM studies is due to the fact that the PETM bears some striking resemblances to the human-caused climate change unfolding today. Most notably, the culprit behind it was a massive injection of heat-trapping greenhouse gases into the atmosphere and oceans, comparable in volume to what our persistent burning of fossil fuels could deliver in coming centuries. The exact knowledge of what went on during the PETM could help us to foresee the future climate change. The response of the oceanic and continental environments to the PETM is different. Many factors might control the response of the environments to the PETM such as paleogeography, paleotopography, paleoenvironment, and paleodepth. To better understand the mechanisms triggering PETM events, two different environments were studied: 1) shallow marine to inner shelf environment (Wadi Nukhul, Sinai; and the Dababiya GSSP, Luxor, Egypt), and 2) terrestrial environments (northwestern India lignite mines) representing wetland, and fluvial environments (Esplugafreda, Spain) both highlighting the climatic changes observed in continental conditions. In the marine realm, the PETM is characterized by negative ö13Ccar and ô13Corg excursions and shifts in Ô15N to ~0%o values above the P/E boundary and persisting along the interval suggesting a bloom and high production of atmospheric N2-fixers. Decrease in carbonate contents could be due to dissolution and/or dilution by increasing detrital input. High Ti, K and Zr and decreased Si contents at the P/E boundary indicate high weathering index (CIA), which coincides with significant kaolinite input and suggests intense chemical weathering under humid conditions at the beginning of the PETM. Two anoxic intervals are observed along the PETM. The lower one may be linked to methane released from the continental shelf with no change in the redox proxies, where the upper anoxic to euxinic conditions are revealed by increasing U, Mo, V, Fe and the presence of small size pyrite framboids (2-5fim). Productivity sensitive elements (Cu, Ni, and Cd) show their maximum concentrated within the upper anoxic interval suggesting high productivity in surface water. The obtained data highlight that intense weathering and subsequent nutrient inputs are crucial parameters in the chain of the PETM events, triggering productivity during the recovery phase. In the terrestrial environments, the establishment of wetland conditions and consequence continental climatic shift towards more humid conditions led to migration of modern mammals northward following the extension of the tropical belts. Relative ages of this mammal event based on bio-chemo- and paleomagnetic stratigraphy support a migration path originating from Asia into Europe and North America, followed by later migration from Asia into India and suggests a barrier to migration that is likely linked to the timing of the India-Asia collision. In contrast, at Esplugafereda, northeastern Spain, the terrestrial environment reacted differently. Two significant S13C shifts with the lower one linked to the PETM and the upper corresponding to the Early Eocene Thermal Maximum (ETM2); 180/160 paleothermometry performed on two different soil carbonate nodule reveal a temperature increase of around 8°C during the PETM. The prominent increase in kaolinite content within the PETM is linked to increased runoff and/or weathering of adjacent and coeval soils. These results demonstrate that the PETM coincides globally with extreme climatic fluctuations and that terrestrial environments are very likely to record such climatic changes. - La transition Paléocène-Eocène (55,8±0,2 Ma) est marquée par un réchauffement extraordinaire communément appelé « Paleocene-Eocene Thermal Maximum » (PETM). Les données géochimiques caractérisant les sédiments marins et continentaux de cette période indiquent que ce réchauffement a été déclenché par une augmentation massive de CO2 lié à la déstabilisation des hydrates de méthane stockés le long des marges océaniques. L'étude des événements PETM constitue donc un bon analogue avec le réchauffement actuel. Le volume de CO2 émis durant le PETM est comparable avec le CO2 lié à l'activité actuelle humaine. La compréhension des causes du réchauffement du PETM peut être cruciale pour prévoir et évaluer les conséquences du réchauffement anthropogénique, en particulier les répercussions d'un tel réchauffement sur les domaines continentaux et océaniques. De nombreux facteurs entrent en ligne de compte dans le cas du PETM, tels que la paléogéographie, la paléotopographie et les paléoenvironnement. Pour mieux comprendre les réponses environnementales aux événements du PETM, 2 types d'environnements ont été choisis : (1) le domaine marin ouvert mais relativement peu profond (Wadi Nukhul. Sinai, Dababiya, Luxor, Egypte), (2) le milieu continental marécageux humide (mines de lignite, Inde) et fluviatile, semi-aride (Esplugafreda, Pyrénées espagnoles). Dans le domaine marin, le PETM est caractérisé par des excursions négatives du ô13Ccar et ô13Corg et un shift persistant des valeurs de 815N à ~ 0 %o indiquant une forte activité des organismes (bactéries) fixant l'azote. La diminution des carbonates observée durant le PETM peut-être due à des phénomènes de dissolution ou une augmentation des apports terrigènes. Des taux élevés en Ti, K et Zr et une diminution des montants de Si, reflétés par des valeurs des indices d'altération (CIA) qui coïncident avec une augmentation significative des apports de kaolinite impliquent une altération chimique accrue, du fait de conditions plus humides au début du PETM. Deux événements anoxiques globaux ont été mis en évidence durant le PETM. Le premier, situé dans la partie inférieur du PETM, serait lié à la libération des hydrates de méthane stockés le long des talus continentaux et ne correspond pas à des variations significatives des éléments sensibles aux changements de conditions redox. Le second est caractérisé par une augmentation des éléments U, Mo, V et Fe et la présence de petit framboids de pyrite dont la taille varie entre 2 et 5pm. Le second épisode anoxique est caractérisé par une forte augmentation des éléments sensibles aux changements de la productivité (Cu, Ni et Co), indiquant une augmentation de la productivité dans les eaux de surface. Les données obtenues mettent en évidence le rôle crucial joué par l'altération et les apports en nutriments qui en découlent. Ces paramètres sont cruciaux pour la succession des événements qui ont conduit au PETM, et plus particulièrement l'augmentation de la productivité dans la phase de récupération. Durant le PETM, le milieu continental est caractérisé par l'établissement de conditions humides qui ont facilité voir provoqué la migration des mammifères modernes qui ont suivi le déplacement de ces ceintures climatiques. L'âge de cette migration est basé sur des arguments chimiostratigraphiques (isotopes stables), biostratigraphiques et paléomagnétiques. Les données bibliographiques ainsi que celles que nous avons récoltées en Inde, montrent que les mammifères modernes ont d'abord migré depuis l'Asie vers l'Europe, puis dans le continent Nord américain. Ces derniers ne sont arrivés en Inde que plus tardivement, suggérant que le temps de leur migration est lié à la collision Inde-Asie. Dans le Nord-Est de l'Espagne (Esplugafreda), la réponse du milieu continental aux événements PETM est assez différente. Comme en Inde, deux excursions signicatives en ô13C ont été observées. La première correspond au PETM et la seconde est corrélée avec l'optimum thermique de l'Eocène précoce (ETM2). Les isotopes stables de l'oxygène mesurés 2 différents types de nodules calcaires provenant de paléosols suggère une augmentation de 10°C pendant le PETM. Une augmentation simultanée des taux de kaolinite indique une intensification de l'altération chimique et/ou de l'érosion de sols adjacents. Ces résultats démontrent que le PETM coïncide globalement avec des variations climatiques extrêmes qui sont très aisément reconnaissables dans les dépôts continentaux.

Social insect genomes exhibit dramatic evolution in gene composition and regulation while preserving regulatory features linked to sociality.

Genomes of eusocial insects code for dramatic examples of phenotypic plasticity and social organization. We compared the genomes of seven ants, the honeybee, and various solitary insects to examine whether eusocial lineages share distinct features of genomic organization. Each ant lineage contains ∼4000 novel genes, but only 64 of these genes are conserved among all seven ants. Many gene families have been expanded in ants, notably those involved in chemical communication (e.g., desaturases and odorant receptors). Alignment of the ant genomes revealed reduced purifying selection compared with Drosophila without significantly reduced synteny. Correspondingly, ant genomes exhibit dramatic divergence of noncoding regulatory elements; however, extant conserved regions are enriched for novel noncoding RNAs and transcription factor-binding sites. Comparison of orthologous gene promoters between eusocial and solitary species revealed significant regulatory evolution in both cis (e.g., Creb) and trans (e.g., fork head) for nearly 2000 genes, many of which exhibit phenotypic plasticity. Our results emphasize that genomic changes can occur remarkably fast in ants, because two recently diverged leaf-cutter ant species exhibit faster accumulation of species-specific genes and greater divergence in regulatory elements compared with other ants or Drosophila. Thus, while the "socio-genomes" of ants and the honeybee are broadly characterized by a pervasive pattern of divergence in gene composition and regulation, they preserve lineage-specific regulatory features linked to eusociality. We propose that changes in gene regulation played a key role in the origins of insect eusociality, whereas changes in gene composition were more relevant for lineage-specific eusocial adaptations.

Effects of grapefruit juice on the pharmacokinetics of the enantiomers of methadone

BACKGROUND AND OBJECTIVES: Cytochrome P450 (CYP) 3A4 is the main CYP isozyme involved in methadone metabolism. We investigated the influence of grapefruit juice, which contains inhibitors of intestinal CYP3A, on the steady-state pharmacokinetics of methadone. METHODS: For 5 days, 8 patients undergoing methadone maintenance treatment received 200 mL water or grapefruit juice 30 minutes before and again together with their daily dose of methadone. Blood sampling for R-, S-, and R,S-methadone plasma determination was performed over a 24-hour period. CYP3A activity was determined by measuring the plasma 1'-hydroxymidazolam/midazolam ratio. RESULTS: A decrease in the midazolam ratio was measured in all patients after grapefruit juice (mean +/- SD before grapefruit juice, 9.3 +/- 5.9; mean +/- SD after grapefruit juice, 3.9 +/- 1.2; P <.05). Grapefruit juice led to a mean 17% increase in the area under the curve extrapolated to 24 hours for both enantiomers of methadone (range, 3% to 29% [P <.005]; range, -4% to 37% [P <.05]; and range, 1% to 32% [P <.01]; for R-, S-, and R,S-methadone, respectively). A similar increase in peak level and decrease in apparent clearance were measured with grapefruit juice, whereas time to peak level, terminal half-life, and apparent volume during the terminal phase of R-, S-, and R,S-methadone were not affected by grapefruit juice. No symptom of overmedication was either detected by the clinical staff or reported by the patients. CONCLUSIONS: Grapefruit juice administration is associated with a modest increase in methadone bioavailability, which is not expected to endanger patients. However, it cannot be excluded that a much stronger effect may occur in some patients, and thus grapefruit juice intake is not recommended during methadone maintenance treatment, in particular in patients initiating such a treatment.

Reduced ejection fraction after myocardial infarction: is it sufficient to justify implantation of a defibrillator?

BACKGROUND: Improved survival after prophylactic implantation of a defibrillator in patients with reduced left ventricular ejection fraction (EF) after myocardial infarction (MI) has been demonstrated in patients who experienced remote MIs in the 1990s. The absolute survival benefit conferred by this recommended strategy must be related to the current risk of arrhythmic death, which is evolving. This study evaluates the mortality rate in survivors of MI with impaired left ventricular function and its relation to pre-hospital discharge baseline characteristics. METHODS: The clinical records of patients who had sustained an acute MI between 1999 and 2000 and had been discharged from the hospital with an EF of < or = 40% were included. Baseline characteristics, drug prescriptions, and invasive procedures were recorded. Bivariate and multivariate analyses were performed using a primary end point of total mortality. RESULTS: One hundred sixty-five patients were included. During a median follow-up period of 30 months (interquartile range, 22 to 36 months) 18 patients died. The 1-year and 2-year mortality rates were 6.7% and 8.6%, respectively. Variables reflecting coronary artery disease and its management (ie, prior MI, acute reperfusion, and complete revascularization) had a greater impact on mortality than variables reflecting mechanical dysfunction (ie, EF and Killip class). CONCLUSIONS: The mortality rate among survivors of MIs with reduced EF was substantially lower than that reported in the 1990s. The strong decrease in the arrhythmic risk implies a proportional increase in the number of patients needed to treat with a prophylactic defibrillator to prevent one adverse event. The risk of an event may even be sufficiently low to limit the detectable benefit of defibrillators in patients with the prognostic features identified in our study. This argues for additional risk stratification prior to the prophylactic implantation of a defibrillator.

Abrupt decrease in serum testosterone levels after an oral glucose load in men: implications for screening for hypogonadism.

OBJECTIVE: This study examines the physiological impact of a glucose load on serum testosterone (T) levels in men with varying glucose tolerance (GT). DESIGN: Cross-sectional study. PATIENTS AND METHODS: 74 men (19-74 years, mean 51·4 ± 1·4 years) underwent a standard 75-g oral glucose tolerance test with blood sampling at 0, 30, 60, 90 and 120 min. Fasting serum glucose, insulin, total T (and calculated free T), LH, SHBG, leptin and cortisol were measured. RESULTS: 57% of the men had normal GT, 30% had impaired GT and 13% had newly diagnosed type 2 diabetes. Glucose ingestion was associated with a 25% decrease in mean T levels (delta = -4·2 ± 0·3 nm, P < 0·0001). T levels remained suppressed at 120 min compared with baseline (13·7 ± 0·6 vs 16·5 ± 0·7 nm, P < 0·0001) and did not differ across GT or BMI. Of the 66 men with normal T levels at baseline, 10 (15%) had levels that decreased to the hypogonadal range (<9·7 nm) at one or more time points. SHBG, LH and cortisol levels were unchanged. Leptin levels decreased from baseline at all time points (P < 0·0001). CONCLUSIONS: Glucose ingestion induces a significant reduction in total and free T levels in men, which is similar across the spectrum of glucose tolerance. This decrease in T appears to be because of a direct testicular defect, but the absence of compensatory changes in LH suggests an additional central component. Men found to have low nonfasting T levels should be re-evaluated in the fasting state.

What if there was a stronger pharmaceutical price competition in Spain? When regulation has a similar effect to collusion

This paper examines statins competition in the Spanish pharmaceutical market, where prices are highly regulated, and simulates a situation in which there is unrestricted price competition. A nested logit demand model is estimated with a panel of monthly data for pharmaceuticals prescribed from 1997 to 2005. The simulation indicates that the regulation of prices is similar in its effects to cooperation among producers, since the regulated prices are close to those that would be observed in a scenario of perfect collusion. Freedom to set prices and a regulatory framework with appropriate incentives would result in a general reduction in prices and may make the current veiled competition in the form of discounts to pharmacists become more visible. The decrease in prices would be partially offset by an increase in consumption but the net effect would be an overall decrease in expenditure. The counterfactual set-up would also lead to important changes in the market shares of both manufacturers and active ingredients, and a reversal of generic drugs. Therefore, pro-competitive regulation would be welfare-enhancing but would imply winners and losers.

Systemic modulation of peripheral eosinophilia (air pouch model) in Schistosoma mansoni infection

Schistosoma mansoni infection induces in their hosts a marked and sustained eosinophilia, which is influenced or modulated by complex mechanisms, that vary according to the phase of infection. To address this phenomenon, we used the air pouch (AP) model in control and infected Swiss webster mice, analyzing the cellular, tissue response and local expression of adhesion molecules [CD18 (beta 2-chain), CD44, ICAM-1 (CD54), L-selectin (CD62L), CD49d (alpha 4-chain), LFA1 (CD11a)]. Infected animals were studied at 3 (pre-oviposition phase), 7 (acute phase), and 14 (chronic phase) weeks after infection (5-6 mice/period of infection). Normal mice were age-matched. Results showed that after egg stimulation, compared with matched controls, the infected mice, at each point of infection, showed a lower eosinophil response in the acute (7 weeks) and chronic phase (14 weeks) of infection. However, when the infected mice were in pre-oviposition phase (3 weeks) their eosinophil response surpassed the control ones. In the AP wall of infected mice, a significant decrease in the expression of ICAM-1 and CD44 in fibroblastic-like cells and a reduction in the number of CD18 and CD11a in migratory cells were observed. The other adhesion molecules were negative or weakly expressed. The results indicated that in the air pouch model, in S. mansoni-infected mice: (1) eosinophil response is strikingly down-regulated, during the acute ovular phase; (2) in the pre-oviposition phase, in contrast, it occurs an up-regulatory modulation of eosinophil response, in which the mechanisms are completely unknown; (3) in the chronic phase of the infection, the down modulation of eosinophil response is less pronounced; 4) Down-regulation of adhesion molecules, specially of ICAM-1 appear to be associated with the lower eosinophil response.