Stéphane Abel , Massimo Marchi CEA, DSV/iBiTeC-S/SB2SM/LBMS, Saclay

1. Simulations de dynamique moléculaires des détergents: Applications à l’étude des micelles directes d’acide gras et de glycolipides Stéphane Abel, Massimo Marchi CEA, DSV/iBiTeC-S/SB2SM/LBMS, Saclay GDR 2748 Protéines membranaires. Les Houches 9 - 13 Mars 2009

2. Avant propos • Les simulations de dynamique moléculaires sont des outils de choix pour étudier la structure des détergents en solution (et notamment des micelles) : plus de 500 articles depuis la fin des années 80 (source web of science). • Deux exemples pour ce GDR de simulations de dynamique moléculaires appliquées à l’étude structurale : • Des micelles d’acide linoléique utilisées dans les études de radiolyse: Examen de l’influence des conditions initiales. • Des micelles de dodecylmaltoside utilisées dans la solubilisation douce des protéines membranaires: Examen de l’influence de la conformation des têtes polaires sur la structure des micelles

3. Les micelles d’acide linoléique • En bref: • L’acide linoléique (AL) est un acide gras -6 avec 18 carbones • Sphériques en solution aqueuse (cmc = 2mM à pH=11.5) de dimension RM≈ 22 - 24 Å • (cryo-TEM ou SANS) • Template pour la formation de nanoparticules métalliques • par radiolyse. • Peu d’autres données structurales pour modéliser les micelles d’AL (et notamment Nagg) • Utilisation d’un modèle géométrique pour le choix de Nagg et de RM ConsidérantquexCH2*Nagg (avec 0 < xCH2 < 17) ,Sh=21 Å2 etL = 23 Å Volume et surface de la micelle : Fraction de surface hydrophile de la micelle: Différentes simulations avec Nagg = 50, 60, 75, 90, meilleurs accords pour Nagg = 60, RM = 19.4 Å et f = 26.6 %

4. Détails des simulations • Deux protocoles de simulation: préformée et self-assemblée • 60 molécules de LIN, 60 Na+ et 9997 mol. d’eau TIP3 pour [LIN]=0.3 M • Système préformée (M60-P): Micelle manuellement construite est placée dans une boite d’eau octaédrique tronquée (a=b=c=75 Å and =β=γ=109.472o) • Système self-assemblée (M60-S): les molécules de LIN et de Na+ sont placées aléatoirement dans une boite cubique d’eau avec a = 69 Å Simulations avec le code ORAC dans les conditions NPT (P= 0. 1 MPa et T= 300 K) avec prise en compte de l’électrostatique et des conditions périodiques. Champs de CHARMM27 avec les paramètres de torsions en cis pour les molécules de LIN. • Association des molécules d’AL en micelles (M60-S) Processus d’agrégation en 2 étapes rapides et lentes.

5. Comparaison entre les micelles préformées et self-agrégées Self-assembled (M60-S)t = 31 ns Preformed (M60-P) t = 10 ns M60-P M60-S f : fraction de surface hydrophile de la micelle sur la surface totale. nhg-OW: nombre moyen d’hydratation.

6. Les alkyls glycosides En bref: Présent dans les membranes cellulaires Biodégradables et non toxiques. S’assemblent en différentes structures (micelles, membranes) Utilisation courante pour l’extraction et la solubilisation des protéines membranaires permettant de garder intact leur activité fonctionnelle. β-Octyl-Glucoside (OG) β-Dodecyl-Maltoside (β-DDM) Tiré de Walian et al. (2004)

7.   Le dodecylmaltoside (DDM) • A deux anomères ( et ) avec des propriétés différentes: • CMC: 1.5 x 10-4 mol.l-1 () et 2.0 x 10-4 mol.l-1 () • Petites micelles sphériques avec la forme  (NDDM75-80) • alors que  forme des micelles oblates(NDDM100-135) • Micelles quasi-sphériques pour  et oblates pour  • Peu d’informations sur la structure des micelles en solution aqueuse (essentiellement SANS et SAXS). • Motivations de l’étude • Construction d’un model moléculaire du détergent DDM pour • simuler les interactions proteine-peptide/glycolipides. • Examen l’influence de la conformation de la tête polaire sur la structure des micelles.

8. DM de micelles avec les deux anomères de DDM • Les simulations • 2 tailles de micelles pré-assemblées avec nombre d’agrégation (NDDM) issus de SAXS et SANS (Dupuy et al. 1997) • Champs de force tirés de CHARMM (chaine alkyl et tête polaire maltose) • Modèle d’eau TIP3 • Conditions périodiques, SPME • Simulations NPT (P = 0.1 MPa et T = 297 K) avec le code ORAC tsim=240 fs tsim=7 ns tsim=14 ns -DDM -DDM

9. Dimensions des micelles <Rg> =26.5 ± 0.1 Å Rgexp = 23.5 ± 1 Å β-DDM Réarrangement des glycolipides -DDM <Rg> =20.5 ± 0.1 Å Rgexp = 18.5 ± 1 Å

10. Paramètres de forme des micelles at acec at ac ac atet <lpla, b, c> -DDM = 7.6 Å lpla, b, c = Micelle a, b, c- Core a, b, c <lpla, b, c> -DDM = 7.2 Å a Dupuy et al., 1997

11. Structure internes des micelles -DDM β-DDM

12. Hydratation des têtes polaires Mol. H2O d’hydratation unique Nwall Exp  8-14 Nwhead Exp  8 Warr, 1986; Cecutti et al. 1991 et Dupuy et al. 1997

13. Diffusion de translation de l’eau interfaciale -DDM β-DDM Bulk TIP3 wb

14. Bibliothèque de paramètres pour les glycolipides à base de glucose (en coll. F.Y Dupradeau, UMR6219, Université d’Amiens) Charges partielles pour 24 glycolipides courants compatibles avec GLYCAM/AMBER http://q4md-forcefieldtools.org/REDDB/index.php

15. Merci pour votre attention !!!