Making sense of AMPA receptor trafficking by modeling molecular mechanisms of synaptic plasticity.

Synaptic plasticity involves a complex molecular machinery with various protein interactions but it is not yet clear how its components give rise to the different aspects of synaptic plasticity. Here we ask whether it is possible to mathematically model synaptic plasticity by making use of known substances only. We present a model of a multistable biochemical reaction system and use it to simulate the plasticity of synaptic transmission in long-term potentiation (LTP) or long-term depression (LTD) after repeated excitation of the synapse. According to our model, we can distinguish between two phases: first, a "viscosity" phase after the first excitation, the effects of which like the activation of NMDA receptors and CaMKII fade out in the absence of further excitations. Second, a "plasticity" phase actuated by an identical subsequent excitation that follows after a short time interval and causes the temporarily altered concentrations of AMPA subunits in the postsynaptic membrane to be stabilized. We show that positive feedback is the crucial element in the core chemical reaction, i.e. the activation of the short-tail AMPA subunit by NEM-sensitive factor, which allows generating multiple stable equilibria. Three stable equilibria are related to LTP, LTD and a third unfixed state called ACTIVE. Our mathematical approach shows that modeling synaptic multistability is possible by making use of known substances like NMDA and AMPA receptors, NEM-sensitive factor, glutamate, CaMKII and brain-derived neurotrophic factor. Furthermore, we could show that the heteromeric combination of short- and long-tail AMPA receptor subunits fulfills the function of a memory tag.

IROme, a New High-Throughput Molecular Tool for the Diagnosis of Inherited Retinal Dystrophies-A Price Comparison with Sanger Sequencing.

The molecular diagnosis of retinal dystrophies (RD) is difficult because of genetic and clinical heterogeneity. Previously, the molecular screening of genes was done one by one, sometimes in a scheme based on the frequency of sequence variants and the number of exons/length of the candidate genes. Payment for these procedures was complicated and the sequential billing of several genes created endless paperwork. We therefore evaluated the costs of generating and sequencing a hybridization-based DNA library enriched for the 64 most frequently mutated genes in RD, called IROme, and compared them to the costs of amplifying and sequencing these genes by the Sanger method. The production cost generated by the high-throughput (HT) sequencing of IROme was established at CHF 2,875.75 per case. Sanger sequencing of the same exons cost CHF 69,399.02. Turnaround time of the analysis was 3 days for IROme. For Sanger sequencing, it could only be estimated, as we never sequenced all 64 genes in one single patient. Sale cost for IROme calculated on the basis of the sale cost of one exon by Sanger sequencing is CHF 8,445.88, which corresponds to the sale price of 40 exons. In conclusion, IROme is cheaper and faster than Sanger sequencing and therefore represents a sound approach for the diagnosis of RD, both scientifically and economically. As a drop in the costs of HT sequencing is anticipated, target resequencing might become the new gold standard in the molecular diagnosis of RD.

Molecular and cellular analysis of genodermatoses

