Dynamic simulation of regulatory networks using SQUAD.

BACKGROUND: The ambition of most molecular biologists is the understanding of the intricate network of molecular interactions that control biological systems. As scientists uncover the components and the connectivity of these networks, it becomes possible to study their dynamical behavior as a whole and discover what is the specific role of each of their components. Since the behavior of a network is by no means intuitive, it becomes necessary to use computational models to understand its behavior and to be able to make predictions about it. Unfortunately, most current computational models describe small networks due to the scarcity of kinetic data available. To overcome this problem, we previously published a methodology to convert a signaling network into a dynamical system, even in the total absence of kinetic information. In this paper we present a software implementation of such methodology. RESULTS: We developed SQUAD, a software for the dynamic simulation of signaling networks using the standardized qualitative dynamical systems approach. SQUAD converts the network into a discrete dynamical system, and it uses a binary decision diagram algorithm to identify all the steady states of the system. Then, the software creates a continuous dynamical system and localizes its steady states which are located near the steady states of the discrete system. The software permits to make simulations on the continuous system, allowing for the modification of several parameters. Importantly, SQUAD includes a framework for perturbing networks in a manner similar to what is performed in experimental laboratory protocols, for example by activating receptors or knocking out molecular components. Using this software we have been able to successfully reproduce the behavior of the regulatory network implicated in T-helper cell differentiation. CONCLUSION: The simulation of regulatory networks aims at predicting the behavior of a whole system when subject to stimuli, such as drugs, or determine the role of specific components within the network. The predictions can then be used to interpret and/or drive laboratory experiments. SQUAD provides a user-friendly graphical interface, accessible to both computational and experimental biologists for the fast qualitative simulation of large regulatory networks for which kinetic data is not necessarily available.

Quantitative integration of geophysical and hydraulic data : from the local towards the regional scale range

