Modelo Matemático Descifra Cómo los Patrones de Glucanos de IgG Cambian con la Edad

La glicosilación de IgG —la unión de cadenas de azúcar a los anticuerpos— no es estática. Cambia con la edad, la enfermedad y la genética, y estos cambios tienen consecuencias diagnósticas y funcionales. Sin embargo, traducir las mediciones brutas de glicanos en conocimiento mecanístico ha seguido siendo difícil. Este estudio da un paso clave al construir y validar un modelo matemático cuantitativo que reconstruye la actividad enzimática subyacente a los perfiles de glicanos de IgG a partir de datos poblacionales a gran escala.

Los investigadores construyeron un modelo cinético basado en reglas de la N-glicosilación de IgG que abarca cuatro compartimentos del Golgi (cisternas cis-, medial-, trans- y la red trans-Golgi). El modelo comienza con el precursor de manosa M5 y simula las acciones secuenciales de siete enzimas: GnT I, GnT II, GnT III, Man II, FucT, GalT y SiaT. Fue calibrado con datos de glicanos medidos por UHPLC de 915 individuos de la isla de Korčula y validado en una cohorte independiente de 890 individuos de la isla de Vis, ambas en Croacia. Cada individuo fue representado por 22 picos cromatográficos de glicanos, y el modelo fue personalizado ajustando seis parámetros individuales de concentración enzimática por persona.

El modelo alcanzó una sólida concordancia con las distribuciones experimentales de glicanos en ambas cohortes. De forma notable, reveló un descenso estadísticamente significativo y consistente en la concentración de GalT relacionado con la edad en ambas poblaciones. GalT es responsable de añadir galactosa a los residuos terminales de GlcNAc —un paso bien conocido que disminuye con el envejecimiento y promueve el desplazamiento hacia formas de glicanos agalactosilados vinculados a la inflamaging. El modelo recuperó eficazmente esta señal biológica únicamente a partir de los perfiles de glicanos, sin medición directa de enzimas.

Más allá de reproducir la biología conocida, el modelo demostró que las concentraciones recuperadas de GalT y otras enzimas podían predecir la edad cronológica con una precisión comparable a la de los biomarcadores de envejecimiento basados en picos de glicanos ya establecidos. Esto posiciona las actividades enzimáticas modeladas como una capa mecánisticamente interpretable entre las variantes genéticas y los productos de glicanos —lo que podría hacer que los estudios de asociación de genoma completo (GWAS) sean más potentes y biológicamente interpretables al dirigirse a la variación a nivel enzimático en lugar de a picos individuales de glicanos compuestos.

Los autores conciben esto como el primer paso de una estrategia en dos etapas: recuperar las actividades enzimáticas mediante modelado matemático y, a continuación, mapear las variantes genéticas sobre dichas actividades. Esto podría mejorar la capacidad de detección de SNPs cuyos efectos se diluyen u oscurecen al analizarse frente a compuestos complejos de picos de glicanos. Las limitaciones incluyen la exclusión del procesamiento temprano en el RE y en el cis-Golgi, la dependencia de poblaciones de ascendencia europea y el uso de concentraciones relativas en lugar de absolutas de glicanos. No obstante, el enfoque es computacionalmente viable, está disponible de forma abierta a través de la plataforma BioUML y es extensible a otros sistemas de glicoproteínas.