Summary Skin is the essential interface between our body and its environment; not only does it prevent water loss and protect us from external insults it also plays an essential role in the central nervous system acting as a major sense organ primarily for touch and pain. The main cell type present in skin, keratinocyte, undergoes a differentiation process leading to the formation of this protecting barrier. This work is intended to contribute to the understanding of how keratinocyte differentiates and skin functions. To do this, we studied two genetic skin diseases: Erythrokeratodermia variabilis and Mal de Meleda. Our approach was to examine the expression and localization of proteins implicated in these two pathologies in normal and diseased tissues and to determine the influence of mutant proteins at the molecular and cellular levels. Connexins are major components of gap junctions, channels allowing direct communication between cells. Our laboratory has identified mutations in both connexin 30.3 (Cx30.3) and 31 (Cx31) to be causally involved in erythrokeratodermia variabilis (EKV), an autosomal dominant disorder of keratinization. In the first chapter, we show a new mutation of Cx31, L209P-Cx31, in 3 EKV patients, extending the field of EKV-causing mutations although the mechanism by which connexin mutations lead to the disease is unclear. In the second chapter, we studied the effect of F137L-Cx30.3 on expression, trafficking and localization of cotransfected Cx31 and Cx30.3 in connexin-deficient HeLa cells. The F137 amino acid, highly conserved in connexin family, is oriented towards the channel pore and F137L mutation in either Cx30.3 or Cx31 lead to EKV. As two genes can lead to EKV when mutated, our hypothesis was that Cx31 and Cx30.3 might cooperate at a molecular level. We were able to demonstrate a physical interaction between Cx31 and Cx30.3. The presence of F137L-Cx30.3 disturbed the trafficking of both connexins, less connexins were integrated into gap junctions and thus, the coupling between cell was diminished. Connexins formed in the presence of F137L-Cx30.3 are degraded at their exit from the endoplasmic reticulum. In conclusion, our results indicate that the genetic heterogeneity of EKV is due to mutations in two interacting proteins. F137L-Cx30.3 has a dominant negative effect and affects Cx31, disturbing cellular communication in epidermal cells. Mal de Meleda is an autosomal recessive inflammatory and a keratotic palmoplantar skin disorder due to mutations in SLURP1 (secreted LY6/PLAUR-related protein 1). SLURP1 belongs to the LY6/PLAUR family of proteins and has the particularity of being secreted instead of being GPI-anchored. The high degree of structural similarity between SLURP1 and the three fingers motif of snake neurotoxins and LYNX 1-C suggests that this protein could interact with the neuronal acetylcholine receptors. In the third chapter, we show that SLURP1 potentiates responses of the a7 nicotinic acetylcholine receptor (nAchR) to acetylcholine. These results identify SLURP1 as a secreted epidermal neuromodulator that is likely to be essential for palmoplantar skin. In the fourth chapter, we show that SLURP1 is expressed in the granular layer of the epidermis but is absent from skin biopsies of Mal de Meleda patients. SLURP1 is also present in secretions such as sweat, tears or saliva. An in vitro analysis on two mutant of SLURP-I demonstrates that W15R-SLURP1 is absent in cells while G86R-SLURP1 is expressed and secreted, suggesting that SLURP1 can lead to the disease by either an absent or an abnormal protein. Finally, in the fifth chapter, we analyse the expression and biological properties of other LY6/PLAUR members, clustered around SLURP] on chromosome 8. Their GPI-anchored or secreted status were analysed in vitro. SLURP1, LYNX1-A and -B are secreted while LYPDC2 and LYNX 1-C are GPI anchored. Three of these proteins are expressed in the epidermis and in cultured keratinocytes. These results suggest that these LY6/PLAUR members may have an important role in skin homeostasis. Résumé Résumé La peau est la barrière essentielle entre notre corps et l'environnement, nous protégeant des agressions extérieures, de la déshydratation et assurant aussi un rôle dans le système nerveux central en tant qu'organe du toucher et de la douleur. Le principal type de cellules présent dans la peau est le kératinocyte qui suit un processus de différenciation aboutissant à la formation de cette barrière protectrice. Ce travail est destiné à comprendre la différenciation des kératinocytes et le fonctionnement de la peau. Pour cela, nous avons étudié deux maladies génodermatoses : l'Erthrokeratodermia Variabilis (EKV) et le Mal de Meleda. Nous avons examiné l'expression et la localisation des protéines impliquées dans ces deux pathologies dans des tissus normaux et malades puis déterminé l'influence des protéines mutantes aux niveaux moléculaires et cellulaires. Les connexines (Cx) sont les composants majeurs des jonctions communicantes, canaux permettant la communication directe entre les cellules. Notre laboratoire a identifié des mutations dans les Cx30.3 et Cx31 comme responsables de l'EKV, génodermatose de transmission autosomique dominante. Dans le ler chapitre, nous décrivons une nouvelle mutation de Cx31, L209-Cx31, et contribuons à l'établissement du catalogue des mutations de Cx31 entraînant cette maladie. Cependant, le mécanisme par lequel les mutations de Cx31 et C3x0.3 provoquent l'EKV est inconnu. Dans le 2ème chapitre, nous étudions les effets de la mutation F137L-Cx30.3 sur l'expression, le trafic et la localisation des Cx31 et Cx30.3 transfectées dans des cellules HeLa, déficientes en connexines. Comme deux gènes peuvent causer une EKV quand ils sont mutés, notre hypothèse était que Cx31 et Cx30.3 pourraient coopérer au niveau moléculaire. Nous avons montré l'existence d'une interaction physique entre ces deux connexines. La présence de la mutation F137L-Cx30.3 perturbe le trafic des deux connexines, moins de connexines sont intégrées dans les jonctions communicantes et donc le couplage entre les cellules est diminué. Les connexons formés en présence de cette mutation sont dégradés à leur sortie du réticulum endoplasmique. En conclusion, nos résultats indiquent que l'hétérogénéité génétique de EKV est due à des mutations dans deux protéines qui interagissent. F137L-Cx30.3 a un effet dominant négatif et affecte Cx31, perturbant la communication entre les cellules épidermiques. Le Mal de Meleda est une maladie récessive de la peau palmoplantaire due à des mutations dans SLURP1. SLURP1 appartient à la famille des protéines contenant un domaine LY6/PLAUR et a la particularité d'être sécrétée. La grande homologie de structure existant entre SLURP1, les neurotoxines de serpent et LYNX1-C suggère que la protéine pourrait interagir avec des récepteurs à acétylcholine (Ach). Dans le 3ème chapitre, nous montrons que SLURP1 module la réponse à l'Ach du récepteur nicotinique α7. Ces résultats identifient SLURP1 comme un neuromodulateur épidermique sécrété, probablement essentiel pour la peau palmoplantaire. Dans le 4ème chapitre, nous montrons que SLURP1 est exprimé dans la couche granuleuse de l'épiderme et qu'il est absent des biopsies des patients. SLURP1 a aussi été détecté dans des sécrétions telles que la sueur, les lamies et la salive. Une analyse in vitro de deux mutants de SLURP1 a montré que W15R-SLURP1 est absent des cellules tandis que G86R-SLURP1 est exprimé et sécrété, suggérant qu'une absence ou une anomalie de SLURP1 peuvent causer la maladie. Finalement, dans le 5ème chapitre, nous analysons l'expression et les propriétés biologiques d'autres membres de la famille LY6/PLAUR localisés autour de SLURP1 sur le chromosome 8. Leur statut de protéines sécrétées ou liées à la membrane par une ancre GPI est analysé in vitro. SLURP1, LYNXI-A et -B sont sécrétées alors que LYPDC2 et LYNX1-C sont liés à la membrane. Trois de ces protéines sont exprimées dans l'épiderme et dans des kératinocytes cultivés. Ces résultats suggèrent que la famille LY6/PLAUR pourrait avoir un rôle important dans l'homéostasie de la peau.

Molecular modeling study of the binding and signaling properties of the TCRpMHC complex

