Si un ser humano, una planta, una mosca de la fruta o un cultivo de bacterias fotosintéticas queda encerrado en una cueva subterránea profunda o en una cámara oscura completamente aislada de la luz solar y de las fluctuaciones de temperatura, ¿pierde el organismo la noción del día y de la noche?

Sorprendentemente, la respuesta es no. Incluso bajo condiciones de oscuridad constante (Constant Darkness, DD), prácticamente todos los seres vivos de la Tierra —desde las cianobacterias unicelulares y los hongos hasta las plantas, los insectos y los mamíferos— mantienen una oscilación autónoma con un período de aproximadamente 24 horas. Los ciclos de sueño y vigilia, la temperatura corporal central, la secreción de cortisol y melatonina, la presión arterial, la vigilancia inmunitaria y la expresión del 10 % al 40 % de todos los genes codificantes de proteínas están gobernados por péndulos moleculares dentro de cada célula.

Sin embargo, desde la perspectiva de la físico-química y la dinámica no lineal, este reloj intracelular plantea tres profundos acertijos científicos:

1. La paradoja de la escala temporal y la estabilidad (Dinámica no lineal): Las reacciones bioquímicas elementales (unión de factores de transcripción, traducción, fosforilación) ocurren en milisegundos o minutos. Además, un bucle simple de retroalimentación negativa de dos variables se amortigua inevitablemente hasta detenerse en un punto de equilibrio. ¿Cómo logra una red molecular, mediante una bifurcación de Hopf, sostener una oscilación de ciclo límite (Limit-Cycle Oscillation) con un período de $86.400\text{ segundos}$ ($24\text{ horas}$)?
2. El dogma de la transcripción genética frente a los relojes puramente proteicos (Biología molecular y estructural): Durante décadas se creyó que todo reloj circadiano requería un bucle de retroalimentación de transcripción-traducción (TTFL) dependiente de ADN y ARNm. Pero en 2005, Takao Kondo y Masahiro Ishiura demolieron este dogma al demostrar que mezclar tan solo tres proteínas purificadas de cianobacterias (KaiA, KaiB y KaiC) con ATP en un tubo de ensayo —sin ADN, sin ARNm y sin ribosomas— genera una oscilación de fosforilación de 24 horas durante días.
3. El misterio biofísico de la compensación de temperatura ($Q_{10} \approx 1{,}0$): Las velocidades de las reacciones químicas y enzimáticas obedecen a la ecuación de Arrhenius, $k(T) = A \exp(-E_a / RT)$, duplicándose o triplicándose por cada aumento de $10^\circ\text{C}$ ($Q_{10} \approx 2\text{–}3$). Si el reloj biológico siguiera la cinética ordinaria de Arrhenius, correría más del doble de rápido en una tarde cálida, volviéndose inútil para medir el tiempo. ¿Cómo logran los relojes biológicos un coeficiente térmico invariante de $Q_{10} \approx 0{,}95\text{–}1{,}05$?

1. De la planta Mimosa pudica (1729) a los mutantes de Drosophila y el Premio Nobel 2017

En 1729, el astrónomo francés Jean-Jacques d’Ortous de Mairan colocó la planta sensitiva Mimosa pudica en una bodega en oscuridad constante (DD) y observó que sus hojas seguían abriéndose cada mañana y cerrándose cada noche sin ver el sol. En 1832, el botánico suizo Augustin Pyramus de Candolle demostró que bajo condiciones constantes este ciclo se repite con un período de libre curso (Free-Running Period) $\tau \approx 22\text{–}23\,\text{h}$.

En ausencia de sincronizadores externos (Zeitgebers), cada especie muestra su propio período endógeno $\tau$: $\tau \approx 23{,}5\,\text{h}$ en el ratón, $\tau \approx 23{,}8\,\text{h}$ en la mosca Drosophila melanogaster y $\tau \approx 24{,}18\,\text{h}$ en el ser humano ($24\text{ h } 11\text{ min}$, medido por Charles Czeisler en 1999). El hecho de que $\tau$ difiera ligeramente de 24,00 horas demuestra que el oscilador es endógeno, término bautizado como circadiano por Franz Halberg en 1959 (del latín circa, «alrededor de», y dies, «día»).

