Exactly Solvable Population Model with Square-Root Growth Noise and Cell-Size Regulation

Resumen Técnico: Modelo de Población Exactamente Soluble con Ruido de Crecimiento de Raíz Cuadrada y Regulación del Tamaño Celular

Planteamiento del Problema
El artículo aborda un problema central no resuelto en la fisiología celular cuantitativa: cómo las fluctuaciones microscópicas en las tasas de crecimiento de células individuales moldean la aptitud (fitness) y las estadísticas de tamaño a nivel de población. Si bien las tecnologías de célula única han permitido una caracterización detallada de las estadísticas de división y las estrategias de control de tamaño, las estadísticas temporales de las tasas de crecimiento de células individuales permanecen insuficientemente restringidas. Los modelos teóricos existentes suelen asumir que la tasa de crecimiento exponencial instantánea $\lambda(t)$ fluctúa independientemente del tamaño celular $s$. Sin embargo, desde una perspectiva microscópica, se espera que las fluctuaciones intrínsecas en la producción de biomasa escalen con la raíz cuadrada del tamaño celular (debido a la naturaleza extensiva del ruido de producción), lo que conduce a un proceso de difusión de tipo Feller. Los autores pretenden cerrar esta brecha introduciendo un modelo de población exactamente soluble que incorpore ruido de crecimiento dependiente del tamaño (de raíz cuadrada) para determinar su impacto en las tasas de crecimiento poblacional, las distribuciones de tamaño en estado estacionario y las fluctuaciones del tamaño de la población.

Metodología
Los autores analizan una población estructurada por tamaño donde las células individuales crecen de acuerdo con una ecuación diferencial estocástica (SDE) que representa una difusión de raíz cuadrada de Feller:
$$\frac{ds}{dt} = k s + \sigma \sqrt{2k \bar{s}_b s} \, \xi(t)$$
donde $k$ es la tasa de crecimiento media, $\bar{s}_b$ es el tamaño medio de nacimiento y $\sigma$ es la intensidad de ruido adimensional, y $\xi(t)$ es un ruido blanco gaussiano. La división y la regulación del tamaño se modelan de forma general: las células se dividen en un tamaño $s_d$ (que puede ser una variable aleatoria dependiente del tamaño de nacimiento $s_b$) y se particionan en dos hijas con tamaños $r s_d$ y $(1-r)s_d$.

El análisis procede a través de varias etapas analíticas:

1. Dinámica de la Masa Poblacional: Los autores aprovechan la aditividad del ruido de raíz cuadrada para demostrar que la masa total de la población $M_t$ obedece la misma SDE que una sola célula, lo que permite obtener soluciones exactas respecto a la masa total y el tamaño de la población.
2. Distribución de Tamaño en Estado Estacionario: Para un caso simplificado de tamaño de división constante, los autores resuelven la ecuación de Fokker-Planck asociada (una ecuación diferencial parcial con términos de deriva y difusión) para derivar la densidad exacta del tamaño celular en estado estacionario.
3. Mapeo de Proceso de Ramificación: El modelo se mapea a un proceso de ramificación de Bellman-Harris. Los autores derivan la distribución del tiempo de generación como el tiempo de primer paso del proceso de Feller y utilizan su transformada de Laplace para calcular los momentos del tamaño de la población.
4. Generalización: Los resultados se extienden para incluir tamaños de división ruidosos y mecanismos generales de regulación de tamaño (por ejemplo, modelos de adder, sizer, timer) desacoplando la dinámica de crecimiento de las estadísticas de división.

Contribuciones Clave y Resultados

1. Tasa de Crecimiento Poblacional Exacta:
   El resultado principal es que la tasa de crecimiento poblacional asintótica $\Lambda$ es exactamente igual a la tasa de crecimiento media de la célula individual ($\Lambda = k$). Esto se mantiene para todas las intensidades de ruido y todos los esquemas de división o regulación de tamaño que mantienen la homeostasis del tamaño.

   • Contraste: Esto contrasta fuertemente con los modelos con ruido de tasa de crecimiento independiente del tamaño, donde $\Lambda$ depende sensiblemente de la varianza y el tiempo de correlación del ruido.
   • Mecanismo: Este "desacoplamiento fuerte" surge específicamente de la estructura de ruido de raíz cuadrada aditiva, la cual no es una propiedad genérica de los modelos de crecimiento estocástico.

2. Distribución de Tamaño Celular en Estado Estacionario:
   Los autores derivan una solución analítica exacta para la distribución de estado estacionario de los tamaños celulares.

   • Estructura de Convolución: La solución para $\sigma > 0$ se obtiene mediante una convolución exponencial de un solo lado de la solución determinista ($\sigma=0$) con un núcleo de ancho $\sigma^2$.
   • Momentos: El tamaño celular medio de la población se desplaza en $-\sigma^2$ respecto al caso determinista, mientras que la varianza se modifica solo al orden de $\sigma^4$.
   • Desacoplo: La forma de la distribución depende de la regulación del tamaño solo a través de la distribución de tamaño de nacimiento del linaje estacionario; el ruido de crecimiento introduce un efecto de suavizado independiente de la división.

3. Fluctuaciones del Tamaño de la Población:
   Se deriva en forma cerrada la distribución del tamaño de población medio reescalado $W = \lim_{t\to\infty} N_t / \langle N_t \rangle$.

   • Tipo de Distribución: $W$ sigue una distribución de Poisson compuesto-exponencial que depende únicamente del parámetro de ruido de crecimiento $\sigma$.
   • Independencia: Esta distribución es invariante al ruido del tamaño de división, a la partición asimétrica o a mecanismos específicos de regulación de tamaño. Estos factores solo afectan al prefactor que relaciona la masa total con el número de células, no la forma de las fluctuaciones.
   • Varianza: La varianza relativa a largo plazo es $\text{Var}(N_t) / \langle N_t \rangle^2 = 2\sigma^2$.
   • Solubilidad: Los autores señalan que, más allá del proceso de Yule (tiempos de generación exponenciales), este es un raro ejemplo de un modelo de ramificación no trivial donde la distribución estacionaria completa del tamaño de población reescalado se conoce en forma cerrada.

Significancia y Reivindicaciones
El artículo identifica el crecimiento exponencial estocástico con ruido de raíz cuadrada como un "punto de referencia exactamente soluble" para estudiar las fluctuaciones de crecimiento en poblaciones con control de tamaño. Los autores afirman que sus resultados revelan una jerarquía de desacoplamientos:

1. La tasa de crecimiento poblacional es independiente del ruido de crecimiento y de los detalles de la división.
2. Las distribuciones de tamaño celular instantáneas dependen de la regulación solo a través de la distribución de tamaño de nacimiento del linaje, modificada por una convolución específica.
3. La forma de las fluctuaciones del tamaño de la población está determinada enteramente por el ruido de crecimiento, independientemente de los mecanismos de división.

Los autores sostienen que estas predicciones difieren cualitativamente de los modelos basados en ruido de tasa de crecimiento independiente del tamaño. Por consiguiente, las firmas derivadas (específicamente la independencia de la tasa de crecimiento poblacional de la intensidad del ruido y la forma específica de las fluctuaciones del tamaño de la población) proporcionan criterios concretos y contrastables para distinguir entre modelos microscópicos de crecimiento competidores utilizando datos de célula única existentes y futuros. El trabajo no propone nuevos protocolos experimentales, sino que ofrece un marco teórico para interpretar datos empíricos.