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.

Evolution and cellular function of monothiol glutaredoxins: involvement in iron-sulphur cluster assembly

A number of bacterial species, mostly proteobacteria, possess monothiol glutaredoxins homologous to the Saccharomyces cerevisiae mitochondrial protein Grx5, which is involved in iron–sulphur cluster synthesis. Phylogenetic profiling is used to predict that bacterial monothiol glutaredoxins also participate in the iron–sulphur cluster (ISC) assembly machinery, because their phylogenetic profiles are similar to the profiles of the bacterial homologues of yeast ISC proteins. High evolutionary cooccurrence is observed between the Grx5 homologues and the homologues of the Yah1 ferredoxin, the scaffold proteins Isa1 and Isa2, the frataxin protein Yfh1 and the Nfu1 protein. This suggests that a specific functional interaction exists between these ISC machinery proteins. Physical interaction analyses using low-definition protein docking predict the formation of strong and specific complexes between Grx5 and several components of the yeast ISC machinery. Two-hybrid analysis has confirmed the in vivo interaction between Grx5 and Isa1. Sequence comparison techniques and cladistics indicate that the other two monothiol glutaredoxins of S. cerevisiae, Grx3 and Grx4, have evolved from the fusion of a thioredoxin gene with a monothiol glutaredoxin gene early in the eukaryotic lineage, leading to differential functional specialization. While bacteria do not contain these chimaeric glutaredoxins, in many eukaryotic species Grx5 and Grx3/4-type monothiol glutaredoxins coexist in the cell.

Diphenyl urea derivatives as inhibitors of transketolase: a structure-based virtual screening.

Transketolase is an enzyme involved in a critical step of the non-oxidative branch of the pentose phosphate pathway whose inhibition could lead to new anticancer drugs. Here, we report new human transketolase inhibitors, based on the phenyl urea scaffold, found by applying structure-based virtual screening. These inhibitors are designed to cover a hot spot in the dimerization interface of the homodimer of the enzyme, providing for the first time compounds with a suggested novel binding mode not based on mimicking the thiamine pyrophosphate cofactor.

Impact of platelet rich plasma and adipose stem cells on lymphangiogenesis in a murine tail lymphedema model.

BACKGROUND: Lymphedema is an underdiagnosed pathology which in industrialized countries mainly affects cancer patients that underwent lymph node dissection and/or radiation. Currently no effective therapy is available so that patients' life quality is compromised by swellings of the concerned body region. This unfortunate condition is associated with body imbalance and subsequent osteochondral deformations and impaired function as well as with an increased risk of potentially life threatening soft tissue infections. METHODS: The effects of PRP and ASC on angiogenesis (anti-CD31 staining), microcirculation (Laser Doppler Imaging), lymphangiogenesis (anti-LYVE1 staining), microvascular architecture (corrosion casting) and wound healing (digital planimetry) are studied in a murine tail lymphedema model. RESULTS: Wounds treated by PRP and ASC healed faster and showed a significantly increased epithelialization mainly from the proximal wound margin. The application of PRP induced a significantly increased lymphangiogenesis while the application of ASC did not induce any significant change in this regard. CONCLUSIONS: PRP and ASC affect lymphangiogenesis and lymphedema development and might represent a promising approach to improve regeneration of lymphatic vessels, restore disrupted lymphatic circulation and treat or prevent lymphedema alone or in combination with currently available lymphedema therapies.

Microtubule Associated Protein 1b (MAP1B) Is a Marker of the Microtubular Cytoskeleton in Podocytes but Is Not Essential for the Function of the Kidney Filtration Barrier in Mice.

