Dans Les Unités Pilotes D’opérations Unitaires De Bioprocédés, Comment La Limitation De La Diffusion De L’oxygène Affecte-T-Elle Les Agrégats Cellulaires ?

La limitation de la diffusion de l’oxygène devient critique lorsque les agrégats cellulaires dépassent la distance sur laquelle l’oxygène peut pénétrer avant d’être consommé. Dans un agrégat sphérique, la concentration en oxygène diminue de la surface vers le centre, car les cellules respirent continuellement. Lorsque le rayon de l’agrégat dépasse une valeur critique, le centre devient appauvri en oxygène ou anoxique, ce qui entraîne un stress métabolique, une réduction de la croissance, une altération de la fonction et potentiellement la mort cellulaire.

La taille critique de l’amas cellulaire correspond au rayon pour lequel la concentration centrale calculée en oxygène atteint zéro ou un seuil de viabilité défini. Elle peut être estimée à l’aide d’un modèle diffusion–réaction en régime stationnaire, puis validée expérimentalement au moyen de mesures de l’oxygène ou de la viabilité dans les agrégats.

Pourquoi l’oxygène devient limitant à l’intérieur des agrégats

Les cellules consomment l’oxygène dans l’ensemble de l’amas

L’oxygène doit d’abord être transféré de la phase gazeuse vers le liquide, puis du liquide vers la surface de l’agrégat, avant de diffuser à travers l’intérieur riche en cellules.

Comme l’oxygène est peu soluble dans les milieux aqueux, une forte densité cellulaire et une respiration rapide peuvent consommer l’oxygène plus vite qu’il n’est renouvelé au centre de l’agrégat.

Le centre subit la limitation la plus importante

Pour un agrégat sphérique présentant une concentration de surface en oxygène (C_s), la concentration en oxygène diminue vers le centre. Le centre est le point le plus éloigné du liquide oxygéné et il est entouré par la plus grande quantité de biomasse consommatrice d’oxygène.

Il en résulte un gradient radial d’oxygène :

Région externe : oxygénée et généralement viable.
Région intermédiaire : limitée en oxygène et soumise à un stress métabolique.
Cœur : fortement appauvri en oxygène ou anoxique lorsque l’agrégat est trop volumineux.

Conséquences pour les cellules et les performances du procédé

La limitation en oxygène peut provoquer :

Une réduction du métabolisme aérobie et de la croissance.
Des modifications de la différenciation ou de la formation du produit.
Une accumulation du stress métabolique.
Des régions centrales nécrotiques ou non viables.
Une plus grande variabilité d’un lot à l’autre dans les cultures de sphéroïdes ou d’agrégats.

À l’échelle pilote, ce problème peut être amplifié par un mauvais transfert global d’oxygène dû à une forte concentration en biomasse, une viscosité accrue du bouillon, une agitation inadéquate ou une insuffisance de barbotage.

Modélisation de la pénétration de l’oxygène dans un agrégat sphérique

Bilan diffusion–réaction en régime stationnaire

Pour un agrégat sphérique de rayon (r_0), on suppose :

Une concentration en oxygène en régime stationnaire.
Une diffusivité effective constante (D) à l’intérieur de l’agrégat.
Une consommation uniforme d’oxygène.
Une géométrie sphérique.
Aucune production d’oxygène à l’intérieur de l’agrégat.
Une concentration en oxygène spécifiée (C_s) à la surface de l’agrégat.

Si la consommation cellulaire d’oxygène est approximée par une cinétique d’ordre zéro (q), l’équation gouvernante est :

[ D\left[\frac{1}{r^2}\frac{d}{dr}\left(r^2\frac{dC}{dr}\right)\right]-q=0 ]

Les conditions aux limites sont :

[ \left.\frac{dC}{dr}\right|_{r=0}=0 ]

[ C(r_0)=C_s ]

La première condition traduit la symétrie au centre : il n’existe aucune direction privilégiée pour le flux d’oxygène à (r=0).

Profil de concentration en oxygène

Le profil d’oxygène obtenu est :

[ C(r)=C_s-\frac{q}{6D}\left(r_0^2-r^2\right) ]

Au centre :

[ C(0)=C_s-\frac{q r_0^2}{6D} ]

Cette équation montre la relation essentielle : la concentration centrale en oxygène diminue avec le carré du rayon de l’agrégat et augmente avec la concentration de surface en oxygène et la diffusivité.