En 1971, Ronald Konopka y Seymour Benzer (Caltech) descubrieron mediante mutagénesis en Drosophila tres alelos mutantes de un único gen en el cromosoma X llamado period ($per$): $per^0$ (arrítmico), $per^S$ (período corto $\tau \approx 19\,\text{h}$) y $per^L$ (período largo $\tau \approx 29\,\text{h}$). En 1984, Jeffrey C. Hall, Michael Rosbash y Michael W. Young clonaron el gen period y descifraron su mecanismo de retroalimentación negativa, hallazgo galardonado con el Premio Nobel de Fisiología o Medicina en 2017.

2. El bucle de retroalimentación transcripción-traducción (TTFL) en eucariotas

  flowchart TD
    subgraph Nucleus["Núcleo celular (Nucleus)"]
        CB["Brazo positivo: Heterodímero CLOCK:BMAL1"]
        EBOX["Elemento potenciador E-box (5'-CACGTG-3')"]
        CB -->|"Unión al ADN y acetilación de histonas (CBP/p300 HAT)"| EBOX
        EBOX -->|"Activa la transcripción (Día subjetivo)"| MRNA_N["ARNm de Per1/2/3 y Cry1/2"]
        EBOX -->|"Impulsa el bucle secundario"| REV_ROR["ARNm de Rev-erbα/β y Rorα/β/γ"]
        INHIB["Complejo represor nuclear PER:CRY:CK1δ/ε (~1,9 MDa)"]
        INHIB -.->|"Inhibición por bloqueo y desplazamiento (Noche tardía)"| CB
        RORE["Elemento RORE en el promotor de Bmal1"]
        RORE -->|"Genera transcripción en antifase de Bmal1"| CB
    end

    subgraph Cytoplasm["Citoplasma (Cytoplasm)"]
        MRNA_C["ARNm citoplasmático de Per / Cry"]
        PROT["Proteínas nacientes PER1/2/3 y CRY1/2"]
        CK1["Caseína Quinasa 1δ/ε (CK1δ/ε)"]
        PHOS["Complejo multifosforilado PER:CRY:CK1"]
        DEG["Ubiquitinación SCF^(β-TrCP / FBXL3) y Proteasoma 26S"]
        REV_PROT["RORα (Activador) vs REV-ERBα (Represor)"]
    end

    MRNA_N -->|"Exportación nuclear"| MRNA_C
    MRNA_C -->|"Traducción ribosómica"| PROT
    PROT -->|"Phosphoswitch y retardo de fase de varias horas (Δτ)"| CK1
    CK1 --> PHOS
    PHOS -->|"Estabilización e importación nuclear (Tarde/Noche)"| INHIB
    PHOS -->|"Degradación citoplasmática"| DEG
    INHIB -->|"Eliminación de represores (Amanecer)"| DEG
    REV_ROR -->|"Traducción y retorno al núcleo"| REV_PROT
    REV_PROT -->|"Unión competitiva"| RORE

En los mamíferos, el heterodímero CLOCK:BMAL1 se une a secuencias E-box (5'-CACGTG-3') para activar la transcripción de Per1/2/3 y Cry1/2. En el citoplasma, las proteínas PER monoméricas son fosforiladas por la Caseína Quinasa $1\delta/\varepsilon$ (CK1$\delta/\varepsilon$) y degradadas por $\text{SCF}^{\beta\text{-TrCP}}$. Tras un retardo de fase ($\Delta\tau$) de 4 a 6 horas, los complejos PER:CRY:CK1$\delta/\varepsilon$ ($\sim 1{,}9\,\text{MDa}$) entran en el núcleo e inhiben a CLOCK:BMAL1 mediante bloqueo (CRY1) y desplazamiento del ADN (fosforilación de CLOCK mediada por PER2:CK1$\delta/\varepsilon$), estabilizados por el bucle secundario REV-ERB$\alpha/\beta$ / ROR$\alpha/\beta/\gamma$ sobre elementos RORE.

