Advances in biliary atresia: from patient care to research

Biliary atresia, the most common cause of liver transplantation in children, remains a challenge for clinicians and investigators. The development of new therapeutic options, besides the typical hepatoportoenterostomy, depends on a greater understanding of its pathogenesis and how it relates to the clinical phenotypes at diagnosis and the rate of disease progression. In this review, we present a perspective of how recent research has advanced the understanding of the disease and has improved clinical care protocols. Molecular and morphological analyses at diagnosis point to the potential contributions of polymorphism in the CFC1 and VEGF genes to the pathogenesis of the disease, and to an association between the degree of bile duct proliferation and long-term outcome. In experimental models, cholangiocytes do not appear to have antigen-presenting properties despite a substantial innate and adaptive immune response that targets the biliary epithelium and produces duct obstruction. Initial clinical trials assessing the efficacy of corticosteroids in decreasing the inflammation and improving outcome do not show a superior effect of corticosteroids as an adjuvant treatment following hepatoportoenterostomy. The best outcome still remains linked to an early diagnosis and surgical treatment. In this regard, the Yellow Alert campaign by the Sociedade Brasileira de Pediatria and the inclusion of the Stool Color Card in the health booklet given to every neonate in Brazil have the potential to decrease the age of diagnosis, shorten the time between diagnosis and surgical treatment, and improve the long-term outcome of children with this devastating disease.

Electron paramagnetic resonance study of lipid and protein membrane components of erythrocytes oxidized with hydrogen peroxide

Electron paramagnetic resonance (EPR) spectroscopy of spin labels was used to monitor membrane dynamic changes in erythrocytes subjected to oxidative stress with hydrogen peroxide (H2O2). The lipid spin label, 5-doxyl stearic acid, responded to dramatic reductions in membrane fluidity, which was correlated with increases in the protein content of the membrane. Membrane rigidity, associated with the binding of hemoglobin (Hb) to the erythrocyte membrane, was also indicated by a spin-labeled maleimide, 5-MSL, covalently bound to the sulfhydryl groups of membrane proteins. At 2% hematocrit, these alterations in membrane occurred at very low concentrations of H2O2 (50 µM) after only 5 min of incubation at 37°C in azide phosphate buffer, pH 7.4. Lipid peroxidation, suggested by oxidative hemolysis and malondialdehyde formation, started at 300 µM H2O2 (for incubation of 3 h), which is a concentration about six times higher than those detected with the probes. Ascorbic acid and α-tocopherol protected the membrane against lipoperoxidation, but did not prevent the binding of proteins to the erythrocyte membrane. Moreover, the antioxidant (+)-catechin, which also failed to prevent the cross-linking of cytoskeletal proteins with Hb, was very effective in protecting erythrocyte ghosts from lipid peroxidation induced by the Fenton reaction. This study also showed that EPR spectroscopy can be useful to assess the molecular dynamics of red blood cell membranes in both the lipid and protein domains and examine oxidation processes in a system that is so vulnerable to oxidation.

Étude structure-fonction des fructose-1,6-bisphosphate aldolases métallo-dépendantes : mécanisme catalytique et développement d’antimicrobiens

