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.

Multispectral Analysis of Paper Production Phases

Paperitehtaissa paperin väriä arvioidaan L*a*b –arvojen perusteella. Tämä lähestymistapa ei ole tarkka ja johtaa siihen, että saman luokkaan kuuluvat paperit saatetaan havaita eri värisinä. Tässä diplomityössä tarkastellaan värien arviointia multispektraalisen kuvantamisen avulla. Menetelmän ajatuksena on verrata standardoidun heijastuman spektriä kohteena olevan paperirullan heijastuman spektriin mittaustapojen avulla. 12:n kirjallisuuteen perustuvan ja yhden tutkimusryhmämme kehittämän mittaustavan joukosta pyrittiin valitsemaan parhaat mittaustavat tutkimuksien perusteella. Mittaustavan tehokkuus tarkoittaa kykä tunnistaa eroavaisuuksia värisävyn, värikylläisyyden ja valoisuuden perusteella. Mittaustapoja tutkittiin sekä keinotekoisilla kuvilla että aidoilla paperinäytteillä. Tuloksena maksimi-minimi-, aritmeettisen keskiarvon minimi- ja geometrisen keskiarvon minimimenetelmä olivat parhaita mittaustapoja.

The ADO Index as a Predictor of Two-Year Mortality in General Practice-Based Chronic Obstructive Pulmonary Disease Cohorts.

BACKGROUND: Existing prediction models for mortality in chronic obstructive pulmonary disease (COPD) patients have not yet been validated in primary care, which is where the majority of patients receive care. OBJECTIVES: Our aim was to validate the ADO (age, dyspnoea, airflow obstruction) index as a predictor of 2-year mortality in 2 general practice-based COPD cohorts. METHODS: Six hundred and forty-six patients with COPD with GOLD (Global Initiative for Chronic Obstructive Lung Disease) stages I-IV were enrolled by their general practitioners and followed for 2 years. The ADO regression equation was used to predict a 2-year risk of all-cause mortality in each patient and this risk was compared with the observed 2-year mortality. Discrimination and calibration were assessed as well as the strength of association between the 15-point ADO score and the observed 2-year all-cause mortality. RESULTS: Fifty-two (8.1%) patients died during the 2-year follow-up period. Discrimination with the ADO index was excellent with an area under the curve of 0.78 [95% confidence interval (CI) 0.71-0.84]. Overall, the predicted and observed risks matched well and visual inspection revealed no important differences between them across 10 risk classes (p = 0.68). The odds ratio for death per point increase according to the ADO index was 1.50 (95% CI 1.31-1.71). CONCLUSIONS: The ADO index showed excellent prediction properties in an out-of-population validation carried out in COPD patients from primary care settings. © 2014 S. Karger AG, Basel.

Determination of Rod and Cone Influence to the Early and Late Dynamic of the Pupillary Light Response.

PURPOSE: This study aims to identify which aspects of the pupil light reflex are most influenced by rods and cones independently by analyzing pupil recordings from different mouse models of photoreceptor deficiency. METHODS: One-month-old wild type (WT), rodless (Rho-/-), coneless (Cnga3-/-), or photoreceptor less (Cnga3-/-; Rho-/- or Gnat1-/-) mice were subjected to brief red and blue light stimuli of increasing intensity. To describe the initial dynamic response to light, the maximal pupillary constriction amplitudes and the derivative curve of the first 3 seconds were determined. To estimate the postillumination phase, the constriction amplitude at 9.5 seconds after light termination was related to the maximal constriction amplitude. RESULTS: Rho-/- mice showed decreased constriction amplitude but more prolonged pupilloconstriction to all blue and red light stimuli compared to wild type mice. Cnga3-/- mice had constriction amplitudes similar to WT however following maximal constriction, the early and rapid dilation to low intensity blue light was decreased. To high intensity blue light, the Cnga3-/- mice demonstrated marked prolongation of the pupillary constriction. Cnga3-/-; Rho-/- mice had no pupil response to red light of low and medium intensity. CONCLUSIONS: From specific gene defective mouse models which selectively voided the rod or cone function, we determined that mouse rod photoreceptors are highly contributing to the pupil response to blue light stimuli but also to low and medium red stimuli. We also observed that cone cells mainly drive the partial rapid dilation of the initial response to low blue light stimuli. Thus photoreceptor dysfunction can be derived from chromatic pupillometry in mouse models.

