Le but de cette thèse est d'étudier des moodèles simples d'agents "auto-propulsés": capables de générer du mouvement en consommant de l'énergie provenant de leur environnement, sans forçage externe. Deux modèles de ce type ont été étudiés lors de cette thèse:-Dans un premier temps un modèle de type "Vicsek" a été étudié, c'est à dire que les particules représentées par un couple (position,vitesse) ont une évolution régie par des règles simples d'alignement et d'auto-propulsion à vitesse constante. Ici, l'alignement est nématique: les particlules s'alignent selon leur grand axe, au contraire d'un alignement polaire il se fait indiféremment tête à queue ou tête à tête. Par rapport aux précédents modèles de ce type la première nouveauté est l'introduction d'une pseudo répulsion (dans l'esprit Vicsek, modélisée par un terme de type couple) donnant une extension spatiale à ces particules. La seconde nouveauté est la présence d'un "taux de retournement" qui rend compte du temps de persistence de la direction de l'auto-propulsion. Dans cette partie nous décrivons divers diagrammes de phases de ce nouveau modèle qui montrent de nouvelles phases non répertoriées précédemment: les arches mais aussi des bandes "smectiques", quelques propriétés de ces structures ont été mesurées. Des équations hydrodynamiques obtenues via la méthode "Boltzmann-Ginzburg-Landau" ayant par ailleurs été dérivées nous effectuons une comparaison: la plupart des phases ainsi que certaines de leurs propriétés sont retrouvés dans le modèle hydrodynamique.-Dans un second temps, nous étudions la bactérie Neisseria Meningitidis qui présente la particularité de générer des "pili", filaments de plusieurs micromètres de long. En dépolymérisant ces structures, à vitesse constantes (~1 µm/s), elle est capable de en générer des forces gigantesques pour le vivant (~100 pN). Cette bactérie a tendance à former des aggrégats sphériques, présentant toutes les propriétés d'un liquide, pour coloniser l'organisme de l'hôte.Des mesures de viscosité et de tension de surface de ces aggrégats ont montré le rôle crucial du nombre de pili. Fort de ces constats nous avons bati un modèle microscopique dont la particularité est l'introduction de potentiels stochastiquement attractifs, c'est à dire que les particles transitent entre un état attractif et un état diffusif. Cette partie retranscrit l'évolution du modèle au cours du temps. Nous avons pu reproduire certaines propriétés des aggrégats, nous avons notamment mis en évidence une variation de la diffusion entre le centre et le bord des aggrégats qui recoupe les données expérimentales. / This thesis purpose is to study simple "self-propelled" agents models: they are able to generate motion by consumming energy comming from their environment, without external forcing. Two models of that kind have been studied:-In the first part a Vicsek-style model has been studied, that is we particles are modeled by a couple (position,velocity) which evolution is dictated by simple rules of alignment and self-propulsion at constant speed. Here the alignment is nematic particles align along their long axis and alignment is not polar, contrarily to a polar alignment particles don't discriminate between head and tail . Compared to previous models of this type, the first novelty is the introduction of a pseudo-repulsion (in the Vicsek-spirit, modelized by a torque-like term) providing spatial extension to these particles. The second addition is a flipping rate which renders the persistence time of the direction of self-propulsion. In this part we describe several phase diagrams of this new model which show new phases not previously classified: arches but also "smectic" bands, some propreties of these structures have been measured. Hydrodynamic equations from the "Boltzmann-Ginzburg-Landau" method have been also developped, comparisons are performed: the hydrodynamic model recovers most phases and some of their propreties.-In the second part we study Neisseria Meningitidis, a bacteria which particularity is to generate pili: filamentous structures several micrometers long. By depolymerizing these structures at constant speed (~1µm/s), it is able to generate gigantic forces for the living word (~ 100pN). This bacteria has a tendancy to form spherical aggregates, showing all propreties of a liquid, in order to colonize the host organism. Viscosity and surface tension measure of these aggregates have shown the crucial role of the pili number. Using these data we've built a microscopic model which particularity is the presence of a stochastically attractive potential, that is to say that particles are transiting between an attractive state and a diffusive one. This part relates the model evolution in time. We've ben able to reproduce some aggregate propreties, in particular we've highlighted a variation of the diffusion between aggregate center and edges which fits experimental data.

