Rôle des facteurs physico-chimiques du micro-environnement intestinal et des boucles inter-hélicales du Domaine I dans l’activité de la toxine insecticide Cry9Ca du bacille de Thuringe

Une fois ingérées par un insecte sensible, les toxines insecticides du bacille de Thuringe doivent être activées par les protéases intestinales de cet insecte. Leur premier domaine, un ensemble de sept hélices-α amphipathiques, est responsable de leur insertion dans la membrane luminale de certaines cellules de l’intestin médian, ce qui crée des pores peu sélectifs. La toxicité et la capacité à former des pores d’une telle toxine, la Cry9Ca, de ses mutants simples R164A et R164K et d’un fragment de 55 kDa résultant d’un clivage protéolytique au niveau de son résidu 164 ont été étudiées à l’aide d’une combinaison de modélisation par homologie, de bioessais, d’expériences de gonflement osmotique avec des vésicules de membrane en bordure en brosse de larves de sphinx du tabac et de mesures électrophysiologiques sur des intestins isolés. Ni les mutations simples ni le clivage protéolytique n’ont altéré la toxicité de la Cry9Ca. Dans une solution à faible force ionique, toutefois, la formation des pores dépend fortement du pH : une augmentation de celui-ci de 6,5 à 10,5 a entraîné une baisse irrégulière et par étapes successives de la perméabilité membranaire. Les quatre préparations de toxine ont néanmoins dépolarisé la membrane apicale d’intestins médians fraîchement isolés baignant dans une solution contenant 122 mM de KCl à pH 10,5. L’activité de la Cry9Ca, et des mutants R164A et R164K, a été grandement stimulée lorsque les expériences ont été effectuées en présence de suc intestinal, de lipides extraits d’un volume équivalent de suc intestinal ou d’un cocktail d’inhibiteurs de protéases solubles dans l’eau. De plus, le rôle des boucles inter-hélicales du Domaine I lors de l’insertion dans la membrane a été étudié avec des mutants doubles de la Cry9Ca dont les mutations introduisaient, neutralisaient ou renversaient une charge électrique. À l’exception de trois d’entres eux, tous ces mutants ont conservé une toxicité et une capacité à former des pores comparables à celles de la toxine parentale. L’ensemble de ces résultats suggère que le micro-environnement de l’intestin médian contribue à minimiser l’influence des charges de surface portées par les résidus des boucles inter-hélicales du Domaine I sur la capacité des toxines du bacille de Thuringe à former des pores. Il indique aussi que, d’une part, selon le site de clivage et les conditions expérimentales utilisées, des protéolyses supplémentaires de la toxine Cry9Ca activée peuvent soit stimuler, soit nuire à son activité et que, d’autre part, le suc intestinal du sphinx du tabac contient probablement un inhibiteur de protéases qui pourrait jouer un rôle important dans l’activité des toxines du bacille de Thuringe.

Calmodulin/KCa3.1 channel interactions as determinant to the KCa3.1 Ca2+ dependent gating : theoretical and experimental analyses

Differentes études ont montré que la sensibilité au Ca2+ du canal KCa3.1, un canal potassique indépendant du voltage, était conférée par la protéine calmoduline (CaM) liée de façon constitutive au canal. Cette liaison impliquerait la région C-lobe de la CaM et un domaine de $\ikca$ directement relié au segment transmembranaire S6 du canal. La CaM pourrait égalment se lier au canal de façon Ca2+ dépendante via une interaction entre un domaine de KCa3.1 du C-terminal (CaMBD2) et la région N-lobe de la CaM. Une étude fut entreprise afin de déterminer la nature des résidus responsables de la liaison entre le domaine CaMBD2 de KCa3.1 et la région N-lobe de la CaM et leur rôle dans le processus d'ouverture du canal par le Ca2+. Une structure 3D du complexe KCa3.1/CaM a d'abord été générée par modélisation par homologie avec le logiciel MODELLER en utilisant comme référence la structure cristalline du complexe SK2.2/CaM (PDB: 1G4Y). Le modèle ainsi obtenu de KCa3.1 plus CaM prévoit que le segment L361-S372 dans KCa3.1 devrait être responsable de la liaison dépendante du Ca2+ du canal avec la région N-lobe de la CaM via les résidus L361 et Q364 de KCa3.1 et E45, E47 et D50 de la CaM. Pour tester ce modèle, les résidus dans le segment L361-S372 ont été mutés en Cys et l'action du MTSET+ (chargé positivement) et MTSACE (neutre) a été mesurée sur l'activité du canal. Des enregistrements en patch clamp en configuration ``inside-out`` ont montré que la liaison du réactif chargé MTSET+ au le mutant Q364C entraîne une forte augmentation du courant, un effet non observé avec le MTSACE. De plus les mutations E45A et E47A dans la CaM, ont empêché l'augmentation du courant initié par MTSET+ sur le mutant Q364C. Une analyse en canal unitaire a confirmé que la liaison MTSET+ à Q364C cause une augmentation de la probabilité d'ouverture de KCa3.1 par une déstabilisation de l'état fermé du canal. Nous concluons que nos résultats sont compatibles avec la formation de liaisons ioniques entre les complexes chargés positivement Cys-MTSET+ à la position 364 de KCa3.1 et les résidus chargés négativement E45 et E47 dans la CaM. Ces données confirment qu'une stabilisation électrostatique des interactions CaM/KCa3.1 peut conduire à une augmentation de la probabilité d'ouverture du canal en conditions de concentrations saturantes de Ca2+.