The activation of the specific immune response against tumor cells is based on the recognition by the CD8+ Cytotoxic Τ Lymphocytes (CTL), of antigenic peptides (p) presented at the surface of the cell by the class I major histocompatibility complex (MHC). The ability of the so-called T-Cell Receptors (TCR) to discriminate between self and non-self peptides constitutes the most important specific control mechanism against infected cells. The TCR/pMHC interaction has been the subject of much attention in cancer therapy since the design of the adoptive transfer approach, in which Τ lymphocytes presenting an interesting response against tumor cells are extracted from the patient, expanded in vitro, and reinfused after immunodepletion, possibly leading to cancer regression. In the last decade, major progress has been achieved by the introduction of engineered lypmhocytes. In the meantime, the understanding of the molecular aspects of the TCRpMHC interaction has become essential to guide in vitro and in vivo studies. In 1996, the determination of the first structure of a TCRpMHC complex by X-ray crystallography revealed the molecular basis of the interaction. Since then, molecular modeling techniques have taken advantage of crystal structures to study the conformational space of the complex, and understand the specificity of the recognition of the pMHC by the TCR. In the meantime, experimental techniques used to determine the sequences of TCR that bind to a pMHC complex have been used intensively, leading to the collection of large repertoires of TCR sequences that are specific for a given pMHC. There is a growing need for computational approaches capable of predicting the molecular interactions that occur upon TCR/pMHC binding without relying on the time consuming resolution of a crystal structure. This work presents new approaches to analyze the molecular principles that govern the recognition of the pMHC by the TCR and the subsequent activation of the T-cell. We first introduce TCRep 3D, a new method to model and study the structural properties of TCR repertoires, based on homology and ab initio modeling. We discuss the methodology in details, and demonstrate that it outperforms state of the art modeling methods in predicting relevant TCR conformations. Two successful applications of TCRep 3D that supported experimental studies on TCR repertoires are presented. Second, we present a rigid body study of TCRpMHC complexes that gives a fair insight on the TCR approach towards pMHC. We show that the binding mode of the TCR is correctly described by long-distance interactions. Finally, the last section is dedicated to a detailed analysis of an experimental hydrogen exchange study, which suggests that some regions of the constant domain of the TCR are subject to conformational changes upon binding to the pMHC. We propose a hypothesis of the structural signaling of TCR molecules leading to the activation of the T-cell. It is based on the analysis of correlated motions in the TCRpMHC structure. - L'activation de la réponse immunitaire spécifique dirigée contre les cellules tumorales est basée sur la reconnaissance par les Lymphocytes Τ Cytotoxiques (CTL), d'un peptide antigénique (p) présenté à la suface de la cellule par le complexe majeur d'histocompatibilité de classe I (MHC). La capacité des récepteurs des lymphocytes (TCR) à distinguer les peptides endogènes des peptides étrangers constitue le mécanisme de contrôle le plus important dirigé contre les cellules infectées. L'interaction entre le TCR et le pMHC est le sujet de beaucoup d'attention dans la thérapie du cancer, depuis la conception de la méthode de transfer adoptif: les lymphocytes capables d'une réponse importante contre les cellules tumorales sont extraits du patient, amplifiés in vitro, et réintroduits après immunosuppression. Il peut en résulter une régression du cancer. Ces dix dernières années, d'importants progrès ont été réalisés grâce à l'introduction de lymphocytes modifiés par génie génétique. En parallèle, la compréhension du TCRpMHC au niveau moléculaire est donc devenue essentielle pour soutenir les études in vitro et in vivo. En 1996, l'obtention de la première structure du complexe TCRpMHC à l'aide de la cristallographie par rayons X a révélé les bases moléculaires de l'interaction. Depuis lors, les techniques de modélisation moléculaire ont exploité les structures expérimentales pour comprendre la spécificité de la reconnaissance du pMHC par le TCR. Dans le même temps, de nouvelles techniques expérimentales permettant de déterminer la séquence de TCR spécifiques envers un pMHC donné, ont été largement exploitées. Ainsi, d'importants répertoires de TCR sont devenus disponibles, et il est plus que jamais nécessaire de développer des approches informatiques capables de prédire les interactions moléculaires qui ont lieu lors de la liaison du TCR au pMHC, et ce sans dépendre systématiquement de la résolution d'une structure cristalline. Ce mémoire présente une nouvelle approche pour analyser les principes moléculaires régissant la reconnaissance du pMHC par le TCR, et l'activation du lymphocyte qui en résulte. Dans un premier temps, nous présentons TCRep 3D, une nouvelle méthode basée sur les modélisations par homologie et ab initio, pour l'étude de propriétés structurales des répertoires de TCR. Le procédé est discuté en détails et comparé à des approches standard. Nous démontrons ainsi que TCRep 3D est le plus performant pour prédire des conformations pertinentes du TCR. Deux applications à des études expérimentales des répertoires TCR sont ensuite présentées. Dans la seconde partie de ce travail nous présentons une étude de complexes TCRpMHC qui donne un aperçu intéressant du mécanisme d'approche du pMHC par le TCR. Finalement, la dernière section se concentre sur l'analyse détaillée d'une étude expérimentale basée sur les échanges deuterium/hydrogène, dont les résultats révèlent que certaines régions clés du domaine constant du TCR sont sujettes à un changement conformationnel lors de la liaison au pMHC. Nous proposons une hypothèse pour la signalisation structurelle des TCR, menant à l'activation du lymphocyte. Celle-ci est basée sur l'analyse des mouvements corrélés observés dans la structure du TCRpMHC.

iLOGP: A Simple, Robust, and Efficient Description of n-Octanol/Water Partition Coefficient for Drug Design Using the GB/SA Approach.

The n-octanol/water partition coefficient (log Po/w) is a key physicochemical parameter for drug discovery, design, and development. Here, we present a physics-based approach that shows a strong linear correlation between the computed solvation free energy in implicit solvents and the experimental log Po/w on a cleansed data set of more than 17,500 molecules. After internal validation by five-fold cross-validation and data randomization, the predictive power of the most interesting multiple linear model, based on two GB/SA parameters solely, was tested on two different external sets of molecules. On the Martel druglike test set, the predictive power of the best model (N = 706, r = 0.64, MAE = 1.18, and RMSE = 1.40) is similar to six well-established empirical methods. On the 17-drug test set, our model outperformed all compared empirical methodologies (N = 17, r = 0.94, MAE = 0.38, and RMSE = 0.52). The physical basis of our original GB/SA approach together with its predictive capacity, computational efficiency (1 to 2 s per molecule), and tridimensional molecular graphics capability lay the foundations for a promising predictor, the implicit log P method (iLOGP), to complement the portfolio of drug design tools developed and provided by the SIB Swiss Institute of Bioinformatics.