L’appréhension des mécanismes d’action biologiques des protéines membranaires nécessite de comprendre les interactions des xénobiotiques avec ces protéines et avec les membranes lipidiques. Les méthodes expérimentales sont parfois coûteuses et ne permettent d’obtenir que des informations partielles sur les interactions xénobiotiques-membrane-protéine. La modélisation moléculaire est une sérieuse alternative. Les simulations de dynamique moléculaire et de dynamique biaisées ont ouvert de nombreuses perspectives en permettant de décrire ces interactions moléculaires à l’échelle atomique. Grâce à des simulations de dynamique moléculaire, nous avons été capables de construire un modèle de transporteur humain ABC : ABCC4/MRP4. Cette protéine a été choisie pour sa présence dans le rein, notamment, et son importance clinique. Nous avons évalué l’influence du cholestérol sur cette protéine. L’étude de domaines spécifiques et l’impact d’un polymorphisme a été reliée à l’activité de transport de cette protéine. Nous avons également étudié l’interaction de xénobiotiques avec ce transporteur humain. Le cycle de transport des transporteurs ABC a été examiné afin de comprendre leur fonctionnement. L’incorporation de cholestérol a montré un impact significatif sur la protéine humaine ABCC4/MRP4 et sur les xénobiotiques étudiés. L’importance de domaines constituant la protéine ABCC4/MRP4 ainsi que l’importance de résidus individuels a clairement été prouvée. Nous avons également pu observer des intermédiaires du cycle de transport d’un transporteur ABC conjointement avec des changements structuraux. / Understanding the biological mechanisms of action of membrane proteins requires the comprehension of the interactions of xenobiotics with these proteins and with lipid membranes. Experimental methods are often demanding and only partially respond to xenobiotic-membrane-protein interactions. In silico molecular modeling is a serious alternative to tackle these issues. Molecular dynamics (MD) and biased dynamics simulations have opened many perspectives by providing an atomistic description of these intermolecular interactions. Using MD simulations, we built a model of the human ABC ABCC4/MRP4 transporter. We explored the influence of cholesterol on this protein as well as the impact of a polymorphism known to shut down the transport activity of this protein. We also studied the interaction of xenobiotics with this human transporter. The transport cycle of the ABC transporters was investigated in an attempt to better understand how it works.Interactions between lipid membranes and xenobiotics were explored by examining their ability to incorporate lipid membranes. Lipid mixtures with cholesterol showed a significant impact on the human protein ABCC4/MRP4 and on the xenobiotics studied. The importance of regions, domains constituting the ABCC4/MRP4 protein as well as the importance of specific residues has been clearly demonstrated. We also observed intermediates in the transport cycle of an ABC transporter in conjunction with structural changes occurring during this cycle.