Modulation de l’adressage membranaire et de la fonction du canal CaV2.3 par les résidus leucine du domaine guanylate kinase impliqués dans la liaison à forte affinité de CaVβ

Les canaux Ca2+ activés par le voltage (CaV) sont des protéines membranaires qui génèrent des courants Ca2+ dans les cellules excitables suite à une dépolarisation membranaire. Ces complexes oligomériques sont classifiés selon les propriétés structurelles de la sous-unité principale qui forme le pore du canal, soit la sous-unité CaVα1. La sous-unité auxiliaire CaVβ module l’expression membranaire et la dépendance au voltage du « gating » de la sous-unité CaVα1 des canaux HVA (« high-voltage-activated ») CaV1 et CaV2. La sous-unité CaVβ est formée par un domaine SH3 (« Src homology-3 ») connecté à un domaine GK (« guanylate kinase-like ») par le biais d’un domaine variable HOOK. Dans le but d’identifier les résidus dans la CaVβ3 qui sont responsables de la densité membranaire du CaV2.3, nous avons produit des mutants de la sous-unité auxiliaire le long de ses domaines fonctionnels. Cela dit, la délétion complète du domaine SH3 ainsi que la délétion du domaine HOOK n’ont pas modifié la densité membranaire de CaV2.3 ni ses propriétés d’activation. Cependant, la délétion de cinq résidus dans le domaine GK interrompt l’expression membranaire et l’expression fonctionnelle de CaV2.3. La mutation de résidus identifiés précédemment comme soutenant une affinité de liaison de l’ordre du nanomolaire dans le domaine GK de CaVβ n’a pas modifié de manière significative l’adressage membranaire de CaV2.3. Toutefois, les mutations de quatre résidus leucine dans les régions α3, α6, β10 et α9 du domaine GK ont grandement réduit l’adressage membranaire du canal CaV2.3. Nos résultats confirment que le domaine GK contient les déterminants moléculaires responsables de la fonction chaperone de CaVβ. Cela dit, l’adressage membranaire induit par CaVβ semble être déterminé par des éléments structuraux qui ne sont pas strictement dépendants d’une liaison à haute affinité de CaVβ sur CaVα1.

Propriétés électrophysiologiques des canaux ioniques formés par la toxine nématicide Cry5Ba du bacille de Thuringe dans les bicouches lipidiques planes

Les toxines Cry sont des protéines synthétisées sous forme de cristaux par la bactérie bacille de Thuringe pendant la sporulation. Elles sont largement utilisées comme agents de lutte biologique, car elles sont toxiques envers plusieurs espèces d’invertébrées, y compris les nématodes. Les toxines Cry5B sont actives contre certaines espèces de nématodes parasites, y compris Ankylostoma ceylanicum un parasite qui infeste le système gastro-intestinal des humains. Jusqu’au présent, le mode d’action des toxines Cry nématicides reste grandement inconnu, sauf que leurs récepteurs spécifiques sont des glycolipides et qu’elles causent des dommages importants aux cellules intestinales. Dans cette étude, on démontre pour la première fois que la toxine nématicide Cry5Ba, membre de la famille des toxines à trois domaines et produite par la bactérie bacille de Thuringe, forme des pores dans les bicouches lipidiques planes en absence de récepteurs. Les pores formés par cette toxine sont de sélectivité cationique, à pH acide ou alcalin. Les conductances des pores formés sous conditions symétriques de 150 mM de KCl varient entre 17 et 330 pS, à pH 6.0 et 9.0. Les niveaux des conductances les plus fréquemment observés diffèrent les uns des autres par environ 17 à 18 pS, ce qui est compatible avec l’existence d’arrangement d’un nombre différent de pores élémentaires similaires, activés de façon synchronisée, ou avec la présence d’oligomères de tailles variables et de différents diamètres de pores.