Evolutionary forces affect multigenomic arbuscular mycorhizal fungi : an experimental approach using Glomus intraradices

Abstract Arbuscular Mycorhizal Fungi (AMF) are important plant symbionts that can improve floristic diversity and ecosystem productivity. These important fungi are obligate biotrophs and form symbioses with roots of the majority of plant species, improving plant nutrient acquisition in exchange of photosynthates. AM fungi are successful both ecologically as they occupy a very large spectrum of environments as well as host range and evolutionarily, as this symbiosis is over 400 million years old. These fungi grow and reproduce clonally by hyphae and multinucleate spores. AMF are coenocytic and recent work has shown that they harbor genetically different nuclei and that AMF populations are genetically diverse. How AMF species diversity is maintained has been addressed theoretically and experimentally at the community level. Much less attention has been drawn to understand how genetic diversity is maintained within populations although closely related individuals are more likely to compete for the same resources and occupy similar niches. How infra-individual genetic diversity is shaped and maintained has received even less attention. In Chapter 2, we show that individuals from a field population may differ in their symbiotic efficiency under reduced phosphate availability: We show there is genetic variation in an AMF field population for fitness-related growth traits in response to different phosphate availability acid host species. Furthermore, AFLP fingerprints of the same individuals growing in contrasting environments diverged suggesting that the composition in nuclei of AMF is dynamical and affected by environmental factors. Thus environmental heterogeneity is likely to play an important role for the maintenance of genetic diversity at the population level. In Chapter 3 we show that single spores do not inherit necessarily the same genetic material. We have found genetic divergences using two different types of molecular marker, as well as phenotypic divergences among single spore lines. Our results stress the importance of considering these organisms as a multilevel hierarchical system and of better knowing their life cycle. They have important consequences for the understanding of AMF genetics, ecology and the development of commercial AMF inocculum. Résumé Les champignons endomycorhiziens arbusculaires (CEA) sont d'importants symbiontes pour les plantes, car ils augmentent la diversité et la productivité des écosystèmes. Ces importants symbiontes sont des biotrophes obligatoires et forment une symbiose avec la plupart des plantes terrestres. Ils améliorent l'acquisition de substances nutritives de leurs hôtes en échange de sucres obtenus par photosynthèse. Ces champignons ont un grand succès écologique, ils colonisent une grande rangée d'environnements ainsi que d'hôtes. Ils ont aussi un succès évolutif certain de part le fait que cette symbiose existe depuis plus de 400 millions d'années. Les CEA sont asexués et croissent clonalement en formant des hyphes et des spores multinuclées. Les CEA sont des coenocytes et des travaux de recherche récents ont montré qu'ils possèdent des noyaux génétiquement différents. D'autres travaux ont aussi révélé que les populations de CEA sont génétiquement diversifiées. Comment la diversité des CEA est maintenue a seulement été adressée par des études théoriques et expérimentalement au niveau des communautés. Très peu d'attention a été portée sur le maintien de la diversité génétique infra et inter populationnelle, or ce sont les individus les plus proches génétiquement qui vont entrer en compétition pour des ressources et niches similaires. La formation et le maintien de la diversité intra-individu des CEA a reçu très peu d'attention. Dans le chapitre 2, nous montrons que des individus CEA d'un même champ différent dans leur efficacité symbiotique lorsque la concentration en phosphoré est réduite. Nous montrons qu'il existe de la variance génétique dans une population de CEA provenant d'un même champ en réponse à différentes concentrations de phosphore, ainsi qu'en réponse à différentes espèces d'hôtes, et ceci pour des traits de croissance vraisemblablement liés au succès reproducteur. De plus grâce à des AFLP nous avons pu montrer que le génome de ces individus subissent des changements lorsqu'ils croissent dans des environnements contrastés. Ceci suggère que les noyaux génétiquement différents des CEA sont des entités dynamiques. Il est fort probable que l'hétérogénéité environnementale joue un rôle dans le maintien de la diversité génétique des populations de CEA. Dans le chapitre 3, nous montrons que toutes les spores d'un même mycélium parental de CEA ne reçoivent pas exactement le même contenu génétique. Nous avons mis en évidence des divergences entre des Lignées monosporales en utilisant deux types de marqueur moléculaires, ainsi que des différences phénotypiques. Nos résutats soulignent l'importance de considézer ces organismes comme dés systëmes hiérarchiques mufti-niveaux, ainsi que de mieux connaître leur cycle de vie. Nos résultats ont d'importantes conséquences pour la compréhension du système génétique des CEA, ainsi que de leur évolution, leur écológie, mais également des conséquences pour la production d' inoccultim commercial.

Two new species of Mycodrosophila (Diptera, Drosophilidae) proposed by molecular and morphological approaches, with a key to American species

