Mutation-induced conformational remodeling explains gain- and loss-of-function phenotypes in Euglena gracilis phytoene synthase

Résumé technique : Remodelage conformationnel induit par des mutations dans la phytoène synthase d'Euglena gracilis

Énoncé du problème
La phytoène synthase (CrtB) catalyse la première étape limitante de la biosynthèse des caroténoïdes, ce qui en fait une cible critique pour l'ingénierie métabolique. Bien que le criblage expérimental chez Euglena gracilis ait identifié des mutations spécifiques (Y197A, E121H et L349F) qui entraînent des phénotypes contrastés de gain de fonction (production de caroténoïdes accrue) et de perte de fonction (production réduite), les mécanismes structurels pilotant ces différences restent obscurs. Les modèles structurels statiques et les observations expérimentales seuls ne parviennent pas à expliquer comment des substitutions de résidus uniques altèrent l'activité enzymatique. L'étude aborde le fossé dans la compréhension de la base structurelle dynamique reliant les mutations spécifiques à l'altération de la productivité en caroténoïdes.

Méthodologie
Les auteurs ont employé des simulations de dynamique moléculaire (MD) de la forme apo (sans substrat) de la CrtB d'E. gracilis pour étudier la dynamique conformationnelle de l'enzyme sauvage et de trois variantes : Y197A, E121H et L349F. L'approche computationnelle comprenait :

• Analyse énergétique : Des calculs MM-GBSA et MM-PBSA pour évaluer les perturbations énergétiques dépendantes des mutations et la stabilité ($\Delta\Delta G$).
• Caractérisation électrostatique : L'analyse PROPKA pour déterminer les décalages des points isoélectriques (pI) et de la charge nette à pH 7,0.
• Métriques structurelles globales : L'analyse de la surface accessible au solvant (SASA) et du rayon de giration ($R_g$) pour évaluer la compacité globale et l'échantillonnage conformationnel.
• Dynamique locale : L'analyse de la fluctuation de la racine carrée moyenne (RMSF) et des $\Delta$RMSF pour cartographier les changements de mobilité spécifiques aux résidus.
• Architecture du site actif : Un examen détaillé de l'accessibilité au solvant, des volumes de la poche de liaison et des réseaux de liaisons hydrogène au sein des motifs riches en Asp (essentiels pour la coordination du Mg²⁺) et du volet hydrophobe (hydrophobic flap).

Résultats clés
Les simulations ont révélé que les mutations induisent un remodelage distinct du paysage énergétique conformationnel, affectant les réseaux d'interactions locales et l'architecture du site actif sans provoquer de changements structurels globaux substantiels :

1. Décalages énergétiques et électrostatiques :

   • E121H a présenté la plus grande redistribution de la dynamique conformationnelle, avec des valeurs de $\Delta\Delta G$ positives indiquant une déstabilisation énergétique par rapport à la protéine native. Elle a également montré le pI prédit le plus élevé et une légère charge positive à pH neutre.
   • Y197A a affiché des changements énergétiques quasi neutres ou favorables.
   • L349F a induit des perturbations énergétiques modestes mais a montré une déstabilisation prononcée dans les prédictions de l'énergie libre de repliement PROPKA.

2. Échantillonnage conformationnel et flexibilité :

   • Variantes de gain de fonction (Y197A, E121H) : Toutes deux ont conservé des volets hydrophobes relativement accessibles. E121H a montré la distribution de $R_g$ la plus large, tandis que Y197A et E121H ont présenté des distributions de SASA plus larges que la protéine native, suggérant un échantillonnage conformationnel accru tout en maintenant une compacité globale.
   • Variante de perte de fonction (L349F) : A étroitement ressemblé à la protéine native dans les distributions globales de SASA et de $R_g$, mais a affiché une constriction localisée dramatique de la poche du volet hydrophobe (réduction de volume d'environ 35 %).

3. Réseaux d'interactions locales et remodelage du site actif :

   • Volet hydrophobe : La mutation L349F a provoqué une réduction significative du volume de la poche du volet et de l'accessibilité au solvant, ce qui est cohérent avec une conformation plus fermée. En revanche, les variantes de gain de fonction ont maintenu une plus grande exposition du volet.
   • Motifs riches en Asp : Tous les mutants ont montré une augmentation de l'occupation des liaisons hydrogène dans la région Asp-riche 1. Notamment, L349F a affiché l'occupation la plus élevée malgré une performance catalytique réduite, suggérant que l'efficacité catalytique repose sur un équilibre entre stabilité et flexibilité plutôt que sur une stabilisation locale maximale.
   • Rigidification dans L349F : L'analyse structurelle du mutant L349F a révélé que le résidu Tyr127 a adopté une conformation plus rigide, formant un réseau d'interaction persistant impliquant des interactions $\pi$-cation et $\pi$-anion avec Arg154 et Glu158. Cela contraste avec la protéine native et les variantes de gain de fonction, où Tyr127 est restée conformationnellement flexible.

Signification et affirmations
L'article affirme fournir un cadre structurel reliant le remodelage conformationnel induit par les mutations autour du centre catalytique aux différences de production de caroténoïdes observées expérimentalement. Les auteurs soutiennent que les phénotypes contrastés ne proviennent pas simplement de perturbations structurelles locales au niveau des sites de mutation, mais de la redistribution de l'échantillonnage conformationnel autour de motifs fonctionnellement importants, spécifiquement le volet hydrophobe et le motif Asp-riche 1 adjacent.

L'étude souligne que le phénotype de perte de fonction dans L349F est caractérisé par une « constriction prononcée » du volet hydrophobe et un « réseau d'interaction rigide », tandis que les variantes de gain de fonction conservent un « volet hydrophobe relativement accessible » et une plus grande flexibilité. Ces découvertes suggèrent que l'organisation structurelle dynamique est un déterminant critique de la fonction de CrtB. Les auteurs concluent que ces informations offrent une base pour l'ingénierie rationnelle des phytoène synthases, bien qu'ils notent explicitement que des simulations futures intégrant les substrats (GGPP) et les ions Mg²⁺ sont nécessaires pour élucider pleinement comment ces changements conformationnels intrinsèques influencent la reconnaissance du substrat et la catalyse.