Time-Resolved Phosphoproteomics-Guided BFS Beam Search Reveals Cell-Type-Specific EGFR Signaling Architectures and SHP2 Inhibitor-Induced Pathway Rewiring

Cette étude présente un cadre computationnel systématique qui intègre la phosphoprotéomique résolue dans le temps avec un algorithme de recherche par faisceau guidé par une recherche en largeur (BFS) pour reconstruire des réseaux de signalisation EGFR spécifiques aux types cellulaires, révélant avec succès comment l'inhibition de SHP2 reconfigure les architectures de voies et entraîne des mécanismes de résistance adaptative.

L'essentiel

Imaginez votre cellule comme une ville animée et high-tech. À l'intérieur de cette ville, des millions de petits messagers (des protéines) courent partout, se passant des notes pour indiquer à la ville comment réagir à des éléments tels que des signaux de croissance ou des médicaments. Lorsqu'un médicament tente d'arrêter un mauvais acteur spécifique dans la ville, les messagers deviennent souvent astucieux : ils réacheminent leurs messages par des rues différentes pour maintenir la ville en fonctionnement de toute façon. C'est ce que les scientifiques appellent la « résistance adaptative ».

Le problème est que nous disposons d'une quantité massive de données sur ces messagers (phosphoprotéomique), mais c'est comme avoir un million de cartes postales éparpillées sans carte. Nous savons ce qui a été envoyé, mais nous ne savons pas exactement quelles rues les messages ont empruntées pour arriver à destination.

Le nouveau « GPS » pour les villes cellulaires
Cet article présente un nouveau programme informatique qui agit comme un GPS intelligent pour déterminer exactement quelles routes ces messages empruntent. Voici comment cela fonctionne, en utilisant des analogies simples :

• La Carte (base de données STRING) : Au lieu de construire une gigantesque bibliothèque de toutes les routes possibles au monde, le programme se connecte à une carte en ligne et en direct (la base de données STRING) qui sait déjà quelles protéines sont amies les unes des autres. Il n'a pas besoin de stocker toute la carte sur son propre disque dur ; il consulte simplement les connexions en temps réel.
• La Stratégie de recherche (BFS + recherche par faisceau) : Imaginez que vous essayez de trouver un chemin depuis la porte de la ville (le début d'un signal) jusqu'au bureau du maire (l'effet final).
  • Le programme utilise une recherche en largeur (BFS) pour examiner toutes les routes possibles simultanément, comme un drone survolant toutes les rues d'un quartier en même temps.
  • Cependant, examiner toutes les routes possibles dans toute la ville prendrait une éternité. Il utilise donc une recherche par faisceau. Imaginez cela comme une lampe torche qui n'éclaire que les 5 routes les plus prometteuses à tout moment, en ignorant les impasses. Cela maintient la recherche focalisée et rapide.
• Filtrage du bruit (MAD et nettoyage) : Toutes les notes trouvées dans la rue ne sont pas importantes. Le programme utilise un filtre statistique (MAD) pour décider quelles notes sont de véritables signaux et lesquelles ne sont que du bruit de fond. Après avoir trouvé toutes les routes possibles, il lance une « équipe de nettoyage » pour éliminer les boucles (où un message tourne en rond) et consulte un annuaire local (Human Protein Atlas) pour s'assurer que les bâtiments sur la route existent réellement dans ce type spécifique de cellule.

Ce qu'ils ont découvert
Les chercheurs ont testé ce GPS sur trois types différents de « villes » (cellules HeLa, MDA-MB-468 et HEK293T). Ils ont constaté que chaque ville possède une disposition unique ; ce qui fonctionne dans l'une ne fonctionne pas dans l'autre.

Ils ont spécifiquement examiné ce qui se passe lorsqu'ils tentent de bloquer un policier de la circulation spécifique nommé SHP2 dans la ville MDA-MB-468 :

1. Le blocus : Lorsqu'ils ont arrêté SHP2, l'ancienne route principale (PTPN11) a été fermée.
2. La déviation : Les messagers ne se sont pas arrêtés ; ils ont immédiatement trouvé de nouveaux raccourcis. Ils ont commencé à utiliser ERBB3 et PIK3CA comme leurs nouvelles entrées principales beaucoup plus souvent.
3. Le rétablissement : Lorsqu'ils ont retiré le médicament (l'ont éliminé par lavage), la ville a commencé lentement à reconstruire l'ancienne route SHP2, et le trafic est revenu vers l'entrée principale d'origine, ERBB2.

La conclusion
Cet article ne se contente pas de dire « les médicaments échouent ». Il fournit une manière systématique et reproductible de dessiner une carte détaillée de la façon dont les cellules réacheminent leurs signaux lorsqu'elles sont attaquées. En transformant des données désordonnées en itinéraires clairs et étape par étape, cet outil aide les scientifiques à comprendre exactement comment les cellules contournent les traitements, ce qui est la première étape pour concevoir de meilleures stratégies afin de les arrêter.

Résumé technique

Résumé technique : Recherche par faisceau BFS guidée par la phosphoprotéomique résolue en temps

Énoncé du problème
La résistance adaptative aux thérapies ciblant les kinases et les phosphatases est souvent pilotée par le réarrangement pharmacologique des réseaux de signalisation intracellulaire. Cependant, le domaine manque actuellement de méthodes computationnelles systématiques capables de quantifier ces changements de voies entre différentes conditions directement à partir de données phosphoprotéomiques résolues en temps. Les approches existantes peinent à transformer des mesures phosphoprotéomiques dynamiques et de haute dimension en hypothèses de signalisation interprétables mécanistiquement et spécifiques au type cellulaire.

Méthodologie
Les auteurs présentent un cadre algorithmique conçu pour reconstruire des voies de signalisation spécifiques au type cellulaire à partir de données phosphoprotéomiques résolues en temps. Le cœur de ce cadre intègre plusieurs stratégies computationnelles distinctes :

• Prétraitement des données : Le système utilise une affectation d'état binaire adaptative aux données basée sur l'écart absolu médian (MAD) pour catégoriser les états d'activation des phosphosites.
• Algorithme de reconstruction de voies : Le cadre utilise une recherche en largeur (BFS) combinée à une recherche par faisceau guidée par le poids des interactions. Cette exploration est ancrée sur des « nœuds actifs » soutenus par l'expérience, identifiés aux limites des zones et aux terminaux au sein des données.
• Intégration de réseaux : La recherche opère sur la base de données d'interactions protéine-protéine STRING (accessible via l'API REST v12.5), permettant l'inclusion de protéines intermédiaires déduites par STRING comme connecteurs intermédiaires, sans nécessiter l'installation d'une base de données locale.
• Validation et nettoyage : Un pipeline de nettoyage de chemins post-énumération assure la génération de routes de signalisation acycliques et biologiquement interprétables. Crucialement, le cadre effectue une validation au niveau des arêtes contre des données d'expression de lignées cellulaires basées sur l'Human Protein Atlas pour garantir la pertinence biologique.