Des progrès significatifs ont été réalisés dans le domaine de l'intégration quantitative des données géophysique et hydrologique l'échelle locale. Cependant, l'extension à de plus grandes échelles des approches correspondantes constitue encore un défi majeur. Il est néanmoins extrêmement important de relever ce défi pour développer des modèles fiables de flux des eaux souterraines et de transport de contaminant. Pour résoudre ce problème, j'ai développé une technique d'intégration des données hydrogéophysiques basée sur une procédure bayésienne de simulation séquentielle en deux étapes. Cette procédure vise des problèmes à plus grande échelle. L'objectif est de simuler la distribution d'un paramètre hydraulique cible à partir, d'une part, de mesures d'un paramètre géophysique pertinent qui couvrent l'espace de manière exhaustive, mais avec une faible résolution (spatiale) et, d'autre part, de mesures locales de très haute résolution des mêmes paramètres géophysique et hydraulique. Pour cela, mon algorithme lie dans un premier temps les données géophysiques de faible et de haute résolution à travers une procédure de réduction déchelle. Les données géophysiques régionales réduites sont ensuite reliées au champ du paramètre hydraulique à haute résolution. J'illustre d'abord l'application de cette nouvelle approche dintégration des données à une base de données synthétiques réaliste. Celle-ci est constituée de mesures de conductivité hydraulique et électrique de haute résolution réalisées dans les mêmes forages ainsi que destimations des conductivités électriques obtenues à partir de mesures de tomographic de résistivité électrique (ERT) sur l'ensemble de l'espace. Ces dernières mesures ont une faible résolution spatiale. La viabilité globale de cette méthode est testée en effectuant les simulations de flux et de transport au travers du modèle original du champ de conductivité hydraulique ainsi que du modèle simulé. Les simulations sont alors comparées. Les résultats obtenus indiquent que la procédure dintégration des données proposée permet d'obtenir des estimations de la conductivité en adéquation avec la structure à grande échelle ainsi que des predictions fiables des caractéristiques de transports sur des distances de moyenne à grande échelle. Les résultats correspondant au scénario de terrain indiquent que l'approche d'intégration des données nouvellement mise au point est capable d'appréhender correctement les hétérogénéitées à petite échelle aussi bien que les tendances à gande échelle du champ hydraulique prévalent. Les résultats montrent également une flexibilté remarquable et une robustesse de cette nouvelle approche dintégration des données. De ce fait, elle est susceptible d'être appliquée à un large éventail de données géophysiques et hydrologiques, à toutes les gammes déchelles. Dans la deuxième partie de ma thèse, j'évalue en détail la viabilité du réechantillonnage geostatique séquentiel comme mécanisme de proposition pour les méthodes Markov Chain Monte Carlo (MCMC) appliquées à des probmes inverses géophysiques et hydrologiques de grande dimension . L'objectif est de permettre une quantification plus précise et plus réaliste des incertitudes associées aux modèles obtenus. En considérant une série dexemples de tomographic radar puits à puits, j'étudie deux classes de stratégies de rééchantillonnage spatial en considérant leur habilité à générer efficacement et précisément des réalisations de la distribution postérieure bayésienne. Les résultats obtenus montrent que, malgré sa popularité, le réechantillonnage séquentiel est plutôt inefficace à générer des échantillons postérieurs indépendants pour des études de cas synthétiques réalistes, notamment pour le cas assez communs et importants où il existe de fortes corrélations spatiales entre le modèle et les paramètres. Pour résoudre ce problème, j'ai développé un nouvelle approche de perturbation basée sur une déformation progressive. Cette approche est flexible en ce qui concerne le nombre de paramètres du modèle et lintensité de la perturbation. Par rapport au rééchantillonage séquentiel, cette nouvelle approche s'avère être très efficace pour diminuer le nombre requis d'itérations pour générer des échantillons indépendants à partir de la distribution postérieure bayésienne. - Significant progress has been made with regard to the quantitative integration of geophysical and hydrological data at the local scale. However, extending corresponding approaches beyond the local scale still represents a major challenge, yet is critically important for the development of reliable groundwater flow and contaminant transport models. To address this issue, I have developed a hydrogeophysical data integration technique based on a two-step Bayesian sequential simulation procedure that is specifically targeted towards larger-scale problems. The objective is to simulate the distribution of a target hydraulic parameter based on spatially exhaustive, but poorly resolved, measurements of a pertinent geophysical parameter and locally highly resolved, but spatially sparse, measurements of the considered geophysical and hydraulic parameters. To this end, my algorithm links the low- and high-resolution geophysical data via a downscaling procedure before relating the downscaled regional-scale geophysical data to the high-resolution hydraulic parameter field. I first illustrate the application of this novel data integration approach to a realistic synthetic database consisting of collocated high-resolution borehole measurements of the hydraulic and electrical conductivities and spatially exhaustive, low-resolution electrical conductivity estimates obtained from electrical resistivity tomography (ERT). The overall viability of this method is tested and verified by performing and comparing flow and transport simulations through the original and simulated hydraulic conductivity fields. The corresponding results indicate that the proposed data integration procedure does indeed allow for obtaining faithful estimates of the larger-scale hydraulic conductivity structure and reliable predictions of the transport characteristics over medium- to regional-scale distances. The approach is then applied to a corresponding field scenario consisting of collocated high- resolution measurements of the electrical conductivity, as measured using a cone penetrometer testing (CPT) system, and the hydraulic conductivity, as estimated from electromagnetic flowmeter and slug test measurements, in combination with spatially exhaustive low-resolution electrical conductivity estimates obtained from surface-based electrical resistivity tomography (ERT). The corresponding results indicate that the newly developed data integration approach is indeed capable of adequately capturing both the small-scale heterogeneity as well as the larger-scale trend of the prevailing hydraulic conductivity field. The results also indicate that this novel data integration approach is remarkably flexible and robust and hence can be expected to be applicable to a wide range of geophysical and hydrological data at all scale ranges. In the second part of my thesis, I evaluate in detail the viability of sequential geostatistical resampling as a proposal mechanism for Markov Chain Monte Carlo (MCMC) methods applied to high-dimensional geophysical and hydrological inverse problems in order to allow for a more accurate and realistic quantification of the uncertainty associated with the thus inferred models. Focusing on a series of pertinent crosshole georadar tomographic examples, I investigated two classes of geostatistical resampling strategies with regard to their ability to efficiently and accurately generate independent realizations from the Bayesian posterior distribution. The corresponding results indicate that, despite its popularity, sequential resampling is rather inefficient at drawing independent posterior samples for realistic synthetic case studies, notably for the practically common and important scenario of pronounced spatial correlation between model parameters. To address this issue, I have developed a new gradual-deformation-based perturbation approach, which is flexible with regard to the number of model parameters as well as the perturbation strength. Compared to sequential resampling, this newly proposed approach was proven to be highly effective in decreasing the number of iterations required for drawing independent samples from the Bayesian posterior distribution.