Performance of 18F-FET versus 18F-FDG-PET for the diagnosis and grading of brain tumors: systematic review and meta-analysis.

BACKGROUND: For the past decade (18)F-fluoro-ethyl-l-tyrosine (FET) and (18)F-fluoro-deoxy-glucose (FDG) positron emission tomography (PET) have been used for the assessment of patients with brain tumor. However, direct comparison studies reported only limited numbers of patients. Our purpose was to compare the diagnostic performance of FET and FDG-PET. METHODS: We examined studies published between January 1995 and January 2015 in the PubMed database. To be included the study should: (i) use FET and FDG-PET for the assessment of patients with isolated brain lesion and (ii) use histology as the gold standard. Analysis was performed on a per patient basis. Study quality was assessed with STARD and QUADAS criteria. RESULTS: Five studies (119 patients) were included. For the diagnosis of brain tumor, FET-PET demonstrated a pooled sensitivity of 0.94 (95% CI: 0.79-0.98) and pooled specificity of 0.88 (95% CI: 0.37-0.99), with an area under the curve of 0.96 (95% CI: 0.94-0.97), a positive likelihood ratio (LR+) of 8.1 (95% CI: 0.8-80.6), and a negative likelihood ratio (LR-) of 0.07 (95% CI: 0.02-0.30), while FDG-PET demonstrated a sensitivity of 0.38 (95% CI: 0.27-0.50) and specificity of 0.86 (95% CI: 0.31-0.99), with an area under the curve of 0.40 (95% CI: 0.36-0.44), an LR+ of 2.7 (95% CI: 0.3-27.8), and an LR- of 0.72 (95% CI: 0.47-1.11). Target-to-background ratios of either FDG or FET, however, allow distinction between low- and high-grade gliomas (P > .11). CONCLUSIONS: For brain tumor diagnosis, FET-PET performed much better than FDG and should be preferred when assessing a new isolated brain tumor. For glioma grading, however, both tracers showed similar performances.

Photoinhibition of Photosystem II. Kinetics, Photoprotection and Mechanism

Photosystem II (PSII) is susceptible to light-induced damage defined as photoinhibition. In natural conditions, plants are capable of repairing the photoinhibited PSII by on-going degradation and re-synthesis of the D1 reaction centre protein of PSII. Photoinhibition is induced by both visible and ultraviolet light and photoinhibition occurs under all light intensities with the same efficiency per photon. In my thesis work, I studied the reaction kinetics and mechanism of photoinhibition of PSII, as well as photoprotection in leaves of higher plants. Action spectroscopy was used to identify photoreceptors of photoinhibition. I found that the action spectrum of photoinhibition in vivo shows resemblance to the absorption spectra of manganese model compounds of the oxygen evolving complex (OEC) suggesting a role for manganese as a photoreceptor of photoinhibition under UV and visible light. In order to study the protective effect of non-photochemical quenching, the action spectrum was measured from leaves of wild type Arabidopsis thaliana and two mutants impaired in nonphotochemical quenching of chlorophyll a excitations. The findings of action spectroscopy and simulations of chlorophyll-based photoinhibition mechanisms suggested that quenching of antenna excitations protects less efficiently than would be expected if antenna chlorophylls were the only photoreceptors of photoinhibition. The reaction kinetics of prolonged photoinhibition was studied in leaves of Cucurbita maxima and Capsicum annuum. The results indicated that photoinhibitory decrease in both the oxygen evolution activity and ratio of variable to maximum fluorescence follows firstorder kinetics in vivo. The persistence of first-order kinetics suggests that already photoinhibited reaction centres do not protect against photoinhibition and that the mechanism of photoinhibition does not have a reversible intermediate. When Cucurbita maxima leaves were photoinhibited with saturating single-turnover flashes and continuous light, the light response curve of photoinhibition was found to be essentially a straight line with both types of illumination, suggesting that similar photoinhibition mechanisms might function during illumination with continuous light and during illumination with short flashes.