An Eigenfunction Approach for Volatility Modeling.

In this paper, we introduce a new approach for volatility modeling in discrete and continuous time. We follow the stochastic volatility literature by assuming that the variance is a function of a state variable. However, instead of assuming that the loading function is ad hoc (e.g., exponential or affine), we assume that it is a linear combination of the eigenfunctions of the conditional expectation (resp. infinitesimal generator) operator associated to the state variable in discrete (resp. continuous) time. Special examples are the popular log-normal and square-root models where the eigenfunctions are the Hermite and Laguerre polynomials respectively. The eigenfunction approach has at least six advantages: i) it is general since any square integrable function may be written as a linear combination of the eigenfunctions; ii) the orthogonality of the eigenfunctions leads to the traditional interpretations of the linear principal components analysis; iii) the implied dynamics of the variance and squared return processes are ARMA and, hence, simple for forecasting and inference purposes; (iv) more importantly, this generates fat tails for the variance and returns processes; v) in contrast to popular models, the variance of the variance is a flexible function of the variance; vi) these models are closed under temporal aggregation.

Probe-level linear model fitting and mixture modeling results in high accuracy detection of differential gene expression

Affiliation: Institut de recherche en immunologie et en cancérologie, Université de Montréal

MODELING HETEROTACHY IN PHYLOGENETICS

Il a été démontré que l’hétérotachie, variation du taux de substitutions au cours du temps et entre les sites, est un phénomène fréquent au sein de données réelles. Échouer à modéliser l’hétérotachie peut potentiellement causer des artéfacts phylogénétiques. Actuellement, plusieurs modèles traitent l’hétérotachie : le modèle à mélange des longueurs de branche (MLB) ainsi que diverses formes du modèle covarion. Dans ce projet, notre but est de trouver un modèle qui prenne efficacement en compte les signaux hétérotaches présents dans les données, et ainsi améliorer l’inférence phylogénétique. Pour parvenir à nos fins, deux études ont été réalisées. Dans la première, nous comparons le modèle MLB avec le modèle covarion et le modèle homogène grâce aux test AIC et BIC, ainsi que par validation croisée. A partir de nos résultats, nous pouvons conclure que le modèle MLB n’est pas nécessaire pour les sites dont les longueurs de branche diffèrent sur l’ensemble de l’arbre, car, dans les données réelles, le signaux hétérotaches qui interfèrent avec l’inférence phylogénétique sont généralement concentrés dans une zone limitée de l’arbre. Dans la seconde étude, nous relaxons l’hypothèse que le modèle covarion est homogène entre les sites, et développons un modèle à mélanges basé sur un processus de Dirichlet. Afin d’évaluer différents modèles hétérogènes, nous définissons plusieurs tests de non-conformité par échantillonnage postérieur prédictif pour étudier divers aspects de l’évolution moléculaire à partir de cartographies stochastiques. Ces tests montrent que le modèle à mélanges covarion utilisé avec une loi gamma est capable de refléter adéquatement les variations de substitutions tant à l’intérieur d’un site qu’entre les sites. Notre recherche permet de décrire de façon détaillée l’hétérotachie dans des données réelles et donne des pistes à suivre pour de futurs modèles hétérotaches. Les tests de non conformité par échantillonnage postérieur prédictif fournissent des outils de diagnostic pour évaluer les modèles en détails. De plus, nos deux études révèlent la non spécificité des modèles hétérogènes et, en conséquence, la présence d’interactions entre différents modèles hétérogènes. Nos études suggèrent fortement que les données contiennent différents caractères hétérogènes qui devraient être pris en compte simultanément dans les analyses phylogénétiques.

Mécanismes moléculaires d’activation du récepteur A des peptides natriurétiques