La Structuration des molécules d’eau à l’interface eau/vapeur d’eau fait l’objet de l’intérêt scientifique depuis des années. La plupart des études sont focalisées sur le bulk d’eau mais des études plus détaillées sur l’eau de surface sont nécessaires. De plus, les interfaces avec les alcanes sont intéressantes d’un point de vue biologique et industriel. Puisque pour des applications biologiques et industrielles les interfaces eau/air et eau/huile possèdent des médiateurs amphiphiles, l’influence d’une monocouche de tensioactif sur la structuration de la surface de l’eau mérite aussi une attention particulière.Dans cette thèse, plusieurs modèles atomistiques d’eau ont été sélectionnés. Des simulations de dynamique moléculaire classique ont été réalisées à 298K pour le bulk d’eau, des systèmes eau/vapeur d’eau et eau/alcane (C5-C9), ainsi que des systèmes eau/DLPC/vapeur d’eau et eau/DLPC/octane (DLPC : dilauroyl phosphatidylcholine). Plusieurs propriétés structurelles, ainsi que les moments dipolaires, la tension de surface et les liaisons hydrogène, du bulk d’eau et des couches de surface d’eau ont été examinées grâce à la fonction de distribution radiale et les diagrammes de Voronoi. L’objectif a été d’estimer l’impact de l’incorporation de la polarisabilité sur les propriétés de l’eau et de sélectionner un modèle optimal (qualité/temps de calcul) pour leur description, ainsi que d’enrichir les données existantes sur la structuration de l’eau à l’interface.Cette étude aborde la structuration de l’eau du bulk et de la surface à l’interface avec la vapeur d’eau ou des alcanes. Un des objectifs de ce travail a été d’évaluer la reproductibilité des données expérimentales en utilisant différents modèles et de s’assurer pour quelles propriétés l’utilisation d’un modèle polarisable est critique. De simples modèles polarisables basés sur les oscillateurs de Drude ont été testés afin de limiter le temps de calcul. Pour le bulk d’eau et les systèmes eau/vapeur d’eau, les modèles TIP4P, SWM4-NDP et COS/G2 ont été les plus performants. Dans la mesure où le modèle TIP4P produit des résultats commensurables avec les modèles polarisables, il a été utilisé pour la simulation des interfaces eau/alcanes (C5-C9). Les molécules à la surface sont organisées de manière plus compacte et moins ordonnée. Cette diminution de l’ordre est principalement due aux liaisons hydrogène qui sont deux fois plus nombreuses dans le bulk qu’en surface. Les analyses de Voronoi ont montré que la coordination tétraédrique n’est pas si claire et que des polyèdres plus complexes sont formés. Les couches de surfaces trouvées s’avèrent être formées de 2 sous-couches, interne et externe, avec des polarités inégales orientées de manières opposées, définissant des zones de charges résiduelles à l’interface.En plus des systèmes avec un contact direct entre l’eau et le fluide apolaire, des interfaces comportant des monocouches de lipides ont été modélisées. La compacité de l’eau de surface, déjà renforcée par la présence d’alcanes, a été augmentée davantage par l’introduction de lipide. Néanmoins, l’orientation de l’eau a été changée et la polarité de la surface inversée, équilibrée par les têtes lipidiques au lieu des sous-couches externes diffuses.Les résultats principaux de cette thèse de doctorat sont les suivants:1. Il a été montré que l’utilisation de modèles d’eau polarisable n’est pas nécessaire pour une évaluation correcte d’un certain nombre de propriétés, mais elle est critique pour les moments dipolaires et la tension de surface.2. Pour la première fois une analyse structurelle a été réalisée en utilisant les diagrammes de Voronoi et un ensemble de modèles d’eau démontrant la différence entre propriétés de l’eau liquide en bulk et en surface.3. Considérant le nombre limité de données existantes, l’étude d’une monocouche de DLPC solide condensée à l’interface eau/vapeur et eau/octane en utilisant différents modèles d’eau a été une contribution originale. / The structuring of water molecules at the water/vapour interface is an object of scientific interest for decades. Most of the existing theoretical studies are focused on bulk water but there is still need of a more detailed research on surface water. In addition, interfaces with alkanes are interesting as being instructive from both biological and industrial perspectives. Since in both bio- and industrial applications water/air and water/oil interfaces are mediated by amphiphiles, the role of a surfactant monolayer on surface water structuring deserves more attention as well.In the present Ph. D. thesis several atomistic water models were chosen and classical molecular dynamics simulations were carried out on bulk water, water/vapour and water/alkane (from pentane to nonane) systems, as well as on water/DLPC/vapour and water/DLPC/octane models, DLPC being dilauroyl phosphatidylcholine. In all cases the temperature was kept at 298 K. Several structural properties of bulk and surface water layers were examined by means of radial distribution functions and Voronoi diagrams. Dipole moments, surface tension and hydrogen bonding were tackled too. The objective was to estimate the impact of accounting for polarisability on the water properties of interest and to select a cost-efficient water model for describing them, as well as to add new data to the existing knowledge about interfacial water structuring.The study addresses the water structuring in bulk and surfacial water at the interface with vapour or alkanes of different chain length. One of the aims of the work was to assess the reproducibility of experimental data using an assortment of polarisable and non-polarisable water models and to check for which properties the utilisation of polarisable models is critical. Simple polarisable models based on Drude oscillators were tested in order to keep the computational costs low. For bulk water and water/vapour systems the models TIP4P, SWM4-NDP and COS/G2 performed the best. Since the TIP4P model produced results commeasurable with the polarisable ones, it was used predominantly further on to simulate water/alkane (C5-C9) interface and to quantify the structural parameters of water obtained from RDFs and Voronoi analyses. The molecules in this layer are organised in a more compact and less ordered manner. The ordering is owed mainly to hydrogen bonds which are twice as many in the bulk compared to the surface. The analysis of the Voronoi diagrams showed that the tetrahedral coordination was blurred and more complex polyhedra were formed. The surface layer was found to consist of two sublayers, inner and outer, with oppositely oriented unequal polarity, defining areas of residual charges at the interface.In addition to the systems with direct contact between water and non-polar fluids, interfaces mediated by lipid monolayers were modelled. The monolayer was meant to seam together the two phases. The compactness of the surfacial water, which was enhanced by the presence of alkanes, was tightened further by the lipid introduction. However, the water orientation was changed and the surfacial polarity was inverted, balanced by the lipid heads instead of the diffuse outer sublayer.The main contributions of the Ph.D. thesis are as follows:1. It is shown that the usage of a polarisable water model is not necessary for correct evaluation of a number of properties, but is critical for characteristics such as dipole moments and surface tension.2. For the first time a structural analysis has been made using Voronoi diagrams and an assortment of water models which demonstrates the difference between bulk and surfacial characteristics of liquid water.3. An original contribution is the study of a solid-condensed DLPC monolayer at the water/vapour interface utilising different water models and at the interface of water/octane, considering the limited experimental data available.

