Simulation de modélisation moléculaire de l’interaction HDP-2P avec l’IP et les phosphoinositides
Le programme de modélisation moléculaire calcule les neuf meilleurs sites de liaison sur la surface moléculaire d’un récepteur (molécule plus grande), qui est ciblé par un ligand (molécule plus petite). Ces sites sont déterminés par la libération des neuf enthalpies les plus élevées, qui désignent les sites présentant les meilleures affinités de liaison. Les affinités de liaison sont exprimées sous forme de valeurs négatives, ce qui indique que la liaison intermoléculaire est un processus exothermique. Dans le cas de l’interaction de HDP-2P avec l’IP, la simulation a prédit quatre sites de liaison (#3, #6, #7 et #8) pour l’IP sur la surface moléculaire de l’HDP-2P dans le centre actif (Tableau 1). De plus, les cinq sites de liaison pour l’IP ont été trouvés en dehors du centre actif de HDP-2P. Les notations des fractions chargées et polaires dans la tête polaire de PI sont illustrées à la figure 1.
Le centre actif du HDP-2P se compose de cinq résidus d’acides aminés clés : Gly30, Gly32, His48, Asp49 et Lys69, qui facilitent l’alignement correct du substrat phospholipidique dans le centre actif pour assurer l’hydrolyse catalytique. Certaines structures arrimées représentatives montrent la liaison de PI à certains résidus d’acides aminés du centre actif, tandis que dans le site de liaison #6, PI s’engage avec les cinq résidus d’acides aminés clés du centre actif par des liaisons ioniques, ionopolaires et hydrogène (Tableau 1). La figure 2 illustre comment PI interagit avec les cinq résidus d’acides aminés clés dans le centre actif de HDP-2P dans le site de liaison #6.
La simulation d’amarrage de l’interaction entre HDP-2P et PI-4-P, qui contient un groupe phosphate à la position 4 du cycle inositol, a prédit uniquement le site de liaison #1 comme l’endroit où PI-4-P se lie au centre actif de HDP-2P. Plus précisément, dans ce site, PI-4-P interagit avec trois résidus d’acides aminés clés du centre actif - Gly30, Gly32 et Lys69 - via des liaisons ioniques, ionopolaires et hydrogène (Tableau 2). La figure 3 illustre l’interaction de PI-4-P avec le centre actif de HDP-2P. Dans les huit sites de liaison restants à la surface moléculaire de HDP-2P, PI-4-P se lie à des emplacements en dehors du centre actif.
Les simulations d’amarrage moléculaire des interactions HDP-2P avec PI-4,5-P₂ (contenant des groupes phosphate aux positions 4 et 5 du cycle inositol) n’ont révélé aucune liaison au centre actif de l’enzyme sur les neuf sites de liaison prédits. Au lieu de cela, PI-4,5-P₂ se lie constamment à une cavité en forme de rainure adjacente au centre actif (Figure 4), formant des liaisons ion-polaires et hydrogène avec les trois résidus HDP-2P, Trp31, Gly33 et Gln34).
Bien que la liaison de PI-4,5-P₂ à une cavité en forme de sillon près du centre actif du PLA2cellulaire ou du venin n’ait pas été documentée auparavant, des caractéristiques structurellement analogues sont bien établies dans les domaines de l’homologie31 (PH) et de l’homologie N-terminale32 (ENTH) de la pleckstrine (PH). Ces domaines de signalisation utilisent des architectures de sillons hydrophobes-polaires-cationiques similaires pour la reconnaissance spécifique de PI-4,5-P2 31,32, malgré leur nature non enzymatique. La présence d’un motif structurel comparable dans HDP-2P soulève des possibilités intrigantes de mimétisme évolutif ou d’évolution convergente pour faciliter la reconnaissance du domaine membranaire par HDP-2P. Si elle est vérifiée par d’autres méthodes, cette étude pourrait représenter le premier cas signalé d’un PLA2 à venin utilisant un tel mécanisme de reconnaissance du phosphoinositide.
Il est intéressant de noter que les simulations d’amarrage moléculaire de CL et PS interagissant avec HDP-2P ont prédit la liaison des deux phospholipides au centre actif de HDP-2P dans huit des neuf sites de liaison, tandis que le seul site de liaison de CL et PS à l’extérieur du centre actif n’était pas situé dans la cavité en forme de sillon (données non présentées). Il convient également de noter que, comme PI, le PC se lie au centre actif du HDP-2P dans quatre des neuf sites de liaison, et que le PI et le PC, comme CL et PS, ne se lient pas à la cavité en forme de rainure (données non présentées). Ces résultats suggèrent que les groupes phosphate sur le cycle inositol de PI ont une affinité de liaison élevée avec la cavité en forme de sillon située à côté du centre actif. Les analyses structurales ont identifié des cavités analogues en forme de sillon avec des taches de surface hydrophobes-polaires-cationiques près du centre catalytique dans les enzymes cellulaires33 et34 PLA2 du venin. La liaison des phosphoinositides à ces cavités peut induire des changements conformationnels physiologiquement significatifs, inhibant potentiellement l’activité enzymatique. Notamment, parmi les phosphoinositides, PI-4,5-P2 présente des caractéristiques de substrat particulièrement médiocres - les isoformes humaines du PLA2 (types cytosoliques, indépendants du calcium et sécrétés) montrent une activité hydrolytique minimale à indétectable envers PI-4,5-P235. Cela suggère que la liaison du phosphoinositide à ces caractéristiques structurelles conservées peut remplir une fonction physiologique distincte (par exemple, l’inhibition allostérique) plutôt que des fonctions catalytiques dans la superfamille du PLA2 , ce qui nécessite une vérification expérimentale future.
Activité hydrolytique de HDP-2P dans les liposomes PC enrichis en PI, PI-4-P, PI-4,5-P2, PS ou CL surveillés à l’aide de la coenzyme A acylée avec des acides gras libres
L’activité hydrolytique de HDP-2P sur des échantillons liposomaux a été évaluée à l’aide d’une solution de coenzyme A. La solution de coenzyme A est incolore mais devient violette lors de l’acylation par les acides gras libres libérés lors de l’hydrolyse des phospholipides catalysée par HDP-2P. L’augmentation de la densité optique qui en résulte, correspondant à l’ampleur de l’acylation de la coenzyme A, a été mesurée par spectrophotométrie à 550 nm. Par conséquent, l’augmentation de la densité optique dans les échantillons liposomaux traités avec de la coenzyme A et du HDP-2P est directement proportionnelle à l’activité hydrolytique du HDP-2P.
La figure 5 présente les courbes d’activité hydrolytique en fonction du rapport molaire HDP-2P/phospholipides dans divers échantillons liposomaux. L’activité hydrolytique du HDP-2P est exprimée en unités arbitraires (u.a.) de densité optique (DO) mesurées à 550 nm. Nous n’avons pas construit de courbe d’étalonnage pour relier les valeurs de l’UA aux concentrations molaires d’acides gras libres, car nos essais expérimentaux ont montré que les acides gras dans les solvants organiques, lorsqu’ils sont ajoutés à des solutions aqueuses contenant de la coenzyme A et de l’acyl-coenzyme A synthétase, ne produisent aucun changement détectable de DO à 550 nm. Cela indique que ces acides gras forment des micelles plutôt que de se lier à la coenzyme A. En revanche, seuls les acides gras générés par la membrane (produits par l’activité du PLA2 ) se lient à la coenzyme A associée à la membrane, ce qui entraîne un changement mesurable de la DO La nature amphipathique de la coenzyme A - ayant à la fois des régions hydrophobes et hydrophiles à sa surface moléculaire - la localise préférentiellement à l’interface membrane-eau, ce qui facilite la liaison des acides gras produits par l’activité du PLA2 à la coenzyme A. Cependant, dans ces conditions d’essai, il n’est pas possible de déterminer la concentration molaire des acides gras. Pour cette raison, les chercheurs utilisant le dosage de la coenzyme A expriment l’activité hydrolytique du PLA2 en unités arbitraires de DO mesurées à 550 nm20.
La plus forte augmentation de la densité optique, correspondant à la plus grande activité hydrolytique de HDP-2P, a été observée dans les liposomes PC enrichis en CL ou PS. Ce résultat s’aligne bien avec le plus grand nombre de sites de liaison pour CL et PS au centre actif de HDP-2P. L’activité hydrolytique de HDP-2P dans les liposomes PC enrichis en PI était 3,0 et 2,5 fois plus faible que dans les liposomes enrichis en CL ou PS, respectivement, mais plus élevée que dans les liposomes PC purs. Cela suggère que l’affinité de liaison au centre actif de HDP-2P de PI est plus élevée que celle de PC, bien que les deux lipides aient le même nombre de sites de liaison au centre actif.
Cet écart apparent peut s’expliquer par des différences structurelles dans les groupes de tête polaire de PI et de PC. Alors que le PC est une molécule neutre, son groupe cationique triméthylammonium s’étend vers l’extérieur dans la solution, tandis que son groupe phosphate anionique est positionné plus près de la zone hydrophobe de la membrane. Cette distribution de charge dans le groupe de tête PC peut créer une répulsion électrostatique contre le HDP-2P de base, réduisant ainsi son affinité de liaison. En revanche, la nature acide de l’IP en fait une espèce lipidique plus attrayante pour l’interaction avec le HDP-2P de base, ce qui pourrait expliquer son activité hydrolytique plus élevée par rapport aux liposomes PC purs.
L’activité hydrolytique de HDP-2P dans les liposomes PC enrichis en PI-4-P, comme le montre la figure 5, suggère une hydrolyse catalytique minimale des phospholipides. Remarquablement, la courbe de densité optique obtenue à partir de liposomes PC enrichis en PI-4,5-P indique une abolition complète de l’activité hydrolytique par HDP-2P dans cet échantillon de liposomes. Ces résultats s’alignent bien avec les simulations d’amarrage, qui prédisent un site de liaison pour PI-4-P au centre actif de HDP-2P, tandis que PI-4,5-P2 ne montre aucune affinité de liaison au centre actif. Cela suggère que les groupes phosphate liés au cycle inositol entravent la liaison du phospholipide au centre actif de HDP-2P. Nous proposons que les groupes phosphate sur le cycle inositol s’étendent vers l’extérieur dans la solution, et en se liant à HDP-2P via des interactions ioniques et iono-polaires intermoléculaires, ils maintiennent HDP-2P à une distance de la surface de la membrane, protégeant ainsi la membrane de l’activité hydrolytique de HDP-2P. Cela est cohérent avec les découvertes récentes révélant que les enzymes PLA2 humaines cytosoliques, indépendantes du calcium et sécrétées ont montré une activité hydrolytique à peine détectable ou nulle vers les substrats PI-4-P et PI-4,5-P2 35.
Modifications de perméabilité induites par HDP-2P dans les membranes liposomales PC enrichies en PI, PI-4-P, PI-4,5-P2, PS ou CL : Évaluation par réaction de substitution de ligand dans une solution de CuSO4 avec NH₃ encapsulé
L’impact de l’activité hydrolytique HDP-2P sur la perméabilité membranaire des liposomes PC enrichis en PI, PI-4-P, PI-4,5-P2, PS ou CL a été étudié à l’aide d’une réaction de substitution de ligand dans un système d’ions complexes, où le ligand de substitution a été encapsulé dans le volume interne des liposomes. Cette méthode simple et fiable est basée sur le déplacement de quatre molécules d’eau agissant comme des ligands dans l’ion complexe [Cu(H2O)6]2+ par NH3
[Cu(H2O)6]2+ (aq) + 4NH3 (aq)
[Cu(H2O)2(NH3)4]2+ (aq) + 4H2O (l)
Au cours du processus de substitution des molécules d’eau par NH3 dans l’ion complexe [Cu(H2O)6]2+, la couleur bleu pâle de [Cu(H2O)6]2+ passe à la couleur bleu foncé de [Cu(H2O)2(NH3)4]2+. Dans un système où NH3 est encapsulé dans le volume interne des liposomes et où [Cu(H2O)6]2+ est présent dans la solution externe, ce changement de couleur -- contrôlé spectrophotométriquement -- indique une augmentation de la perméabilité de la membrane liposomale à NH3.
Pour établir un tel système, nous avons soniquisé une dispersion phospholipidique dans un tampon Tris-HCl contenant NH3, formant des liposomes sonicés avec NH₃ à l’intérieur et à l’extérieur des liposomes. Le NH3 externe a été éliminé par dialyse, après quoi la solution de CuSO4 a été ajoutée à l’échantillon de liposomes contenant du NH3 dans son volume interne. Nous avons ensuite traité les liposomes avec du HDP-2P et surveillé les changements de densité optique par spectrophotométrie.
Pour déterminer la longueur d’onde optimale pour mesurer les changements de densité optique, nous avons enregistré les spectres d’absorbance de [Cu(H2O)2(NH3)4]2+ dans une solution CuSO4 0,025 M, NH3 0,100 M et [Cu(H2O)6]2+ dans une solution CuSO4 0,025 M, en utilisant la lumière visible dans la gamme de 500 nm à 900 nm. Comme le montre la figure 6, l’irradiation à 600 nm a produit une absorbance maximale pour le [Cu(H2O)2(NH3)4]2+ et une absorbance minimale pour le [Cu(H2O)6]2+, ce qui nous a amenés à sélectionner 600 nm comme longueur d’onde optimale pour nos expériences de perméabilité membranaire.
La figure 7 montre que la perméabilité membranaire à NH₃ est la plus élevée dans les liposomes PC enrichis en CL et PS anioniques. Cette observation est cohérente avec nos données précédemment rapportées montrant que HDP-2P, un PLA2 du groupe IIA, interagit efficacement avec le phospholipide anionique palmitoyloléoyl-glycérophosphoglycérol, conduisant à l’hydrolyse du substrat - un effet non observé avec le phospholipide neutre palmitoyloléoylglycérophosphocholine36. En général, la plupart des enzymes PLA2 du groupe IIA présentent une activité hydrolytique significativement plus élevée envers les phospholipides anioniques 33,37,38. La perméabilité des liposomes PC+PI est inférieure à celle des liposomes PC+CL et PC+PS mais reste supérieure à celle des liposomes PC purs. Ces résultats concordent bien avec les résultats des simulations d’amarrage et des analyses d’activité hydrolytique HDP-2P. En d’autres termes, les échantillons lyposomaux qui présentent une faible perméabilité devraient être structurellement stables en raison de l’incapacité du HDP-2P à hydrolyser les phospholipides et en raison d’une plus faible affinité du HDP-2 à se lier au PI et aux phosphoinositides, comme prédit par les données d’amarrage moléculaire (Figure 2, Figure 3 et Figure 4). Par conséquent, l’activité hydrolytique enzymatique par HDP-2P (Figure 5) et la perméabilité membranaire, comme le montre la Figure 7, devraient être positivement corrélées.
Cependant, dans les liposomes PC+PI-4-P et PC+PI-4,5-P2 , la perméabilité n’a pas augmenté lors du traitement par HDP-2P. Bien que les résultats de perméabilité dans les liposomes PC+PI-4,5-P2 soient bien corrélés avec les simulations d’amarrage et les données d’activité hydrolytique HDP-2P pour les mêmes liposomes, l’absence d’augmentation de la perméabilité dans les liposomes PC+PI-4-P ne correspond pas entièrement aux résultats de l’activité hydrolytique, car un faible degré d’activité hydrolytique a été enregistré dans ces liposomes. Ces résultats suggèrent qu’une légère augmentation localisée de l’activité hydrolytique dans les liposomes PC+PI-4-P ne perturbe pas l’empilement bicouche des phospholipides à un degré qui modifierait la perméabilité de la membrane à NH3.