Protease-activated receptor 2 signalling promotes dendritic cell antigen transport and T-cell activation in vivo.

Deficiency of protease-activated receptor-2 (PAR2) modulates inflammation in several models of inflammatory and autoimmune disease, although the underlying mechanism(s) are not understood. PAR2 is expressed on endothelial and immune cells, and is implicated in dendritic cell (DC) differentiation. We investigated in vivo the impact of PAR2 activation on DCs and T cells in PAR2 wild-type (WT) and knockout (KO) mice using a specific PAR2 agonist peptide (AP2). PAR2 activation significantly increased the frequency of mature CD11c(high) DCs in draining lymph nodes 24 hr after AP2 administration. Furthermore, these DCs exhibited increased expression of major histocompatibility complex (MHC) class II and CD86. A significant increase in activated (CD44(+) CD62(-)) CD4(+) and CD8(+) T-cell frequencies was also observed in draining lymph nodes 48 hr after AP2 injection. No detectable change in DC or T-cell activation profiles was observed in the spleen. The influence of PAR2 signalling on antigen transport to draining lymph nodes was assessed in the context of delayed-type hypersensitivity. PAR2 WT mice that were sensitized by skin-painting with fluorescein isothiocyanate (FITC) to induce delayed-type hypersensitivity possessed elevated proportion of FITC(+) DCs in draining lymph nodes 24 hr after FITC painting when compared with PAR2 KO mice (0.95% versus 0.47% of total lymph node cells). Collectively, these results demonstrate that PAR2 signalling promotes DC trafficking to the lymph nodes and subsequent T-cell activation, and thus provides an explanation for the pro-inflammatory effect of PAR2 in animal models of inflammation.

Transport and metabolism of [214-C]abscisic acid in maize root

The tips of intact maize (cv. LG 11) roots, maintained vertically, were pretreated with a droplet of buffer solution or a bead of anion exchange resin, both containing [214-C]abscisic acid (ABA). A significant basipetal ABA movement was observed and two metabolites of ABA (possibly phaseic acid and dihydrophaseic acid) were found. ABA pretreatment enhanced the gravireaction of 10 mm apical root segments kept both in the dark and in the light. The possibility that ABA could be one of the endogenous growth inhibitors produced or released by the cap cells is discussed.

Aldosterone paradox: differential regulation of ion transport in distal nephron.

The mechanisms through which aldosterone promotes apparently opposite effects like salt reabsorption and K(+) secretion remain poorly understood. The identification, localization, and physiological analysis of ion transport systems in distal nephron have revealed an intricate network of interactions between several players, revealing the complex mechanism behind the aldosterone paradox. We review the mechanisms involved in differential regulation of ion transport that allow the fine tuning of salt and K(+) balance.

Purine transport and the cell cycle.

Background: Alliance evolutions, i.e. ruptures and resolutions over the course of psychotherapy, have been shown to be important descriptive features in different forms of psychotherapy, and in particular in psychodynamic psychotherapy. This case study of a client presenting elements of adjustment disorder undergoing short-term dynamic psychotherapy is drawn from a systematic naturalistic study and aims at illustrating, on a session-by-session-level, the processes of alliance ruptures and resolutions, by comparing both the client's and the therapist's perspectives. Method: Two episodes of alliance evolution were more fully studied, in relation to the evolution of transference, as well as the client's defensive functioning and core conflictual theme. These concepts were measured by means of valid, reliable observer-rater methods, based on session transcripts: the Defense Mechanisms Rating Scales (DMRS) for defensive functioning and the Core Conflictual Relationship Theme (CCRT) for the conflicts. Alliance was measured after each session using the Helping Alliance questionnaire (HAq-II). Results: The results indicated that these episodes of alliance rupture and resolutions may be understood as key moments of the whole therapeutic process reflecting the client's main relationship stakes. Illustrations are provided based on the client's in-session processes and related to the alliance development over the course of the entire therapy.

Biometal muscle to restore atrial transport function in a permanent atrial fibrillation animal model: a potential tool in the treatment of end-stage heart failure.