Thèse numérisée par la Direction des bibliothèques de l'Université de Montréal.

Nous avons présenté COMMA, une méthode pour décrire et comparer les architectures dynamiques de différentes protéines. Il extrait propriétés dynamiques de ensembles conformationnels pour identifier les voies de communication, des chaînes de résidus liés par des interactions stables qui se déplacent ensemble, et cliques indépendants, des groupes de résidus qui fluctuent de manière concertée. Il fournit une description de l'infostery d'un complexe de protéines qui va au-delà des mesures au-delà de classiques de la façon dont une protéine se déplace ou change de forme. Nous avons montré l'efficacité de notre approche pour fournir des idées mécanistiques sur les effets des mutations délétères en identifiant les résidus qui jouent un rôle clé dans la propagation de ces effets. En outre COMMA révèle un lien entre les clusters de coévoluant résidus et les réseaux de corrélations dynamiques. Il permet de comparer les différents types de communication se produisant entre les résidus et de hiérarchiser les différentes régions d'une protéine en fonction de l'efficacité de leur communication. En outre, nous avons présenté une approche pour exploiter les séquences et les dynamiques structurelles pour prédire un paysage mutationnel. La discussion des exemples, a révélé l'interprétation physique sur la façon dont l'étude de la conservation apporte des idées importantes sur la sensibilité des positions conservées à des mutations. Notre méthode proposée, peut détecter des régions de protéines qui sont sujettes à des troubles ou des réarrangements conformationnels substantiels. De plus, il nous a permis de proposer des mutations qui régulent la stabilité des bobines enroulées désordonnées. / We presented COMMA, a method to describe and compare the dynamical architectures of different proteins or different variants of the same protein. COMMA extracts dynamical properties from conformational ensembles to identify communication pathways, chains of residues linked by stable interactions that move together, and independent cliques, clusters of residues that fluctuate in a concerted way. It provides a description of the infostery of a protein or protein complex that goes beyond the notions of chain, domain and secondary structure element/motif, and beyond classical measures of how a protein moves and/or changes its shape. We showed the efficiency of our approach in providing mechanistic insights on the effects of deleterious mutations by pinpointing residues playing key roles in the propagation of these effects. In addition COMMA reveals a link between clusters of coevolving residues and networks of dynamical correlations. It enables to contrast the different types of communication occurring between residues and to hierarchise the different regions of a protein depending on their communication efficiency. Furthermore, we presented an approach to exploit both the sequences and structural dynamics to predict a mutational landscape. The discussion of examples, revealed physical interpretation on how the study of conservation brings significant insights on the sensitivity of conserved positions to mutations. Our proposed method, can detect protein regions that are prone to disorder or substantial conformational rearrangements. Moreover, it enabled us to suggest mutations that regulate the stability of the disordered coiled-coils.