Les fructose-1,6-bisphosphate aldolases (FBPA) sont des enzymes glycolytiques (EC 4.1.2.13) qui catalysent la transformation réversible du fructose-1,6-bisphosphate (FBP) en deux trioses-phosphates, le glycéraldéhyde-3-phosphate (G3P) et le dihydroxyacétone phosphate (DHAP). Il existe deux classes de FBPA qui diffèrent au niveau de leur mécanisme catalytique. Les classes I passent par la formation d’un intermédiaire covalent de type iminium alors que les classes II, métallodépendantes, utilisent généralement un zinc catalytique. Contrairement au mécanisme des classes I qui a été très étudié, de nombreuses interrogations subsistent au sujet de celui des classes II. Nous avons donc entrepris une analyse détaillée de leur mécanisme réactionnel en nous basant principalement sur la résolution de structures cristallographiques. De nombreux complexes à haute résolution furent obtenus et ont permis de détailler le rôle de plusieurs résidus du site actif de l’enzyme. Nous avons ainsi corrigé l’identification du résidu responsable de l’abstraction du proton de l’O4 du FBP, une étape cruciale du mécanisme. Ce rôle, faussement attribué à l’Asp82 (chez Helicobacter pylori), est en fait rempli par l’His180, un des résidus coordonant le zinc. L’Asp82 n’en demeure pas moins essentiel car il oriente, active et stabilise les substrats. Enfin, notre étude met en évidence le caractère dynamique de notre enzyme dont la catalyse nécessite la relocalisation du zinc et de nombreux résidus. La dynamique de la protéine ne permet pas d’étudier tous les aspects du mécanisme uniquement par l’approche cristallographique. En particulier, le résidu effectuant le transfert stéréospécifique du proton pro(S) sur le carbone 3 (C3) du DHAP est situé sur une boucle qui n’est visible dans aucune de nos structures. Nous avons donc développé un protocole de dynamique moléculaire afin d’étudier sa dynamique. Validé par l’étude d’inhibiteurs de la classe I, l’application de notre protocole aux FBPA de classe II a confirmé l’identification du résidu responsable de cette abstraction chez Escherichia coli (Glu182) mais pointe vers un résidu diffèrent chez H. pylori (Glu149 au lieu de Glu142). Nos validations expérimentales confirment ces observations et seront consolidées dans le futur. Les FBPA de classe II sont absentes du protéome humain mais sont retrouvées chez de nombreux pathogènes, pouvant même s'y révéler essentielles. Elles apparaissent donc comme étant une cible idéale pour le développement de nouveaux agents anti-microbiens. L’obtention de nouveaux analogues des substrats pour ces enzymes a donc un double intérêt, obtenir de nouveaux outils d’étude du mécanisme mais aussi développer des molécules à visée pharmacologique. En collaboration avec un groupe de chimistes, nous avons optimisé le seul inhibiteur connu des FBPA de classe II. Les composés obtenus, à la fois plus spécifiques et plus puissants, permettent d’envisager une utilisation pharmacologique. En somme, c’est par l’utilisation de techniques complémentaires que de nouveaux détails moléculaires de la catalyse des FBPA de classe II ont pu être étudiés. Ces techniques permettront d’approfondir la compréhension fine du mécanisme catalytique de l’enzyme et offrent aussi de nouvelles perspectives thérapeutiques.

Optimisation et approximation adiabatique

L'approximation adiabatique en mécanique quantique stipule que si un système quantique évolue assez lentement, alors il demeurera dans le même état propre. Récemment, une faille dans l'application de l'approximation adiabatique a été découverte. Les limites du théorème seront expliquées lors de sa dérivation. Ce mémoire à pour but d'optimiser la probabilité de se maintenir dans le même état propre connaissant le système initial, final et le temps d'évolution total. Cette contrainte sur le temps empêche le système d'être assez lent pour être adiabatique. Pour solutionner ce problème, une méthode variationnelle est utilisée. Cette méthode suppose connaître l'évolution optimale et y ajoute une petite variation. Par après, nous insérons cette variation dans l'équation de la probabilité d'être adiabatique et développons en série. Puisque la série est développée autour d'un optimum, le terme d'ordre un doit nécessairement être nul. Ceci devrait nous donner un critère sur l'évolution la plus adiabatique possible et permettre de la déterminer. Les systèmes quantiques dépendants du temps sont très complexes. Ainsi, nous commencerons par les systèmes ayant des énergies propres indépendantes du temps. Puis, les systèmes sans contrainte et avec des fonctions d'onde initiale et finale libres seront étudiés.

Étude structurale conformationnelle des toxines de l’anthrax par cryo-microscopie et dynamique moléculaire

Les toxines de l’anthrax font partie de la famille des toxines A-B dans laquelle la moitié B se fixe à la membrane de la cellule permettant par la suite la translocation de la moitié A. Dans le cas de l’anthrax, la moitié B est représentée par le Protective Antigen (PA) et la moitié A par les deux protéines Edema Factor (EF) et Lethal Factor (LF). Après le recrutement par les récepteurs cellulaires (CMG2 et TEM8), PA s’organise en heptamère. Il peut fixer jusqu'à 3 ligands (EF et LF) avant d'être endocyté. Les modèles actuels de PA suggèrent que la baisse de pH à l’intérieur des endosomes permet un changement de conformation de la forme pré-pore vers la forme pore et que les ligands EF et LF passeraient au travers le pore pour entrer dans le cytoplasme. Cependant, le diamètre du pore est environ dix fois inférieur à celui des ligands (10 Å contre 100 Å). Un processus de folding/unfolding a été proposé mais demeure controversé. Afin d'identifier le processus de passage des facteurs EF et LF dans le cytoplasme, nous avons déterminé par cryo-microscopie électronique combinée avec l’analyse d’image les structures tridimensionnelles des complexes formés par PA et LF aux étapes prépore et pore. Par la suite, une étude complémentaire par dynamique moléculaire nous a permis de modéliser à haute résolution les différentes interactions qui ont lieu au sein du complexe. La structure 3D du complexe prépore combiné à 3 LF a été déterminée à une résolution de 14 Å. Nous avons aussi calculé une structure préliminaire du complexe pore également combiné à 3 LF Celles-ci n’ont jamais été résolues auparavant et leur connaissance permet d’envisager l’étude en profondeur du mécanisme infectieux de l’Anthrax in vivo.

