A hierarchical approach to model parameter optimization for developmental systems.

In the context of Systems Biology, computer simulations of gene regulatory networks provide a powerful tool to validate hypotheses and to explore possible system behaviors. Nevertheless, modeling a system poses some challenges of its own: especially the step of model calibration is often difficult due to insufficient data. For example when considering developmental systems, mostly qualitative data describing the developmental trajectory is available while common calibration techniques rely on high-resolution quantitative data. Focusing on the calibration of differential equation models for developmental systems, this study investigates different approaches to utilize the available data to overcome these difficulties. More specifically, the fact that developmental processes are hierarchically organized is exploited to increase convergence rates of the calibration process as well as to save computation time. Using a gene regulatory network model for stem cell homeostasis in Arabidopsis thaliana the performance of the different investigated approaches is evaluated, documenting considerable gains provided by the proposed hierarchical approach.

Critical role for the lactate transporter, monocarboxylate transporter 1 (mct1), in the regeneration of peripheral nerves

Axons, and particularly regenerating axons, have high metabolic needs in order to maintain critical functions such as axon transport and membrane depolarization. Though some of the required energy likely comes form extracellular glucose and ATP generated in the soma, we and others hypothesize that some of the energy may be supplied by lactate. Unlike glucose that requires glycolytic enzymes to produce pyruvate, lactate can be converted directly to pyruvate by lactate dehydrogenase and transported into mitochondria for oxidative metabolism. In order to be transported into or out of cells, lactate requires specific monocarboxylate transporters (MCTs), the most abundant of which is MCT1. If MCT1 and lactate are critical for nerve function and regeneration, we hypothesize that MCT1 heterozygote null mice, which appear phenotypically normal despite having approximately 40% MCT1 as compared to wildtype littermate mice, would have reduced capacity for repair following nerve injury. To investigate this, adult MCT1 heterozygote null mice or wild-type mice underwent unilateral sciatic nerve crush in the proximal thigh. We found that regeneration of the sciatic nerve, as measured by recovery of compound muscle action potentials (CMAP) in the lateral plantar muscles following proximal sciatic nerve stimulation, was delayed from a median of 21 days in wildtype mice to 38.5 days in MCT1 heterozygote mice. In fact, half of the MCT1 heterozygote null mice had no recovery of CMAP by the endpoint of the study at 42 days, while all of the wild-type mice had recovered. In addition, the maximal amplitude of CMAP recovery in MCT1 heterozygote mull mice was reduced from a mean of 3 mV to 0.5 mV. As would be expected, the denervated gastrocnemius muscle of MCT1 heterozygote null mice remained atrophic at 42 days compared to wild-type mice. Our experiments show that lactate supplied through MCT1 is necessary for nerve regeneration. Experiments are underway to determine whether loss of MCT1 prevents nerve regrowth directly due to reduced energy supply to axons or indirectly by dysfunctional Schwann cells normally dependent on lactate supply through MCT1.

Evidence from glut2-null mice that glucose is a critical physiological regulator of feeding.

A role for glucose in the control of feeding has been proposed, but its precise physiological importance is unknown. Here, we evaluated feeding behavior in glut2-null mice, which express a transgenic glucose transporter in their beta-cells to rescue insulin secretion (ripglut1;glut2-/- mice). We showed that in the absence of GLUT2, daily food intake was increased and feeding initiation and termination following a fasting period were abnormal. This was accompanied by suppressed regulation of hypothalamic orexigenic and anorexigenic neuropeptides expression during the fast-to-refed transition. In these conditions, however, there was normal regulation of the circulating levels of insulin, leptin, or glucose but a loss of regulation of plasma ghrelin concentrations. To evaluate whether the abnormal feeding behavior was due to suppressed glucose sensing, we evaluated feeding in response to intraperitoneal or intracerebroventricular glucose or 2-deoxy-D-glucose injections. We showed that in GLUT2-null mice, feeding was no longer inhibited by glucose or activated by 2-deoxy-D-glucose injections and the regulation of hypothalamic neuropeptide expression by intracerebroventricular glucose administration was lost. Together, these data demonstrate that absence of GLUT2 suppressed the function of central glucose sensors, which control feeding probably by regulating the hypothalamic melanocortin pathway. Furthermore, inactivation of these glucose sensors causes overeating.

