A transport-only null model for apparent heterogeneity in diffusively dosed organoid arrays

Resumen Técnico: Un modelo nulo de transporte puro para la heterogeneidad aparente en arreglos de organoides

Planteamiento del problema
Los organoides de cáncer de hígado derivados de pacientes se utilizan cada vez más para estudiar la biología tumoral; sin embargo, los organoides nominalmente idénticos sometidos al mismo protocolo de dosificación suelen exhibir una variación fenotípica sustancial. Aunque esta variación se atribuye frecuentemente a la heterogeneidad biológica o al error de medición, un tercer mecanismo suele pasarse por alto: la exposición intracelular desigual impulsada por el transporte espacial. En un reservorio difusible compartido, la posición de un organoide respecto a la fuente de dosificación, los límites de la cámara y los organoides vecinos dicta el material que llega y entra en cada compartimento. En consecuencia, una entrega espacialmente heterogénea puede generar fenotipos heterogéneos incluso cuando todos los parámetros biológicos son idénticos. El artículo aborda la falta de un "nulo de transporte puro" utilizable para cuantificar cuánta variación entre organoides puede ser generada únicamente por la geometría y los protocolos de entrega.

Metodología
Los autores construyen un modelo nulo de transporte a fenotipo que acopla la difusión en masa, la cinética de superficie y la dinámica del estado intracelular.

1. Transporte en masa y cinética de superficie: El modelo trata el dominio extracelular como una cámara acotada que contiene $N$ compartimentos pequeños y bien separados (organoides). Acopla la difusión en masa con la adsorción dependiente de encuentros, la residencia superficial reversible y la internalización productiva.

   • Teoría de la renovación: El modelo emplea ecuaciones de renovación para describir los ciclos repetidos de adsorción, residencia, desorción y re-adsorción. Distingue entre "reinicio" (donde una molécula desorbida olvida su ubicación y reingresa a un fluido bien mezclado) y "búsqueda continua" (donde se retiene la memoria espacial).
   • Condiciones de contorno: El intercambio superficial se modela como un filtro reactivo dependiente de la frecuencia en lugar de un sumidero estático, incorporando un formalismo de tiempo local de frontera para manejar la dinámica de encuentros no markoviana.

2. Reducción dimensional mediante asíntotas emparejadas:

   • Para dominios bidimensionales, las asíntotas emparejadas reducen el problema del dominio perforado a un sistema de funciones de Green de $N \times N$ evaluado en los centros de los organoides.
   • La reducción expresa el flujo productivo en términos de fuerzas monopolo y capacitancias dependientes de la forma (en 3D) o interacciones logarítmicas (en 2D).
   • Se demuestra que el error de esta reducción espacial es $O(\varepsilon^2)$ para objetivos circulares, donde $\varepsilon$ es la relación entre el radio del organoide y el tamaño de la cámara.

3. Dinámica del estado intracelular: El flujo de internalización productiva impulsa un modelo de estado celular de cinco dimensiones (progenitor, diferenciado, maduro, apoptótico/estresado y puntuación de enfermedad). El modelo asume parámetros biológicos idénticos en todos los organoides; por lo tanto, cualquier variación en el fenotipo de salida surge estrictamente de los flujos productivos dependientes de la geometría.

4. Verificación numérica:

   • Espacial: El modelo reducido se verifica frente a referencias radiales exactas y soluciones de elementos finitos completos (elementos P2) del problema del dominio perforado para $N=2, 5, 10$ objetivos.
   • Temporal: La reconstrucción en el dominio del tiempo se valida mediante inversión de Laplace numérica. Los autores comparan un sustituto Gamma de dos momentos contra la inversión exacta, observando que, si bien el sustituto es altamente preciso para leyes de residencia markovianas (exponenciales), puede fallar para leyes no markovianas (de cola larga).

Contribuciones clave

