Estimating parameters for generalized mass action models using constraint propagation.

As modern molecular biology moves towards the analysis of biological systems as opposed to their individual components, the need for appropriate mathematical and computational techniques for understanding the dynamics and structure of such systems is becoming more pressing. For example, the modeling of biochemical systems using ordinary differential equations (ODEs) based on high-throughput, time-dense profiles is becoming more common-place, which is necessitating the development of improved techniques to estimate model parameters from such data. Due to the high dimensionality of this estimation problem, straight-forward optimization strategies rarely produce correct parameter values, and hence current methods tend to utilize genetic/evolutionary algorithms to perform non-linear parameter fitting. Here, we describe a completely deterministic approach, which is based on interval analysis. This allows us to examine entire sets of parameters, and thus to exhaust the global search within a finite number of steps. In particular, we show how our method may be applied to a generic class of ODEs used for modeling biochemical systems called Generalized Mass Action Models (GMAs). In addition, we show that for GMAs our method is amenable to the technique in interval arithmetic called constraint propagation, which allows great improvement of its efficiency. To illustrate the applicability of our method we apply it to some networks of biochemical reactions appearing in the literature, showing in particular that, in addition to estimating system parameters in the absence of noise, our method may also be used to recover the topology of these networks.

Fish oil supplementation does not alter energy efficiency in healthy males

BACKGROUND AND AIMS: Fish oil (FO) supplementation prevents the development of obesity and insulin resistance, and upregulate the expression of UCP3 in skeletal muscle in rodents. This may represent indirect evidence that FO promotes fat oxidation and/or alter energy efficiency. The aim of this study was to evaluate whether such effects can be observed in humans. The metabolic effects of FO were assessed during exercise in order to obtain a direct measurement of energy efficiency. METHODS: Eight healthy male volunteers were studied with and without supplementation with 7.2 g/day FO (including 1.1 g/day eicosopentaenoic acid and 0.7 g/day decosahexaenoic acid) during 14 days. Their VO(2 max) was measured on cycle ergometer. Thereafter, energy metabolism (substrate oxidation, energy expenditure and energy efficiency) was assessed during a 30 min cycling exercise at 50% VO(2 max) performed 2 h 30 after a standardized, high carbohydrate breakfast. RESULTS: VO(2 max) was 38.6+/-2.2 after FO and 38.4+/-2.0 (mL x kg(-1) x min(-1)) in control conditions (NS). Basal plasma glucose, insulin and NEFA concentrations, and energy metabolism were similar with FO and in controls. During exercise, the increases in plasma NEFA concentrations, energy expenditure, glucose and lipid oxidation, and the decreases in glycaemia and insulinemia were not altered by FO intake. Energy efficiency was 22.4+/-0.6% after FO vs 21.8+/-0.7% in controls. In order to ascertain that the absence of effects of FO was not due to consumption of a carbohydrate meal immediately before exercise, 4 of the 8 subjects were re-studied in fasting conditions, FO also failed to alter energy efficiency in this subset of studies. CONCLUSION: FO supplementation did not significantly alter energy metabolism and energy efficiency during exercise in healthy humans.

Improved SNR efficiency in gradient echo coronary MRA with high temporal resolution using parallel imaging.

Contemporary coronary magnetic resonance angiography techniques suffer from signal-to-noise ratio (SNR) constraints. We propose a method to enhance SNR in gradient echo coronary magnetic resonance angiography by using sensitivity encoding (SENSE). While the use of sensitivity encoding to improve SNR seems counterintuitive, it can be exploited by reducing the number of radiofrequency excitations during the acquisition window while lowering the signal readout bandwidth, therefore improving the radiofrequency receive to radiofrequency transmit duty cycle. Under certain conditions, this leads to improved SNR. The use of sensitivity encoding for improved SNR in three-dimensional coronary magnetic resonance angiography is investigated using numerical simulations and an in vitro and an in vivo study. A maximum 55% SNR enhancement for coronary magnetic resonance angiography was found both in vitro and in vivo, which is well consistent with the numerical simulations. This method is most suitable for spoiled gradient echo coronary magnetic resonance angiography in which a high temporal and spatial resolution is required.

