Una proteína de 100 aminoácidos tiene ~10⁴⁷ conformaciones posibles. Si las probara al azar
a 10¹³ por segundo, tardaría más que la edad del universo en encontrar la correcta. Sin embargo,
en la célula se pliega en milisegundos. Cómo lo logra — la paradoja de Levinthal —
es uno de los problemas centrales de la biología molecular y la base de AlphaFold.
La Función de Energía
E = E_bond + E_LJ + E_HP
Energía total = bonds (enlaces covalentes) + Lennard-Jones (interacciones de Van der Waals) + HP (núcleo hidrofóbico)
Cada término captura una fuerza diferente:
- E_bond — Energía elástica de los enlaces péptídicos. Penaliza elongaciones o compresiones excesivas de la cadena principal.
- E_LJ — Potencial de Lennard-Jones entre residuos no covalentes. Repulsión a corta distancia (núcleos electrónicos), atracción a media distancia (Van der Waals).
- E_HP — Modelo HP (hidrofóbico-polar). Los residuos hidrofóbicos pagan un costo energético por estar expuestos al agua; buscan formar un núcleo hidrofóbico compacto en el interior. Es la fuerza dominante en el plegamiento real.
El Paisaje Energético
La idea clave de la teoría moderna del plegamiento es el paisaje energético (energy landscape). La conformación de la proteína es un punto en un espacio de altísima dimensionalidad. La energía libre define un "relieve" sobre ese espacio:
Embudo de plegamiento
El paisaje no es plano: tiene forma de embudo. Hay muchas conformaciones de alta energía (extendidas) y pocas de baja energía (plegadas). La proteína "cae" cuesta abajo en el embudo, convergiendo hacia el estado nativo.
Mínimos locales
El paisaje tiene valles locales donde la proteína puede quedar atrapada: conformaciones incorrectas pero estables. Los chaperones moleculares in vivo ayudan a la proteína a escapar de estos trampas y seguir hacia el mínimo global.
Estado nativo
El mínimo global de energía libre. En la célula es la conformación biológicamente activa. La hipótesis de Anfinsen (Nobel 1972): toda la información necesaria para alcanzarlo está en la secuencia de aminoácidos.
Velocidad y cooperatividad
El plegamiento real es cooperativo: no hay pasos intermedios estables. La proteína pasa directamente de desnaturalizada a nativa en milisegundos, como una transición de fase de primer orden.
Los Tres Métodos de Búsqueda
Dinámica Molecular (MD)
Integra las ecuaciones de Newton F = ma para cada átomo. Calcula fuerzas de E_bond, E_LJ y E_HP en cada paso. Físicamente riguroso pero lento: la energía se disipa gradualmente.
Monte Carlo (MC)
Propone conformaciones aleatorias. Acepta si baja la energía; acepta con probabilidad e^(−ΔE/kT) si sube. El criterio de Metropolis. Explora el espacio de manera estocástica sin seguir trayectorias físicas.
Recocido Simulado (SA)
Como MC pero con temperatura decreciente: al principio acepta movimientos malos (exploración global); al final los rechaza (convergencia local). Escapa de mínimos locales mejor que MD puro.
AlphaFold2 (DeepMind, 2021) resolvió el problema del plegamiento de proteínas. Usando redes neuronales entrenadas con la base de datos PDB, predice la estructura 3D a partir de la secuencia con precisión atómica. En el CASP14, superó a todos los grupos humanos por un margen histórico. Desde 2022, AlphaFold ha predicho la estructura de ~200 millones de proteínas, incluyendo todo el proteoma humano. La biología estructural cambió para siempre.
Malplegamiento y Enfermedades
Cuando la proteína no alcanza su estado nativo y se agrega en estructuras anómalas (hojas β cruzadas, fibrillas amiloides), causa enfermedades degenerativas:
- Alzheimer — Agregados de péptido β-amiloide y proteína tau hiperfosforilada.
- Parkinson — Fibrillas de α-sinucleína (cuerpos de Lewy).
- Priones (Creutzfeldt-Jakob) — La proteína PrP normal se convierte en la forma mal plegada PrP^Sc, que actúa como molde para plegar más proteínas incorrectamente. Transmisible sin ADN ni ARN.
- Diabetes tipo II — Depósitos de amilina en células beta del páncreas.
Experimentos Guiados
Experimento 1 — Comparar métodos: MD vs. SA
- Genera una cadena aleatoria. Activa Dinámica Molecular: la proteína oscila y se pliega lentamente, siguiendo trayectorias físicas.
- Resetea y activa Recocido Simulado (SA): la temperatura inicial alta permite saltos conformacionales grandes. Observa cómo al bajar T, la exploración se vuelve más fina.
- SA generalmente encuentra conformaciones de menor energía que MD en el mismo tiempo: explora mejor el espacio porque no está limitado por barreras físicas.
Experimento 2 — Núcleo hidrofóbico
- Observa los residuos hidrofóbicos (naranja) y polares (azul) durante el plegamiento.
- En la conformación final, los residuos hidrofóbicos deberían estar en el interior (núcleo), lejos del "agua" virtual.
- Este efecto hidrofóbico es la fuerza dominante: crea el colapso inicial que guía el embudo de plegamiento hacia la estructura nativa.
- Una proteína con más hidrofóbicos tiende a tener un núcleo más compacto y una energía final más baja.
Experimento 3 — Mínimos locales y trampas
- Con Monte Carlo a temperatura baja (T pequeño), la proteína queda atrapada en el primer mínimo que encuentra.
- Aumenta T: los movimientos malos se aceptan con mayor probabilidad y la proteína escapa de trampas locales.
- Reduce T muy lentamente (SA manual): llega a conformaciones de menor energía que con T fija.
- Esto replica exactamente el rol de los chaperones in vivo: aportan energía para deshacer conformaciones incorrectas y reiniciar la búsqueda.
Conexiones