Podocytes are essential for the function of the kidney glomerular filter. A highly differentiated cytoskeleton is requisite for their integrity. Although much knowledge has been gained on the organization of cortical actin networks in podocyte's foot processes, less is known about the molecular organization of the microtubular cytoskeleton in primary processes and the cell body. To gain an insight into the organization of the microtubular cytoskeleton of the podocyte, we systematically analyzed the expression of microtubule associated proteins (Maps), a family of microtubules interacting proteins with known functions as regulator, scaffold and guidance proteins. We identified microtubule associated protein 1b (MAP1B) to be specifically enriched in podocytes in human and rodent kidney. Using immunogold labeling in electron microscopy, we were able to demonstrate an enrichment of MAP1B in primary processes. A similar association of MAP1B with the microtubule cytoskeleton was detected in cultured podocytes. Subcellular distribution of MAP1B HC and LC1 was analyzed using a double fluorescent reporter MAP1B fusion protein. Subsequently we analyzed mice constitutively depleted of MAP1B. Interestingly, MAP1B KO was not associated with any functional or structural alterations pointing towards a redundancy of MAP proteins in podocytes. In summary, we established MAP1B as a specific marker protein of the podocyte microtubular cytoskeleton.

Macroporous Scaffolds for Bone Engineering. Studies on Cell Culture and Ectopic Bone Formation

Bone engineering is a rapidly developing area of reconstructive medicine where bone inducing factors and/or cells are combined with a scaffold material to regenerate the structure and function of the original tissue. The aim of this study was to compare the suitability of different macroporous scaffold types for bone engineering applications. The two scaffold categories studied were a) the mechanically strong and stable titanium fiber meshes and b) the elastic and biodegradable porous polymers. Furthermore, bioactive modifications were applied to these basic scaffold types, and their effect on the osteogenic responses was evaluated in cell culture and ectopic bone formation studies. The osteogenic phenotype of cultured cell-scaffold constructs was heightened with a sol-gel derived titania coating, but not with a mixed titania-silica coating. The latter coating also resulted in delayed ectopic bone formation in bone marrow stromal cell seeded scaffolds. However, the better bone contact in early implantation times and more even bone tissue distribution at later times indicated enhanced osteoconductivity of both the coated scaffold types. Overall, the most promising bone engineering results were obtained with titania coated fiber meshes. Elastic and biodegradable poly(ε-caprolactone/D,L-lactide) based scaffolds were also developed in this study. The degradation rates of the scaffolds in vitro were governed by the hydrophilicity of the polymer matrix, and the porous architecture was controlled by the amount and type of porogen used. A continuous phase macroporosity was obtained using a novel CaCl2 • 6H2O porogen. Dynamic culture conditions increased cell invasion, but decreased cell numbers and osteogenicity, within the scaffolds. Osteogenic differentiation in static cultures and ectopic bone formation in cell seeded scaffolds were enhanced in composites, with 30 wt-% of bioactive glass filler.

Bcl-2 family proteins and cytoskeleton changes involved in DM-1 cytotoxic effect on melanoma cells

Melanoma is one of the most aggressive types of skin cancer and its incidence rate is still increasing. All existing treatments are minimally effective. Consequently, new therapeutic agents for melanoma treatment should be developed. The DM-1 compound is a curcumin analog that possesses several curcumin characteristics, such as antiproliferative, antitumor, and anti-metastatic properties. The aim of this study was to evaluate the different signaling pathways involved in the cytotoxic effect of DM-1 on melanoma cells. The apoptotic process and cytoskeletal changes were evaluated by immunoblotting and immunofluorescence, respectively, in melanoma cells. After DM-1 treatment, SK-MEL-5 melanoma cells showed actin filament disorganization with spicule formation throughout the cytoskeleton and significant reduction of focal adhesion as well as they were present only at cell extremities, conferring a poor connection between the cell and the substrate. Besides this, there was significant filopodium retraction and loss of typical cytoskeleton scaffold. These modifications contributed to cell detachment followed by cell death. Furthermore, DM-1-induced apoptosis was triggered by multiple Bcl-2 proteins involved in both the extrinsic and the intrinsic apoptotic pathways. SK-MEL-5 cells showed a death mechanism mainly by Bcl-2/Bax ratio decrease, whereas A375 cells presented apoptosis induction by Mcl-1 and Bcl-xL downregulation. In SK-MEL-5 and A375 melanoma cells, there was a significant increase in the active form of caspase 9, and the inactive form of the effector caspase 3 was decreased in both cell lines. Expression of cleaved poly ADP ribose polymerase was increased after DM-1 treatment in these melanoma cell lines, demonstrating that the apoptotic process occurred. Altogether, these data elucidate the cellular and molecular mechanisms involved in the cytotoxicity induced by the antitumor agent DM-1 in melanoma cells.