Computational analysis of the genetic and environmental contributions to disease-related human phenotypes: From predicting adverse lipid response of HIV patients receiving antiretroviral therapy to genome-wide association studies of cardiovascular and lipid related disorders

Genetic variants influence the risk to develop certain diseases or give rise to differences in drug response. Recent progresses in cost-effective, high-throughput genome-wide techniques, such as microarrays measuring Single Nucleotide Polymorphisms (SNPs), have facilitated genotyping of large clinical and population cohorts. Combining the massive genotypic data with measurements of phenotypic traits allows for the determination of genetic differences that explain, at least in part, the phenotypic variations within a population. So far, models combining the most significant variants can only explain a small fraction of the variance, indicating the limitations of current models. In particular, researchers have only begun to address the possibility of interactions between genotypes and the environment. Elucidating the contributions of such interactions is a difficult task because of the large number of genetic as well as possible environmental factors.In this thesis, I worked on several projects within this context. My first and main project was the identification of possible SNP-environment interactions, where the phenotypes were serum lipid levels of patients from the Swiss HIV Cohort Study (SHCS) treated with antiretroviral therapy. Here the genotypes consisted of a limited set of SNPs in candidate genes relevant for lipid transport and metabolism. The environmental variables were the specific combinations of drugs given to each patient over the treatment period. My work explored bioinformatic and statistical approaches to relate patients' lipid responses to these SNPs, drugs and, importantly, their interactions. The goal of this project was to improve our understanding and to explore the possibility of predicting dyslipidemia, a well-known adverse drug reaction of antiretroviral therapy. Specifically, I quantified how much of the variance in lipid profiles could be explained by the host genetic variants, the administered drugs and SNP-drug interactions and assessed the predictive power of these features on lipid responses. Using cross-validation stratified by patients, we could not validate our hypothesis that models that select a subset of SNP-drug interactions in a principled way have better predictive power than the control models using "random" subsets. Nevertheless, all models tested containing SNP and/or drug terms, exhibited significant predictive power (as compared to a random predictor) and explained a sizable proportion of variance, in the patient stratified cross-validation context. Importantly, the model containing stepwise selected SNP terms showed higher capacity to predict triglyceride levels than a model containing randomly selected SNPs. Dyslipidemia is a complex trait for which many factors remain to be discovered, thus missing from the data, and possibly explaining the limitations of our analysis. In particular, the interactions of drugs with SNPs selected from the set of candidate genes likely have small effect sizes which we were unable to detect in a sample of the present size (<800 patients).In the second part of my thesis, I performed genome-wide association studies within the Cohorte Lausannoise (CoLaus). I have been involved in several international projects to identify SNPs that are associated with various traits, such as serum calcium, body mass index, two-hour glucose levels, as well as metabolic syndrome and its components. These phenotypes are all related to major human health issues, such as cardiovascular disease. I applied statistical methods to detect new variants associated with these phenotypes, contributing to the identification of new genetic loci that may lead to new insights into the genetic basis of these traits. This kind of research will lead to a better understanding of the mechanisms underlying these pathologies, a better evaluation of disease risk, the identification of new therapeutic leads and may ultimately lead to the realization of "personalized" medicine.

A random linear functional approach to efficiency bounds

This paper suggests a method for obtaining efficiency bounds in models containing either only infinite-dimensional parameters or both finite- and infinite-dimensional parameters (semiparametric models). The method is based on a theory of random linear functionals applied to the gradient of the log-likelihood functional and is illustrated by computing the lower bound for Cox's regression model

Efficiency and significance of multiple vocal signals in sibling competition