Les générateurs thermoélectriques convertissent directement l'énergie thermique en énergie électrique. Ils pourraient devenir de plus en plus utiles à des fins de récupération d'énergie et font l'objet de recherches actives. Cependant, les meilleurs matériaux thermoélectriques sont rares et polluants.Le Silicium et le Germanium seraient des matériaux intéressants si leur efficacité thermoélectrique était améliorée. Pour ce faire, la nanostructuration est une voie possible, par exemple en introduisant des faces rugueuses ou de nouvelles interfaces semi-transparentes.Récemment, des nanofils polyphasés (composés d'une alternance de phases cubiques et hexagonales de Si et Ge) ont été fabriqués, mais la caractérisation expérimentale de nanostructures aussi complexes comprenant des matériaux exotiques peut être difficile. Dans cette thèse, nous étudions en détail le transport thermique dans des nanostructures avec des simulations numériques. Une méthode Monte Carlo originale a été développée, avec une description "full band" des matériaux. Elle inclut des modèles pour les faces rugueuses et les interfaces entre matériaux. Des simulations de Dynamique Moléculaire sont également effectuées pour caractériser les propriétés des interfaces.Nous confirmons que les phases hexagonales de Si et Ge ont une conductivité thermique inférieure à celle des phases cubiques correspondantes. Le modèle "full band" montre que le flux thermique est fortement anisotrope. Des modèles semi-analytiques habituels n'ont pas pu reproduire la conductivité thermique des nanostructures simulées avec des faces rugueuses.De plus, ces faces ont tendance à concentrer le flux de chaleur dans la direction principale de la nanostructure. Enfin, certaines interfaces polyphasées peuvent avoir une conductance thermique presque aussi faible que les interfaces Si/Ge, et pourrait ainsi améliorer significativement l'efficacité thermoélectrique des nanofils polyphasés. La méthode Monte Carlo présentée peut facilement être utilisée pour étudier une large gamme de matériaux, et elle est capable de modéliser des nanostructures arbitrairement complexes. A l'avenir, les simulations en Dynamique Moléculaire seront utilisées pour paramétrer un modèle plus réaliste d'interfaces. / Thermoelectric generators are able to directly convert heat into electrical energy. They could have a great potential in terms of energy harvesting, but unfortunately, the best thermoelectric materials are rare and pollutant.Silicon and Germanium would be attractive materials if their thermoelectric efficiency were improved. For this purpose, nanostructuring is a possible route, for instance via the introduction of rough boundaries or interfaces between materials.Recently, polytype nanowires (composed of a sequence of cubic and hexagonal phases of Si and Ge) have been fabricated, but the experimental characterization of such complex nanostructures with exotic materials is challenging.In this thesis, we study the details of thermal transport in nanostructures with numerical simulations. An original Monte Carlo method is developed, with a full band emph{ab initio} description of materials. It includes models for the rough boundaries and the solid-solid interfaces. Molecular Dynamics simulations are also performed to characterize the properties of interfaces.We confirm that the hexagonal phases of Si and Ge have lower thermal conductivity than their cubic counterparts. The full band model shows a strong anisotropy in the heat flux.Usual semi-analytical models failed to reproduce the thermal conductivity of simulated nanostructures with rough boundaries. Besides, those boundaries tend to focus the heat flux in the main direction of the nanostructure. Finally, some polytype interfaces can have an interfacial conductance almost as low as Si/Ge interfaces, and thus could improve significantly the thermoelectric efficiency of polytype nanowires. The presented Monte Carlo method could easily be used with a wide range of materials,and it can model arbitrarily complex nanostructures. In the future, the results from Molecular Dynamics simulation will be used to parametrize a more realistic model of solid-solid interfaces.