Le récepteur A des peptides natriurétiques (NPRA) fait partie de la famille des guanylates cyclases membranaires. L’activation du NPRA par ses agonistes naturels, ANP et BNP, induit une production de GMPc qui est responsable de leur rôle dans l’homéostasie cardiovasculaire, l’inhibition de l’hypertrophie et de la fibrose cardiaques et la régulation de la lipolyse. Le NPRA est un homodimère non covalent composé d’un domaine extracellulaire de liaison du ligand (ECD), d’un unique domaine transmembranaire (TM), d’un domaine d’homologie aux kinases et d’un domaine guanylate cyclase. Bien que le NPRA ait un rôle physiologique important, les mécanismes moléculaires régissant son processus d’activation restent inconnus. Nous avons donc analysé les premières étapes du processus d’activation du NPRA. Nous avons d'abord étudié le rôle de la dimérisation des ECD dans l’activation du récepteur. Nous avons utilisé les techniques de liaison de radioligand, de FRET et de modélisation moléculaire, pour caractériser la liaison à l’ECD des agonistes naturels, d’un superagoniste et d’un antagoniste. L’ANP se lie à un dimère d’ECD préformé et la dimérisation spontanée est l’étape limitante du processus de liaison. De plus, comme le démontrent nos études de FRET, tous les peptides, incluant l’antagoniste, stabilisent le récepteur sous sa forme dimérique. Cependant, l’antagoniste A71915 stabilise le dimère d’ECD dans une conformation différente de celle induite par l’ANP. La dimérisation du NPRA semble donc nécessaire, mais non suffisante à l’activation du récepteur. L’état d’activation du NPRA dépend plutôt de l’orientation des sous unités dans le dimère. Nous avons ensuite étudié le mécanisme moléculaire de transduction du signal à travers la membrane. Plusieurs études ont suggéré que l’activation du NPRA implique un changement de conformation du domaine juxtamembranaire (JM). Cependant, les études de cristallographie de l’ECD soluble de NPRA n’ont pas permis de documenter la structure du JM et le changement de conformation impliqué dans la transduction du signal reste inconnu. Pour analyser ce changement de conformation, nous avons d’abord séquentiellement substitué les neuf acides aminés du JM par une cystéine. En étudiant la capacité des mutants à former des dimères covalents de façon constitutive ou induite par l’ANP, nous avons pu évaluer la proximité relative des résidus du JM, avant et après activation du NPRA. Ces résultats ont démontré la proximité élevée de certains résidus spécifiques et sont en contradiction avec les données cristallographiques. Nous avons également démontré que le domaine intracellulaire impose une contrainte conformationnelle au JM à l’état de base, qui est levée après liaison de l’ANP. En introduisant de 1 à 5 alanines dans l’hélice-α transmembranaire, nous avons montré qu’une rotation des TM de 40° induit une activation constitutive du NPRA. Le signal d’activation pourrait donc être transmis à travers la membrane par un mécanisme de rotation des TM. En utilisant nos données expérimentales, nous avons généré le premier modèle moléculaire illustrant la conformation active du NPRA, où les domaines JM et TM sont représentés. Dans son ensemble, cette étude apporte une meilleure compréhension des mécanismes moléculaires régissant les premières étapes du processus complexe d’activation du NPRA.

Timing verification in transaction modeling

Les systèmes Matériels/Logiciels deviennent indispensables dans tous les aspects de la vie quotidienne. La présence croissante de ces systèmes dans les différents produits et services incite à trouver des méthodes pour les développer efficacement. Mais une conception efficace de ces systèmes est limitée par plusieurs facteurs, certains d'entre eux sont: la complexité croissante des applications, une augmentation de la densité d'intégration, la nature hétérogène des produits et services, la diminution de temps d’accès au marché. Une modélisation transactionnelle (TLM) est considérée comme un paradigme prometteur permettant de gérer la complexité de conception et fournissant des moyens d’exploration et de validation d'alternatives de conception à des niveaux d’abstraction élevés. Cette recherche propose une méthodologie d’expression de temps dans TLM basée sur une analyse de contraintes temporelles. Nous proposons d'utiliser une combinaison de deux paradigmes de développement pour accélérer la conception: le TLM d'une part et une méthodologie d’expression de temps entre différentes transactions d’autre part. Cette synergie nous permet de combiner dans un seul environnement des méthodes de simulation performantes et des méthodes analytiques formelles. Nous avons proposé un nouvel algorithme de vérification temporelle basé sur la procédure de linéarisation des contraintes de type min/max et une technique d'optimisation afin d'améliorer l'efficacité de l'algorithme. Nous avons complété la description mathématique de tous les types de contraintes présentées dans la littérature. Nous avons développé des méthodes d'exploration et raffinement de système de communication qui nous a permis d'utiliser les algorithmes de vérification temporelle à différents niveaux TLM. Comme il existe plusieurs définitions du TLM, dans le cadre de notre recherche, nous avons défini une méthodologie de spécification et simulation pour des systèmes Matériel/Logiciel basée sur le paradigme de TLM. Dans cette méthodologie plusieurs concepts de modélisation peuvent être considérés séparément. Basée sur l'utilisation des technologies modernes de génie logiciel telles que XML, XSLT, XSD, la programmation orientée objet et plusieurs autres fournies par l’environnement .Net, la méthodologie proposée présente une approche qui rend possible une réutilisation des modèles intermédiaires afin de faire face à la contrainte de temps d’accès au marché. Elle fournit une approche générale dans la modélisation du système qui sépare les différents aspects de conception tels que des modèles de calculs utilisés pour décrire le système à des niveaux d’abstraction multiples. En conséquence, dans le modèle du système nous pouvons clairement identifier la fonctionnalité du système sans les détails reliés aux plateformes de développement et ceci mènera à améliorer la "portabilité" du modèle d'application.