NLRP12 is a neutrophil-specific, negative regulator of in vitro cell migration but does not modulate LPS- or infection-induced NF-κB or ERK signalling.

NOD-like receptors (NLR) are a family of cytosolic pattern recognition receptors that include many key drivers of innate immune responses. NLRP12 is an emerging member of the NLR family that is closely related to the well-known inflammasome scaffold, NLRP3. Since its discovery, various functions have been proposed for NLRP12, including the positive regulation of dendritic cell (DC) and neutrophil migration and the inhibition of NF-κB and ERK signalling in DC and macrophages. We show here that NLRP12 is poorly expressed in murine macrophages and DC, but is strongly expressed in neutrophils. Using myeloid cells from WT and Nlrp12(-/)(-) mice, we show that, contrary to previous reports, NLRP12 does not suppress LPS- or infection-induced NF-κB or ERK activation in myeloid cells, and is not required for DC migration in vitro. Surprisingly, we found that Nlrp12 deficiency caused increased rather than decreased neutrophil migration towards the chemokine CXCL1 and the neutrophil parasite Leishmania major, revealing NLRP12 as a negative regulator of directed neutrophil migration under these conditions.

Synthesis and structure - Activity relationship study of potent cytotoxic analogues of the marine alkaloid Lamellarin D

The marine alkaloid, Lamellarin D (Lam-D), has shown potent cytotoxicity in numerous cancer cell lines, and was recently identified as a potent topoisomerase I inhibitor. A library of open lactone analogs of Lam-D was prepared from a methyl 5,6-dihydropyrrolo[2,1-a]isoquinoline-3- carboxylate scaffold (1) by introducing various aryl groups through sequential and regioselective bromination, followed by Pd(0)-catalyzed Suzuki cross-coupling chemistry. The compounds were obtained in a 24-44% overall yield, and tested in a panel of three human tumor cell lines, MDA-MB- 231 (breast), A-549 (lung), and HT-29 (colon), to evaluate their cytotoxic potential. From these data the SAR study concluded that more than 75% of the open-chain Lam-D analogs tested showed cytotoxicity in a low micromolar GI50 range.

Synthesis and structure - Activity relationship study of potent cytotoxic analogues of the marine alkaloid Lamellarin D

The marine alkaloid, Lamellarin D (Lam-D), has shown potent cytotoxicity in numerous cancer cell lines, and was recently identified as a potent topoisomerase I inhibitor. A library of open lactone analogs of Lam-D was prepared from a methyl 5,6-dihydropyrrolo[2,1-a]isoquinoline-3- carboxylate scaffold (1) by introducing various aryl groups through sequential and regioselective bromination, followed by Pd(0)-catalyzed Suzuki cross-coupling chemistry. The compounds were obtained in a 24-44% overall yield, and tested in a panel of three human tumor cell lines, MDA-MB- 231 (breast), A-549 (lung), and HT-29 (colon), to evaluate their cytotoxic potential. From these data the SAR study concluded that more than 75% of the open-chain Lam-D analogs tested showed cytotoxicity in a low micromolar GI50 range.

Total synthesis of lamellarin D and analogs library

Larnellarins are a group of marine natural products isolated from the prosobranch mollusc Lamellaria sp., the ascidian Didemnum sp., and the sponge Dendrilla Cactos. Several of them exhibit interesting biological activities. Natural as well as synthetic lamellarins should be excellent candidates for the development of new drugs due to their unique skeletal structure and their important biological activities especially as antitumor agents. Lamelarin O has been recently characterized as a topoisomerase 1-targeted anti tumor agent. A variety of synthetic approaches have been developed for this family of alkaloids. Herein we describe a new route to the synthesis of Lamellarin D, from a methyl 2-pyrrolecarboxylate. Transformation of the starting material into the scaffold, a substituted 5,6-dihydropyrrolo (2,l ­a)isoquinoline (5,6-DHPl), was afforded by N-alkylation followed by intramolecular Heck cyclization. From this scaffold the synthetic strategy is based on two sequential regioselective bromination!Suzuki cross-coupling reactions which permitted the introduction of differently substituted aryl groups on positions 1 and 2 followed by oxidation, deprotection, and lactonization.