ABSTRACT There are approximately 130 species of MycodrosophilaOldenberg, 1914 worldwide, although only nine species were recorded in American countries so far, three of which are exclusively Nearctic, five exclusively Neotropical and one found in both biogeographic regions (Mycodrosophila projectans). Such a small number of American species is likely a consequence of collecting bias, which favors the capture of frugivorous drosophilids, and to the general absence of Neotropical Mycodrosophila studies in the last 50 years. Here, we describe two commonly sampled species of Mycodrosophila from the Amazonian and Pampa Brazilian biomes, which share morphological similarities with Mycodrosophila neoprojectans and M. projectans, respectively. We compared sequences of the mitochondrial gene cytochrome oxidase subunit I (COI), external morphology characteristics and male terminalia among these species. Based on a DNA barcoding approach coupled to morphological differences, we proposed the delimitation of two new species, Mycodrosophila hofmanni sp. nov. and Mycodrosophila valentae sp. nov. An updated key to identifying Neotropical and Nearctic Mycodrosophila species is also provided.

Molecular determinants of desensitization in an ENaC/degenerin channel.

Epithelial Na(+) channel (ENaC)/degenerin family members are involved in mechanosensation, blood pressure control, pain sensation, and the expression of fear. Several of these channel types display a form of desensitization that allows the channel to limit Na(+) influx during prolonged stimulation. We used site-directed mutagenesis and chemical modification, functional analysis, and molecular dynamics simulations to investigate the role of the lower palm domain of the acid-sensing ion channel 1, a member of the ENaC/degenerin family. The lower palm domains of this trimeric channel are arranged around a central vestibule, at ∼20 Å above the plasma membrane and are covalently linked to the transmembrane channel parts. We show that the lower palm domains approach one another during desensitization. Residues in the palm co-determine the pH dependence of desensitization, its kinetics, and the stability of the desensitized state. Mutations of palm residues impair desensitization by preventing the closing movement of the palm. Overexpression of desensitization-impaired channel mutants in central neurons allowed--in contrast to overexpression of wild type--a sustained signaling response to rapid pH fluctuations. We identify and describe here the function of an important regulatory domain that most likely has a conserved role in ENaC/degenerin channels.

Pyrroloquinoline quinone biosynthesis gene pqqC, a novel molecular marker for studying the phylogeny and diversity of phosphate-solubilizing pseudomonads.

Many root-colonizing pseudomonads are able to promote plant growth by increasing phosphate availability in soil through solubilization of poorly soluble rock phosphates. The major mechanism of phosphate solubilization by pseudomonads is the secretion of gluconic acid, which requires the enzyme glucose dehydrogenase and its cofactor pyrroloquinoline quinone (PQQ). The main aim of this study was to evaluate whether a PQQ biosynthetic gene is suitable to study the phylogeny of phosphate-solubilizing pseudomonads. To this end, two new primers, which specifically amplify the pqqC gene of the Pseudomonas genus, were designed. pqqC fragments were amplified and sequenced from a Pseudomonas strain collection and from a natural wheat rhizosphere population using cultivation-dependent and cultivation-independent approaches. Phylogenetic trees based on pqqC sequences were compared to trees obtained with the two concatenated housekeeping genes rpoD and gyrB. For both pqqC and rpoD-gyrB, similar main phylogenetic clusters were found. However, in the pqqC but not in the rpoD-gyrB tree, the group of fluorescent pseudomonads producing the antifungal compounds 2,4-diacetylphloroglucinol and pyoluteorin was located outside the Pseudomonas fluorescens group. pqqC sequences from isolated pseudomonads were differently distributed among the identified phylogenetic groups than pqqC sequences derived from the cultivation-independent approach. Comparing pqqC phylogeny and phosphate solubilization activity, we identified one phylogenetic group with high solubilization activity. In summary, we demonstrate that the gene pqqC is a novel molecular marker that can be used complementary to housekeeping genes for studying the diversity and evolution of plant-beneficial pseudomonads.

New Findings in a Global Approach to Dissect the Whole Phenotype of PLA2G6 Gene Mutations.

Mutations in PLA2G6 gene have variable phenotypic outcome including infantile neuroaxonal dystrophy, atypical neuroaxonal dystrophy, idiopathic neurodegeneration with brain iron accumulation and Karak syndrome. The cause of this phenotypic variation is so far unknown which impairs both genetic diagnosis and appropriate family counseling. We report detailed clinical, electrophysiological, neuroimaging, histologic, biochemical and genetic characterization of 11 patients, from 6 consanguineous families, who were followed for a period of up to 17 years. Cerebellar atrophy was constant and the earliest feature of the disease preceding brain iron accumulation, leading to the provisional diagnosis of a recessive progressive ataxia in these patients. Ultrastructural characterization of patients' muscle biopsies revealed focal accumulation of granular and membranous material possibly resulting from defective membrane homeostasis caused by disrupted PLA2G6 function. Enzyme studies in one of these muscle biopsies provided evidence for a relatively low mitochondrial content, which is compatible with the structural mitochondrial alterations seen by electron microscopy. Genetic characterization of 11 patients led to the identification of six underlying PLA2G6 gene mutations, five of which are novel. Importantly, by combining clinical and genetic data we have observed that while the phenotype of neurodegeneration associated with PLA2G6 mutations is variable in this cohort of patients belonging to the same ethnic background, it is partially influenced by the genotype, considering the age at onset and the functional disability criteria. Molecular testing for PLA2G6 mutations is, therefore, indicated in childhood-onset ataxia syndromes, if neuroimaging shows cerebellar atrophy with or without evidence of iron accumulation.

Molecular modeling of peptide-MHC class I complexes : applications to peptide based immunotherapy