Animals can compete for resources by displaying various acoustic signals that may differentially affect the outcome of competition. We propose the hypothesis that the most efficient signal to deter opponents should be the one that most honestly reveals motivation to compete. We tested this hypothesis in the barn owl (Tyto alba) in which nestlings produce more calls of longer duration than siblings to compete for priority access to the indivisible prey item their parents will deliver next. Because nestlings increase call rate to a larger extent than call duration when they become hungrier, call rate would signal more accurately hunger level. This leads us to propose three predictions. First, a high number of calls should be more efficient in deterring siblings to compete than long calls. Second, the rate at which an individual calls should be more sensitive to variation in the intensity of the sibling vocal competition than the duration of its calls. Third, call rate should influence competitors' vocalization for a longer period of time than call duration. To test these three predictions we performed playback experiments by broadcasting to singleton nestlings calls of varying durations and at different rates. According to the first prediction, singleton nestlings became less vocal to a larger extent when we broadcasted more calls compared to longer calls. In line with the second prediction, nestlings reduced vocalization rate to a larger extent than call duration when we broadcasted more or longer calls. Finally, call rate had a longer influence on opponent's vocal behavior than call duration. Young animals thus actively and differentially use multiple signaling components to compete with their siblings over parental resources.

The potential effects of social interactions on reproductive efficiency of stallions

The reproductive efficiency of stabled domestic stallions is often lower than what could be expected from observations in feral herds. In the wild, stallions typically live with mares in harem bands, with other stallions in bachelor bands, or occasionally in mixed-sex transitional bands. We, therefore, argue that permanent contact with mares may increase reproductive efficiency of stallions suffering from low libido and/or fertility. We also provide a summary of our present knowledge of natural conditions, management, and husbandry of domestic stallions, and of intra- and intersexual behavioral interactions in horses.

Genomic determinants of the efficiency of internal ribosomal entry sites of viral and cellular origin

Variation in cellular gene expression levels has been shown to be inherited. Expression is controlled at transcriptional and post-transcriptional levels. Internal ribosome entry sites (IRES) are used by viruses to bypass inhibition of cap-dependent translation, and by eukaryotic cells to control translation under conditions when protein synthesis is inhibited. We aimed at identifying genomic determinants of variability in IRES-mediated translation of viral [Encephalomyocarditis virus (EMCV)] and cellular IRES [X-linked inhibitor-of-apoptosis (XIAP) and c-myc]. Bicistronic lentiviral constructs expressing two fluorescent reporters were used to transduce laboratory and B lymphoblastoid cell lines [15 CEPH pedigrees (n = 205) and 50 unrelated individuals]. IRES efficiency varied according to cell type and among individuals. Control of IRES activity has a significant genetic component (h(2) of 0.47 and 0.36 for EMCV and XIAP, respectively). Quantitative linkage analysis identified a suggestive locus (LOD 2.35) on chromosome 18q21.2, and genome-wide association analysis revealed of a cluster of SNPs on chromosome 3, intronic to the FHIT gene, marginally associated (P = 5.9E-7) with XIAP IRES function. This study illustrates the in vitro generation of intermediate phenotypes by using cell lines for the evaluation of genetic determinants of control of elements such as IRES.

Guideline adaptation: an approach to enhance efficiency in guideline development and improve utilisation.