Evaluation of a modified Broström-Gould procedure for treatment of chronic lateral ankle instability: A retrospective study with critical analysis of outcome scoring.

BACKGROUND: Chronic lateral ankle instability accounts for 20% of the ankle injuries. This study evaluates functional outcome of the modified Broström-Gould technique using suture anchors, with 4 different clinical scores. METHODS: A consecutive series of 41 patients were included with a minimum follow-up of one year. The function was assessed using 4 clinical scores including: the AOFAS for hind foot; the FAAM; the CAIT and the CAIS. RESULTS: Out of 41 patients; 27 patients were very satisfied, 11 satisfied and 3 were not satisfied. Ankle mobility returned to normal in 93% of patients. At follow-up the AOFAS was 89/100 (37-100), the FAAM 85/100% (35-100%), the CAIT 20/30 (5-30), and the CAIS 74/100% (27-100%). CONCLUSION: Outcome of modified Broström-Gould procedure is good with high satisfaction rate in terms of ankle mobility. The disparity in outcome of scores, signals towards the need of a standard evaluation system.

Seismic attenuation and velocity dispersion in heterogeneous partially saturated porous rocks

Using a numerical approach, we explore wave-induced fluid flow effects in partially saturated porous rocks in which the gas-water saturation patterns are governed by mesoscopic heterogeneities associated with the dry frame properties. The link between the dry frame properties and the gas saturation is defined by the assumption of capillary pressure equilibrium, which in the presence of heterogeneity implies that neighbouring regions can exhibit different levels of saturation. To determine the equivalent attenuation and phase velocity of the synthetic rock samples considered in this study, we apply a numerical upscaling procedure, which permits to take into account mesoscopic heterogeneities associated with the dry frame properties as well as spatially continuous variations of the pore fluid properties. The multiscale nature of the fluid saturation is taken into account by locally computing the physical properties of an effective fluid, which are then used for the larger-scale simulations. We consider two sets of numerical experiments to analyse such effects in heterogeneous partially saturated porous media, where the saturation field is determined by variations in porosity and clay content, respectively. In both cases we also evaluate the seismic responses of corresponding binary, patchy-type saturation patterns. Our results indicate that significant attenuation and modest velocity dispersion effects take place in this kind of media for both binary patchy-type and spatially continuous gas saturation patterns and in particular in the presence of relatively small amounts of gas. The numerical experiments also show that the nature of the gas distribution patterns is a critical parameter controlling the seismic responses of these environments, since attenuation and velocity dispersion effects are much more significant and occur over a broader saturation range for binary patchy-type gas-water distributions. This analysis therefore suggests that the physical mechanisms governing partial saturation should be accounted for when analysing seismic data in a poroelastic framework. In this context, heterogeneities associated with the dry frame properties, which do not play important roles in wave-induced fluid flow processes per se, should be taken into account since they may determine the kind of gas distribution pattern taking place in the porous rock.

Respiratory variability in mechanically ventilated patients. Critical Care 2011

Introduction: Increased respiratory pattern variability is associated with improved oxygenation. Pressure support (PS) is a widely used partial-assist mechanical ventilation (MV) mode, in which each breathing cycle is initiated by flow or pressure variation at the airway due to patient inspiratory effort. Neurally adjusted ventilatory assist (NAVA) is relatively new and uses the electrical activity of the diaphragm (Eadi) to deliver ventilatory support proportional to the patient's inspiratory demand. We hypothesize that respiratory variability should be greater with NAVA compared with PS.Methods: Twenty-two patients underwent 20 minutes of PS followed by 20 minutes of NAVA. Flow and Eadi curves were used to obtain tidal volume (Vt) and ∫Eadi for 300 to 400 breaths in each patient. Patient-specific cumulative distribution functions (CDF) show the percentage Vt and ∫Eadi within a clinically defined (±10%) variability band for each patient. Values are normalized to patient-specific medians for direct comparison. Variability in Vt (outcome) is thus expressed in terms of variability in ∫Eadi (demand) on the same plot.Results: Variability in Vt relative to variability in ∫Eadi is significantly greater for NAVA than PS (P = 0.00012). Hence, greater variability in outcome Vt is obtained for a given demand in ∫Eadi, under NAVA, as illustrated in Figure 1 for a typical patient. A Fisher 2 × 2 contingency analysis showed that 45% of patients under NAVA had a Vt variability in equal proportion to ∫Eadi variability, versus 0% for PS (P < 0.05).Conclusions: NAVA yields greater variability in tidal volume, relative to ∫Eadi demand, and a better match between Vt and ∫Eadi. These results indicate that NAVA could achieve improved oxygenation compared with PS when sufficient underlying variability in ∫Eadi is present, due to its ability to achieve higher tidal volume variability from a given variability in ∫Eadi.