Summary The specific CD8+ T cell immune response against tumors relies on the recognition by the T cell receptor (TCR) on cytotoxic T lymphocytes (CTL) of antigenic peptides bound to the class I major histocompatibility complex (MHC) molecule. Such tumor associated antigenic peptides are the focus of tumor immunotherapy with peptide vaccines. The strategy for obtaining an improved immune response often involves the design of modified tumor associated antigenic peptides. Such modifications aim at creating higher affinity and/or degradation resistant peptides and require precise structures of the peptide-MHC class I complex. In addition, the modified peptide must be cross-recognized by CTLs specific for the parental peptide, i.e. preserve the structure of the epitope. Detailed structural information on the modified peptide in complex with MHC is necessary for such predictions. In this thesis, the main focus is the development of theoretical in silico methods for prediction of both structure and cross-reactivity of peptide-MHC class I complexes. Applications of these methods in the context of immunotherapy are also presented. First, a theoretical method for structure prediction of peptide-MHC class I complexes is developed and validated. The approach is based on a molecular dynamics protocol to sample the conformational space of the peptide in its MHC environment. The sampled conformers are evaluated using conformational free energy calculations. The method, which is evaluated for its ability to reproduce 41 X-ray crystallographic structures of different peptide-MHC class I complexes, shows an overall prediction success of 83%. Importantly, in the clinically highly relevant subset of peptide-HLAA*0201 complexes, the prediction success is 100%. Based on these structure predictions, a theoretical approach for prediction of cross-reactivity is developed and validated. This method involves the generation of quantitative structure-activity relationships using three-dimensional molecular descriptors and a genetic neural network. The generated relationships are highly predictive as proved by high cross-validated correlation coefficients (0.78-0.79). Together, the here developed theoretical methods open the door for efficient rational design of improved peptides to be used in immunotherapy. Résumé La réponse immunitaire spécifique contre des tumeurs dépend de la reconnaissance par les récepteurs des cellules T CD8+ de peptides antigéniques présentés par les complexes majeurs d'histocompatibilité (CMH) de classe I. Ces peptides sont utilisés comme cible dans l'immunothérapie par vaccins peptidiques. Afin d'augmenter la réponse immunitaire, les peptides sont modifiés de façon à améliorer l'affinité et/ou la résistance à la dégradation. Ceci nécessite de connaître la structure tridimensionnelle des complexes peptide-CMH. De plus, les peptides modifiés doivent être reconnus par des cellules T spécifiques du peptide natif. La structure de l'épitope doit donc être préservée et des structures détaillées des complexes peptide-CMH sont nécessaires. Dans cette thèse, le thème central est le développement des méthodes computationnelles de prédiction des structures des complexes peptide-CMH classe I et de la reconnaissance croisée. Des applications de ces méthodes de prédiction à l'immunothérapie sont également présentées. Premièrement, une méthode théorique de prédiction des structures des complexes peptide-CMH classe I est développée et validée. Cette méthode est basée sur un échantillonnage de l'espace conformationnel du peptide dans le contexte du récepteur CMH classe I par dynamique moléculaire. Les conformations sont évaluées par leurs énergies libres conformationnelles. La méthode est validée par sa capacité à reproduire 41 structures des complexes peptide-CMH classe I obtenues par cristallographie aux rayons X. Le succès prédictif général est de 83%. Pour le sous-groupe HLA-A*0201 de complexes de grande importance pour l'immunothérapie, ce succès est de 100%. Deuxièmement, à partir de ces structures prédites in silico, une méthode théorique de prédiction de la reconnaissance croisée est développée et validée. Celle-ci consiste à générer des relations structure-activité quantitatives en utilisant des descripteurs moléculaires tridimensionnels et un réseau de neurones couplé à un algorithme génétique. Les relations générées montrent une capacité de prédiction remarquable avec des valeurs de coefficients de corrélation de validation croisée élevées (0.78-0.79). Les méthodes théoriques développées dans le cadre de cette thèse ouvrent la voie du design de vaccins peptidiques améliorés.

Molecular and Functional Characterization of Neuroblastoma-Initiating Cells

Neuroblastoma (NB) is the most common extracranial malignant tumor in young children and arises at any site of the sympathetic nervous system. The disease exhibits a remarkable phenotypic diversity ranging from spontaneous regression to fatal disease. Poor outcome results from a rapidly progressive, metastatic and drug-resistant disease. Recent studies have suggested that solid tumors may arise from a minor population of cancer stem cells (CSCs) with stem cell markers and typical properties such as self-renewal ability, asymmetric division and drug resistance. In this model, CSCs possess the exclusive ability to initiate and maintain the tumor, and to produce distant metastases. Tumor cell subpopulations with stem-like phenotypes have indeed been identified in several cancer including leukemia, breast, brain and colon cancers. CSC hypothesis still needs to be validated in the other cancers including NB.NB originates from neural crest-derived malignant sympatho-adrenal cells. We have identified rare cells that express markers in conformity with neural crest stem cells and their derived lineages within primary NB tissue and cell lines, leading us to postulate the existence of CSCs in NB tumors.In the absence of specific markers to isolate CSCs, we adapted to NB tumor cells the sphere functional assay, based on the ability of stem cells to grow as spheres in non-adherent conditions. By serial passages of spheres from bone marrow NB metastases, a subset of cells was gradually selected and its specific gene expression profile identified by micro-array time-course analysis. The differentially expressed genes in spheres are enriched in genes implicated in development including CD133, ABC-transporters, WNT and NOTCH genes, identified in others solid cancers as CSCs markers, and other new markers, all referred by us as the Neurosphere Expression Profile (NEP). We confirmed the presence of a cell subpopulation expressing a combination of the NEP markers within a few primary NB samples.The tumorigenic potential of NB spheres was assayed by in vivo tumor growth analyses using orthotopic (adrenal glands) implantations of tumor cells into immune-compromised mice. Tumors derived from the sphere cells were significantly more frequent and were detected earlier compared to whole tumor cells. However, NB cells expressing the neurosphere-associated genes and isolated from the bulk tumors did not recapitulate the CSC-like phenotype in the orthotopic model. In addition, the NB sphere cells lost their higher tumorigenic potential when implanted in a subcutaneous heterotopic in vivo model.These results highlighted the complex behavior of CSC functions and led us to consider the stem-like NB cells as a dynamic and heterogeneous cell population influenced by microenvironment signals.Our approach identified for the first time candidate genes that may be associated with NB self-renewal and tumorigenicity and therefore would establish specific functional targets for more effective therapies in aggressive NB.