BACKGROUND: Developing and updating high-quality guidelines requires substantial time and resources. To reduce duplication of effort and enhance efficiency, we developed a process for guideline adaptation and assessed initial perceptions of its feasibility and usefulness. METHODS: Based on preliminary developments and empirical studies, a series of meetings with guideline experts were organised to define a process for guideline adaptation (ADAPTE) and to develop a manual and a toolkit made available on a website (http://www.adapte.org). Potential users, guideline developers and implementers, were invited to register and to complete a questionnaire evaluating their perception about the proposed process. RESULTS: The ADAPTE process consists of three phases (set-up, adaptation, finalisation), 9 modules and 24 steps. The adaptation phase involves identifying specific clinical questions, searching for, retrieving and assessing available guidelines, and preparing the draft adapted guideline. Among 330 registered individuals (46 countries), 144 completed the questionnaire. A majority found the ADAPTE process clear (78%), comprehensive (69%) and feasible (60%), and the manual useful (79%). However, 21% found the ADAPTE process complex. 44% feared that they will not find appropriate and high-quality source guidelines. DISCUSSION: A comprehensive framework for guideline adaptation has been developed to meet the challenges of timely guideline development and implementation. The ADAPTE process generated important interest among guideline developers and implementers. The majority perceived the ADAPTE process to be feasible, useful and leading to improved methodological rigour and guideline quality. However, some de novo development might be needed if no high quality guideline exists for a given topic.

Efficiency of coronary artery disease screening guidelines in asymptomatic and atypical chest pain type 2 diabetic patients

RESUME Objectifs. Évaluer la prévalence de maladie coronarienne chez les patients diabétiques de type 2 asymptomatiques ou avec angor atypique selon les recommandations américaines de l'American Diabetes Association et de l'American College of Cardiology. Méthodes. Cent cinquante-quatre patients diabétiques de type 2 asymptomatiques ou avec angor atypique et présentant au minimum 2 facteurs de risque cardio-vasculaires additionnels ont été dépistés par échocardiographie de stress (71%, n=109), scintigraphie myocardique de perfusion (26%, n=40) ou l'association des 2 examens (3%, n=5). Résultats. L'échocardiographie de stress s'est révélée positive chez 16 patients (14%) et 14 ont eu une coronarographie révélant des sténoses significatives chez 12 (86%). La scintigraphie myocardique de perfusion était positive chez 16 patients (36%). Huit patients ont eu une coronarographie et 4 (50%) présentaient des sténoses significatives. Au total, 31 patients (20%) ont montré des signes d'ischémie lors de l'examen non-invasif et 15 (10%) ont présenté des sténoses significatives à la coronarographie. Les facteurs prédictifs indépendants de la maladie coronarienne étaient le tabagisme (OR 6.5, p=0.05), la microalbuminurie (OR 3.9, p=0.03), ainsi que les souffles fémoraux (OR 17.1, p=0.008). Conclusions. En suivant les recommandations américaines, un patient sur cinq présentait une ischémie lors des examens non-invasifs, tandis que 1 sur 10 avait des sténoses significatives à la coronarographie. L'analyse multivariée suggère que des marqueurs des complications micro- et macro-vasculaires en combinaison avec des facteurs de risque cardio-vasculaire classiques pourraient améliorer le pouvoir diagnostic de ces recommandations. SUMMARY Aims. We evaluated the prevalence of coronary artery disease in asymptomatic and atypical chest pain type 2 diabetic patients according to the American Diabetes Association and American College of Cardiology guidelines. Methods. Asymptomatic or atypical chest pain type 2 diabetic patients (n=154), with at least two additional cardiovascular risk factors, were screened for coronary artery disease using stress echocardiography (71%, n=109), myocardial perfusion imaging (26%, n=40) or both (3%, n=5). Results. Stress echocardiography was positive in 16 patients (14%) and 14 had a coronary angiography, revealing significant stenoses in 12 (86%). Myocardial perfusion imaging was positive in 16 patients (36%). Eight patients underwent angiography and 4 (50%) presented significant stenoses. Overall, 31 patients (20%) demonstrated signs of ischemia on non-invasive tests and 15 (10%) presented significant stenoses on coronary angiography. Independent predictors of coronary artery disease were smoking (OR 6.5, p=0.05), microalbuminuria (OR 3.9, p=0.03) and femoral murmur (OR 17.1, p=0.008). Conclusions. Following the guidelines, one in five diabetic patient presented ischemia on noninvasive tests, while one in ten presented significant coronary stenoses. Multivariate analysis suggests that adding markers of micro- and macro-vascular complications to classical cardiovascular risk factors may enhance the diagnostic efficiency of the guidelines.

Effect of water velocity on the uptake of polychlorinatedd biphenyls (PCBs) by silicone rubber (SR) and low-density polyethylene (LDPE) passive samplers: an assessment of the efficiency of performance reference compounds (PRCs) in river-like flow conditions

One aim of this study is to determine the impact of water velocity on the uptake of indicator polychlorinated biphenyls (iPCBs) by silicone rubber (SR) and low-density polyethylene (LDPE) passive samplers. A second aim is to assess the efficiency of performance reference compounds (PRCs) to correct for the impact of water velocity. SR and LDPE samplers were spiked with 11 or 12 PRCs and exposed for 6 weeks to four different velocities (in the range of 1.6 to 37.7 cm s−1) in river-like flow conditions using a channel system supplied with river water. A relationship between velocity and the uptakewas found for each iPCB and enables to determine expected changes in the uptake due to velocity variations. For both samplers, velocity increases from 2 to 10 cm s−1, 30 cm s−1 (interpolated data) and 100 cm s−1 (extrapolated data) lead to increases of the uptake which do not exceed a factor of 2, 3 and 4.5, respectively. Results also showed that the influence of velocity decreased with increasing the octanol-water coefficient partition (log Kow) of iPCBs when SR is used whereas the opposite effect was observed for LDPE. Time-weighted average (TWA) concentrations of iPCBs in water were calculated from iPCB uptake and PRC release. These calculations were performed using either a single PRC or all the PRCs. The efficiency of PRCs to correct the impact of velocity was assessed by comparing the TWA concentrations obtained at the four tested velocities. For SR, a good agreement was found among the four TWA concentrations with both methods (average RSD b 10%). Also for LDPE, PRCs offered a good correction of the impact of water velocity (average RSD of about 10 to 20%). These results contribute to the process of acceptance of passive sampling in routine regulatory monitoring programs.

Justice and Court Administrations, Their Workings and Efficiency in Switzerland : Aspects of Sentencing And Its Outcome In Swiss Cantons

The comparison of the operations of the administration of justice among cantons shows on one side large differences in the three major types of sentencing, in the use of pre-trial detention and the unsuspended prison sanction. When combined, one finds however very weak relationships when considering absolute, percentage or weighted results. On the other side, the outcome of these different policies is much paradoxical as there are no differences when comparing recidivism rates among cantons, despite strong differences in the use of pre-trial detention and the sentencing with prison sanctions. The paradoxical outcome of crime policies in terms of recidivism - e.g. the absence of differences of the outcome based on sanctions in the domain of less severe delinquency - suggests the need for more empirically informed crime policies. The role of justice administrators could be to participate in the dissemination of those findings as well as the dissemination of best practices among cantons with regard to outcomes and the use of resources - especially with consideration to the use of the prison sanction as it is the most costly and the most inefficient of all sanctions. Furthermore, the observance of the principle of equality before the law would be most likely be promoted.

Development and application of computational methodologies in qualitative modeling

La présente thèse s'intitule "Développent et Application des Méthodologies Computationnelles pour la Modélisation Qualitative". Elle comprend tous les différents projets que j'ai entrepris en tant que doctorante. Plutôt qu'une mise en oeuvre systématique d'un cadre défini a priori, cette thèse devrait être considérée comme une exploration des méthodes qui peuvent nous aider à déduire le plan de processus regulatoires et de signalisation. Cette exploration a été mue par des questions biologiques concrètes, plutôt que par des investigations théoriques. Bien que tous les projets aient inclus des systèmes divergents (réseaux régulateurs de gènes du cycle cellulaire, réseaux de signalisation de cellules pulmonaires) ainsi que des organismes (levure à fission, levure bourgeonnante, rat, humain), nos objectifs étaient complémentaires et cohérents. Le projet principal de la thèse est la modélisation du réseau de l'initiation de septation (SIN) du S.pombe. La cytokinèse dans la levure à fission est contrôlée par le SIN, un réseau signalant de protéines kinases qui utilise le corps à pôle-fuseau comme échafaudage. Afin de décrire le comportement qualitatif du système et prédire des comportements mutants inconnus, nous avons décidé d'adopter l'approche de la modélisation booléenne. Dans cette thèse, nous présentons la construction d'un modèle booléen étendu du SIN, comprenant la plupart des composantes et des régulateurs du SIN en tant que noeuds individuels et testable expérimentalement. Ce modèle utilise des niveaux d'activité du CDK comme noeuds de contrôle pour la simulation d'évènements du SIN à différents stades du cycle cellulaire. Ce modèle a été optimisé en utilisant des expériences d'un seul "knock-out" avec des effets phénotypiques connus comme set d'entraînement. Il a permis de prédire correctement un set d'évaluation de "knock-out" doubles. De plus, le modèle a fait des prédictions in silico qui ont été validées in vivo, permettant d'obtenir de nouvelles idées de la régulation et l'organisation hiérarchique du SIN. Un autre projet concernant le cycle cellulaire qui fait partie de cette thèse a été la construction d'un modèle qualitatif et minimal de la réciprocité des cyclines dans la S.cerevisiae. Les protéines Clb dans la levure bourgeonnante présentent une activation et une dégradation caractéristique et séquentielle durant le cycle cellulaire, qu'on appelle communément les vagues des Clbs. Cet évènement est coordonné avec la courbe d'activation inverse du Sic1, qui a un rôle inhibitoire dans le système. Pour l'identification des modèles qualitatifs minimaux qui peuvent expliquer ce phénomène, nous avons sélectionné des expériences bien définies et construit tous les modèles minimaux possibles qui, une fois simulés, reproduisent les résultats attendus. Les modèles ont été filtrés en utilisant des simulations ODE qualitatives et standardisées; seules celles qui reproduisaient le phénotype des vagues ont été gardées. L'ensemble des modèles minimaux peut être utilisé pour suggérer des relations regulatoires entre les molécules participant qui peuvent ensuite être testées expérimentalement. Enfin, durant mon doctorat, j'ai participé au SBV Improver Challenge. Le but était de déduire des réseaux spécifiques à des espèces (humain et rat) en utilisant des données de phosphoprotéines, d'expressions des gènes et des cytokines, ainsi qu'un réseau de référence, qui était mis à disposition comme donnée préalable. Notre solution pour ce concours a pris la troisième place. L'approche utilisée est expliquée en détail dans le dernier chapitre de la thèse. -- The present dissertation is entitled "Development and Application of Computational Methodologies in Qualitative Modeling". It encompasses the diverse projects that were undertaken during my time as a PhD student. Instead of a systematic implementation of a framework defined a priori, this thesis should be considered as an exploration of the methods that can help us infer the blueprint of regulatory and signaling processes. This exploration was driven by concrete biological questions, rather than theoretical investigation. Even though the projects involved divergent systems (gene regulatory networks of cell cycle, signaling networks in lung cells), as well as organisms (fission yeast, budding yeast, rat, human), our goals were complementary and coherent. The main project of the thesis is the modeling of the Septation Initiation Network (SIN) in S.pombe. Cytokinesis in fission yeast is controlled by the SIN, a protein kinase signaling network that uses the spindle pole body as scaffold. In order to describe the qualitative behavior of the system and predict unknown mutant behaviors we decided to adopt a Boolean modeling approach. In this thesis, we report the construction of an extended, Boolean model of the SIN, comprising most SIN components and regulators as individual, experimentally testable nodes. The model uses CDK activity levels as control nodes for the simulation of SIN related events in different stages of the cell cycle. The model was optimized using single knock-out experiments of known phenotypic effect as a training set, and was able to correctly predict a double knock-out test set. Moreover, the model has made in silico predictions that have been validated in vivo, providing new insights into the regulation and hierarchical organization of the SIN. Another cell cycle related project that is part of this thesis was to create a qualitative, minimal model of cyclin interplay in S.cerevisiae. CLB proteins in budding yeast present a characteristic, sequential activation and decay during the cell cycle, commonly referred to as Clb waves. This event is coordinated with the inverse activation curve of Sic1, which has an inhibitory role in the system. To generate minimal qualitative models that can explain this phenomenon, we selected well-defined experiments and constructed all possible minimal models that, when simulated, reproduce the expected results. The models were filtered using standardized qualitative ODE simulations; only the ones reproducing the wave-like phenotype were kept. The set of minimal models can be used to suggest regulatory relations among the participating molecules, which will subsequently be tested experimentally. Finally, during my PhD I participated in the SBV Improver Challenge. The goal was to infer species-specific (human and rat) networks, using phosphoprotein, gene expression and cytokine data and a reference network provided as prior knowledge. Our solution to the challenge was selected as in the final chapter of the thesis.