Simulations numériques de la dynamique des protéines : translation de ligands, flexibilité et dynamique des boucles

La flexibilité est une caractéristique intrinsèque des protéines qui doivent, dès le mo- ment de leur synthèse, passer d’un état de chaîne linéaire à un état de structure tridimen- sionnelle repliée et enzymatiquement active. Certaines protéines restent flexibles une fois repliées et subissent des changements de conformation de grande amplitude lors de leur cycle enzymatique. D’autres contiennent des segments si flexibles que leur structure ne peut être résolue par des méthodes expérimentales. Dans cette thèse, nous présentons notre application de méthodes in silico d’analyse de la flexibilité des protéines : • À l’aide des méthodes de dynamique moléculaire dirigée et d’échantillonnage pa- rapluie, nous avons caractérisé les trajectoires de liaison de l’inhibiteur Z-pro- prolinal à la protéine Prolyl oligopeptidase et identifié la trajectoire la plus pro- bable. Nos simulations ont aussi identifié un mode probable de recrutement des ligands utilisant une boucle flexible de 19 acides aminés à l’interface des deux domaines de la protéine. • En utilisant les méthodes de dynamique moléculaire traditionnelle et dirigée, nous avons examiné la stabilité de la protéine SAV1866 dans sa forme fermée insérée dans une membrane lipidique et étudié un des modes d’ouverture possibles par la séparation de ses domaines liant le nucléotide. • Nous avons adapté auproblème de la prédiction de la structure des longues boucles flexibles la méthode d’activation et de relaxation ART-nouveau précédemment uti- lisée dans l’étude du repliement et de l’agrégation de protéines. Appliqué au replie- ment de boucles de 8 à 20 acides aminés, la méthode démontre une dépendance quadratique du temps d’exécution sur la longueur des boucles, rendant possible l’étude de boucles encore plus longues.

Étude en dynamique moléculaire par approximation des liaisons fortes de l'influence des défauts ponctuels dans la relaxation du silicium amorphe

Mémoire numérisé par la Division de la gestion de documents et des archives de l'Université de Montréal

Structural rules for the formation of backbone-backbone interactions between closely packed RNA double helices

Les interactions entre les squelettes sucre-phosphate de nucléotides jouent un rôle important dans la stabilisation des structures tertiaires de larges molécules d’ARN. Elles sont régies par des règles particulières qui gouverne leur formation mais qui jusque là demeure quasiment inconnues. Un élément structural d’ARN pour lequel les interactions sucre-phosphate sont importantes est le motif d’empaquetage de deux doubles hélices d’ARN le long du sillon mineur. Ce motif se trouve à divers endroits dans la structure du ribosome. Il consiste en deux doubles hélices interagissant de manière à ce que le squelette sucre-phosphate de l’une se niche dans le sillon mineur de l’autre et vice versa. La surface de contact entre les deux hélices est majoritairement formée par les riboses et implique au total douze nucléotides. La présente thèse a pour but d’analyser la structure interne de ce motif et sa dépendance de stabilité résultant de l’association optimale ou non des hélices, selon leurs séquences nucléotidiques. Il est démontré dans cette thèse qu’un positionnement approprié des riboses leur permet de former des contacts inter-hélices, par l’entremise d’un choix particulier de l’identité des pairs de bases impliquées. Pour différentes pairs de bases participant à ce contact inter-hélices, l’identité optimale peut être du type Watson-Crick, GC/CG, or certaines pairs de bases non Watson-Crick. Le choix adéquat de paires de bases fournit une interaction inter-hélice stable. Dans quelques cas du motif, l’identité de certaines paires de bases ne correspond pas à la structure la plus stable, ce qui pourrait refléter le fait que ces motifs devraient avoir une liberté de formation et de déformation lors du fonctionnement du ribosome.