Estimation du rayon critique

Le rayon critique idéalisé est obtenu en fixant la concentration centrale à zéro :

[ C(0)=0 ]

Par conséquent :

[ r_c=\sqrt{\frac{6DC_s}{q}} ]

Si (r_0<r_c), le modèle prévoit que l’oxygène reste disponible dans l’ensemble de l’agrégat.

Si (r_0>r_c), la concentration calculée deviendrait négative au centre. Physiquement, la concentration en oxygène ne peut pas être négative ; un cœur appauvri en oxygène ou anoxique se forme à la place.

Définir une taille critique des amas plus pratique

Utiliser un seuil de viabilité plutôt que le seul critère d’oxygène nul

Le critère d’oxygène nul est utile pour identifier le début de l’anoxie, mais de nombreuses cellules sont soumises à un stress avant que l’oxygène n’atteigne zéro.

Un seuil plus pratique est une concentration centrale minimale acceptable, (C_{\min}), représentant une exigence de viabilité, de productivité ou de phénotype :

[ C(0)=C_{\min} ]

Le rayon admissible correspondant est :

[ r_{\text{crit}}= \sqrt{\frac{6D(C_s-C_{\min})}{q}} ]

Cette définition est généralement plus utile pour l’exploitation d’une unité pilote, car l’objectif du procédé n’est pas simplement d’éviter l’anoxie : il est de maintenir des performances cellulaires acceptables dans tout l’agrégat.

Distinguer le rayon du diamètre

Les équations utilisent le rayon, alors que les mesures des agrégats sont souvent exprimées en diamètre.

Ainsi :

[ d_{\text{crit}}=2r_{\text{crit}} ]

Confondre le rayon et le diamètre entraîne une erreur d’un facteur deux dans la taille critique rapportée.

Comment déterminer les tailles critiques dans une unité pilote

Mesurer les paramètres qui contrôlent le transport interne de l’oxygène

Le calcul nécessite des estimations de :

Diffusivité effective, (D) : diffusivité de l’oxygène dans l’agrégat, en tenant compte de l’effet du compactage cellulaire et de la matrice extracellulaire.
Concentration de surface en oxygène, (C_s) : concentration en oxygène à la limite de l’agrégat.
Consommation spécifique ou volumique d’oxygène, (q) : vitesse de consommation d’oxygène de l’agrégat dans les conditions de culture concernées.
Concentration minimale acceptable, (C_{\min}) : concentration associée à une viabilité ou une fonction acceptables.

Ces valeurs doivent être obtenues dans des conditions représentatives de température, de milieu, de densité cellulaire et d’état physiologique.

Tenir compte du transfert global d’oxygène

La concentration à la surface de l’agrégat n’est pas nécessairement égale à la concentration d’oxygène dissous dans le milieu, définie par le contrôleur du bioréacteur.

Une différence de concentration peut exister, car l’oxygène doit traverser un film liquide externe avant de pénétrer dans l’agrégat. Le transfert d’oxygène à l’échelle du réacteur et la diffusion à l’échelle de l’agrégat doivent donc être considérés comme des résistances couplées.

L’unité pilote doit caractériser ou contrôler :

L’intensité de l’agitation.
Le débit de barbotage et la distribution des bulles.
Le coefficient volumétrique de transfert d’oxygène, (k_La).
La consigne d’oxygène dissous.
La composition gazeuse et la pression de fonctionnement.
La viscosité du bouillon et la concentration en biomasse.

Augmenter l’agitation, améliorer le barbotage ou fonctionner sous pression élevée peut accroître la disponibilité de l’oxygène à la surface de l’agrégat. Toutefois, ces actions n’éliminent pas la limitation interne de diffusion si l’amas lui-même est trop volumineux.

Comparer les prédictions du modèle aux mesures

Une détermination pratique combine calcul et observation :

Mesurer la distribution des tailles d’agrégats au cours du temps.
Estimer la consommation d’oxygène dans des conditions de culture représentatives.
Estimer ou mesurer la diffusivité effective de l’oxygène.
Déterminer les conditions d’oxygène dissous à la surface ou dans le milieu.
Calculer la concentration centrale et le rayon critique prévus.
Mesurer l’oxygène interne, la viabilité ou l’activité métabolique.
Comparer le début observé du stress du cœur de l’agrégat avec la prédiction du modèle.