Interaction study of moxifloxacin with Cu(II) ion using square-wave voltammetry and its application in the determination in tablets

This work presents an electroanalytical method for the determination of moxifloxacin (MOXI) in tablets by its interaction with Cu(II) ion and subsequent electrochemical reduction at hanging mercury drop electrode (HMDE). A well-defined reduction peak at -0.21 V vs. Ag/AgCl in Phosphate buffer 0.04 mol L-1 pH 8.0 was observed for the complex reduction MOXI-Cu(II), using square-wave voltammetry (SWV). Using a 10 s of accumulation time at -0.40 V was found a limit detection of 3.60x10-8 mol l-1. The obtained results have shown good agreement with those obtained by spectrophotometric method.

Metabolic evaluation of dairy cows submitted to three different strategies to decrease the effects of negative energy balance in early postpartum

In early lactation dairy cattle suffer metabolic alterations caused by negative energy balance, which predisposes to fatty liver and ketosis. The aim of this study was to evaluate the metabolic condition of high yielding dairy cows subjected to three treatments for preventing severe lipomobilization and ketosis in early lactation. Fifty four multiparous Holstein cows yielding >30 L/day were divided into four groups: control (CN= no treatment), glucose precursor (PG= propylene-glycol), hepatic protector (Mp= Mercepton®), and energy supplement with salts of linolenic and linoleic faty acids (Mg-E= Megalac-E®). Treatments were administrated randomly at moment of calving until 8 weeks postpartum. Blood samples were collected on days 1, 7, 14, 21, 28, 35, 42 and 49 postpartum. Body condition score (BCS) was evaluated at the same periods and milk yield was recorded at 2nd, 4th, 5th, 6th, 7th, and 8th weeks of lactation. Concentrations of non-esterified fatty acids (NEFA), albumin, AST, ß-hydroxybutyrate (BHBA), cholesterol, glucose, total protein, urea and triglycerides were analyzed in blood samples. Cut-off points for subclinical ketosis were defined when BHBA >1.4 mmol/L and NEFA >0.7 mmol/L. General occurrence of subclinical ketosis was 24% during the period. An ascendant curve of cholesterol and glucose was observed from the 1st to the 8th week of lactation, while any tendency was observed with BHBA and NEFA, although differences among treatments were detected (p<0.05). BCS decreased from a mean of 3.85 at 1st week to 2.53 at 8th week of lactation (p=0.001). Milk yield was higher in the Mg-E group compared with the other treatment groups (p<0.05) Compared with the CN group, the treatments with Mp and PG did not show significant differences in blood biochemistry and milk yield. Cows receiving PG and Mg-E showed higher values of BHBA and NEFA (P<0.05), indicating accentuated lipomobilization. Supplementation with Mg-E also resulted in significant higher concentrations of cholesterol, BHBA, urea, AST and lower values of glycemia. This performance may be explained by the highest milk yield observed with this treatment. Treatments with PG and Mp did not improve milk yield, compared with control cows, but did not show metabolic evidence of ketosis, fat mobilization or fatty liver. These results suggest that treatment with Mg-E improves milk production but induces a higher negative energy balance leading to moderated lipomobilization and ketone bodies production, increasing the risk of fatty liver.

Algorithm to determine the intersection curves between bezier surfaces by the solution of multivariable polynomial system and the differential marching method