Cette thèse porte sur l'étude des effets quantiques en dynamique moléculaire (DM). La DM est une méthode numérique qui permet l'étude des propriétés de la matière condensée. Cependant, la méthode étant basée sur la mécanique classique, les effets quantiques associés à la dynamique des noyaux, tels que l'énergie de point zéro ou l'effet tunnel, ne sont pas pris en compte. Ces effets quantiques nucléaires peuvent cependant jouer un rôle majeur, en particulier aux basses températures et/ou dans les systèmes contenant des atomes légers comme l'hydrogène. La dynamique moléculaire par intégrales de chemins (PIMD) est souvent utilisée, dans ce cas, pour tenir compte de la nature quantique des noyaux. Cette approche fournit des résultats quantiques exacts, mais son coût en temps de calcul élevé limite son domaine d'application. La méthode du bain thermique quantique (QTB) a été proposée comme une alternative à la PIMD. L'approche QTB est particulièrement intéressante car son coût en temps de calcul est équivalent à celui de la DM standard permettant ainsi l'étude de systèmes complexes et de plus grande taille.La première partie de cette thèse est consacrée à l'étude de la méthode QTB. Nous avons étudié le comportement de la méthode sur différents systèmes modèles afin d'étudier ses limites. En particulier, le problème du "zero point energy leakage" est étudié en détail et nous montrons que l'augmentation du coefficient de friction du QTB permet de limiter ce problème. Nous avons également développé une combinaison de la méthode QTB avec la méthode PIMD. Cette méthode combinée QTB-PIMD permet de réduire le coût en temps de calcul des simulations PIMD standards.Dans une deuxième partie, nous avons utilisé ces méthodes pour étudier la conduction de l'hydrogène dans des matériaux pérovskites. Nous nous intéressons d'abord à l'impact des effets quantiques sur la diffusion de l'hydrogène dans BaZrO3, un matériau d'électrolyte potentiel pour piles à hydrogène. L'hydrogène étant l'élément le plus léger, un impact important des effets quantiques est attendu. Nous trouvons que les effets quantiques sont effectivement importants à basse température, mais leur impact sur la diffusion reste faible aux températures de fonctionnement typiques des piles à hydrogène. Enfin, nous avons étudié les mécanismes de diffusion de l'hydrogène dans GdBaCo2O5.5. Nous mettons en évidence une diffusion anisotrope dans ce matériau et deux mécanismes principaux de diffusion. / This thesis deals with the study of quantum effects in molecular dynamics (MD). MD is a powerful numerical method to investigate the properties of condensed matter systems. However, since the method is based on classical mechanics, quantum effects associated with the dynamics of the nuclei, such as zero-point energy or tunnelling, are not taken into account. These nuclear quantum effects can, however, play a major role in particular at low temperatures and/or in systems containing light atoms such as hydrogen. In these cases, a standard way to account for the quantum nature of the nuclei is to use path integral molecular dynamics (PIMD). This method provides exact quantum results however its high computational cost limits its range of applicability. The quantum thermal bath (QTB) method has been proposed as an alternative to PIMD. The QTB method is particularly appealing because of its computational cost that is equivalent to standard MD thus allowing to study large and complex systems.The first part of this thesis is devoted to the study of the QTB method. We have studied the behavior of the method in selected model systems in order to investigate its limitations. We have focused, in particular, on the zero-point energy leakage problem and found that increasing the friction coefficient of the QTB can significantly limit this problem. We also have developed another way to use the QTB method by combining it with PIMD. This combined QTB-PIMD method allows, in particular, to decrease the computational cost of standard PIMD simulations.In a second part, we have used these methods to study hydrogen conduction in perovskite materials. We have first investigated the impact of quantum effects on the diffusion of hydrogen in BaZrO3, a potential electrolyte material for hydrogen fuel cells. Since hydrogen is the lightest element, we expect quantum effects to have a significant impact on its dynamics. We find that quantum effects are indeed significant at low temperatures although their impact on the diffusion remains low at the typical working temperatures of hydrogen fuel cells. Finally, we have investigated the diffusion mechanisms of hydrogen in GdBaCo2O5.5. We evidence that the diffusion is anisotropic in this material and two main diffusion mechanisms.

