EADock : design of a new molecular docking algorithm and some of its applications

3 Summary 3. 1 English The pharmaceutical industry has been facing several challenges during the last years, and the optimization of their drug discovery pipeline is believed to be the only viable solution. High-throughput techniques do participate actively to this optimization, especially when complemented by computational approaches aiming at rationalizing the enormous amount of information that they can produce. In siiico techniques, such as virtual screening or rational drug design, are now routinely used to guide drug discovery. Both heavily rely on the prediction of the molecular interaction (docking) occurring between drug-like molecules and a therapeutically relevant target. Several softwares are available to this end, but despite the very promising picture drawn in most benchmarks, they still hold several hidden weaknesses. As pointed out in several recent reviews, the docking problem is far from being solved, and there is now a need for methods able to identify binding modes with a high accuracy, which is essential to reliably compute the binding free energy of the ligand. This quantity is directly linked to its affinity and can be related to its biological activity. Accurate docking algorithms are thus critical for both the discovery and the rational optimization of new drugs. In this thesis, a new docking software aiming at this goal is presented, EADock. It uses a hybrid evolutionary algorithm with two fitness functions, in combination with a sophisticated management of the diversity. EADock is interfaced with .the CHARMM package for energy calculations and coordinate handling. A validation was carried out on 37 crystallized protein-ligand complexes featuring 11 different proteins. The search space was defined as a sphere of 15 R around the center of mass of the ligand position in the crystal structure, and conversely to other benchmarks, our algorithms was fed with optimized ligand positions up to 10 A root mean square deviation 2MSD) from the crystal structure. This validation illustrates the efficiency of our sampling heuristic, as correct binding modes, defined by a RMSD to the crystal structure lower than 2 A, were identified and ranked first for 68% of the complexes. The success rate increases to 78% when considering the five best-ranked clusters, and 92% when all clusters present in the last generation are taken into account. Most failures in this benchmark could be explained by the presence of crystal contacts in the experimental structure. EADock has been used to understand molecular interactions involved in the regulation of the Na,K ATPase, and in the activation of the nuclear hormone peroxisome proliferatoractivated receptors a (PPARa). It also helped to understand the action of common pollutants (phthalates) on PPARy, and the impact of biotransformations of the anticancer drug Imatinib (Gleevec®) on its binding mode to the Bcr-Abl tyrosine kinase. Finally, a fragment-based rational drug design approach using EADock was developed, and led to the successful design of new peptidic ligands for the a5ß1 integrin, and for the human PPARa. In both cases, the designed peptides presented activities comparable to that of well-established ligands such as the anticancer drug Cilengitide and Wy14,643, respectively. 3.2 French Les récentes difficultés de l'industrie pharmaceutique ne semblent pouvoir se résoudre que par l'optimisation de leur processus de développement de médicaments. Cette dernière implique de plus en plus. de techniques dites "haut-débit", particulièrement efficaces lorsqu'elles sont couplées aux outils informatiques permettant de gérer la masse de données produite. Désormais, les approches in silico telles que le criblage virtuel ou la conception rationnelle de nouvelles molécules sont utilisées couramment. Toutes deux reposent sur la capacité à prédire les détails de l'interaction moléculaire entre une molécule ressemblant à un principe actif (PA) et une protéine cible ayant un intérêt thérapeutique. Les comparatifs de logiciels s'attaquant à cette prédiction sont flatteurs, mais plusieurs problèmes subsistent. La littérature récente tend à remettre en cause leur fiabilité, affirmant l'émergence .d'un besoin pour des approches plus précises du mode d'interaction. Cette précision est essentielle au calcul de l'énergie libre de liaison, qui est directement liée à l'affinité du PA potentiel pour la protéine cible, et indirectement liée à son activité biologique. Une prédiction précise est d'une importance toute particulière pour la découverte et l'optimisation de nouvelles molécules actives. Cette thèse présente un nouveau logiciel, EADock, mettant en avant une telle précision. Cet algorithme évolutionnaire hybride utilise deux pressions de sélections, combinées à une gestion de la diversité sophistiquée. EADock repose sur CHARMM pour les calculs d'énergie et la gestion des coordonnées atomiques. Sa validation a été effectuée sur 37 complexes protéine-ligand cristallisés, incluant 11 protéines différentes. L'espace de recherche a été étendu à une sphère de 151 de rayon autour du centre de masse du ligand cristallisé, et contrairement aux comparatifs habituels, l'algorithme est parti de solutions optimisées présentant un RMSD jusqu'à 10 R par rapport à la structure cristalline. Cette validation a permis de mettre en évidence l'efficacité de notre heuristique de recherche car des modes d'interactions présentant un RMSD inférieur à 2 R par rapport à la structure cristalline ont été classés premier pour 68% des complexes. Lorsque les cinq meilleures solutions sont prises en compte, le taux de succès grimpe à 78%, et 92% lorsque la totalité de la dernière génération est prise en compte. La plupart des erreurs de prédiction sont imputables à la présence de contacts cristallins. Depuis, EADock a été utilisé pour comprendre les mécanismes moléculaires impliqués dans la régulation de la Na,K ATPase et dans l'activation du peroxisome proliferatoractivated receptor a (PPARa). Il a également permis de décrire l'interaction de polluants couramment rencontrés sur PPARy, ainsi que l'influence de la métabolisation de l'Imatinib (PA anticancéreux) sur la fixation à la kinase Bcr-Abl. Une approche basée sur la prédiction des interactions de fragments moléculaires avec protéine cible est également proposée. Elle a permis la découverte de nouveaux ligands peptidiques de PPARa et de l'intégrine a5ß1. Dans les deux cas, l'activité de ces nouveaux peptides est comparable à celles de ligands bien établis, comme le Wy14,643 pour le premier, et le Cilengitide (PA anticancéreux) pour la seconde.

