Caractérisation de structures explorées dans les simulations de dynamique moléculaire.

L’objectif de cette thèse est d’analyser et prédire les conformations d’un système moléculaire en combinant la théorie des graphes et la chimie computationnelle.Dans le cadre des simulations de dynamique moléculaire, une molécule peut avoir une ou plusieurs conformations au cours du temps. Dans les trajectoires de simulation de dynamique moléculaire, on peut avoir des trajectoires n’explorant qu’une seule conformation ou des trajectoires explorant plusieurs conformations, donc plusieurs transitions entre conformations sont observées. L’exploration de ces conformations dépend du temps de la simulation et de l'énergie (température) fixée dans le système. Pour avoir une bonne exploration des conformations d’un système moléculaire, il faut générer et analyser plusieurs trajectoires à différentes énergies. Notre objectif est de proposer un algorithme universel qui permet d’analyser la dynamique conformationnelle de ces trajectoires d’une façon rapide et automatique. Les trajectoires fournissent les positions cartésiennes des atomes du système moléculaire à des intervalles de temps réguliers. Chaque intervalle contenant un ensemble de positions est appelé image. L’algorithme utilise des règles de géométrie (distances, angles, etc.) sur les positions pour trouver les liaisons (liaisons covalentes, liaisons hydrogène et interactions électrostatiques) créées entre les atomes, permettant par la suite d’obtenir le graphe mixte qui modélise une conformation. Nous ne considérons un changement conformationnel que s’il y a un changement dans les liaisons calculées à partir des positions données. L’algorithme permet de donner l’ensemble des conformations explorées sur une ou plusieurs trajectoires, la durée d’exploration de chaque conformation, ainsi que le graphe de transitions qui contient tous les changements conformationnels observés.Les conformations se caractérisent par une énergie appelée énergie potentielle. Cette énergie est représentée par une courbe appelée surface d’énergie potentielle. En chimie théorique et computationnelle, certains s’intéressent à trouver des points particuliers sur cette surface. Il s'agit des minima qui représentent les conformations les plus stables et des maxima ou états de transition qui représentent les points de passage d'une conformation à une autre. En effet, d'une part, la conformation la plus stable est celle de plus basse énergie. D'autres part, pour aller d’une conformation à une autre il faut une énergie supplémentaire, le point maximum représente l'état de transition. Les méthodes développées pour calculer ces points nécessitent une connaissance de l’énergie potentielle ce qui est coûteux en temps et en calculs. Notre objectif est de proposer une méthode alternative en utilisant des mesures ah doc basées sur des propriétés des graphes qu’on a utilisées dans le premier algorithme et sans faire appel à la géométrie ni aux calculs moléculaires. Ces mesures permettent de générer des conformations avec un classement énergétique ainsi de définir le coût énergétique de chaque transition permise. Les conformations possibles avec les transitions représentent respectivement les sommets et les arcs de ce qu’on appelle le “graphe des possibles”. Les hypothèses utilisées dans le modèle proposé est que seules les liaisons hydrogène peuvent changer entre les conformations et que le nombre de liaisons hydrogène présentes dans le système permet de déterminer son coût énergétique.L’algorithme d'analyser des trajectoires a été testé sur trois types de systèmes moléculaires en phase gazeuse de taille et de complexité croissantes. Bien que la complexité théorique de l’algorithme est exponentielle (tests d’isomorphisme) les résultats ont montré que l’algorithme est rapide (quelques secondes). De plus, cet algorithme peut être facilement adapté et appliqué à d’autres systèmes. Pour la prédiction conformationnelle, le modèle proposé a été testé sur des peptides isolés.