Interaction between cell membranes and protein inclusions in the large-deformation regime

Resumen Técnico: Interacción entre membranas celulares e inclusiones proteicas en el régimen de gran deformación

Planteamiento del problema
Las membranas biológicas son superficies dinámicas cuya morfología y función están críticamente gobernadas por las inclusiones proteicas (IP). Si bien la interacción entre las IP y las membranas lipídicas ha sido modelada extensamente en el régimen de pequeña deformación (SD, por sus siglas en inglés) mediante expansiones perturbativas alrededor de configuraciones planas, muchos procesos biológicamente relevantes involucran deformaciones significativas de la membrana. Estos incluyen escenarios donde las IP se mueven a altas velocidades, induciendo invaginaciones significativas, o donde múltiples IP se posicionan en estrecha proximidad, creando gradientes pronunciados que no pueden describirse mediante teorías linealizadas. Este trabajo aborda la brecha en la comprensión de las interacciones membrana-proteína en el régimen de gran deformación (LD), centrándose específicamente en la ruptura de las aproximaciones SD y la emergencia de comportamientos elásticos complejos.

Metodología
Los autores emplean un enfoque híbrido que combina simulaciones numéricas de elementos finitos (FE) con derivaciones analíticas aproximadas.

• Marco numérico: El estudio utiliza el software de elementos finitos IRENE para resolver la ecuación de forma no lineal completa derivada de la energía libre de Helfrich. La membrana se modela como una variedad bidimensional embebida en tres dimensiones utilizando la calada de Monge ($z(x)$). La ecuación gobernante incluye términos para la curvatura media ($H$), la curvatura gaussiana ($K$), la rigidez de flexión ($\kappa$) y la tensión superficial ($\sigma$).
• Condiciones de contorno: Las simulaciones imponen restricciones específicas, incluyendo fronteras de membrana ancladas en un radio $R$, ángulos de contacto fijos ($\alpha$) en la interfaz de la proteína (modelada como infinitamente rígida) y condiciones de no deslizamiento para los flujos de membrana.
• Enfoque analítico: Para obtener una visión física, los autores derivan una solución analítica aproximada para el régimen LD asumiendo una curvatura media cero ($H=0$). Esto conduce a un perfil de catenoide, lo que permite expresiones explícitas de la altura y la fuerza de la membrana como funciones del ángulo de contacto y el desplazamiento.
• Dinámica de flujo: En presencia de flujos de membrana, el estudio resuelve ecuaciones diferenciales parciales (PDE) acopladas que describen la conservación de la masa, el equilibrio de cantidad de movimiento (incluyendo términos convectivos y tensiones viscosas), y la ecuación de forma, mientras que se desprecia el fluido externo para evitar la paradoja de Stokes dentro del dominio computacional.

Contribuciones clave y resultados

1. Limitaciones de la teoría linealizada:
   Al comparar las soluciones de FE con la ecuación linealizada (válida para gradientes pequeños), los autores demuestran que la aproximación lineal falla cuando el ángulo de contacto supera $\tan \alpha \approx 0.4$ o cuando los desplazamientos verticales son grandes. En estos regímenes LD, se debe resolver la ecuación de forma no lineal completa para capturar la geometría correcta de la membrana.

2. Solución analítica aproximada (Curvatura media cero):
   Los autores proponen una solución analítica basada en la suposición de curvatura media cero, lo que produce un perfil de catenoide. Aunque es una aproximación para el caso general, captura con éxito las características cualitativas del régimen LD. Predice que el desplazamiento vertical espontáneo ($h$) de una proteína con un ángulo de contacto fijo escala logarítmicamente con el tamaño del dominio ($h \sim \ln R$) para valores grandes de $R$.

3. Relación fuerza-desplazamiento:
   El estudio calcula la fuerza ejercida por la membrana sobre una IP. Un hallazgo clave es que esta fuerza exhibe un comportamiento no monotónico con respecto al desplazamiento vertical de la proteína. La fuerza aumenta hasta un máximo y luego disminuye, una característica que parece ser independiente de la geometría específica de la proteína. Esta no monotonicidad sugiere el potencial inicio de inestabilidades dinámicas relevantes para la formación y retracción de túbulos. La magnitud de estas fuerzas se estima en el rango de 1–10 pN, consistente con las fuerzas ejercidas por motores moleculares.

4. Interacción mediada por la membrana entre dos IP:
   Para un sistema de dos IP que interactúan, la energía potencial mediada por la membrana no sigue un simple decaimiento de ley de potencia como se observa en los regímenes SD. En su lugar, el potencial exhibe un decaimiento de sub-ley de potencia, disminuyendo más lentamente que una ley de potencia con la distancia entre proteínas. Esto refleja la naturaleza compleja del medio elástico. Además, el potencial de interacción es una función impar de la orientación de la proteína:

• Las IP con orientaciones idénticas (misma "carga") se repelen.
• Las IP con orientaciones opuestas (cargas opuestas) se atraen.
  Esto confirma la analogía entre la orientación de la IP y la carga eléctrica en el régimen LD. El sistema favorece energéticamente una configuración antisimétrica, lo que potencialmente conduce a patrones de tipo antiferromagnético de proteínas alternantes.

5. Deformaciones inducidas por el flujo y velocidad característica:
   En presencia de flujos de membrana, los autores identifican una escala de velocidad característica, $v^* = \kappa / (L\eta)$, que corresponde a un número de Scriven-Love de unidad. Esta velocidad separa dos regímenes distintos:

• $v < v^*$: El flujo tiene un efecto insignificante en la forma de la membrana; la rigidez de flexión domina.
• $v > v^*$: Surgen deformaciones inducidas por el flujo, con una longitud de onda característica $\lambda \sim \kappa/(v\eta)$.
  El estudio sugiere que por encima de este umbral, la competencia entre las fuerzas viscosas y de flexión conduce a cambios significativos de forma.

Significancia y afirmaciones
El artículo afirma proporcionar predicciones cuantitativas para sistemas de membrana-proteína en escenarios biológicamente relevantes que involucran grandes deformaciones, yendo más allá de las limitaciones de las teorías perturbativas SD. Los resultados tienen implicaciones directas en la comprensión de:

• Clasificación y agrupamiento de proteínas: Las fuerzas no monotónicas y las reglas específicas de atracción/repulsión basadas en la orientación sugieren mecanismos para cómo las proteínas podrían organizarse en cúmulos o patrones en la membrana.
• Tráfico de membranas: La identificación de una velocidad crítica para las deformaciones inducidas por el flujo implica que el movimiento de las proteínas puede impulsar activamente los cambios de forma de la membrana, influyendo potencialmente en la formación y el transporte de vesículas.
• Ejemplos biológicos: Los autores ilustran la relevancia de sus hallazgos utilizando la bacteriorrodopsina (BR) en vesículas gigantes unilamelares (GUV). Estiman que la velocidad de difusión lateral de las moléculas de BR es suficiente para inducir deformaciones impulsadas por el flujo con longitudes de onda comparables al espaciamiento entre proteínas, lo que sugiere que las interacciones mediadas por el flujo podrían influir significativamente en el patrón de la BR.

En conjunto, el trabajo establece un marco para analizar la compleja interacción no lineal entre la elasticidad de la membrana, la geometría de la proteína y los flujos hidrodinámicos, ofreciendo una descripción más realista de la mecánica de la membrana en entornos celulares dinámicos.