The design and characterization of receptor-selective APRIL variants.

A proliferation-inducing ligand (APRIL), a member of the TNF ligand superfamily with an important role in humoral immunity, is also implicated in several cancers as a prosurvival factor. APRIL binds two different TNF receptors, B cell maturation antigen (BCMA) and transmembrane activator and cylclophilin ligand interactor (TACI), and also interacts independently with heparan sulfate proteoglycans. Because APRIL shares binding of the TNF receptors with B cell activation factor, separating the precise signaling pathways activated by either ligand in a given context has proven quite difficult. In this study, we have used the protein design algorithm FoldX to successfully generate a BCMA-specific variant of APRIL, APRIL-R206E, and two TACI-selective variants, D132F and D132Y. These APRIL variants show selective activity toward their receptors in several in vitro assays. Moreover, we have used these ligands to show that BCMA and TACI have a distinct role in APRIL-induced B cell stimulation. We conclude that these ligands are useful tools for studying APRIL biology in the context of individual receptor activation.

Glucagon-like peptide-1 and control of insulin secretion.

Although glucose is the major regulator of insulin secretion by pancreatic beta cells, its action is modulated by several neural and hormonal stimuli. In particular, hormones secreted by intestinal endocrine cells stimulate glucose-induced insulin secretion very potently after nutrient absorption. These hormones, called gluco-incretins or insulinotropic hormones, are major regulators of postprandial glucose homeostasis. The main gluco-incretins are GIP (gastric inhibitory polypeptide or glucose-dependent insulinotropic polypeptide) and GLP-1 (glucagon-like polypeptide-1). The secretion of GIP, a 42 amino acid polypeptide secreted by duodenal K cells, is triggered by fat and glucose. GIP stimulation of insulin secretion depends on the presence of specific beta-cell receptors and requires glucose at a concentration at least equal to or higher than the normoglycaemic level of approximately 5 mM. GIP accounts for about 50% of incretin activity, and the rest may be due to GLP-1 which is produced by proteolytic processing of the preproglucagon molecule in intestinal L cells. GLP-1 is the most potent gluco-incretin characterized so far. As with GIP, its stimulatory action requires a specific membrane receptor and normal or elevated glucose concentrations. Contrary to GIP, the incretin effect of GLP-1 is maintained in non-insulin-dependent diabetic patients. This peptide or agonists of its beta-cell receptor could provide new therapeutic tools for the treatment of Type II diabetic hyperglycaemia.

N-myc controls proliferation, morphogenesis, and patterning of the inner ear