Résultats clés
Le cadre a été évalué sur cinq ensembles de données phosphoprotéomiques publiés couvrant trois types cellulaires : HeLa, MDA-MB-468 et EGFR Flp-In HEK293T. L'analyse a démontré les résultats spécifiques suivants :

• Spécificité du type cellulaire : La méthode a capturé avec succès des architectures de signalisation EGFR distinctes, uniques à chaque type cellulaire.
• Effets de l'inhibition de SHP2 : Dans les cellules MDA-MB-468 soumises à trois conditions pharmacologiques, l'inhibition de SHP2 a été montrée comme abolissant les voies médiées par PTPN11. Cette inhibition a induit un déplacement dans la distribution des effecteurs de premier saut, augmentant l'engagement avec ERBB3 (passant de 21,5 % à 25,2 % des chemins) et PIK3CA (passant de 9,2 % à 14,3 % des chemins).
• Réversibilité et récupération : Après lavage de l'inhibiteur de SHP2, le système a révélé une récupération partielle des voies médiées par PTPN11. Parallèlement, ERBB2 est réapparu comme l'effecteur de premier saut dominant, représentant 30,3 % des chemins.

Signification et revendications
L'article positionne ce cadre comme une approche systématique et reproductible pour convertir les mesures phosphoprotéomiques résolues en temps en hypothèses de signalisation interprétables mécanistiquement. Les auteurs revendiquent l'applicabilité directe de cette méthode à la modélisation de la résistance aux médicaments et à la conception de thérapies combinées. En quantifiant le réarrangement des voies induit par les médicaments et en identifiant des déplacements spécifiques dans la distribution des effecteurs, le cadre vise à combler le fossé entre les données phosphoprotéomiques brutes et des informations exploitables sur les mécanismes de résistance adaptative.