BACKGROUND: Half of the patients with end-stage heart failure suffer from persistent atrial fibrillation (AF). Atrial kick (AK) accounts for 10-15% of the ejection fraction. A device restoring AK should significantly improve cardiac output (CO) and possibly delay ventricular assist device (VAD) implantation. This study has been designed to assess the mechanical effects of a motorless pump on the right chambers of the heart in an animal model. METHODS: Atripump is a dome-shaped biometal actuator electrically driven by a pacemaker-like control unit. In eight sheep, the device was sutured onto the right atrium (RA). AF was simulated with rapid atrial pacing. RA ejection fraction (EF) was assessed with intracardiac ultrasound (ICUS) in baseline, AF and assisted-AF status. In two animals, the pump was left in place for 4 weeks and then explanted. Histology examination was carried out. The mean values for single measurement per animal with +/-SD were analysed. RESULTS: The contraction rate of the device was 60 per min. RA EF was 41% in baseline, 7% in AF and 21% in assisted-AF conditions. CO was 7+/-0.5 l min(-1) in baseline, 6.2+/-0.5 l min(-1) in AF and 6.7+/-0.5 l min(-1) in assisted-AF status (p<0.01). Histology of the atrium in the chronic group showed chronic tissue inflammation and no sign of tissue necrosis. CONCLUSIONS: The artificial muscle restores the AK and improves CO. In patients with end-stage cardiac failure and permanent AF, if implanted on both sides, it would improve CO and possibly delay or even avoid complex surgical treatment such as VAD implantation.

Numerical simulation of left ventricular assist device implantations: comparing the ascending and the descending aorta cannulations.

In this work we present numerical simulations of continuous flow left ventricle assist device implantation with the aim of comparing difference in flow rates and pressure patterns depending on the location of the anastomosis and the rotational speed of the device. Despite the fact that the descending aorta anastomosis approach is less invasive, since it does not require a sternotomy and a cardiopulmonary bypass, its benefits are still controversial. Moreover, the device rotational speed should be correctly chosen to avoid anomalous flow rates and pressure distribution in specific location of the cardiovascular tree. With the aim of assessing the differences between these two approaches and device rotational speed in terms of flow rate and pressure waveforms, we set up numerical simulations of network of one-dimensional models where we account for the presence of an outflow cannula anastomosed to different locations of the aorta. Then, we use the resulting network to compare the results of the two different cannulations for several stages of heart failure and different rotational speed of the device. The inflow boundary data for the heart and the cannulas are obtained from a lumped parameters model of the entire circulatory system with an assist device, which is validated with clinical data. The results show that ascending and descending aorta cannulations lead to similar waveforms and mean flow rate in all the considered cases. Moreover, regardless of the anastomosis region, the rotational speed of the device has an important impact on wave profiles; this effect is more pronounced at high RPM.

Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study.

The identification of genetically homogeneous groups of individuals is a long standing issue in population genetics. A recent Bayesian algorithm implemented in the software STRUCTURE allows the identification of such groups. However, the ability of this algorithm to detect the true number of clusters (K) in a sample of individuals when patterns of dispersal among populations are not homogeneous has not been tested. The goal of this study is to carry out such tests, using various dispersal scenarios from data generated with an individual-based model. We found that in most cases the estimated 'log probability of data' does not provide a correct estimation of the number of clusters, K. However, using an ad hoc statistic DeltaK based on the rate of change in the log probability of data between successive K values, we found that STRUCTURE accurately detects the uppermost hierarchical level of structure for the scenarios we tested. As might be expected, the results are sensitive to the type of genetic marker used (AFLP vs. microsatellite), the number of loci scored, the number of populations sampled, and the number of individuals typed in each sample.

Simulation of amoxicillin pharmacokinetics in humans for the prevention of streptococcal endocarditis in rats.

