Etude par dynamique moléculaire de formation d'excimères et de mélanges de lipides oxydés dans les membranes lipides modèles

Molecular dynamics study of pyrene excimer formation and oxidation in lipid bilayer models

• Auteur / Autrice : Pierre Ayoub
• Direction : Fabrice Thalmann
• Type : Thèse de doctorat
• Discipline(s) : Physique
• Date : Soutenance le 16/12/2015
• Etablissement(s) : Strasbourg
• Ecole(s) doctorale(s) : École doctorale Physique et chimie-physique (Strasbourg ; 1994-....)
• Partenaire(s) de recherche : Laboratoire : Institut Charles Sadron (Strasbourg ; 1985-....)
• Jury : Président / Présidente : Joachim Wittmer
• Jury : Rapporteurs / Rapporteuses : Luca Monticelli, Olivier Sandre

Mots clés

Mots clés contrôlés

Modèles mathématiques

Dynamique moléculaire

Membranes lipidiques

Lipides des membranes

Excimères

Mots clés libres

Modèle membrane

Dynamique moléculaire

DOPC

POPC

Gros-grains

Formation d' excimère

Pyrène

Oxidation

Coefficient de diffusion

Model membrane

Molecular dynamics

DOPC

POPC

Coarse-grain

Excimer formation

Pyrene

Oxidation

Diffusion coefficient

Résumé

Nous proposons une nouvelle approche pour déterminer le coefficient de diffusion dans des membranes lipidiques se basant sur la formation d'excimères. Alors que les autres modèles statistiques considèrent le système comme un ensemble de points sur un réseau, nous utilisons un modèle à gros grain afin d'étudier des bicouches lipidiques simulées à l'aide du champs de force Martini. Nous déterminons le taux de réaction dépendant du temps à partir des probabilités de survie obtenues a posteriori à l'aide des trajectoires numeriques des bicouches symétriques de DOPC (1,2-Dioleoyl-sn-glycero-3-phosphocholine) et POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) simulées à 283 K et 293 K respectivement. Les dynamiques de collision sont obtenues en distinguant virtuellement les molécules simulées. Les sondes fluorescentes sont supposées semblables aux lipides, et par conséquent, ne modifient pas la dynamique. Nous obtenons une expression générale pour la probabilité de survie en combinant approximation des paires indépendantes et propriétés d'échelle, mais aucune hypothèse n'est faite pour le taux de formation d'excimère. En superposant les intensités d'émission de fluorescence normalisées, déterminées numériquement, aux courbes de titrations expérimentales, nous obtenons deux ensembles de résultats pour le coefficient de diffusion latéral, selon que l'association entre feuillets est autorisée ou pas. Nous utilisons un rayon de capture de 0.5 nm, la distance à partir de laquelle les deux sondes réagissent pour former un excimère. En comparant la dynamique Martini aux expériences de fluorescence, il est possible d'estimer le facteur d'accélération.

We propose a novel approach to extract the lateral diffusion coefficient in lipid bilayers using excimer formation. In contrast to previous statistical models that modeled the system as points undergoing jumps from site to site on a lattice, we use coarse-grained molecular dynamics to study lipid bilayers simulated using the Martini force field. We derive time dependent reaction rates from survival probabilities obtained a posteriori from numerically generated trajectories of symmetric DOPC (1,2-Dioleoyl-sn-glycero-3-phosphocholine) and POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) bilayers at 283K and 293K respectively. Collision dynamics are determined by virtually relabeling the simulated molecules. The fluorescent probes are assumed to behave like ordinary membrane lipids and therefore the dynamics remain unaffected. We derive a generalized expression for the survival probability combining independent pairs and size scaling assumptions, but no assumption is made regarding the kinetic rate of the excimer formation process. By fitting the numerically determined normalized fluorescence emission intensities to experimental titration curves, we obtain two sets of results for the lateral diffusion coefficients depending whether interleaflet excimer association is allowed or not. We use a capture radius of 0.5 nm, the distance at which the probes react to form excimers. By relating Martini dynamics to real fluorescence experiments, we estimate the numerical Martini acceleration factor. We also study mixtures of oxidized-non oxidized DOPC and POPC bilayers using a hydroperoxidized model of these lipids for different concentrations of the oxidized component (3.1%, 25% and 50%). Using pair correlation functions, we extract structural information on the systems and determine whether the two components are prone to mixing or not. Finally, we calculate the thermodynamic mixing parameters within the framework of the virial expansion.