Myc family members play crucial roles in regulating cell proliferation, size, and differentiation during organogenesis. Both N-myc and c-myc are expressed throughout inner ear development. To address their function in the mouse inner ear, we generated mice with conditional deletions in either N-myc or c-myc. Loss of c-myc in the inner ear causes no apparent defects, whereas inactivation of N-myc results in reduced growth caused by a lack of proliferation. Reciprocally, the misexpression of N-myc in the inner ear increases proliferation. Morphogenesis of the inner ear in N-myc mouse mutants is severely disturbed, including loss of the lateral canal, fusion of the cochlea with the sacculus and utriculus, and stunted outgrowth of the cochlea. Mutant cochleas are characterized by an increased number of cells exiting the cell cycle that express the cyclin-dependent kinase inhibitor p27Kip1 and lack cyclin D1, both of which control the postmitotic state of hair cells. Analysis of different molecular markers in N-myc mutant ears reveals the development of a rudimentary organ of Corti containing hair cells and the underlying supporting cells. Differentiated cells, however, fail to form the highly ordered structure characteristic for the organ of Corti but appear as rows or clusters with an excess number of hair cells. The Kölliker's organ, a transient structure neighboring the organ of Corti and a potential source of ectopic hair cells, is absent in the mutant ears. Collectively, our data suggest that N-myc regulates growth, morphogenesis, and pattern formation during the development of the inner ear.

Human high-density lipoprotein particles prevent activation of the JNK pathway induced by human oxidised low-density lipoprotein particles in pancreatic beta cells.

AIMS/HYPOTHESIS: We explored the potential adverse effects of pro-atherogenic oxidised LDL-cholesterol particles on beta cell function. MATERIALS AND METHODS: Isolated human and rat islets and different insulin-secreting cell lines were incubated with human oxidised LDL with or without HDL particles. The insulin level was monitored by ELISA, real-time PCR and a rat insulin promoter construct linked to luciferase gene reporter. Cell apoptosis was determined by scoring cells displaying pycnotic nuclei. RESULTS: Prolonged incubation with human oxidised LDL particles led to a reduction in preproinsulin expression levels, whereas the insulin level was preserved in the presence of native LDL-cholesterol. The loss of insulin production occurred at the transcriptional levels and was associated with an increase in activator protein-1 transcriptional activity. The rise in activator protein-1 activity resulted from activation of c-Jun N-terminal kinases (JNK, now known as mitogen-activated protein kinase 8 [MAPK8]) due to a subsequent decrease in islet-brain 1 (IB1; now known as MAPK8 interacting protein 1) levels. Consistent with the pro-apoptotic role of the JNK pathway, oxidised LDL also induced a twofold increase in the rate of beta cell apoptosis. Treatment of the cells with JNK inhibitor peptides or HDL countered the effects mediated by oxidised LDL. CONCLUSIONS/INTERPRETATION: These data provide strong evidence that oxidised LDL particles exert deleterious effects in the progression of beta cell failure in diabetes and that these effects can be countered by HDL particles.

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.

Quantitative modelling of transcription factor binding sites

Abstract One of the most important issues in molecular biology is to understand regulatory mechanisms that control gene expression. Gene expression is often regulated by proteins, called transcription factors which bind to short (5 to 20 base pairs),degenerate segments of DNA. Experimental efforts towards understanding the sequence specificity of transcription factors is laborious and expensive, but can be substantially accelerated with the use of computational predictions. This thesis describes the use of algorithms and resources for transcriptionfactor binding site analysis in addressing quantitative modelling, where probabilitic models are built to represent binding properties of a transcription factor and can be used to find new functional binding sites in genomes. Initially, an open-access database(HTPSELEX) was created, holding high quality binding sequences for two eukaryotic families of transcription factors namely CTF/NF1 and LEFT/TCF. The binding sequences were elucidated using a recently described experimental procedure called HTP-SELEX, that allows generation of large number (> 1000) of binding sites using mass sequencing technology. For each HTP-SELEX experiments we also provide accurate primary experimental information about the protein material used, details of the wet lab protocol, an archive of sequencing trace files, and assembled clone sequences of binding sequences. The database also offers reasonably large SELEX libraries obtained with conventional low-throughput protocols.The database is available at http://wwwisrec.isb-sib.ch/htpselex/ and and ftp://ftp.isrec.isb-sib.ch/pub/databases/htpselex. The Expectation-Maximisation(EM) algorithm is one the frequently used methods to estimate probabilistic models to represent the sequence specificity of transcription factors. We present computer simulations in order to estimate the precision of EM estimated models as a function of data set parameters(like length of initial sequences, number of initial sequences, percentage of nonbinding sequences). We observed a remarkable robustness of the EM algorithm with regard to length of training sequences and the degree of contamination. The HTPSELEX database and the benchmarked results of the EM algorithm formed part of the foundation for the subsequent project, where a statistical framework called hidden Markov model has been developed to represent sequence specificity of the transcription factors CTF/NF1 and LEF1/TCF using the HTP-SELEX experiment data. The hidden Markov model framework is capable of both predicting and classifying CTF/NF1 and LEF1/TCF binding sites. A covariance analysis of the binding sites revealed non-independent base preferences at different nucleotide positions, providing insight into the binding mechanism. We next tested the LEF1/TCF model by computing binding scores for a set of LEF1/TCF binding sequences for which relative affinities were determined experimentally using non-linear regression. The predicted and experimentally determined binding affinities were in good correlation.