Probabilistic evidential assessment of gunshot residue particle evidence (Part II) : Bayesian parameter estimation for experimental count data.

Part I of this series of articles focused on the construction of graphical probabilistic inference procedures, at various levels of detail, for assessing the evidential value of gunshot residue (GSR) particle evidence. The proposed models - in the form of Bayesian networks - address the issues of background presence of GSR particles, analytical performance (i.e., the efficiency of evidence searching and analysis procedures) and contamination. The use and practical implementation of Bayesian networks for case pre-assessment is also discussed. This paper, Part II, concentrates on Bayesian parameter estimation. This topic complements Part I in that it offers means for producing estimates useable for the numerical specification of the proposed probabilistic graphical models. Bayesian estimation procedures are given a primary focus of attention because they allow the scientist to combine (his/her) prior knowledge about the problem of interest with newly acquired experimental data. The present paper also considers further topics such as the sensitivity of the likelihood ratio due to uncertainty in parameters and the study of likelihood ratio values obtained for members of particular populations (e.g., individuals with or without exposure to GSR).

A repeated GGA motif is critical for the activity and stability of the riboregulator RsmY of Pseudomonas fluorescens.

The riboregulator RsmY of Pseudomonas fluorescens strain CHA0 is an example of small regulatory RNAs belonging to the global Rsm/Csr regulatory systems controlling diverse cellular processes such as glycogen accumulation, motility, or formation of extracellular products in various bacteria. By binding multiple molecules of the small regulatory protein RsmA, RsmY relieves the negative effect of RsmA on the translation of several target genes involved in the biocontrol properties of strain CHA0. RsmY and functionally related riboregulators have repeated GGA motifs predicted to be exposed in single-stranded regions, notably in the loops of hairpins. The secondary structure of RsmY was corroborated by in vivo cleavage with lead acetate. RsmY mutants lacking three or five (out of six) of the GGA motifs showed reduced ability to derepress the expression of target genes in vivo and failed to bind the RsmA protein efficiently in vitro. The absence of GGA motifs in RsmY mutants resulted in reduced abundance of these transcripts and in a shorter half-life (< or = 6 min as compared with 27 min for wild type RsmY). These results suggest that both the interaction of RsmY with RsmA and the stability of RsmY strongly depend on the GGA repeats and that the ability of RsmY to interact with small regulatory proteins such as RsmA may protect this RNA from degradation.

S-system parameter estimation for noisy metabolic profiles using newton-flow analysis.

Biochemical systems are commonly modelled by systems of ordinary differential equations (ODEs). A particular class of such models called S-systems have recently gained popularity in biochemical system modelling. The parameters of an S-system are usually estimated from time-course profiles. However, finding these estimates is a difficult computational problem. Moreover, although several methods have been recently proposed to solve this problem for ideal profiles, relatively little progress has been reported for noisy profiles. We describe a special feature of a Newton-flow optimisation problem associated with S-system parameter estimation. This enables us to significantly reduce the search space, and also lends itself to parameter estimation for noisy data. We illustrate the applicability of our method by applying it to noisy time-course data synthetically produced from previously published 4- and 30-dimensional S-systems. In addition, we propose an extension of our method that allows the detection of network topologies for small S-systems. We introduce a new method for estimating S-system parameters from time-course profiles. We show that the performance of this method compares favorably with competing methods for ideal profiles, and that it also allows the determination of parameters for noisy profiles.