En humanos, la mutación Ser662Gly en PER2 (o Thr44Ala en CSNK1D/CK1$\delta$) impide la fosforilación estabilizadora del dominio FASP, acelera la degradación de PER2 vía Ser478, acorta el período a $\tau \approx 21{,}6\,\text{h}$ y causa el Síndrome Familiar de Fase del Sueño Avanzada (FASPS): los pacientes se duermen a las 19:00 y se despiertan completamente alerta a las 3:00 de la madrugada. Por el contrario, la mutación $\text{CRY1}\,\Delta 11$ alarga $\tau$ y provoca el síndrome de fase retrasada (FDSPS).

3. Biología matemática de los osciladores genéticos: Modelo de Goodwin–Goldbeter y bifurcación de Hopf

En un sistema bidimensional de retroalimentación negativa sin retardo ($\dot{M} = f(P) - \gamma_M M,\; \dot{P} = k_s M - \gamma_P P$ con $f'(P) < 0$), la divergencia del campo vectorial es estrictamente negativa en todo el plano positivo:

Por el criterio negativo de Bendixson, es matemáticamente imposible que exista una órbita cerrada (ciclo límite). Para cruzar una bifurcación de Hopf, Brian Goodwin (1965) y Albert Goldbeter (1995) formularon el modelo de tres variables ($M$, proteína citoplasmática $P_C$, proteína nuclear $P_N$):

La oscilación de ciclo límite requiere tres ingredientes matemáticos: (1) retroalimentación negativa, (2) un retardo de fase $\Delta\tau$ suficientemente largo (desfase de $180^\circ$) y (3) fuerte no linealidad / ultrasensibilidad. Con degradación lineal de primer orden, el teorema de Griffith (1968) exige un coeficiente de Hill $n > \sec^3(\pi/3) = 8$. Los relojes reales superan esta barrera mediante la degradación enzimática saturada de Michaelis-Menten (ultrasensibilidad de orden cero, $K_d \ll P_C$, Goldbeter 1995) y la titulación molecular / secuestro estequiométrico 1:1 de CLOCK:BMAL1 por PER:CRY (Kim & Forger, 2012), que genera un umbral no lineal abrupto incluso con $n = 1$.

4. El asombroso descubrimiento de Takao Kondo y Masahiro Ishiura (2005): El reloj KaiABC en un tubo de ensayo

En 2005, Masato Nakajima, Takao Kondo y colaboradores (Science, 2005) demostraron en la cianobacteria Synechococcus elongatus PCC 7942 que mezclar únicamente tres proteínas purificadas (KaiA, KaiB y KaiC) más ATP en un tubo de ensayo —sin ADN ni ARNm— reconstituye una oscilación de fosforilación de 24 horas compensada térmicamente (oscilador postraduccional, PTO).

  flowchart TD
    subgraph DayPhase["Día subjetivo: Fase de fosforilación ordenada"]
        U["Hexámero U-KaiC (No fosforilado: Ser431 / Thr432)"]
        T["Hexámero T-KaiC (Solo Fosfo-Thr432: pThr432)"]
        ST["Hexámero ST-KaiC (Doble fosforilación: pSer431 / pThr432)"]
        U -->|"El dímero KaiA se une al A-loop del anillo CII y estimula la autoquinasa"| T
        T -->|"Fosforilación secuencial de Ser431"| ST
    end

    subgraph NightPhase["Noche subjetiva: Cambio de plegamiento, secuestro de KaiA y desfosforilación"]
        B_SWITCH["ATPasa ultralenta del anillo CI + Cambio metamórfico de KaiB (gsKaiB → fsKaiB)"]
        S["Hexámero S-KaiC (Solo Fosfo-Ser431: pSer431)"]
        SEQ["Secuestro estequiométrico de KaiA (Complejo nocturno KaiA12-KaiB6-KaiC6) e intercambio de monómeros"]
        ST -->|"Acoplamiento alostérico del anillo CII al anillo CI"| B_SWITCH
        B_SWITCH -->|"Inhibe a KaiA y revierte la fosforilación de Thr432 (regeneración de ATP)"| S
        S -->|"fsKaiB captura y secuestra los dímeros libres de KaiA"| SEQ
        SEQ -->|"Desfosforilación de Ser431 y disociación al amanecer"| U
    end

El hexámero de doble anillo KaiC (anillos CI y CII) es la ATPasa más lenta de la naturaleza, hidrolizando tan solo $\sim 15\text{ moléculas de ATP por monómero al día}$. Durante el día, KaiA estimula la autofosforilación ordenada de Thr432 (T) y Ser431 (S) ($\text{U} \to \text{T} \to \text{ST}$). Al fosforilarse Ser431, el acoplamiento alostérico activa la hidrólisis de ATP en el anillo CI, donde la proteína metamórfica KaiB cambia radicalmente su plegamiento tridimensional (gsKaiB $\to$ fsKaiB tipo tiorredoxina), se une al anillo CI y secuestra todo el KaiA libre en el complejo nocturno $\text{KaiA}_{12}\text{KaiB}_6\text{KaiC}_6$. Sin KaiA, KaiC se desfosforila regenerando ATP a partir de ADP ($\text{ST} \to \text{S} \to \text{U}$) y sincroniza los hexámeros mediante intercambio de monómeros (monomer shuffling).

5. El misterio de la compensación de temperatura ($Q_{10} \approx 1{,}0$)

Mientras que la ecuación de Arrhenius $k(T) = A \exp(-E_a / RT)$ impone un $Q_{10} = \frac{k(T+10)}{k(T)} \approx 2\text{–}3$ a las reacciones ordinarias, los relojes circadianos logran $Q_{10} \approx 1{,}0$ mediante dos principios:

1. Equilibrio antagónico de Peter Ruoff (1992): En términos de los coeficientes de control del período $C_i^\tau = \frac{\partial \ln \tau}{\partial \ln k_i}$ ($\sum_i C_i^\tau = -1$) y las energías de activación $E_{a,i}$: $$\frac{d \ln \tau}{d T} = \sum_{i=1}^{N} \frac{\partial \ln \tau}{\partial \ln k_i} \frac{d \ln k_i}{d T} = \frac{1}{R T^2} \sum_{i=1}^{N} C_i^{\tau} E_{a,i} \approx 0$$ Las reacciones que acortan el período ($C_i^\tau < 0$) se compensan exactamente con las reacciones antagónicas que lo alargan ($C_j^\tau > 0$).
2. Compensación intrínseca en las enzimas CK1$\delta/\varepsilon$ y KaiC: Terauchi y Kondo (2007) demostraron que la actividad ATPasa basal de KaiC por sí sola tiene $Q_{10} \approx 1{,}0$. Asimismo, Yasushi Isojima, Yusuke Shinohara y Hiroki R. Ueda (2009, 2017) descubrieron que la fosforilación de PER2 por CK1$\delta/\varepsilon$ in vitro está compensada térmicamente ($Q_{10} \approx 1{,}0$), ya que la aceleración catalítica a mayor temperatura es contrarrestada por una mayor afinidad por el producto fosforilado (inhibición por producto mediada por el bolsillo aniónico K224).

6. El núcleo supraquiasmático (SCN), la melanopsina y la cronomedicina

Para sincronizar el reloj humano ($\tau \approx 24{,}18\,\text{h}$) con el día solar de 24,00 horas (fotoentrenamiento), el 1–2 % de las células ganglionares de la retina (ipRGCs) expresan melanopsina (OPN4), sensible a la luz azul ($\lambda_{\text{max}} \approx 480\,\text{nm}$). A través del tracto retinohipotalámico (RHT), liberan glutamato y PACAP en el núcleo supraquiasmático (SCN), fosforilando a CREB e induciendo la transcripción aguda de Per1/2 según la curva de respuesta de fase (PRC): retraso de fase al anochecer y avance de fase al amanecer.

En el SCN, unas 20.000 neuronas acopladas mediante VIP, AVP y GABA forman una red de osciladores de Kuramoto. El jet lag se produce por la desincronización interna entre el SCN (reajustado rápidamente por la luz) y los relojes periféricos (hígado, tracto gastrointestinal) sincronizados por el horario de las comidas. En cronomedicina, administrar estatinas de vida corta antes de dormir, aplicar quimioterapia cronomodulada (5-FU y oxaliplatino según Francis Lévi) y programar las infusiones de inmunoterapia (inhibidores de PD-1 como pembrolizumab o nivolumab) por la mañana en lugar de por la tarde casi duplica la supervivencia global en pacientes con melanoma y cáncer de pulmón (Qian et al., Lancet Oncology, 2021).