L'étude de la dynamique intrinsèque de l'ADN-B permet de caractériser l'espace conformationnel exploré par cette macromolécule. Cette dynamique, qui dépend de la séquence, est un facteur clé dans les mécanismes d'interaction avec les protéines, l'ADN étant plus ou moins prédisposé à s'adapter à son partenaire. Lors de cette thèse, nous avons sondé la dynamique de l'ADN à deux échelles de temps, en nous basant sur l'étude de quatre dodécamères par Résonance Magnétique Nucléaire (RMN). Nous nous sommes d'abord intéressés à la dynamique des groupements phosphates, qui correspondent à des mouvements rapides (picoseconde-nanoseconde). Nous avons ainsi confirmé l'effet de la séquence dinucléotidique sur cette dynamique, qui peut être prédit et quantifié. Nous avons également mis en lumière pour la première fois l'étroite inter-dépendance qui existe entre déplacements chimiques du phosphore, distances internucléotides et constantes de couplage dipolaire résiduel. L'interprétation de ces observables RMN en termes de conformations des phosphates, de paramètres hélicoïdaux et de taille des sillons, montre qu'en fait ces couplages reflètent la mécanique intrinsèque de l'ADN en solution. En interfaçant ces résultats avec l'effet de séquence observé sur la dynamique des phosphates, il est aussi possible de saisir de quelle façon la conformation moyenne de la double hélice et l'espace conformationnel associé sont modulés par la séquence au niveau dinucléotidique. Enfin, des dynamiques moléculaires réalisées avec les très récents champs de force CHARMM36 et Parmbsc0ezOLI, confrontées aux données expérimentales, ont permis d'apprécier le réalisme croissant des ADN simulés et ont aidé à préciser des éléments de la dynamique qui échappent à l'expérience. Le deuxième volet de cette thèse a porté sur les mouvements de l'ADN se produisant à l'échelle de la milliseconde, encore très peu étudiés. Nous avons mis au point des expériences de dispersion-relaxation qui ont apporté la preuve de l’existence d’un échange conformationnel d'un type totalement nouveau. Cet échange ne semble apparaitre que sur un type particulier de séquence, riche en A:T. Certaines régions de l’ADN, probablement spécifiques, peuvent ainsi localement évoluer vers une forme très faiblement peuplée, dont la structure détaillée reste à caractériser. L'ensemble de ces résultats offre un panorama des capacités dynamiques de l'ADN, dépendantes de la séquence, et ouvre ainsi de nombreuses perspectives vers une meilleure compréhension des mécanismes qui guident la formation des complexes ADN-protéines. / The study of B-DNA intrinsic dynamics enables to characterize the conformational landscape explored by this macromolecule. Indeed, binding of DNA to proteins is modulated by subtle sequence-dependent variations inherent to the dynamics of free DNA, which facilitate or disfavor the structural fit with cognate partners.In this thesis, the DNA dynamics was investigated at two time-scales, on the basis of the Nuclear Magnetic Resonance (NMR) study of four dodecamers. First, we examined the fast dynamics (pico-nanosecond) of phosphate linkages. We confirmed that the dinucleotide sequence modulates the backbone dynamics, an effect that can be quantified and predicted. Then, our experimental data enabled to establish that phosphorus chemical shifts, internucleotide distances and residual dipolar couplings constants are closely correlated. The translation of the NMR observables in terms of phosphate conformations, helicoidal parameters and minor groove dimension, allowed the structural interpretation of the couplings and led to the first coherent description of the intrinsic DNA mechanics in solution. Owing our knowledge of the effect of the sequence on the backbone behavior, it is now possible to understand how the DNA shape and the associated conformational landscape are modulated at the dinucleotide level. Finally, the performance of molecular dynamics (MD) simulations with the recent force-fields Parmbsc0εζOLI and CHARMM36 was tested extensively against our NMR data. We found impressive progress towards a realistic representation of DNA, despite residual shortcomings. This advance allowed to reveal new aspects of the DNA dynamics, which cannot be assessed from experiments.The second part of this thesis focused on slow motions in B-DNA, which are still largely under-investigated. Using and developing sophisticated relaxation-dispersion NMR experiments, we demonstrated the existence of a new conformational exchange at the millisecond time-scale, which seems to only occur in a particular type of sequence, A:T rich. Thus, in addition to the familiar structural patterns that are the signature of the B double helix, some short DNA regions, likely specific, are able to explore another conformational state, weakly populated, whose detailed structure still needs to be characterized.Overall, these results provide original insights on the DNA dynamic repertoire, sequence-dependent, and open the way towards a better understanding of the mechanisms underlying the formation of DNA-protein complexes.