Des sondes sensibles à l’oxygène, des indicateurs fluorescents de l’oxygène, des microélectrodes et des colorations spatiales de viabilité peuvent aider à déterminer si un cœur appauvri en oxygène prévu par le modèle se forme effectivement.

Utiliser des expériences dépendant de la taille

Les cultures peuvent être conduites ou échantillonnées sur une plage de tailles d’agrégats, en maintenant les autres conditions aussi constantes que possible.

La représentation de la concentration centrale en oxygène, de la viabilité ou de la productivité en fonction du rayon de l’agrégat aide à identifier la taille à partir de laquelle les performances commencent à diminuer. Ce seuil empirique peut ensuite être comparé à l’estimation diffusion–réaction.

Quels facteurs contrôlent la taille critique des amas cellulaires ?

Augmenter la concentration de surface en oxygène

Une valeur plus élevée de (C_s) augmente la quantité d’oxygène disponible pour l’agrégat et le rayon critique prévu.

En pratique, cela peut impliquer une consigne d’oxygène dissous plus élevée, un enrichissement en oxygène, une amélioration du transfert gaz–liquide ou une pression accrue dans le réacteur. L’avantage est limité si le transfert d’oxygène vers le liquide devient le principal goulot d’étranglement.

Augmenter la diffusivité effective

Une valeur plus élevée de (D) permet à l’oxygène de pénétrer plus rapidement dans l’agrégat.

Le compactage des agrégats, une matrice extracellulaire dense et une porosité interne réduite peuvent diminuer la diffusivité effective. Le contrôle de l’agrégation cellulaire et de la structure de la matrice peut donc être aussi important que la modification de l’agitation du réacteur.

Réduire la vitesse de consommation d’oxygène

Une valeur plus faible de (q) augmente la taille admissible de l’amas.

Une alimentation contrôlée peut modérer la croissance de la biomasse et la demande respiratoire. Le fonctionnement en fed-batch peut également aider à éviter une viscosité excessive et une accumulation de biomasse qui réduisent le (k_La) global du réacteur.

Contrôler la taille des agrégats

La stratégie la plus directe consiste à maintenir le rayon des agrégats en dessous de la valeur critique calculée ou validée expérimentalement.

En pratique, cela peut impliquer la gestion de la densité d’inoculation, de l’agitation, du temps de culture, de la géométrie des supports ou de la désagrégation des agrégats. Les petits supports ou les couches cellulaires actives minces réduisent les distances de diffusion et contribuent à maintenir une plus grande proportion de cellules viables.

Comprendre les compromis

Une agitation plus forte n’est pas une solution universelle

L’agitation peut améliorer le transfert global d’oxygène et réduire les gradients de concentration dans le liquide, mais une agitation excessive peut endommager les agrégats fragiles ou les cellules sensibles au cisaillement.

Elle peut également modifier la taille et la morphologie des agrégats, ce qui rend difficile la distinction entre la réponse biologique et la réponse hydrodynamique.

Une concentration plus élevée en oxygène dissous peut avoir des coûts biologiques

Augmenter la disponibilité de l’oxygène peut prévenir l’hypoxie, mais des conditions fortement oxygénées peuvent modifier le comportement cellulaire ou accroître le stress oxydatif dans les systèmes sensibles à l’oxygène.

L’objectif doit donc être une plage d’oxygène adéquate, et non la concentration maximale possible en oxygène dissous.

Le modèle d’ordre zéro est une approximation

Le modèle d’ordre zéro suppose une vitesse constante de consommation d’oxygène partout où l’oxygène est présent. En réalité, les cellules présentent souvent une consommation d’oxygène qui dépend de la concentration locale en oxygène, avec un comportement généralement proche du premier ordre lorsque la concentration en oxygène est faible.

Par conséquent, l’équation simple du rayon critique doit être utilisée principalement comme outil de présélection et de conception. Un modèle de consommation du premier ordre ou de type Monod peut être nécessaire lorsque la cinétique en faible teneur en oxygène, les changements métaboliques ou une forte hétérogénéité spatiale sont importants.

Le régime stationnaire ne décrit pas nécessairement tous les agrégats

Les agrégats peuvent croître, se compacter, se fragmenter ou modifier leur état métabolique au cours du temps.