The pharmacokinetic determinants of successful antibiotic prophylaxis of endocarditis are not precisely known. Differences in half-lives of antibiotics between animals and humans preclude extrapolation of animal results to human situations. To overcome this limitation, we have mimicked in rats the amoxicillin kinetics in humans following a 3-g oral dose (as often used for prophylaxis of endocarditis) by delivering the drug through a computerized pump. Rats with catheter-induced vegetations were challenged with either of two strains of antibiotic-tolerant viridans group streptococci. Antibiotics were given either through the pump (to simulate the whole kinetic profile during prophylaxis in humans) or as an intravenous bolus which imitated only the peak level of amoxicillin (18 mg/liter) in human serum. Prophylaxis by intravenous bolus was inoculum dependent and afforded a limited protection only in rats challenged with the minimum inoculum size infecting > or = 90% of untreated controls. In contrast, simulation of kinetics in humans significantly protected animals challenged with 10 to 100 times the inoculum of either of the test organisms infecting > or = 90% of untreated controls. Thus, simulation of the profiles of amoxicillin prophylaxis in human serum was more efficacious than mere imitation of the transient peak level in rats. This confirms previous studies suggesting that the duration for which the serum amoxicillin level remained detectable (not only the magnitude of the peak) was an important parameter in successful prophylaxis of endocarditis. The results also suggest that single-dose prophylaxis with 3 g of amoxicillin in humans might be more effective than predicted by conventional animal models in which only peak levels of antibiotic in human serum were stimulated.

NF-kappaB inhibits sodium transport via down-regulation of SGK1 in renal collecting duct principal cells.

Tubulointerstitial inflammation is a common feature of renal diseases. We have investigated the relationship between inflammation and Na(+) transport in the collecting duct (CD) using the mCCD(cl1) and mpkCDD(cl4) principal cell models. Lipopolysaccharide (LPS) decreased basal and aldosterone-stimulated amiloride-sensitive transepithelial current in a time-dependent manner. This effect was associated with a decrease in serum and glucocorticoid-regulated kinase 1 (SGK1) mRNA and protein levels followed by a decrease in epithelial sodium channel (ENaC) alpha-subunit mRNA levels. The LPS-induced decrease in SGK1 expression was confirmed in isolated rat CD. This decreased expression of either SGK1 or the ENaC alpha-subunit was not due to enhanced degradation of mRNA. In contrast, LPS inhibited transcriptional activity of the SGK1 promoter measured by luciferase-reporter gene assay. The effect of LPS was not mediated by inhibition of mineralocorticoid or glucocorticoid receptor, because expression of both receptors was unchanged and blockade of either receptor by spironolactone or RU486, respectively, did not prevent the down-regulation of SGK1. The effect of LPS was mediated by the canonical NF-kappaB pathway, as overexpression of a constitutively active mutant, IKKbeta (inhibitor of nuclear factor kappaB kinase-beta) decreased SGK1 mRNA levels, and knockdown of p65 NF-kappaB subunit by small interfering RNA increased SGK1 mRNA levels. Chromatin immunoprecipitation showed that LPS increased p65 binding to two NF-kappaB sites along the SGK1 promoter. In conclusion, we show that activation of the NF-kappaB pathway down-regulates SGK1 expression, which might lead to decreased ENaC alpha-subunit expression, ultimately resulting in decreased Na(+) transport.

Imaging and quantifying salt-tracer transport in a riparian groundwater system by means of 3D ERT monitoring

Determining groundwater flow paths of infiltrated river water is necessary for studying biochemical processes in the riparian zone, but their characterization is complicated by strong temporal and spatial heterogeneity. We investigated to what extent repeat 3D surface electrical resistance tomography (ERT) can be used to monitor transport of a salt-tracer plume under close to natural gradient conditions. The aim is to estimate groundwater flow velocities and pathways at a site located within a riparian groundwater system adjacent to the perialpine Thur River in northeastern Switzerland. Our ERT time-lapse images provide constraints on the plume's shape, flow direction, and velocity. These images allow the movement of the plume to be followed for 35 m. Although the hydraulic gradient is only 1.43 parts per thousand, the ERT time-lapse images demonstrate that the plume's center of mass and its front propagate with velocities of 2x10(-4) m/s and 5x10(-4) m/s, respectively. These velocities are compatible with groundwater resistivity monitoring data in two observation wells 5 m from the injection well. Five additional sensors in the 5-30 m distance range did not detect the plume. Comparison of the ERT time-lapse images with a groundwater transport model and time-lapse inversions of synthetic ERT data indicate that the movement of the plume can be described for the first 6 h after injection by a uniform transport model. Subsurface heterogeneity causes a change of the plume's direction and velocity at later times. Our results demonstrate the effectiveness of using time-lapse 3D surface ERT to monitor flow pathways in a challenging perialpine environment over larger scales than is practically possible with crosshole 3D ERT.