Clinical and Molecular Characterization of a Rare KRAS p.Val14Ile Variant in Noonan Syndrome: Integrative Phenotypic Analysis and Molecular Dynamics Simulation

Resumen Técnico: Caracterización Clínica y Molecular de una Variante Rara de KRAS p.Val14Ile en el Síndrome de Noonan

Planteamiento del Problema
El síndrome de Noonan (SN) es una rasopatía genéticamente heterogénea causada por la desregulación de la vía de señalización RAS/mitogen-activated protein kinase (RAS/MAPK). Aunque PTPN11 es el gen causante más común, las variantes de KRAS son raras (<5% de los casos) y a menudo están insuficientemente caracterizadas en cuanto a sus relaciones específicas entre genotipo y fenotipo, así como sus consecuencias moleculares. La literatura existente se centra predominantemente en variantes recurrentes o de "hotspot", dejando la comprensión del impacto funcional de sitios no canónicos, como la variante de sentido erróneo de novo* rara *NM_004985.4:c.40G>A (p.Val14Ile), incompleta. Este estudio aborda la necesidad de caracterizar sistemáticamente el espectro clínico y los mecanismos moleculares de esta variante específica para explicar su distintividad fenotípica.

Metodología
El estudio empleó un enfoque integrativo que combina el análisis de casos clínicos, la revisión sistemática de la literatura y el modelado computacional avanzado:

1. Caracterización Clínica y Genética: Se identificó un paciente varón de 15 años diagnosticado con SN en el Hospital de la Gente de Dongying. Se realizó una secuenciación del exoma completo (WES) en el probando y sus padres, seguida de una secuenciación de Sanger para confirmar la herencia. La variante fue clasificada utilizando las guías ACMG y respaldada por herramientas de predicción in silico (PolyPhen-2, SIFT, MutationTaster, REVEL) y análisis de conservación (ConSurf).
2. Revisión Sistemática de la Literatura y Metaanálisis: Una búsqueda exhaustiva (PubMed, Web of Science, Embase, CNKI, ClinVar, HGMD) identificó 57 pacientes con SN portadores de variantes de KRAS. Estos se dividieron en un grupo de estudio (p.Val14Ile, N=10) y un grupo de control (otras variantes de KRAS, N=47). Se realizó un análisis clínico comparativo utilizando términos de la Ontología del Fenotipo Humano (HPO). La significancia estadística se evaluó mediante pruebas de chi-cuadrado, la prueba exacta de Fisher y regresión logística binaria para calcular los Odds Ratios (OR) con intervalos de confianza (IC) del 95%. El análisis agrupado utilizó modelos de efectos aleatorios para evaluar las tendencias fenotípicas generales.
3. Simulaciones de Dinámica Molecular (MD) y Acoplamiento: Para investigar los mecanismos estructurales, se realizaron experimentos de acoplamiento molecular y simulaciones de MD de 100 ns utilizando GROMAS 2024.4. Se utilizaron el campo de fuerza AMBER14SB y GAFF2 para la proteína y los ligandos (GTP), respectivamente. Las simulaciones compararon los complejos de KRAS–GTP de tipo silvestre (WT) y el mutante p.Val14Ile–GTP. Las métricas clave incluyeron la desviación de la raíz cuadrática media (RMSD), la fluctuación de la raíz cuadrática media (RMSF), el radio de giro (Rg), el área de superficie accesible al solvente (SASA), el recuento de puentes de hidrógeno y los paisajes de energía libre (FEL). Las energías libres de unión se calcularon utilizando el método de mecánica molecular/superficie de Gauss generalizada (MM/GBSA).

Result últimos Resultados Clave