Un modèle transitoire peut être nécessaire lorsque la taille des agrégats évolue rapidement ou lorsque l’exposition à l’oxygène varie fortement pendant le mélange et le changement d’échelle.

L’oxygène dissous global ne prouve pas une suffisance interne

Un bioréacteur peut afficher une valeur acceptable au niveau du contrôleur alors que les centres des agrégats restent limités en oxygène.

La mesure globale décrit le liquide à l’emplacement de la sonde ; elle ne mesure pas directement la concentration en oxygène à l’intérieur des amas cellulaires.

Comment appliquer cela à votre unité pilote

L’approche la plus fiable consiste à combiner l’estimation diffusion–réaction avec des mesures du transfert d’oxygène et de la viabilité à l’échelle pilote.

Si votre objectif principal est de prévoir le début de l’anoxie : utilisez (r_c=\sqrt{6DC_s/q}) comme première estimation et identifiez le moment où le centre de l’agrégat atteint zéro oxygène.
Si votre objectif principal est de maintenir la viabilité cellulaire : remplacez l’oxygène nul par un (C_{\min}) défini biologiquement et calculez le rayon maximal admissible à partir de (r_{\text{crit}}=\sqrt{6D(C_s-C_{\min})/q}).
Si votre objectif principal est l’exploitation du réacteur : mesurez ou estimez le (k_La), contrôlez l’agitation et le barbotage, puis vérifiez que le transfert global d’oxygène maintient la concentration de surface requise pour les agrégats.
Si votre objectif principal est le changement d’échelle : suivez les distributions de taille des agrégats, la viscosité, la consommation d’oxygène et la viabilité interne, plutôt que de vous fier uniquement à la consigne d’oxygène dissous.
Si votre objectif principal est l’optimisation de la conception : comparez des supports plus petits ou des couches cellulaires plus minces avec des consignes d’oxygène et une agitation plus élevées, tout en vérifiant les dommages dus au cisaillement et le stress oxydatif.

En considérant la taille des agrégats, la consommation cellulaire d’oxygène, la diffusivité interne et le transfert d’oxygène du réacteur comme un système couplé, les unités pilotes peuvent définir des limites de fonctionnement qui protègent à la fois la viabilité cellulaire et les performances du procédé.

Tableau récapitulatif :

Facteur clé	Effet sur les agrégats ou le procédé	Contrôle pratique
Rayon de l’agrégat	Les amas plus volumineux créent des gradients d’oxygène plus prononcés et des cœurs anoxiques	Contrôler l’agrégation, le temps de culture, l’agitation ou la désagrégation
Concentration de surface en oxygène ($C_s$)	Une disponibilité accrue de l’oxygène augmente le rayon critique	Optimiser l’oxygène dissous, la composition gazeuse, la pression et le $k_La$
Diffusivité effective ($D$)	Une diffusivité plus faible limite la pénétration de l’oxygène	Gérer le compactage, la matrice extracellulaire et la porosité
Vitesse de consommation d’oxygène ($q$)	Une respiration plus élevée réduit la concentration centrale en oxygène	Contrôler la croissance de la biomasse, l’alimentation et les conditions de culture
Rayon critique	Estimé par $r_c=\sqrt{6DC_s/q}$ pour une concentration centrale nulle en oxygène	L’utiliser comme estimation initiale de présélection
Seuil de viabilité	Une limite pratique utilise $r_{crit}=\sqrt{6D(C_s-C_{min})/q}$	Définir $C_{min}$ à partir de données de viabilité, de productivité ou de phénotype
Validation à l’échelle pilote	Confirme si les limites d’oxygène modélisées correspondent aux performances biologiques	Mesurer la taille des agrégats, l’oxygène, la viabilité et l’activité métabolique

Vous devez maintenir des agrégats cellulaires viables et productifs à l’échelle pilote ? Contactez LABPARK pour discuter des unités pilotes d’opérations unitaires éducatives et professionnelles destinées aux applications de bioprocédés et de biotechnologie. LABPARK aide les universités, les instituts de recherche et les entreprises à évaluer le transfert d’oxygène, les limites de taille des agrégats, le $k_La$, l’agitation, le barbotage et les conditions de changement d’échelle grâce à des solutions pratiques d’unités pilotes. Utilisez des données de procédé représentatives pour identifier les contraintes de diffusion, protéger la viabilité cellulaire et sélectionner les conditions de fonctionnement avec une plus grande confiance.