Études des premières étapes d’oligomérisation de la région Non-Aβ Component de l’a-synucléine et de Aβ1-40

Les protéines amyloïdes sont retrouvées sous forme de fibres dans de nombreuses maladies neurodégénératives. En tentant d’élucider le mécanisme de fibrillation, les chercheurs ont découvert que cette réaction se fait par un phénomène de nucléation passant par des oligomères. Il semblerait que ces espèces soient la principale cause de la toxicité observée dans les cellules des patients atteints d’amyloïdose. C’est pourquoi un intérêt particulier est donc porté aux premières étapes d’oligomérisation. Dans ce mémoire, nous nous intéressons à une séquence d’acide aminé fortement hydrophobe de l’α-synucléine appelée composante non β -amyloïde (Non-Amyloid β Component ou NAC). Cette dernière est retrouvée sous forme de fibres dans les corps et les neurites de Lewy des patients atteints de la maladie de Parkinson. De plus, elle constitue une composante minoritaire des fibres impliquées dans la maladie d’Alzheimer. Nous avons observé les changements structuraux qui ont lieu pour le monomère, le dimère et le trimère de la séquence NAC de l’α-synucléine. Nous nous sommes aussi intéressés aux conséquences structurelles observées dans des oligomères hétérogènes qui impliqueraient, Aβ1−40. Pour cela nous utilisons des dynamiques moléculaires, d’échange de répliques couplées au potentiel gros-grain, OPEP. Nous constatons une disparition des hélices α au profit des feuillets β , ainsi que le polymorphisme caractéristique des fibres amyloïdes. Certaines régions se sont démarquées par leurs capacités à former des feuillets β . La disparition de ces régions lorsque NAC est combinée à Aβ laisse entrevoir l’importance de l’emplacement des résidus hydrophobes dans des structures susceptibles de former des fibres amyloïdes.

Caratérisation de la surface d’énergie potentielle des matériaux complexes et son application sur la cinétique du SiO2/Si

Dans ce rapport de mémoire, nous avons utilisé les méthodes numériques telles que la dynamique moléculaire (code de Lammps) et ART-cinétique. Ce dernier est un algorithme de Monte Carlo cinétique hors réseau avec construction du catalogue d'événements à la volée qui incorpore exactement tous les effets élastiques. Dans la première partie, nous avons comparé et évalué des divers algorithmes de la recherche du minimum global sur une surface d'énergie potentielle des matériaux complexes. Ces divers algorithmes choisis sont essentiellement ceux qui utilisent le principe Bell-Evans-Polanyi pour explorer la surface d'énergie potentielle. Cette étude nous a permis de comprendre d'une part, les étapes nécessaires pour un matériau complexe d'échapper d'un minimum local vers un autre et d'autre part de contrôler les recherches pour vite trouver le minimum global. En plus, ces travaux nous ont amené à comprendre la force de ces méthodes sur la cinétique de l'évolution structurale de ces matériaux complexes. Dans la deuxième partie, nous avons mis en place un outil de simulation (le potentiel ReaxFF couplé avec ART-cinétique) capable d'étudier les étapes et les processus d'oxydation du silicium pendant des temps long comparable expérimentalement. Pour valider le système mis en place, nous avons effectué des tests sur les premières étapes d'oxydation du silicium. Les résultats obtenus sont en accord avec la littérature. Cet outil va être utilisé pour comprendre les vrais processus de l'oxydation et les transitions possibles des atomes d'oxygène à la surface du silicium associée avec les énergies de barrière, des questions qui sont des défis pour l'industrie micro-électronique.

Développement d'une sonde fluorescente bioactivable pour l'étude du rôle in vitro et in vivo des protéases dégradant l’apolipoprotéine A-I.