Intraspecific genetic variation in arbuscular mycorrhizal fungi and its consequences for molecular biology, ecology, and development of inoculum

It has been known for some time that different arbuscular mycorrhizal fungal (AMF) taxa confer differences in plant growth. Although genetic variation within AMF species has been given less attention, it could potentially be an ecologically important source of variation. Ongoing studies on variability in AMF genes within Glomus intraradices indicate that at least for some genes, such as the BiP gene, sequence variability can be high, even in coding regions. This suggests that genetic variation within an AMF may not be selectively neutral. This clearly needs to be investigated in more detail for other coding regions of AMF genomes. Similarly, studies on AMF population genetics indicate high genetic variation in AMF populations, and a considerable amount of variation seen in phenotypes in the population can be attributed to genetic differences among the fungi. The existence of high within-species genetic variation could have important consequences for how investigations on AMF gene expression and function are conducted. Furthermore, studies of within-species genetic variability and how it affects variation in plant growth will help to identify at what level of precision ecological studies should be conducted to identify AMF in plant roots in the field. A population genetic approach to studying AMF genetic variability can also be useful for inoculum development. By knowing the amount of genetic variability in an AMF population, the maximum and minimum numbers of spores that will contain a given amount of genetic diversity can be estimated. This could be particularly useful for developing inoculum with high adaptability to different environments.

Surface collective oscillations of metal clusters and spheres: Random-phase-approximation sum-rules approach

Using the once and thrice energy-weighted moments of the random-phase-approximation strength function, we have derived compact expressions for the average energy of surface collective oscillations of clusters and spheres of metal atoms. The L=0 volume mode has also been studied. We have carried out quantal and semiclassical calculations for Na and Ag systems in the spherical-jellium approximation. We present a rather thorough discussion of surface diffuseness and quantal size effects on the resonance energies.

Cheminformatic studies of molecular properties and their influence on bioactivity in the context of drug discovery