The DNA damage response to adeno-associated virus : centrosome overduplication and induction of cell death independently of the spindle checkpoint

Summary: Adeno-associated virus type 2 (AAV2) is a small virus containing single-stranded DNA of approximately 4.7kb in size. Both ends of the viral genome are flanked with inverted terminal repeat sequences (ITRs), which serve as primers for viral replication. Previous work in our laboratory has shown that AAV2 DNA with ultraviolet radiation-generated crosslinks (UV-AAV2) provokes a DNA damage response in the host cell by mimicking a stalled replication fork. Infection of cells with UV-AAV2 leads to a p53-and Chk1-mediated cell cycle arrest at the G2/M border of the cell cycle. However, tumour cells lacking the tumour suppressor protein p53 cannot sustain this arrest and enter a prolonged impaired mitosis, the outcome of which is cell death. The aim of my thesis was to investigate how UV-inactivated AAV2 kilts p53-deficient cancer cells. I found that the UV-AAV2-induced DNA damage signalling induces centriole overduplication in infected cells. The virus is able to uncouple the centriole duplication cycle from the cell cycle, leading to amplified centrosome numbers. Chk1 colocalises with centrosomes in the infected cells and the centrosome overduplication is dependent on the presence of Chk1, as well as on the activities of ATR and Cdk kinases and on the G2 arrest. The UV-AAV2-induced DNA damage signalling inhibits the degradation of cyclin B 1 and securin by the anaphase promoting complex, suggesting that the spindle checkpoint is activated in these mitotic cells. Interference with the spindle checkpoint components Mad2 and BubR1 revealed that the UV-AAV2-provoked mitotic catastrophe occurs independently of spindle checkpoint function, This work shows that, in the p53 deficient cells, UV-AAV2 triggers mitotic catastrophe associated with a dramatic Chk1-dependent overduplication of centrioles and the consequent formation of multiple spindle poles in mitosis. Résumé Le virus associé à l'adénovirus type 2 (AAV2) est un petit virus contenant un simple brin d'ADN d'environ 4.7kb. Des expériences antérieures dans notre laboratoire ont montré que les liens intramoléculaires sur l'ADN de AAV2 provoqués paz l'irradiation aux ultraviolets (UV) ressemblent à une fourche de réplication bloquée, ce qui provoque une réponse aux dommages à l'ADN dans la cellule hôte. L'infection des cellules avec UV-AAV2 résulte en un arrêt du cycle cellulaire à la transition G2/M entraîné par les protéines ATR et Chk1. Cependant, les cellules tumorales auxquelles il manque le suppresseur de tumeur p53 ne peuvent pas tenir cet arrêt et entrent dans une mitose anormale et prolongée qui se terminera par la mort cellulaire. Le but de ma thèse était d'étudier comment l'AAV2 inactivé par l'irradiation UV tue les cellules cancéreuses n'ayant pas p53. Je montre ici que le signal de dommages à l'ADN induit par UV-AAV2 génère une surduplication des centrioles dans les cellules infectées. Le virus est capable de dissocier le cycle de duplication du centriole du cycle cellulaire ce qui crée un nombre amplifié de centrosomes. Chk1 est co-localisé avec le centrosome dans les cellules infectées et la swduplication du centrosome est dépendante de la présence de Chk1, de l'activité des kinases ATR et Cdk et de l'arrêt en G2 de la cellule. Le signal d'ADN endommagé induit par UV-AAV2 réprime la dégradation des protéines cycline B1 et securine par le complexe promoteur de l'anaphase (APC), ce qui suggère que le point de contrôle du fuseau mitotique est activé dans ces cellules en mitose. L'étude d'interférence avec des éléments du point de contrôle du fuseau mitotique, Mad2 et BubR1, a révélé que la catastrophe mitotique provoquée paz UV-AAV2 survient indépendamment du point de contrôle du fuseau mitotique. Ce travail montre que dans les cellules déficientes en p53, UV-AAV2 induit une catastrophe mitotique associée à une surduplication des centrioles dépendant de Chk1 et ayant pour conséquence dramatique la formation de multiples fuseaux mitotiques dans la cellule en mitose.