The determination of the intersection curve between Bézier Surfaces may be seen as the composition of two separated problems: determining initial points and tracing the intersection curve from these points. The Bézier Surface is represented by a parametric function (polynomial with two variables) that maps a point in the tridimensional space from the bidimensional parametric space. In this article, it is proposed an algorithm to determine the initial points of the intersection curve of Bézier Surfaces, based on the solution of polynomial systems with the Projected Polyhedral Method, followed by a method for tracing the intersection curves (Marching Method with differential equations). In order to allow the use of the Projected Polyhedral Method, the equations of the system must be represented in terms of the Bernstein basis, and towards this goal it is proposed a robust and reliable algorithm to exactly transform a multivariable polynomial in terms of power basis to a polynomial written in terms of Bernstein basis .

A study of the facial aging - a multidisciplinary approach

This paper describes a mathematical and graphical model for face aging. It considers the possibility of predicting the aging process by offering an initial quantification of this process as it applies to the face. It is concerned with physical measurements and a general law of time dependence. After measuring and normalizing a photograph of a person, one could predict, with a known amount of error, the appearance of that person at a different age. The technique described has served its purpose successfully, with a representative amount of patient data behaving sufficiently near the general aging curve of each parameter. That model uses a warping technique to emulate the aging changes on the face of women. Frequently the warping methods are based on the interpolation between images or general mathematical functions to calculate the pixel attributes. The implemented process considers the age features of selected parts of a face such as the face outline and the shape of the lips. These age features were obtained by measuring the facial regions of women that have been photographed throughout their lives. The present work is first concerned with discussing a methodology to define the aging parameters that can be measured, and second with representing the age effects graphically.

Viscosity of gums in vitro and their ability to reduce postprandial hyperglycemia in normal subjects

Experiments were carried out in vitro with three viscous polysaccharides (guar gum, pectin, and carboxymethylcellulose (CMC)) of similar initial viscosity submitted to conditions that mimic events occurring in the stomach and duodenum, and their viscosity in these situations was compared to their actions on postprandial hyperglycemia in normal human subjects. Guar gum showed greater viscosity than the other gums during acidification and/or alkalinization and also showed larger effects on plasma glucose levels (35% reduction in maximum rise in plasma glucose) and on the total area under the curve of plasma glucose (control: 20,314 ± 1007 mg dl-1 180 min-1 vs guar gum: 18,277 ± 699 mg dl-1 180 min-1, P<0.01). Pectin, which showed a marked reduction in viscosity at 37oC and after events mimicking those that occur in the stomach and duodenum, did not have a significant effect on postprandial hyperglycemia. The performance of viscosity and the glycemia response to CMC were at an intermediate level between guar gum and pectin. In conclusion, these data suggest that temperature, the process of acidification, alkalinization and exposure to intestinal ions induce different viscosity changes in gums having similar initial viscosity, establishing a direct relationship between a minor decrease of gum viscosity in vitro and a reduction of postprandial hyperglycemia

Determination of anaerobic threshold in rats using the lactate minimum test

The break point of the curve of blood lactate vs exercise load has been called anaerobic threshold (AT) and is considered to be an important indicator of endurance exercise capacity in human subjects. There are few studies of AT determination in animals. We describe a protocol for AT determination by the "lactate minimum test" in rats during swimming exercise. The test is based on the premise that during an incremental exercise test, and after a bout of maximal exercise, blood lactate decreases to a minimum and then increases again. This minimum value indicates the intensity of the AT. Adult male (90 days) Wistar rats adapted to swimming for 2 weeks were used. The initial state of lactic acidosis was obtained by making the animals jump into the water and swim while carrying a load equivalent to 50% of body weight for 6 min (30-s exercise interrupted by a 30-s rest). After a 9-min rest, blood was collected and the incremental swimming test was started. The test consisted of swimming while supporting loads of 4.5, 5.0, 5.5, 6.0 and 7.0% of body weight. Each exercise load lasted 5 min and was followed by a 30-s rest during which blood samples were taken. The blood lactate minimum was determined from a zero-gradient tangent to a spline function fitting the blood lactate vs workload curve. AT was estimated to be 4.95 ± 0.10% of body weight while interpolated blood lactate was 7.17 ± 0.16 mmol/l. These results suggest the application of AT determination in animal studies concerning metabolism during exercise.