1. Formulación del nulo de transporte puro: El artículo formula una receta directa para mapear los centros de los organoides medidos y la geometría de dosificación hacia la captación productiva y la variación fenotípica bajo una biología idéntica. Esto proporciona una referencia específica para el diseño contra la cual se puede probar la heterogeneidad observada.
2. Cuantificación de la dispersión espacial: A través de 2000 arreglos aleatorios de diez organoides con dosificación localizada, el modelo predice un coeficiente de variación (CV) de maduración de transporte puro mediano de 0.623 y un CV de la puntuación de enfermedad de 0.835. Una auditoría de diseño de un solo factor muestra que este valor mediano puede oscilar entre 0.27 y 0.86 dependiendo de la geometría.
3. Clarificaciones mecanísticas:
   • Búsqueda conservadora vs. con pérdida: En una cámara conservadora (paredes reflectantes, sin pérdida en masa), la probabilidad total de captación productiva eventual es 1. La desorción retrasa la entrega y redistribuye la captación, pero no reduce la captación eventual total a menos que exista un canal de pérdida competidor.
   • Forma de la ley de residencia: Las leyes de residencia con medias iguales pero formas distintas (por ejemplo, Gamma vs. Exponencial vs. Cola larga) producen fenotipos transitorios distintos. El transitorio depende de la forma de la ley, no solo de su media.
   • Normalización de la competencia: La captación bruta confunde la competencia con la distancia a la fuente. Se deriva una relación normalizada por la competencia que cancela el forzamiento de la fuente, aislando la interacción objetivo-objetivo.
   • Identificabilidad de la mezcla: La mezcla parcial (interpolación entre reinicio y búsqueda continua) es identificable a partir de las fracciones de captación final solo cuando la retención es heterogénea; la captación resuelta en el tiempo es genéricamente informativa incluso bajo una retención homogénea.
4. Implicaciones de diseño:
   • Geometría de dosificación: La dosificación distribuida (por ejemplo, cuatro puertos) reduce la dispersión base aproximadamente cinco veces en comparación con un único puerto localizado tras igualar la exposición media del arreglo.
   • Límites del tamaño de lote: Aumentar el tamaño del lote promedia únicamente la variación de proceso independiente dentro del lote. Los efectos de línea y lote compartidos producen un resultado generalmente no degenerado de probabilidad de aceptación de lote grande, lo que significa que los lotes más grandes no necesariamente convergen a un resultado determinista si existen efectos compartidos.

Resultados

• Variabilidad base: Para un puerto localizado, el CV de maduración mediano es 0.623. Este cae a 0.124 para una fuente espacialmente uniforme (tras igualar la exposición media), demostrando una reducción de aproximadamente cinco veces en la dispersión.
• Sensibilidad: La variabilidad es altamente sensible al número de objetivos, al tamaño de los objetivos, a la relación de aspecto de la cámara y a la ubicación de la dosificación (por ejemplo, esquina vs. centro). Es mayormente insensible a las tasas de pérdida en masa en el límite conservador.
• Convergencia: El modelo reducido logra una precisión de $O(\varepsilon^2)$ para el flujo productivo. A separaciones de aproximadamente cuatro radios de los objetivos, el error en el CV del fenotipo es inferior a $10^{-3}$, aunque los errores de flujo individual siguen siendo mayores ($0.5\% - 2\%$).
• Aceptación de lote: La probabilidad de aceptación del lote no tiende a 0 ni a 1 a medida que aumenta el tamaño del lote si existen efectos compartidos (de línea o de lote). En su lugar, converge a un límite interior determinado por la distribución de los efectos compartidos.

Significancia y Reivindicaciones
El artículo afirma proporcionar un modelo nulo específico para la geometría en lugar de un gemelo digital calibrado. Su principal significancia es establecer una línea base: antes de atribuir la heterogeneidad observada a la variación biológica, se debe cuantificar la geometría medida y el protocolo de entrega para ver cuánto de esa variación pueden explicar.

• Una medición que exceda la predicción del modelo identifica estrictamente la variación biológica o cinética que el transporte por sí solo no puede explicar.
• No exceder la predicción no prueba la homogeneidad biológica; simplemente indica que los datos son indistinguibles del ruido inducido por el transporte.
• El marco ofrece recomendaciones de diseño específicas y probables: la dosificación distribuida reduce significativamente la heterogeneidad aparente, y aumentar el tamaño del lote tiene rendimientos decrecientes si no se abordan los efectos de proceso compartidos.

Los autores declaran explícitamente que los ingredientes matemáticos (teoría de la renovación, asíntotas emparejadas, reducciones de Dirichlet-a-Neumann) son heredados de la literatura existente. La novedad reside en su composición en un modelo de transporte a fenotipo verificado para arreglos de organoides de cáncer de hígado y en la derivación de las consecuencias específicas respecto al igualamiento del forzamiento, la interpretación de la retención y los límites de aceptación de lote.