Intra- and inter-hemispheric interactions in somatosensory processing of pain : dynamical causal modeling analysis of fMRI data

La douleur est une expérience perceptive comportant de nombreuses dimensions. Ces dimensions de douleur sont inter-reliées et recrutent des réseaux neuronaux qui traitent les informations correspondantes. L’élucidation de l'architecture fonctionnelle qui supporte les différents aspects perceptifs de l'expérience est donc une étape fondamentale pour notre compréhension du rôle fonctionnel des différentes régions de la matrice cérébrale de la douleur dans les circuits corticaux qui sous tendent l'expérience subjective de la douleur. Parmi les diverses régions du cerveau impliquées dans le traitement de l'information nociceptive, le cortex somatosensoriel primaire et secondaire (S1 et S2) sont les principales régions généralement associées au traitement de l'aspect sensori-discriminatif de la douleur. Toutefois, l'organisation fonctionnelle dans ces régions somato-sensorielles n’est pas complètement claire et relativement peu d'études ont examiné directement l'intégration de l'information entre les régions somatiques sensorielles. Ainsi, plusieurs questions demeurent concernant la relation hiérarchique entre S1 et S2, ainsi que le rôle fonctionnel des connexions inter-hémisphériques des régions somatiques sensorielles homologues. De même, le traitement en série ou en parallèle au sein du système somatosensoriel constitue un autre élément de questionnement qui nécessite un examen plus approfondi. Le but de la présente étude était de tester un certain nombre d'hypothèses sur la causalité dans les interactions fonctionnelle entre S1 et S2, alors que les sujets recevaient des chocs électriques douloureux. Nous avons mis en place une méthode de modélisation de la connectivité, qui utilise une description de causalité de la dynamique du système, afin d'étudier les interactions entre les sites d'activation définie par un ensemble de données provenant d'une étude d'imagerie fonctionnelle. Notre paradigme est constitué de 3 session expérimentales en utilisant des chocs électriques à trois différents niveaux d’intensité, soit modérément douloureux (niveau 3), soit légèrement douloureux (niveau 2), soit complètement non douloureux (niveau 1). Par conséquent, notre paradigme nous a permis d'étudier comment l'intensité du stimulus est codé dans notre réseau d'intérêt, et comment la connectivité des différentes régions est modulée dans les conditions de stimulation différentes. Nos résultats sont en faveur du mode sériel de traitement de l’information somatosensorielle nociceptive avec un apport prédominant de la voie thalamocorticale vers S1 controlatérale au site de stimulation. Nos résultats impliquent que l'information se propage de S1 controlatéral à travers notre réseau d'intérêt composé des cortex S1 bilatéraux et S2. Notre analyse indique que la connexion S1→S2 est renforcée par la douleur, ce qui suggère que S2 est plus élevé dans la hiérarchie du traitement de la douleur que S1, conformément aux conclusions précédentes neurophysiologiques et de magnétoencéphalographie. Enfin, notre analyse fournit des preuves de l'entrée de l'information somatosensorielle dans l'hémisphère controlatéral au côté de stimulation, avec des connexions inter-hémisphériques responsable du transfert de l'information à l'hémisphère ipsilatéral.

Modeling protein evolution using secondary structures

L’évolution des protéines est un domaine important de la recherche en bioinformatique et catalyse l'intérêt de trouver des outils d'alignement qui peuvent être utilisés de manière fiable et modéliser avec précision l'évolution d'une famille de protéines. TM-Align (Zhang and Skolnick, 2005) est considéré comme l'outil idéal pour une telle tâche, en termes de rapidité et de précision. Par conséquent, dans cette étude, TM-Align a été utilisé comme point de référence pour faciliter la détection des autres outils d'alignement qui sont en mesure de préciser l'évolution des protéines. En parallèle, nous avons élargi l'actuel outil d'exploration de structures secondaires de protéines, Helix Explorer (Marrakchi, 2006), afin qu'il puisse également être utilisé comme un outil pour la modélisation de l'évolution des protéines.

Effect of modeling methods on the body and head-neck-trunck moments of inertia calculations in individuals of different morphology

Thèse numérisée par la Division de la gestion de documents et des archives de l'Université de Montréal