Building complexity in O2-binding copper complexes : site-selective metalation and intermolecular O2-binding at dicopper and heterometallic complexes derived from an unsymmetric ligand

A novel unsymmetric dinucleating ligand (LN3N4) combining a tridentate and a tetradentate binding sites linked through a m-xylyl spacer was synthesized as ligand scaffold for preparing homo- and dimetallic complexes, where the two metal ions are bound in two different coordination environments. Site-selective binding of different metal ions is demonstrated. LN3N4 is able to discriminate between CuI and a complementary metal (M′ = CuI, ZnII, FeII, CuII, or GaIII) so that pure heterodimetallic complexes with a general formula [CuIM′(LN3N4)]n+ are synthesized. Reaction of the dicopper(I) complex [CuI 2(LN3N4)]2+ with O2 leads to the formation of two different copper-dioxygen (Cu2O2) intermolecular species (O and TP) between two copper atoms located in the same site from different complex molecules. Taking advantage of this feature, reaction of the heterodimetallic complexes [CuM′(LN3N4)]n+ with O2 at low temperature is used as a tool to determine the final position of the CuI center in the system because only one of the two Cu2O2 species is formed

Cartilage Appearance Using an Environmental Scanning Electron Microscope.

Because of technical principles, samples to be observed with electron microscopy need to be fixed in a chemical process and exposed to vacuum conditions that can produce some changes in the morphology of the specimen. The aim of this work was to obtain high-resolution images of the fresh articular cartilage surface with an environmental scanning electron microscope (ESEM), which is an instrument that permits examination of biological specimens without fixation methods in a 10 Torr chamber pressure, thus minimizing the risk of creating artifacts in the structure. Samples from weight-bearing areas of femoral condyles of New Zealand white rabbits were collected and photographed using an ESEM. Images were analyzed using a categorization based in the Jurvelin classification system modified by Hong and Henderson. Appearance of the observed elevations and depressions as described in the classification were observed, but no fractures or splits of cartilage surface, thought to be artifacts, were detected. The ESEM is a useful tool to obtain images of fresh articular cartilage surface appearance without either employing fixation methods or exposing the specimen to extreme vacuum conditions, reducing the risk of introducing artifacts within the specimen. For all these reasons it could become a useful tool for quality control of the preservation process of osteochondral allografting in a bank of musculoskeletal tissues.

Study of ligand-based virtual screening tools in computer-aided drug design

Virtual screening is a central technique in drug discovery today. Millions of molecules can be tested in silico with the aim to only select the most promising and test them experimentally. The topic of this thesis is ligand-based virtual screening tools which take existing active molecules as starting point for finding new drug candidates. One goal of this thesis was to build a model that gives the probability that two molecules are biologically similar as function of one or more chemical similarity scores. Another important goal was to evaluate how well different ligand-based virtual screening tools are able to distinguish active molecules from inactives. One more criterion set for the virtual screening tools was their applicability in scaffold-hopping, i.e. finding new active chemotypes. In the first part of the work, a link was defined between the abstract chemical similarity score given by a screening tool and the probability that the two molecules are biologically similar. These results help to decide objectively which virtual screening hits to test experimentally. The work also resulted in a new type of data fusion method when using two or more tools. In the second part, five ligand-based virtual screening tools were evaluated and their performance was found to be generally poor. Three reasons for this were proposed: false negatives in the benchmark sets, active molecules that do not share the binding mode, and activity cliffs. In the third part of the study, a novel visualization and quantification method is presented for evaluation of the scaffold-hopping ability of virtual screening tools.