En 2016, 216 millions de cas et 445000 décès liés au paludisme ont été reportés par l'Organisation Mondiale de la Santé (OMS), 90 % de ces cas ont eu lieu sur le continent africain, très majoritairement causés par Plasmodium falciparum. Chez les femmes enceintes le paludisme peut avoir de lourdes conséquences en termes de morbidité et de mortalité aussi bien pour la mère que pour l'enfant à venir. Outre l'utilisation de moustiquaires imprégnées, l'OMS recommande l'utilisation d'un traitement préventif intermittent par sulfadoxine-pyriméthamine (SP) pour les femmes enceintes des zones endémiques d'Afrique sub-saharienne. La résistance du parasite à la SP est un réel problème en l'absence d'alternative thérapeutique et le niveau de résistance est déjà très élevé dans certaines régions d'Afrique. Cette résistance est causée par différents polymorphismes nucléotidiques simples (SNPs) sur les gènes codant pour DHPS et DHFR. Ces mutations modifient la structure 3D de ces enzymes, diminuant leur liaison à la SP. Ce travail s'est attaché à faire un état des lieux récent de la prévalence des mutations des gènes pfdhfr et pfdhps à travers 7 sites d'Afrique Centrale. Un octuple mutant combinant trois mutations sur le gène pfdhfr et cinq mutations sur le gène pfdhps (CirnI + vagKgs) a été découvert à Yaoundé en 2015 ; nous l'avons retrouvé au Nigeria et au Cameroun mais pour la première fois à un haut niveau de prévalence à Maroua (52,2%). Sur plusieurs sites d'Afrique Centrale, nous décrivons l'augmentation de la prévalence de la mutation K540E, typique des parasites résistants de l'Afrique de l'Est et sur laquelle se base l'implémentation de la SP par l'OMS. D'autre part, à partir d'un modèle d'homologie de pfDHPS, nous avons utilisé des techniques de dynamique moléculaire pour mieux comprendre les modifications de la structure 3D liées à différents haplotypes d'intérêts retrouvés en Afrique Centrale. Cette partie de l'étude propose une autre approche pour évaluer et visualiser l'effet structurel des mutations, élément supplémentaire à la compréhension de leur impact sur la résistance à la SP. / Malaria was responsible for 445 000 death and 216 million new cases worldwide in 2016, 90% of these cases occurred in Africa. Pregnant women are particularly susceptible to malaria, resulting in medical consequences on both mother and child. In the endemic countries of sub Saharan Africa, in addition to the use of insecticidal nets, the World Health Organization (WHO) recommends a preventive treatment against Plasmodium falciparum malaria in pregnancy using Sulfadoxine Pyrimethamine (SP). Resistance of the parasite to SP is a significant problem in the absence of alternative treatment and the level of resistance is very high in some regions of Africa. This resistance is mediated by several single nucleotide polymorphisms (SNPs) in the genes coding for DHPS and DHFR, changing the 3D structure of the two target proteins and decreasing drug binding. This study offers an updated view of the epidemiology of these mutations throughout Central Africa. We identified for the first time the octuple mutant combining 3 mutations on pfdhfr gene and 5 mutations on pfdhps gene (CirnI + vagKgs), discovered in 2015 in Yaoundé, as the major haplotype in Maroua, northern Cameroon (52,2%). On several Central-African sites, we describe the increase of the prevalence of the typically East-African K540E mutation on which is based the implementation of SP by the WHO. Knowledge of the prevalence of resistance markers is crucial for the adaptation of SP recommendation in Central Africa. This study also explores the changes occurring in the 3D structure of the most common mutated pfDHPS haplotypes identified in our epidemiology study, using molecular dynamics on a homology model of the protein. It offers a new vision of the effect of the mutations on the structure, a further step to understanding the impact of mutations in resistance.

Mémoire numérisé par la Direction des bibliothèques de l'Université de Montréal.