An analysis of the electrocardiogram QT interval

No reports testing the efficacy of the use of the QT/RR ratio <1/2 for detecting a normal QTc interval were found in the literature. The objective of the present study was to determine if a QT/RR ratio <=1/2 can be considered to be equal to the normal QTc and to compare the QT and QTc measured and calculated clinically and by a computerized electrocardiograph. Ratios (140 QT/RR) of 28 successive electrocardiograms obtained from 28 consecutive patients in a tertiary level teaching hospital were analyzed clinically by 5 independent observers and by a computerized electrocardiograph. The QT/RR ratio provided 56% sensitivity and 78% specificity, with an area under the receiver operator characteristic curve of 75.8% (95%CI: 0.68 to 0.84). The divergence in QT and QTc interval measurements between clinical and computerized evaluation were 0.01 ± 0.03 s (95%CI: 0.04-0.02) and 0.01 ± 0.04 s (95%CI: -0.05-0.03), respectively. The QT and QTc values measured clinically and by a computerized electrocardiograph were similar. The QT/RR ratio <=1/2 was not a satisfactory index for QTc evaluation because it could not predict a normal QTc value.

Automated grading of left ventricular segmental wall motion by an artificial neural network using color kinesis images

The present study describes an auxiliary tool in the diagnosis of left ventricular (LV) segmental wall motion (WM) abnormalities based on color-coded echocardiographic WM images. An artificial neural network (ANN) was developed and validated for grading LV segmental WM using data from color kinesis (CK) images, a technique developed to display the timing and magnitude of global and regional WM in real time. We evaluated 21 normal subjects and 20 patients with LVWM abnormalities revealed by two-dimensional echocardiography. CK images were obtained in two sets of viewing planes. A method was developed to analyze CK images, providing quantitation of fractional area change in each of the 16 LV segments. Two experienced observers analyzed LVWM from two-dimensional images and scored them as: 1) normal, 2) mild hypokinesia, 3) moderate hypokinesia, 4) severe hypokinesia, 5) akinesia, and 6) dyskinesia. Based on expert analysis of 10 normal subjects and 10 patients, we trained a multilayer perceptron ANN using a back-propagation algorithm to provide automated grading of LVWM, and this ANN was then tested in the remaining subjects. Excellent concordance between expert and ANN analysis was shown by ROC curve analysis, with measured area under the curve of 0.975. An excellent correlation was also obtained for global LV segmental WM index by expert and ANN analysis (R² = 0.99). In conclusion, ANN showed high accuracy for automated semi-quantitative grading of WM based on CK images. This technique can be an important aid, improving diagnostic accuracy and reducing inter-observer variability in scoring segmental LVWM.

Toward including the effect of manufacturing processes in the pre-estimated losses of the switched reluctance motor

The estimation of losses plays a key role in the process of building any electrical machine. How to estimate those losses while designing any machine; by obtaining the characteristic of the electrical steel from the catalogue and calculate the losses. However, this way is inaccurate since the electrical steel performs several manufacturing processes during the process of building any machine, which affects directly the magnetic property of the electrical steel and accordingly the characteristic of the electrical steel will be affected. That means the B–H curve of the steel that was obtained from the catalogue will be changed. Moreover, during loading and rotating the machine, some important changes occur to the B–H characteristic of the electrical steel such as the stress on the laminated iron. Accordingly, the pre-estimated losses are completely far from the actual losses because they were estimated based on the data of the electrical steel obtained from the catalogue. So in order to estimate the losses precisely significant factors of the manufacturing processes must be included. The paper introduces the systematic estimation of the losses including the effect of one of the manufacturing factors. Similarly, any other manufacturing factor can be included in the pre-designed losses estimations.