Les effets bénéfiques des lipoprotéines de haute densité (HDL) contre l'athérosclérose ont été attribués, en grande partie, à leur composante protéique majeure, l'apolipoprotéine A-I (apoA-I). Cependant, il y a des indications que l'apoA-I peut être dégradée par des protéases localisées dans les plaques athérosclérotiques humaines, ce qui pourrait réduire l'efficacité des thérapies basées sur les HDL sous certaines conditions. Nous décrivons ici le développement et l'utilisation d'une nouvelle sonde bioactivatable fluorescente dans le proche infrarouge, apoA-I-Cy5.5, pour l'évaluation des activités protéolytiques spécifiques qui dégradent l'apoA-I in vitro, in vivo et ex vivo. La fluorescence basale de la sonde est inhibée par la saturation du fluorophore Cy5.5 sur la protéine apoA-I, et la fluorescence émise par le Cy5.5 (proche infrarouge) est révélée après clivage de la sonde. La protéolyse in vitro de l'apoA-I par des protéases a montré une augmentation de la fluorescence allant jusqu'à 11 fois (n=5, P ≤ 0.05). En utilisant notre nouvelle sonde apoA-I-Cy5.5 nous avons pu quantifier les activités protéolytiques d'une grande variété de protéases, incluant des sérines (chymase), des cystéines (cathepsine S), et des métalloprotéases (MMP-12). En outre, nous avons pu détecter l'activation de la sonde apoA-I-Cy5.5 sur des sections d'aorte de souris athérosclérotiques par zymographie in situ et avons observé qu'en présence d'inhibiteurs de protéases à large spectre, la sonde pourrait être protégée des activités protéolytiques des protéases (-54%, n=6, P ≤ 0,001). L'infusion in vivo de la sonde apoA-I-Cy5.5 dans les souris athérosclérotiques (Ldlr -/--Tg (apoB humaine)) a résulté en utilisant un système d'imagerie moléculaire combinant la fluorescence moléculaire tomographique et la résonance magnétique,en un signal de fluorescence dans l'aorte plus important que celui dans les aortes des souris de type sauvage C57Bl/6J (CTL). Les mesures in vivo ont été confirmées par l'imagerie ex vivo de l'aorte qui a indiqué une augmentation de 5 fois du signal fluorescent dans l'aorte des souris ATX (n=5) par rapport à l'aorte des souris (n=3) CTL (P ≤ 0,05). L'utilisation de cette sonde pourrait conduire à une meilleure compréhension des mécanismes moléculaires qui sous-tendent le développement et la progression de l'athérosclérose et l'amélioration des stratégies thérapeutiques à base de HDL.

Measurement and simulation of anisotropy in the Infrared and Raman spectra of beta-FeSi2 single crystals

In this paper we show that the orthorhombic phase of FeSi2 (stable at room temperature) displays a sizable anisotropy in the infrared spectra, with minor effects in the Raman data too. This fact is not trivial at all, since the crystal structure corresponds to a moderate distortion of the fluorite symmetry. Our analysis is carried out on small single crystals grown by flux transport, through polarization-resolved far-infrared reflectivity and Raman measurements. Their interpretation has been obtained by means of the simulated spectra with tight-binding molecular dynamics.

Diffusive transport properties in monovalent and divalent metal-ion halide melts: A computer simulation study

Self- and cross-velocity correlation functions and related transport coefficients of molten salts are studied by molecular-dynamics simulation. Six representative systems are considered, i.e., NaCl and KCl alkali halides, CuCl and CuBr noble-metal halides, and SrCl2 and ZnCl2 divalent metal-ion halides. Computer simulation results are compared with experimental self-diffusion coefficients and electrical conductivities. Special attention is paid to dynamic cross correlations and their dependence on the Coulomb interactions as well as on the size and mass differences between anions and cations.

Collective dynamic properties in simple crystals: Influence of the structural disorder

Collective dynamic properties in Lennard-Jones crystals are investigated by molecular dynamics simulation. The study is focused on properties such as the dynamic structure factors, the longitudinal and transverse currents and the density of states. The influence on these properties of the structural disorder is analyzed by comparing the results for one-component crystals with those for liquids and supercooled liquids at analogous conditions. The effects of species-disorder on the collective properties of binary crystals are also discussed.

Heat transfer between nanoparticles: Thermal conductance for near-field interactions

We analyze the heat transfer between two nanoparticles separated by a distance lying in the near-field domain in which energy interchange is due to the Coulomb interactions. The thermal conductance is computed by assuming that the particles have charge distributions characterized by fluctuating multipole moments in equilibrium with heat baths at two different temperatures. This quantity follows from the fluctuation-dissipation theorem for the fluctuations of the multipolar moments. We compare the behavior of the conductance as a function of the distance between the particles with the result obtained by means of molecular dynamics simulations. The formalism proposed enables us to provide a comprehensive explanation of the marked growth of the conductance when decreasing the distance between the nanoparticles.