El problema del plegamiento del genoma en la cápside viral se sitúa en la intersección de la física, la química y la virología. Trabajos relacionados en el archivo exploran aspectos afines: cómo los genes bacterianos se ensamblan en gotas vivas, cómo funcionan los transistores iónicos en nanoporos, cómo se manifiestan los efectos cuánticos en sistemas biológicos. Este artículo se distingue porque incluye por primera vez de manera sistemática en la dinámica molecular, simultáneamente, iones explícitos, agua explícita, sitios de carga discretos en la cápside y electrostática coulombiana completa. Anteriormente, estos cálculos usaban o bien un solvente implícito, o bien modelos simplificados de polímero, o condiciones iónicas fijas. Aquí, el apantallamiento iónico surge naturalmente de la dinámica de los iones, lo que permite verificar las predicciones de las teorías de campo medio a nivel molecular.
Los autores construyeron un modelo de grano grueso en el que el ARN se representa como una cadena lineal de 1122 perlas de tres tipos (fosfato, azúcar, azúcar) basado en el campo de fuerza MARTINI-2. La cápside se modeló como una doble capa esférica con un radio interno de aproximadamente 12 nm, en cuya superficie interna se fijaron 180 sitios cargados positivamente que imitan los dominios N-terminales de las proteínas. El agua se describió con un modelo polarizable de tres perlas, y los iones de sodio y cloro con partículas estándar de MARTINI. Las interacciones electrostáticas se calcularon con el método PPPM con una permitividad dieléctrica de 2.5. Las simulaciones se realizaron con el paquete LAMMPS a una temperatura de 310 K con termostato de Langevin. Para validar el modelo, se midió la longitud de persistencia de la cadena libre a diferentes concentraciones de sal (10–100 mM) y se comparó con datos experimentales de ARN monocatenario. Luego, el polímero se colocó dentro de la cápside y se analizaron los perfiles radiales y angulares de densidad de monómeros después de alcanzar el equilibrio.
El modelo reprodujo la dependencia experimental de la longitud de persistencia del ARN monocatenario con la concentración de sal: al aumentar la sal de 10 a 100 mM, la longitud de persistencia disminuye de ~2.5 a ~1.0 nm, lo que concuerda con los datos de Chen et al. para poli(rU). Dentro de la cápside, aumentar la concentración de sal de 10 a 100 mM desplaza sistemáticamente la distribución radial del genoma desde la pared hacia el centro, como predicen las teorías de campo medio. Sin embargo, el tiempo de equilibrio aumenta siete veces: de 64 ns a 10 mM a 470 ns a 100 mM. La variación independiente de la carga de la cápside y de la densidad de sitios cargados mostró que la conexión electrostática total gobierna la organización global, mientras que los detalles geométricos modulan el contacto local del genoma con la pared. Los perfiles angulares revelaron un doble pico alrededor de cada sitio cargado con 180 sitios, que desaparece al aumentar el número de sitios a 1470, cuando la distancia entre ellos se vuelve menor que la longitud de persistencia de la cadena. Representar los dominios N-terminales con una sola perla (que sobresale) divide el máximo único en dos, lo que indica que el genoma se enrolla alrededor de la protuberancia, en lugar de adherirse directamente.
Este trabajo se distingue de los estudios vecinos: a diferencia de los trabajos sobre genes bacterianos o transistores iónicos, aquí se verifican directamente las predicciones de las teorías de campo medio para el plegamiento de un polielectrolito en una capa esférica cargada, con consideración explícita de los detalles moleculares. Los trabajos vecinos usan modelos continuos, solvente implícito o condiciones iónicas fijas. En este estudio se muestra por primera vez que los iones discretos y los sitios de carga sobresalientes producen desviaciones cuantitativas del campo medio, especialmente en la distribución angular del genoma. Esto significa que las aproximaciones continuas son válidas para una imagen radial gruesa, pero no capturan la organización local, importante para comprender la especificidad del empaquetamiento.
El desarrollo natural es incluir la estructura secundaria del ARN (apareamiento de bases) y la topología ramificada, lo que podría cambiar el papel relativo de la carga total y la distribución de cargas. También es importante ampliar el modelo a iones multivalentes (Mg2+ y otros), que afectan más fuertemente el plegamiento del ARN. Los dominios N-terminales largos y flexibles, que van más allá de una o dos perlas, podrían añadir nuevos efectos conformacionales. Finalmente, pasar de la organización de equilibrio dentro de una cápside ya formada a la cinética de ensamblaje permitiría comparar directamente con experimentos de autoensamblaje viral.
Los resultados son importantes para la virología (comprensión del empaquetamiento del genoma), la física de polielectrolitos (comportamiento de cadenas cargadas en volumen confinado) y la nanotecnología (diseño de cápsides para la entrega de ácidos nucleicos terapéuticos).
Los próximos pasos son incluir la estructura secundaria del ARN, modelar el ensamblaje de la cápside alrededor del genoma y probar los efectos de los iones multivalentes.
En la física de polielectrolitos sigue abierta la cuestión del papel exacto de los iones discretos y las correlaciones en geometría confinada. Este trabajo muestra que estos efectos ya son notables con sal monovalente y cargas moderadas, lo que plantea la necesidad de desarrollar teorías más precisas que vayan más allá del campo medio. Para la virología, el problema clave no resuelto es la relación entre la secuencia del ARN, su estructura secundaria y la eficiencia de empaquetamiento en la cápside; este modelo crea una plataforma para estudiarla.
Esta sección está dirigida al autor del trabajo y a sus colegas. Los trabajos vecinos los encontró la búsqueda por significado en nuestro archivo, no por coincidencia de palabras. No evaluamos: mostramos el campo alrededor.
Qué hace el trabajo
Los autores construyeron un modelo de grano grueso del ARN de cadena sencilla dentro de la cápside viral, incluyendo explícitamente iones, agua y electrostática. El modelo reprodujo la dependencia experimental de la longitud de persistencia del ARN con la concentración de sal y mostró que aumentar la sal desplaza el genoma hacia el interior de la cápside, mientras que el tiempo para alcanzar el equilibrio aumenta siete veces.
Dónde está su fuerza
La fortaleza del trabajo radica en la inclusión sistemática de detalles moleculares (cargas discretas, iones explícitos, agua polarizable) en el modelo electrostático del empaquetamiento del genoma, lo que permite contrastar las predicciones de las teorías de campo medio a nivel molecular.
Hacia dónde puede avanzar
Vale la pena probar
Por qué puede importar
Comprender cómo las condiciones iónicas y la discreción de las cargas controlan el empaquetamiento del genoma puede ayudar a diseñar estrategias para interferir con el ensamblaje de los virus o para estabilizar las cápsides en aplicaciones biotecnológicas.
Trabajos cercanos en nuestro archivo
El archivo está vacío alrededor de este trabajo: el análisis más cercano queda en 0.51, frente a un 0.64 habitual. Parece que aquí al lado casi nadie ha perforado.
El tiempo de equilibrio en el modelo aumenta siete veces al incrementar la concentración de sal de 10 a 100 mM, lo que significa que en un ambiente más salado el genoma 'tarda más en encontrar' su plegamiento final. Para comparar, a 10 mM el equilibrio se alcanza en 64 nanosegundos, mientras que a 100 mM toma casi medio microsegundo.