Bifurcations and multistability in inducible three-gene toggle switch networks

Problème
Les systèmes biologiques reposent sur des réseaux de régulation génique capables de multistabilité pour conduire des phénotypes cellulaires distincts, tels que la différenciation et les décisions de destin de développement. Bien que les commutateurs de type "toggle switch" bistables à deux gènes soient bien caractérisés, l'extension de ces concepts à des réseaux d'ordre supérieur (par exemple, des systèmes à trois gènes) pour comprendre la tristabilité et des dynamiques plus complexes reste un défi. Une limitation primaire des modèles théoriques existants est le traitement des paramètres de régulation (constantes de dissociation, taux de transcription) comme des variables abstraites déconnectées des mécanismes biologiques spécifiques qui les contrôlent in vivo. Plus précisément, le rôle des molécules effectrices (inducteurs) dans la modulation de l'activité des facteurs de transcription allostériques est souvent omis ou simplifié. Cette étude comble l'écart entre la topologie abstraite du réseau et les mécanismes biologiques ajustables expérimentalement en examinant comment les concentrations d'effecteurs pilotent les changements dynamiques dans les commutateurs à trois gènes.

Méthodologie
Les auteurs construisent un cadre mathématique pour un réseau de répression mutuelle à trois gènes (un commutateur de type "toggle switch") en utilisant des modèles de fonctions de Hill dérivés de principes de mécanique statistique.

• Modèle de base : Le système est défini par des équations différentielles décrivant la production et la dégradation de trois répresseurs ($R_1, R_2, R_3$). Le modèle incorpore une interprétation thermodynamique de l'induction via le modèle de Monod-Wyman-Changeux (MWC). Ici, la concentration totale d'un facteur de transcription est mise à l'échelle par une fonction de probabilité, $p_{act}(c)$, qui dépend de la concentration d'un effecteur allostérique ($c$). Cette probabilité détermine la fraction de répresseurs actifs.
• Symétrie et perturbation : L'étude suppose initialement un système symétrique où tous les répresseurs se lient avec une affinité égale. Elle introduit ensuite systématiquement des perturbations pour briser cette symétrie, analysant comment les biais dans les forces de répression affectent les seuils de bifurcation.
• Extensions par auto-activation : Le modèle est étendu pour inclure l'auto-activation, où une protéine peut agir à la fois comme répresseur d'autres gènes et comme activateur de sa propre expression. Deux interprétations biologiques distinctes du contrôle par l'effecteur sont testées :
  • Cas I : La liaison de l'effecteur rend la protéine inactive, quel que soit sa fonction (répresseur ou activateur).
  • Cas II : La liaison de l'effeteur modifie directement la conformation de la protéine, changeant sa fonction d'un répresseur actif en un activateur actif.
• Modes de liaison : Les liaisons compétitives (liaison répresseur/activateur mutuellement exclusive) et non exclusives (liaison simultanée) sont toutes deux analysées.
• Analyse : Les auteurs effectuent une analyse de bifurcation, une analyse de stabilité linéaire (via les valeurs propres de la matrice Jacobienne) et des simulations numériques pour cartographier les points fixes, les points de selle et les régimes dynamiques (monostabilité, bistabilité, tristabilité, quadristabilité) en fonction de la concentration de l'inducteur et des paramètres du modèle.

Contributions clés et résultats

1. Changements dynamiques pilotés par l'inducteur : Dans le commutateur à trois gènes de base, l'augmentation de la concentration d'un inducteur unique (ciblant un seul répresseur) conduit le système à travers des régimes dynamiques distincts. À de faibles concentrations d'inducteur, le système est tristable (trois états stables, chacun dominé par un répresseur). À mesure que la concentration de l'inducteur augmente, le système subit des bifurcations nœud-selle, passant d'abord à un régime tristable restreint (chemins d'accès moins nombreux entre les états) pour finalement s'effondrer en un commutateur bistable entre les deux répresseurs non induits.
2. Sensibilité à la coopérativité et à la symétrie : La concentration d'inducteur seuil requise pour passer de la tristabilité à la bistabilité est relativement insensible aux augmentations du coefficient de Hill (coopérativité) une fois que $n \geq 4$. Cependant, le système est hautement sensible à la rupture de symétrie. De petites perturbations dans les forces de répression ($\alpha$) peuvent éliminer entièrement le régime tristable ou rétrécir la fenêtre des dynamiques complexes, suggérant que la "proximité" d'un système biologique par rapport à la symétrie peut être déduite de sa réponse dynamique à l'induction.
3. Impact de l'auto-activation : L'introduction de l'auto-activation ne change pas nécessairement le nombre d'états stables dans le régime de liaison compétitive, mais approfondit les puits de potentiel des états stables, rendant les transitions entre eux moins probables. Cependant, avec une auto-activation suffisamment forte, le système peut présenter une quadristabilité (quatre états stables), incluant un état où les trois gènes sont exprimés.
4. Dynamiques distinctes selon les mécanismes d'effecteur : L'étude démontre que l'interprétation biologique du contrôle des facteurs de transcription par les effecteurs conduit à des paysages dynamiques fondamentalement différents :
   • Dans le Cas I (l'effeteur inactive la protéine), la réponse du système est principalement une suppression de l'activité.
   • Dans le Cas II (l'effeteur change la fonction), le système présente des seuils de bifurcation et des configurations d'états stables différents. Par exemple, dans les scénarios de liaison non exclusive, le Cas II permet des états stables exprimant le gène ciblé par l'effeteur uniquement à des concentrations élevées de l'effeteur, alors que le Cas I permet cela à de faibles concentrations.
5. Complexité de la liaison non exclusive : Lorsque les répresseurs et les activateurs se lient de manière non exclusive, l'espace des paramètres s'étend pour permettre des états stables avec plusieurs gènes exprimés simultanément. L'existence de ces états dépend fortement de l'équilibre entre la force d'activation et la force de répression, ainsi que du mécanisme spécifique de contrôle par l'effeteur.

Signification
L'article affirme que modéliser le contrôle par l'effeteur à travers une interprétation thermodynamique concrète (modèle MWC) transforme non seulement les tendances phénotypiques du réglage dynamique, mais aussi l'ensemble des régimes dynamiques disponibles. En reliant l'analyse de bifurcation directement aux paramètres accessibles expérimentalement (concentrations d'effecteurs), l'étude offre un cadre pour :

• Déchiffrer les mécanismes de régulation : Distinguer entre différents mécanismes biologiques (par exemple, si un effeteur inactive une protéine ou change sa fonction) sur la base de la réponse entrée-sortie du système et du comportement de bifurcation.
• Réglage expérimental : Fournir une base théorique pour utiliser les concentrations d'inducteurs comme des "boutons de réglage" pour sonder la topologie sous-jacente des réseaux de régulation génique et inférer des paramètres tels que la coopérativité et la symétrie.
• Comprendre la multistabilité : Éclairer comment de simples changements dans l'architecture de régulation (ajout de l'auto-activation, changement de l'exclusivité de liaison) et les mécanismes d'effeteurs peuvent générer des comportements de multistabilité complexes (tristabilité, quadristabilité) pertinents pour les décisions de destin cellulaire dans les systèmes naturels et synthétiques.

Les auteurs concluent que leur approche complète les travaux existants en définissant rigoureusement le fonctionnement des effecteurs, permettant ainsi l'identification des paramètres physiques pilotant l'expression dans les données expérimentales et offrant une structure ajustable pour encoder un comportement multistable inductible.