Nestmate recognition and levels of aggression are not altered by changes in genetic diversity in a socially polymorphic ant

The ability to distinguish nestmates from foreign individuals is central to the functioning of insect societies. In ants, workers from multiple-queen colonies are often less aggressive than workers from single-queen ones. In line with this observation, it has been hypothesized that workers from multiple-queen colonies have less precise recognition abilities than workers from single-queen ones because their colonies contain genetically more diverse individuals, which results in a broader template of recognition cues. Here, we assessed the impact of social structure ( queen number) variation on nestmate recognition and aggression in a large population of the socially polymorphic ant Formica selysi. We staged unilateral aggression tests on the nest surface. Workers from single-and multiple-queen colonies had good nestmate recognition ability and did not differ significantly in their level of aggression towards foreign, immobilized workers ( cue-bearers). In particular, workers from multiple-queen colonies efficiently recognized non-nestmates despite the higher genetic diversity in their colony. Cue-bearers from single- and multiple-queen colonies elicited similar reactions. However, the level of aggression was higher between than within social forms, suggesting that workers detect a signal that is specific to the colony social structure. Finally, the level of aggression was not correlated with the genetic distance between colonies. Overall, we found no evidence for the hypothesis that the presence of multiple breeders in the same colony decreases recognition abilities and found no simple relationship between genetic diversity and aggression level. (c) 2007 The Association for the Study of Animal Behaviou

Fabrication and characterization of silicon position sensitive particle detectors

Position sensitive particle detectors are needed in high energy physics research. This thesis describes the development of fabrication processes and characterization techniques of silicon microstrip detectors used in the work for searching elementary particles in the European center for nuclear research, CERN. The detectors give an electrical signal along the particles trajectory after a collision in the particle accelerator. The trajectories give information about the nature of the particle in the struggle to reveal the structure of the matter and the universe. Detectors made of semiconductors have a better position resolution than conventional wire chamber detectors. Silicon semiconductor is overwhelmingly used as a detector material because of its cheapness and standard usage in integrated circuit industry. After a short spread sheet analysis of the basic building block of radiation detectors, the pn junction, the operation of a silicon radiation detector is discussed in general. The microstrip detector is then introduced and the detailed structure of a double-sided ac-coupled strip detector revealed. The fabrication aspects of strip detectors are discussedstarting from the process development and general principles ending up to the description of the double-sided ac-coupled strip detector process. Recombination and generation lifetime measurements in radiation detectors are discussed shortly. The results of electrical tests, ie. measuring the leakage currents and bias resistors, are displayed. The beam test setups and the results, the signal to noise ratio and the position accuracy, are then described. It was found out in earlier research that a heavy irradiation changes the properties of radiation detectors dramatically. A scanning electron microscope method was developed to measure the electric potential and field inside irradiated detectorsto see how a high radiation fluence changes them. The method and the most important results are discussed shortly.

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.

Design and implementation of FPGA-based LQ control of active magnetic bearings

The need for high performance, high precision, and energy saving in rotating machinery demands an alternative solution to traditional bearings. Because of the contactless operation principle, the rotating machines employing active magnetic bearings (AMBs) provide many advantages over the traditional ones. The advantages such as contamination-free operation, low maintenance costs, high rotational speeds, low parasitic losses, programmable stiffness and damping, and vibration insulation come at expense of high cost, and complex technical solution. All these properties make the use of AMBs appropriate primarily for specific and highly demanding applications. High performance and high precision control requires model-based control methods and accurate models of the flexible rotor. In turn, complex models lead to high-order controllers and feature considerable computational burden. Fortunately, in the last few years the advancements in signal processing devices provide new perspective on the real-time control of AMBs. The design and the real-time digital implementation of the high-order LQ controllers, which focus on fast execution times, are the subjects of this work. In particular, the control design and implementation in the field programmable gate array (FPGA) circuits are investigated. The optimal design is guided by the physical constraints of the system for selecting the optimal weighting matrices. The plant model is complemented by augmenting appropriate disturbance models. The compensation of the force-field nonlinearities is proposed for decreasing the uncertainty of the actuator. A disturbance-observer-based unbalance compensation for canceling the magnetic force vibrations or vibrations in the measured positions is presented. The theoretical studies are verified by the practical experiments utilizing a custom-built laboratory test rig. The test rig uses a prototyping control platform developed in the scope of this work. To sum up, the work makes a step in the direction of an embedded single-chip FPGA-based controller of AMBs.