Figure 1 : Notations des fractions chargées et polaires de la tête polaire du phosphatidylinositol. Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 2 : Interaction de PI avec le centre actif de HDP-2P prédite pour le site de liaison #6 par un logiciel d’amarrage moléculaire. PI est donné en représentation bâton et HDP-2P est donné en représentations linéaires (schéma à gauche) et en surface moléculaire (schéma à droite). Pour les représentations en bâton, l’émeraude représente le carbone, l’orange le phosphore, le rouge l’oxygène et le blanc l’hydrogène. Pour la surface moléculaire et les lignes, le vert représente le carbone, le bleu l’azote, le rouge l’oxygène, le blanc l’hydrogène et le jaune le soufre. Les lignes brisées jaunes identifient les liaisons intermoléculaires. Les flèches pointent vers les résidus d’acides aminés au centre actif. Abréviations : PI = Phosphatidylinositol ; HDP-2P = sous-unité de base de la phospholipase de venin de vipère. Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 3 : Interaction de PI-4-P avec le centre actif de HDP-2P prédite pour le site de liaison #1 par un logiciel d’amarrage moléculaire. Le PI-4-P est donné en représentation en bâton et le HDP-2P est donné en lignes (schéma de gauche) et en surface moléculaire (schéma de droite). Pour les bâtonnets, l’émeraude représente le carbone, l’orange est le phosphore, le rouge est l’oxygène et le blanc est l’hydrogène. Pour la surface moléculaire et les lignes, le vert représente le carbone, le bleu l’azote, le rouge l’oxygène, le blanc l’hydrogène et le jaune le soufre. Les lignes brisées jaunes identifient les liaisons intermoléculaires. Les flèches pointent vers les résidus d’acides aminés au centre actif. Abréviations : PI-4-P = Phosphatidylinositol-4-phosphate ; HDP-2P = sous-unité de base de la phospholipase de venin de vipère. Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 4 : Interaction de PI-4,5-P2 avec la surface moléculaire de HDP-2P prédite pour le site de liaison #6 par un logiciel d’amarrage moléculaire. PI-4,5-P2 se lie à une grotte en forme de rainure située sous le centre actif du HDP-2P. Les cinq résidus d’acides aminés clés au centre actif sont donnés dans la notation de trois lettres. PI-4,5-P2 est donné en représentation bâton et HDP-2P est donné en représentations linéaires (schéma à gauche) et en surface moléculaire (schéma à droite). Pour les bâtonnets, l’émeraude représente le carbone, l’orange est le phosphore, le rouge est l’oxygène et le blanc est l’hydrogène. Pour la surface moléculaire et les lignes, le vert représente le carbone, le bleu l’azote, le rouge l’oxygène, le blanc l’hydrogène et le jaune le soufre. Abréviations : PI-4,5-P2 = phosphatidylinositol-4,5-diphosphate ; HDP-2P = sous-unité de base de la phospholipase de venin de vipère. Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 5 : Courbes d’activité hydrolytique à partir d’échantillons de liposomes à divers rapports molaires HDP-2P/phospholipides. Liposomes en PC+CL (vert foncé), PC+PS (orange), PC+PI (bleu foncé), PC+PI-4-P (bleu), PC+PI-4,5-P2 (vert) et PC (violet). Les valeurs de densité optique en unités arbitraires représentent les moyennes des lectures triples à 550 nm. Les barres d’erreur indiquent l’écart-type des lectures en trois exemplaires. La signification statistique des différences de valeurs de densité optique entre les échantillons de liposomes a été évaluée à l’aide de l’ANOVA, toutes les valeurs p étant bien inférieures à 0,05. Abréviations : PC = phosphatidylcholine ; CL = cardiolipin cardiaque bovin ; PS = phosphatidylsérine ; PI = Phosphatidylinositol ; PI-4-P = Phosphatidylinositol-4-phosphate; PI-4,5-P2 = Phosphatidylinositol-4,5-diphosphate ; a.u. = unités arbitraires ; HDP-2P = sous-unité de base de la phospholipase de venin de vipère. Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 6 : Spectres d’absorbance de solutions aqueuses de complexes de cuivre. Complexes de cuivre [Cu(H2O)6]2+ (ligne orange) et [Cu(H2O)2(NH3)4]2+ (ligne bleue), irradiés par la lumière visible dans la gamme de 500 nm à 900 nm. Le complexe de cuivre [Cu(H2O)6]2+ a été réalisé dans une solution de 0,025 M CuSO4 et le complexe de cuivre [Cu(H2O)2(NH3)4]2+ a été réalisé dans une solution de 0,025 M CuSO4 contenant 0,100 M NH3. Les valeurs d’absorbance dans cette figure pour chaque point de données sont les moyennes des lectures d’absorbance effectuées en trois exemplaires. Les barres d’erreur représentent l’écart-type des lectures en trois exemplaires. Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 7 : Courbes de perméabilité membranaire à partir d’échantillons de liposomes. Liposomes en PC+CL (violet), PC+PS (bleu), PC+PI (bleu foncé), PC+PI-4-P (vert foncé), PC+PI-4,5-P2 (orange) et PC (vert) à différents rapports molaires HDP-2P/phospholipides. Les valeurs de perméabilité de la membrane sont exprimées en valeurs de densité optique en unités arbitraires (a.u.) mesurées à 600 nm. Chaque point de données représente la moyenne des lectures en trois exemplaires. Les barres d’erreur indiquent l’écart-type des lectures en trois exemplaires. La signification statistique des différences de valeurs de densité optique entre les échantillons de liposomes a été évaluée à l’aide de l’ANOVA, toutes les valeurs p étant bien inférieures à 0,05. Veuillez cliquer ici pour voir une version agrandie de cette figure.
Tableau 1 : Énergies de liaison d’affinité prédites pour les complexes moléculaires candidats de HDP-2P avec PI. Des simulations de modélisation moléculaire ont été utilisées pour prédire les types de liaisons et les affinités de site de liaison, ainsi que les fractions chargées et polaires du groupe de tête de phosphatidylinositol et des résidus d’acides aminés HDP2P qui participent à leurs interactions intermoléculaires. Les groupes chargés et polaires du phosphatidylinositol sont mis en évidence à la figure 1. Notez que Pb dans NHpbδ+ désigne une liaison peptidique. Abréviations : PI = Phosphatidylinositol ; HDP-2P = sous-unité de base de la phospholipase de venin de vipère. Veuillez cliquer ici pour télécharger ce tableau.
Tableau 2 : Énergies de liaison d’affinité prédites pour les complexes moléculaires candidats de HDP-2P avec PI-4-P. Des simulations de modélisation moléculaire ont été utilisées pour prédire les types de liaisons et d’affinités entre les sites de liaison, ainsi que les fractions chargées et polaires du groupe de tête du phosphatidyl-inositol-4-phosphate (PI-4-P) et des résidus d’acides aminés HDP2P qui participent à leurs interactions intermoléculaires. Notez que Pb dans NHpbδ+ désigne une liaison peptidique. Veuillez cliquer ici pour télécharger ce tableau.
Fichier supplémentaire 1 : Flux de travail d’amarrage moléculaire. Vue d’ensemble de la procédure d’amarrage moléculaire à l’aide d’AutoDock, y compris la construction des ligands, la préparation des récepteurs, la configuration de la boîte de grille et l’exécution de simulations d’amarrage pour prédire les interactions PI-HDP-2P. Veuillez cliquer ici pour télécharger ce fichier.