Molecular shape has long been known to be an important property for the process of molecular recognition. Previous studies postulated the existence of a drug-like shape space that could be used to artificially bias the composition of screening libraries, with the aim to increase the chance of success in Hit Identification. In this work, it was analysed to which extend this assumption holds true. Normalized Principal Moments of Inertia Ratios (NPRs) have been used to describe the molecular shape of small molecules. It was investigated, whether active molecules of diverse targets are located in preferred subspaces of the NPR shape space. Results illustrated a significantly stronger clustering than could be expected by chance, with parts of the space unlikely to be occupied by active compounds. Furthermore, a strong enrichment of elongated, rather flat shapes could be observed, while globular compounds were highly underrepresented. This was confirmed for a wide range of small molecule datasets from different origins. Active compounds exhibited a high overlap in their shape distributions across different targets, making a purely shape­ based discrimination very difficult. An additional perspective was provided by comparing the shapes of protein binding pockets with those of their respective ligands. Although more globular than their ligands, it was observed that binding sites shapes exhibited a similarly skewed distribution in shape space: spherical shapes were highly underrepresented. This was different for unoccupied binding pockets of smaller size. These were on the contrary identified to possess a more globular shape. The relation between shape complementarity and exhibited bioactivity was analysed; a moderate correlation between bioactivity and parameters including pocket coverage, distance in shape space, and others could be identified, which reflects the importance of shape complementarity. However, this also suggests that other aspects are of relevance for molecular recognition. A subsequent analysis assessed if and how shape and volume information retrieved from pocket or respective reference ligands could be used as a pre-filter in a virtual screening approach. ln Lead Optimization compounds need to get optimized with respect to a variety of pararneters. Here, the availability of past success stories is very valuable, as they can guide medicinal chemists during their analogue synthesis plans. However, although of tremendous interest for the public domain, so far only large corporations had the ability to mine historical knowledge in their proprietary databases. With the aim to provide such information, the SwissBioisostere database was developed and released during this thesis. This database contains information on 21,293,355 performed substructural exchanges, corresponding to 5,586,462 unique replacements that have been measured in 35,039 assays against 1,948 molecular targets representing 30 target classes, and on their impact on bioactivity . A user-friendly interface was developed that provides facile access to these data and is accessible at http//www.swissbioisostere.ch. The ChEMBL database was used as primary data source of bioactivity information. Matched molecular pairs have been identified in the extracted and cleaned data. Success-based scores were developed and integrated into the database to allow re-ranking of proposed replacements by their past outcomes. It was analysed to which degree these scores correlate with chemical similarity of the underlying fragments. An unexpectedly weak relationship was detected and further investigated. Use cases of this database were envisioned, and functionalities implemented accordingly: replacement outcomes are aggregatable at the assay level, and it was shawn that an aggregation at the target or target class level could also be performed, but should be accompanied by a careful case-by-case assessment. It was furthermore observed that replacement success depends on the activity of the starting compound A within a matched molecular pair A-B. With increasing potency the probability to lose bioactivity through any substructural exchange was significantly higher than in low affine binders. A potential existence of a publication bias could be refuted. Furthermore, often performed medicinal chemistry strategies for structure-activity-relationship exploration were analysed using the acquired data. Finally, data originating from pharmaceutical companies were compared with those reported in the literature. It could be seen that industrial medicinal chemistry can access replacement information not available in the public domain. In contrast, a large amount of often-performed replacements within companies could also be identified in literature data. Preferences for particular replacements differed between these two sources. The value of combining different endpoints in an evaluation of molecular replacements was investigated. The performed studies highlighted furthermore that there seem to exist no universal substructural replacement that always retains bioactivity irrespective of the biological environment. A generalization of bioisosteric replacements seems therefore not possible. - La forme tridimensionnelle des molécules a depuis longtemps été reconnue comme une propriété importante pour le processus de reconnaissance moléculaire. Des études antérieures ont postulé que les médicaments occupent préférentiellement un sous-ensemble de l'espace des formes des molécules. Ce sous-ensemble pourrait être utilisé pour biaiser la composition de chimiothèques à cribler, dans le but d'augmenter les chances d'identifier des Hits. L'analyse et la validation de cette assertion fait l'objet de cette première partie. Les Ratios de Moments Principaux d'Inertie Normalisés (RPN) ont été utilisés pour décrire la forme tridimensionnelle de petites molécules de type médicament. Il a été étudié si les molécules actives sur des cibles différentes se co-localisaient dans des sous-espaces privilégiés de l'espace des formes. Les résultats montrent des regroupements de molécules incompatibles avec une répartition aléatoire, avec certaines parties de l'espace peu susceptibles d'être occupées par des composés actifs. Par ailleurs, un fort enrichissement en formes allongées et plutôt plates a pu être observé, tandis que les composés globulaires étaient fortement sous-représentés. Cela a été confirmé pour un large ensemble de compilations de molécules d'origines différentes. Les distributions de forme des molécules actives sur des cibles différentes se recoupent largement, rendant une discrimination fondée uniquement sur la forme très difficile. Une perspective supplémentaire a été ajoutée par la comparaison des formes des ligands avec celles de leurs sites de liaison (poches) dans leurs protéines respectives. Bien que plus globulaires que leurs ligands, il a été observé que les formes des poches présentent une distribution dans l'espace des formes avec le même type d'asymétrie que celle observée pour les ligands: les formes sphériques sont fortement sous­ représentées. Un résultat différent a été obtenu pour les poches de plus petite taille et cristallisées sans ligand: elles possédaient une forme plus globulaire. La relation entre complémentarité de forme et bioactivité a été également analysée; une corrélation modérée entre bioactivité et des paramètres tels que remplissage de poche, distance dans l'espace des formes, ainsi que d'autres, a pu être identifiée. Ceci reflète l'importance de la complémentarité des formes, mais aussi l'implication d'autres facteurs. Une analyse ultérieure a évalué si et comment la forme et le volume d'une poche ou de ses ligands de référence pouvaient être utilisés comme un pré-filtre dans une approche de criblage virtuel. Durant l'optimisation d'un Lead, de nombreux paramètres doivent être optimisés simultanément. Dans ce contexte, la disponibilité d'exemples d'optimisations réussies est précieuse, car ils peuvent orienter les chimistes médicinaux dans leurs plans de synthèse par analogie. Cependant, bien que d'un extrême intérêt pour les chercheurs dans le domaine public, seules les grandes sociétés pharmaceutiques avaient jusqu'à présent la capacité d'exploiter de telles connaissances au sein de leurs bases de données internes. Dans le but de remédier à cette limitation, la base de données SwissBioisostere a été élaborée et publiée dans le domaine public au cours de cette thèse. Cette base de données contient des informations sur 21 293 355 échanges sous-structuraux observés, correspondant à 5 586 462 remplacements uniques mesurés dans 35 039 tests contre 1948 cibles représentant 30 familles, ainsi que sur leur impact sur la bioactivité. Une interface a été développée pour permettre un accès facile à ces données, accessible à http:/ /www.swissbioisostere.ch. La base de données ChEMBL a été utilisée comme source de données de bioactivité. Une version modifiée de l'algorithme de Hussain et Rea a été implémentée pour identifier les Matched Molecular Pairs (MMP) dans les données préparées au préalable. Des scores de succès ont été développés et intégrés dans la base de données pour permettre un reclassement des remplacements proposés selon leurs résultats précédemment observés. La corrélation entre ces scores et la similarité chimique des fragments correspondants a été étudiée. Une corrélation plus faible qu'attendue a été détectée et analysée. Différents cas d'utilisation de cette base de données ont été envisagés, et les fonctionnalités correspondantes implémentées: l'agrégation des résultats de remplacement est effectuée au niveau de chaque test, et il a été montré qu'elle pourrait également être effectuée au niveau de la cible ou de la classe de cible, sous réserve d'une analyse au cas par cas. Il a en outre été constaté que le succès d'un remplacement dépend de l'activité du composé A au sein d'une paire A-B. Il a été montré que la probabilité de perdre la bioactivité à la suite d'un remplacement moléculaire quelconque est plus importante au sein des molécules les plus actives que chez les molécules de plus faible activité. L'existence potentielle d'un biais lié au processus de publication par articles a pu être réfutée. En outre, les stratégies fréquentes de chimie médicinale pour l'exploration des relations structure-activité ont été analysées à l'aide des données acquises. Enfin, les données provenant des compagnies pharmaceutiques ont été comparées à celles reportées dans la littérature. Il a pu être constaté que les chimistes médicinaux dans l'industrie peuvent accéder à des remplacements qui ne sont pas disponibles dans le domaine public. Par contre, un grand nombre de remplacements fréquemment observés dans les données de l'industrie ont également pu être identifiés dans les données de la littérature. Les préférences pour certains remplacements particuliers diffèrent entre ces deux sources. L'intérêt d'évaluer les remplacements moléculaires simultanément selon plusieurs paramètres (bioactivité et stabilité métabolique par ex.) a aussi été étudié. Les études réalisées ont souligné qu'il semble n'exister aucun remplacement sous-structural universel qui conserve toujours la bioactivité quel que soit le contexte biologique. Une généralisation des remplacements bioisostériques ne semble donc pas possible.