Role of rhinovirus load in the upper respiratory tract and severity of symptoms in lung transplant recipients.

BACKGROUND: Rhinovirus is the most common cause of respiratory viral infections and leads to frequent respiratory symptoms in lung transplant recipients. However, it remains unknown whether the rhinovirus load correlates with the severity of symptoms. OBJECTIVES: This study aimed to better characterize the pathogenesis of rhinoviral infection and the way in which viral load correlates with symptoms. STUDY DESIGN: We assessed rhinovirus load in positive upper respiratory specimens of patients enrolled prospectively in a cohort of 116 lung transplant recipients. Rhinovirus load was quantified according to a validated in-house, real-time, reverse transcription polymerase chain reaction in pooled nasopharyngeal and pharyngeal swabs. Symptoms were recorded in a standardised case report form completed at each screening/routine follow-up visit, or during any emergency visit occurring during the 3-year study. RESULTS: Rhinovirus infections were very frequent, including in asymptomatic patients not seeking a specific medical consultation. Rhinovirus load ranged between 4.1 and 8.3 log copies/ml according to the type of visit and clinical presentation. Patients with highest symptom scores tended to have higher viral loads, particularly those presenting systemic symptoms. When considering symptoms individually, rhinovirus load was significantly higher in the presence of symptoms such as sore throat, fever, sputum production, cough, and fatigue. There was no association between tacrolimus levels and rhinovirus load. CONCLUSIONS: Rhinovirus infections are very frequent in lung transplant recipients and rhinoviral load in the upper respiratory tract is relatively high even in asymptomatic patients. Patients with the highest symptom scores tend to have a higher rhinovirus load.

Insertion of a malE B-Galactosidase fusion protein into the envelope of Escherichia coli disrupts biogenesis of outer membrane proteins and processing of inner membrane proteins

The synthesis of a membrane-bound MalE ,B-galactosidase hybrid protein, when induced by growth of Escherichia coli on maltose, leads to inhibition of cell division and eventually a reduced rate of mass increase. In addition, the relative rate of synthesis of outer membrane proteins, but not that of inner membrane proteins, was reduced by about 50%o. Kinetic experiments demonstrated that this reduction coincided with the period of maximum synthesis of the hybrid protein (and another maltose-inducible protein, LamB). The accumulation of this abnormal protein in the envelope therefore appeared specifically to inhibit the synthesis, the assembly of outer membrane proteins, or both, indicating that the hybrid protein blocks some export site or causes the sequestration of some limiting factor(s) involved in the export process. Since the MalE protein is normally located in the periplasm, the results also suggest that the synthesis of periplasmic and outer membrane proteins may involve some steps in common. The reduced rate of synthesis of outer membrane proteins was also accompanied by the accumulation in the envelope of at least one outer membrane protein and at least two inner membrane proteins as higher-molecular-weight forms, indicating that processing (removal of the N-terminal signal sequence) was also disrupted by the presence of the hybrid protein. These results may indicate that the assembly of these membrane proteins is blocked at a relatively late step rather than at the level of primary recognition of some site by the signal sequence. In addition, the results suggest that some step common to the biogenesis of quite different kinds of envelope protein is blocked by the presence of the hybrid protein.

Selection of internal control genes for real-time quantitative PCR in ovary and uterus of sows across pregnancy

Reproductive traits play a key role in pig production in order to reduce costs and increase economic returns. Among others, gene expression analyses represent a useful approach to study genetic mechanisms underlying reproductive traits in pigs. The application of reverse-transcription quantitative PCR requires the selection of appropriate reference genes, whose expression levels should not be affected by the experimental conditions, especially when comparing gene expression across different physiological stages.