• Fenotipo Clínico: El grupo p.Val14Ile exhibió perfiles fenotípicos ampliamente similares a otras variantes de KRAS en la mayoría de los sistemas orgánicos (rasgos faciales, restricción del crecimiento, cardiovascular, esquelético). Se observó una diferencia estadísticamente significativa en las anomalías linfáticas y hematológicas, las cuales fueron significativamente menos frecuentes en el grupo p.Val14Ile en comparación con el grupo de control (25.00% vs. 86.67%; P = 0.013; OR = 0.05). El análisis agrupado de HPO indicó una tendencia general hacia menores probabilidades de los fenotipos evaluados en los portadores de p.Val14Ile (Pooled OR = 0.63, 95% CI: 0.42–0.96) bajo el modelo de efectos aleatorios, sin heterogeneidad significativa (I² = 0%).
• Clasificación Genética: La variante KRAS p.Val14Ile fue clasificada como patogénica según los criterios ACMG (PS2_VeryStrong, PS3, PS4, PM1, PP2, PM2_Supporting). Se localiza en la región del bucle P altamente conservada del dominio G. El análisis de empalme (splicing) in silico sugirió que la variante puede interrumpir los motivos potenciadores de empalme exónico (ESE), afectando potencialmente el empalme del pre-ARNm además de la sustitución de aminoácidos.
• Dinámica Molecular y Acoplamiento:
  • Afinidad de Unión: El acoplamiento molecular predijo una energía de unión de GTP más débil para el mutante (−7.829 kcal·mol⁻¹) en comparación con el WT (−8.387 kcal·mol⁻¹).
  • Estabilidad Estructural: Los análisis de RMSD y Rg mostraron que el complejo mutante exhibió mayores fluctuaciones conformacionales (notablemente a los 35 ns y 80 ns) en comparación con el estable complejo WT, aunque la estructura global permaneció intacta (SASA ~90 nm²).
  • Flexibilidad Local: El análisis de RMSF reveló una mayor flexibilidad local en la región de los residuos 30–40 (que abarca el Switch I) en el mutante.
  • Red de Interacción: El mutante mantuvo menos puentes de hidrógeno con GTP (6–8) en comparación con el WT (7–9). El análisis de energía de interacción a nivel de residuo mostró contribuciones alteradas, particularmente una disminución en las contribuciones de Tyr32 y Gly13 y un aumento en Gly15.
  • Energía Libre: Los cálculos de MM/GBSA indicaron una energía libre de unión menos favorable para el mutante (−34.27 kcal·mol⁻¹) en comparación con el WT (−48.64 kcal·mol⁻¹).
  • Paisaje Conformacional: El análisis del paisaje de energía libre mostró que el mutante muestreó una distribución más dispersa de estados de baja energía en comparación con los grupos concentrados del WT, lo que sugiere una alteración en el conjunto conformacional en lugar de un reordenamiento estructural grueso.

Contribuciones Clave

1. Expansión Clínica: El estudio expande el espectro clínico del SN asociado a KRAS al proporcionar una comparación sistemática de la variante p.Val14Ile frente a otras variantes de KRAS, destacando una reducción específica en las anomalías linfáticas y hematológicas.
2. Perspectiva Mecanística: Proporciona una explicación estructural y dinámica de la patogenicidad de la variante, demostrando que la sustitución p.Val14Ile no interrumpe el plegamiento general de la proteína, sino que altera la dinámica conformacional local, las redes de interacción y el conjunto conformacional del complejo KRAS–GTP.
3. Integración Metodológica: El artículo integra datos clínicos con simulaciones moleculares de alta resolución para cerrar la brecha entre un hallazgo genético raro y su potencial mecanismo molecular, centrándose específicamente en cómo los cambios sutiles en el bucle P afectan la dinámica del Switch I y las interacciones con efectores.

Significancia y Reivindicaciones
Los autores afirman que este estudio proporciona una explicación mecánica potencial para las características fenotípicas de la variante rara p.Val14Ile. Proponen un modelo donde la variante induce perturbaciones dinámicas localizadas en la región de unión a nucleótidos y en el Switch I, lo que conduce a una redistribución del conjunto conformacional de KRAS y a la alteración de las redes de interacción, en lugar de un simple efecto de ganancia o pérdida de función.

Los autores declaran explícitamente que estos hallazgos sugieren que las alteraciones específicas de la variante en la dinámica conformacional de KRAS pueden contribuir a la heterogeneidad clínica de las rasopatías. Enfatizan que, si bien los datos respaldan una tendencia general de menores probabilidades fenotípicas para p.Val14Ile, el tamaño limitado de la muestra y los amplios intervalos de confianza para los fenotipos individuales requieren una interpretación cautelosa. El estudio concluye que la variante p.Val14Ile produce un paisaje conformacional funcional distinto, ofreciendo un marco estructural para futuras investigaciones experimentales sobre los mecanismos de la enfermedad, reconociendo al mismo tiempo que los hallazgos computacionales requieren validación mediante ensayos bioquímicos y celulares directos.