Vasopressin-dependent coupling between sodium transport and water flow in a mouse cortical collecting duct cell line.

Water balance is achieved through the ability of the kidney to control water reabsorption in the connecting tubule and the collecting duct. In a mouse cortical collecting duct cell line (mCCD(c11)), physiological concentrations of arginine vasopressin increased both electrogenic, amiloride-sensitive, epithelial sodium channel (ENaC)-mediated sodium transport measured by the short-circuit current (Isc) method and water flow (Jv apical to basal) measured by gravimetry with similar activation coefficient K(1/2) (6 and 12 pM, respectively). Jv increased linearly according to the osmotic gradient across the monolayer. A small but highly significant Jv was also measured under isoosmotic conditions. To test the coupling between sodium reabsorption and water flow, mCCD(c11) cells were treated for 24 h under isoosmotic condition with either diluent, amiloride, vasopressin or vasopressin and amiloride. Isc, Jv, and net chemical sodium fluxes were measured across the same monolayers. Around 30% of baseline and 50% of vasopressin-induced water flow is coupled to an amiloride-sensitive, ENaC-mediated, electrogenic sodium transport, whereas the remaining flow is coupled to an amiloride-insensitive, nonelectrogenic sodium transport mediated by an unknown electroneutral transporter. The mCCD(c11) cell line is a first example of a mammalian tight epithelium allowing quantitative study of the coupling between sodium and water transport. Our data are consistent with the 'near isoosmotic' fluid transport model.

Adipose tissue integrity as a prerequisite for systemic energy balance: a critical role for peroxisome proliferator-activated receptor gamma.

Peroxisome proliferator-activated receptor gamma (PPARgamma) is an essential regulator of adipocyte differentiation, maintenance, and survival. Deregulations of its functions are associated with metabolic diseases. We show here that deletion of one PPARgamma allele not only affected lipid storage but, more surprisingly, also the expression of genes involved in glucose uptake and utilization, the pentose phosphate pathway, fatty acid synthesis, lipolysis, and glycerol export as well as in IR/IGF-1 signaling. These deregulations led to reduced circulating adiponectin levels and an energy crisis in the WAT, reflected in a decrease to nearly half of its intracellular ATP content. In addition, there was a decrease in the metabolic rate and physical activity of the PPARgamma(+/-) mice, which was abolished by thiazolidinedione treatment, thereby linking regulation of the metabolic rate and physical activity to PPARgamma. It is likely that the PPARgamma(+/-) phenotype was due to the observed WAT dysfunction, since the gene expression profiles associated with metabolic pathways were not affected either in the liver or the skeletal muscle. These findings highlight novel roles of PPARgamma in the adipose tissue and underscore the multifaceted action of this receptor in the functional fine tuning of a tissue that is crucial for maintaining the organism in good health.

APRIL is critical for plasmablast survival in the bone marrow and poorly expressed by early-life bone marrow stromal cells.

The persistence of serum IgG antibodies elicited in human infants is much shorter than when such responses are elicited later in life. The reasons for this rapid waning of antigen-specific antibodies elicited in infancy are yet unknown. We have recently shown that adoptively transferred tetanus toxoid (TT)-specific plasmablasts (PBs) efficiently reach the bone marrow (BM) of infant mice. However, TT-specific PBs fail to persist in the early-life BM, suggesting that they fail to receive the molecular signals that support their survival/differentiation. Using a proliferation-inducing ligand (APRIL)- and B-cell activating factor (BAFF) B-lymphocyte stimulator (BLyS)-deficient mice, we demonstrate here that APRIL is a critical factor for the establishment of the adult BM reservoir of anti-TT IgG-secreting cells. Through in vitro analyses of PB/plasma cell (PC) survival/differentiation, we show that APRIL induces the expression of Bcl-X(L) by a preferential binding to heparan sulfate proteoglycans at the surface of CD138(+) cells. Last, we identify BM-resident macrophages as the main cells that provide survival signals to PBs and show that this function is slowly acquired in early life, in parallel to a progressive acquisition of APRIL expression. Altogether, this identifies APRIL as a critical signal for PB survival that is poorly expressed in the early-life BM compartment.