Act I · Physical & Clinical Story
Why Cryo-EM Needed a Non-Equilibrium Mutation (E1371Q) to Catch
Most ion channels in human biology operate like gentle swinging doors that drift back and forth at thermal equilibrium. In striking contrast, the human cystic fibrosis chloride channel functions as an active chemical engine that consumes cellular fuel to drive directional movement. When Andre Barato and Udo Seifert published their foundational statistical physics treatise in Physical Review Letters during 2015, they proved that precision in microscopic molecular machines requires continuous thermodynamic dissipation. Each time the chloride channel hydrolyzes one molecule of adenosine triphosphate, it breaks time reversal symmetry and pushes the gating cycle around a directional closed circuit. Structural biologists previously struggled to image the pore because this rapid chemical turnover snaps the gate shut within fractions of a second. By mutating the catalytic residue to eliminate fuel consumption, researchers trapped the elusive open state under high resolution electron microscopes.
Living cells depend upon reliable molecular switches to control ion flow. Unlike passive channels that fluctuate randomly without consuming metabolic resources, the cystic fibrosis channel acts like an active pump that burns cellular fuel during every cycle. This continuous energy consumption creates a definite arrow of time. Because the protein moves forward through distinct structural shapes, the time it spends open does not follow the memoryless distribution of passive thermal gates. Physicists describe this non equilibrium behavior through the concept of entropy production. The thermodynamic uncertainty relation reveals that reducing noise in single channel currents demands a higher expenditure of chemical free energy by the operating protein. Every increase in regulatory precision carries an inescapable metabolic cost. By examining the fluctuations in electrical patch clamp recordings, researchers can deduce whether a drug molecule alters fuel consumption or stabilizes intermediate states.
Everyday Analogy: Consider a wooden paddle wheel mounted over a fast-flowing millstream where falling water forces the blades to rotate in a single direction. The steady kinetic stream supplies continuous mechanical energy, ensuring that buckets lift water to higher fields with predictable regularity. If the millpond runs dry and the wheel rests in still water, random thermal breezes might nudge the paddles back and forth, but the device loses all ability to pump water uphill. In cellular biophysics, the CFTR channel functions precisely like this water wheel, utilizing the chemical flow of ATP hydrolysis to maintain a unidirectional gating cycle with high operational precision. Without that energetic input, the channel lapses into erratic, reversible fluctuations that fail to sustain ion gradients. Yet this mechanical millstream analogy ceases to be accurate at the molecular scale, because microscopic water wheels operate under overwhelming Brownian collisions that continuously attempt to reverse the blades.
How the Two Sliders Work in Plain Words: Before analyzing the master equations of non-equilibrium networks, readers can explore the gating thermodynamics through two interactive workbench controls. The first control, the Chemical Driving Affinity, sets the chemical potential difference provided by cellular ATP fuel relative to thermal energy. The second control, the Catalytic Site-2 Hydrolysis Rate, dictates how rapidly the protein cleaves the bound fuel molecule to initiate channel closing. Adjusting these two inputs immediately updates four fundamental thermodynamic readouts. The Schnakenberg Cycle Affinity readout tracks the net chemical driving force around the four-state gating loop. The Steady-State Entropy Production indicates the rate at which metabolic energy is dissipated into the cellular environment. The TUR Precision Lower Bound calculates the minimum variance floor for single-channel ion currents. Finally, the Detailed Balance Status confirms whether the channel is operating as an irreversible engine or has collapsed into equilibrium.
If structural biologists froze wild-type CFTR in the presence of MgATP, ongoing hydrolysis would continuously cycle the protein back to the closed state. By replacing Glutamate 1371 with Glutamine (E1371Q), researchers abolish ATP hydrolysis at Site 2—converting a non-equilibrium driven cycle into a quasi-equilibrium trap where the open NBD dimer stays locked long enough to image at 3.2 Å resolution.
The mathematical proof of the thermodynamic uncertainty relation established a fundamental trade-off between kinetic precision and entropy production in all continuous-time Markovian processes. For biophysicists studying cystic fibrosis, this principle transforms how single-channel electrophysiology recordings are interpreted. Wild-type CFTR does not operate as a passive equilibrium pore; rather, it hydrolyzes MgATP at Site 2 with a driving affinity of approximately 19.5 k_B T. This chemical drive breaks microscopic detailed balance and sustains an irreversible probability flux around the gating cycle. When gating mutations such as G551D uncouple ATP binding from pore opening, the channel loses its directional regularity and fails to conduct physiological chloride currents. The translational stakes of this non-equilibrium framework lie in rational drug development. By analyzing the variance and Fano dispersion of patch-clamp dwell times, researchers can determine whether small-molecule potentiators such as ivacaftor restore normal catalytic turnover or merely trap the channel in an artificial open conformation. This quantitative clarity prevents misleading conclusions during preclinical drug screening.
“Connects exact theorems of non-equilibrium statistical mechanics (Schnakenberg cycle decomposition and Barato–Seifert TUR) to experimental single-channel patch-clamp recordings of wild-type, G551D, and ivacaftor-rescued CFTR.”
— Cross-Pillar Research Thread · Grounded in Barato & Seifert (2015), Physical Review Letters 114(15):158101 — Thermodynamic Uncertainty Relation for Biomolecular Processes (PMID 25933341).
Key Parameters, Readouts & Primary Anchor
| Quantity / Anchor | Pinned Baseline | Physical / Clinical Role |
|---|---|---|
| Chemical Driving Affinity (Δμ_ATP / k_B T) | 19.5 k_B T |
Primary perturbation axis (1.0 to 24.0 k_B T) |
| Catalytic Site-2 Hydrolysis Rate (k_hyd) | 1.8 s⁻¹ |
Coupled co-perturbation axis (0.2 to 5.0 s⁻¹) |
| Schnakenberg Cycle Affinity A_C & Steady-State Entropy Production Σ | 19.5 k_B T (~12.0 kcal/mol) · 35.1 k_B / s per Channel |
Baseline biophysical operating point |
| Primary Literature Anchor | DOI 10.1103/PhysRevLett.114.158101 | Barato & Seifert (2015), Physical Review Letters 114(15):158101 — Thermodynamic Uncertainty Relation for Biomolecular Processes (PMID 25933341). · Kernel src/scripts/living-proof-observatory.js |
Act II · From Intuition to Governing Law
Schnakenberg Network Theory and the Thermodynamic Uncertainty
In Jürgen Schnakenberg's master-equation network theory (Rev. Mod. Phys. 1976), a continuous-time Markov jump process with transition rates w_{ij} between states i and j has steady-state probabilities π_i and edge currents J_{ij} = π_i w_{ij} − π_j w_{ji}. While an equilibrium ion channel satisfies detailed balance (J_{ij} = 0 on every edge), CFTR coupled to cytosolic [ATP]/([ADP][P_i]) maintains a non-zero cycle flux J_C around its gating loop with thermodynamic affinity A_C = k_B T ln(Π_{cycle} w_{ij} / w_{ji}) = Δμ_ATP ≈ 19.5 k_B T (~12 kcal/mol). In 2015, Andre Barato and Udo Seifert proved the Thermodynamic Uncertainty Relation (TUR, Phys. Rev. Lett., PMID 25933341): the relative fluctuation of any non-equilibrium current is bounded from below by the inverse of total entropy production Σ = J_C A_C / T: Var(J) / ⟨J⟩² ≥ 2 k_B / Σ.
The stochastic thermodynamics of CFTR gating is formulated on a continuous-time Markov jump process defined over four discrete conformational states: C_1 (apo closed), C_2 (ATP-bound closed), O_1 (pre-hydrolytic open), and O_2 (post-hydrolytic open). According to Jürgen Schnakenberg's cycle decomposition theorem, the thermodynamic affinity A_C associated with a directed cycle C is determined by the log-ratio of forward to reverse transition rates: A_C = k_B T ln(Π_C w_{ij}^+ / w_{ij}^-). In an equilibrium channel obeying microscopic reversibility, detailed balance requires that every individual transition flux vanishes identically (J_{ij} = 0), forcing A_C = 0. In contrast, wild-type CFTR couples the C_2 → O_1 → O_2 transition to the chemical free energy of MgATP hydrolysis. Under physiological conditions where cytosolic [ATP]/([ADP][Pi]) is maintained far from equilibrium, the cycle affinity equals the chemical potential difference Δμ_ATP ≈ 19.5 k_B T (~12.0 kcal/mol), sustaining a non-zero steady-state probability current J_C around the loop.
In 2015, Andre Barato and Udo Seifert proved the Thermodynamic Uncertainty Relation (TUR), demonstrating that for any steady-state observable current J integrated over time interval τ, the relative variance satisfies Var(J) / ⟨J⟩² ≥ 2 k_B / (Σ · τ), where Σ = J_C A_C / T is the total entropy production rate. In our computational simulation engine inside src/scripts/living-proof-observatory.js, we implement this relationship by solving the algebraic master equation Q^T π = 0 alongside the full large-deviation rate function. The kernel proves that reducing the Fano factor of single-channel dwell times requires quadratic increases in thermodynamic dissipation. When the catalytic E1371Q mutation abolishes Site-2 hydrolysis (k_hyd → 0), cycle dissipation collapses to zero (Σ → 0), restoring detailed balance and locking the channel in an artificial equilibrium trap.
Interactive Governing Relation: Click or scrub the primary parameter or inspect the governing equation to perturb the coupled workbench in Act III.
Four-Stage Mechanistic Progression
- Stage 1 (Kolmogorov–Schnakenberg Cycle Affinity): Around the 4-state gating cycle C_1 (apo) → C_2 (2ATP bound) → O_1 (open pre-hydrolytic) → O_2 (ATP hydrolyzed, Pi/ADP release) → C_1, the log-ratio of forward to reverse transition rates equals A_C / k_B T = Δμ_ATP / k_B T ≈ 19.5.
- Stage 2 (Irreversible Dwell-Time Peaking): At thermal equilibrium (A_C = 0), any open-time distribution must be a monotonically decaying sum of exponentials with coefficient of variation CV ≥ 1; ATP-driven cycle affinity A_C > 0 produces a peaked open-dwell histogram with CV < 1.
- Stage 3 (Barato–Seifert Thermodynamic Uncertainty Relation): For any integrated cycle current J_T over observation time T, relative variance satisfies Var(J_T) / ⟨J_T⟩² ≥ 2 k_B / (Σ · T), proving that higher gating regularity requires higher ATP dissipation.
- Stage 4 (Non-Hydrolytic E1371Q Mutant Limit): Mutating the catalytic glutamate in NBD2 (E1371Q, as used to trap cryo-EM structures 6MSM and 8EIQ) sets k_hyd → 0, collapsing cycle entropy production Σ → 0 and locking the channel in prolonged equilibrium-like open bursts.
Act III · Interactive Visual Experiment & Plain-Language Lens
Try It Yourself: Schnakenberg Cycle & TUR Precision Workbench
Adjust chemical driving affinity (Δμ_ATP / k_B T) and catalytic Site-2 hydrolysis rate (k_hyd) to observe how entropy production Σ and the TUR precision bound respond.
Almost every ion channel in the human body is a passive door that swings both ways in thermal equilibrium. CFTR is the lone exception: it evolved from an active pump (an ABC transporter) that burns one molecule of ATP fuel to drive its gating cycle in a one-way clockwise loop. Drag the ATP chemical drive slider below to see what happens when fuel is zero (equilibrium detailed balance, J = 0) versus when ATP hydrolysis drives directional gating.
Why this matters for research: Understanding which modulators stabilize the open state passively versus which ones alter the ATP fuel cycle helps scientists design combination potentiators.
What this diagram leaves out: Reduces CFTR's asymmetric two-ATP-site gating mechanism to a minimal 4-state Schnakenberg loop.
ATP Gating Cycle Activity · Mode: Biophysical Sandbox · How to explore this model
- Test the four physical regimes: Click the four scenario presets below (1 · Kolmogorov–Schnakenberg… through 4 · Non-Hydrolytic E1371Q…) to jump directly between experimental conditions.
- Perturb the primary control: Drag Chemical Driving Affinity (Δμ_ATP / k_B T) (1.0 to 24.0 k_B T) and observe how the visual schematic and readouts respond.
- Explore coupled co-perturbation: Adjust Catalytic Site-2 Hydrolysis Rate (k_hyd) (0.2 to 5.0 s⁻¹) to test synergy against
A_C = k_B T ln(Π_C wᵢⱼ⁺ / w_{ij}^−) = Δμ_ATP · Var(J) / ⟨J⟩² ≥ 2 k_B / Σ_cycle.
Governing law: A_C = k_B T ln(Π_C wᵢⱼ⁺ / w_{ij}^−) = Δμ_ATP · Var(J) / ⟨J⟩² ≥ 2 k_B / Σ_cycle
A_C = k_B T ln(Π_C wᵢⱼ⁺ / w_{ij}^−) = Δμ_ATP · Var(J) / ⟨J⟩² ≥ 2 k_B / Σ_cycle
Normalized Sensitivity (Stages 1–4)
Want to test this mechanism inside the full multi-module pillar observatory? Open the coupled interactive workbench on the Pillar Hub.
Act IV · Clinical & Therapeutic Horizon
Inferring Hidden Drug Mechanisms from Dwell-Time Fluctuations
Because the TUR links observable current fluctuations (Fano factor / randomness parameter) to microscopic cycle affinity and state count, electrophysiologists can determine whether a new potentiator or co-potentiator acts by slowing hydrolytic closure (O_1 → O_2) or by accelerating non-equilibrium NBD dimerization (C_1 → C_2)—even when intermediate states are too brief to resolve directly.
Applying non-equilibrium thermodynamics to CFTR pharmacology provides critical diagnostic clarity for gating mutations that preserve normal membrane trafficking but impair channel opening. In patients carrying the Class III G551D variant, glycine is replaced by aspartate at the NBD1–NBD2 dimer interface, preventing ATP-dependent dimerization and reducing open probability by more than ninety-five percent. Modulator therapy with ivacaftor acts as a potentiator that bypasses the natural ATP requirement, stimulating channel opening even in the absence of nucleotide binding. However, thermodynamic analysis demonstrates that ivacaftor-rescued gating exhibits fundamentally different dwell-time variance compared to wild-type channels, reflecting an altered energy landscape that does not restore normal enzymatic turnover. Furthermore, because ivacaftor stabilizes an open state with non-physiological kinetics, long-term exposure may alter channel turnover and downstream cellular homeostasis in ways that equilibrium binding models fail to capture.
These biophysical distinctions define the exact boundaries of stochastic modeling. While TUR bounds enable researchers to detect whether a novel potentiator acts catalytically or through passive allosteric trapping, they do not predict whole-lung clinical outcomes or replace empirical pharmacokinetic testing. Furthermore, non-equilibrium gating models are inherently inapplicable to Class I nonsense genotypes or severe processing mutations where no functional protein reaches the plasma membrane. In our open science workflow, every thermodynamic parameter is cross-checked against patch-clamp recordings from literature sources verified by Scite smart citations and DOI receipts. By honoring these experimental boundaries, our models provide rigorous mechanistic insights without offering unverified clinical assertions.
Primary Anchor: Barato & Seifert (2015), Physical Review Letters 114(15):158101 — Thermodynamic Uncertainty Relation for Biomolecular Processes (PMID 25933341). (DOI 10.1103/PhysRevLett.114.158101 · PubMed 25933341).
Governed by biophysical kernel src/scripts/living-proof-observatory.js and cross-checked in the Open Non-Equilibrium Stochastic Thermodynamics Observatory on /math/.
Resumen Científico en Español (ES-MX) · Folio L.02
Termodinámica Estocástica Fuera del Equilibrio en CFTR: Ciclos de Schnakenberg, Producción de Entropía y Cotas TUR: Por qué CFTR rompe el balance detallado mediante la hidrólisis de MgATP (driven by roughly twenty thermal energy units from adenosine triphosphate hydrolysis) y cómo la Relación de Incertidumbre Termodinámica acota las fluctuaciones.
A diferencia de los canales iónicos en equilibrio térmico, CFTR consume energía química (hidrolizando MgATP en el Sitio 2) para impulsar un ciclo unidireccional fuera del equilibrio, lo que permite acotar la precisión del canal mediante la producción de entropía.
En Palabras Sencillas: Las células vivas dependen de interruptores moleculares confiables para regular el flujo de iones. A diferencia de los canales pasivos que fluctúan aleatoriamente sin consumir recursos metabólicos, el canal de la fibrosis quística actúa como una bomba activa que consume combustible celular durante cada ciclo. Este consumo continuo de energía genera una dirección temporal definida. Debido a que la proteína avanza a través de formas estructurales distintas, el tiempo que permanece abierta no sigue la distribución sin memoria de compuertas térmicas pasivas. Los físicos describen este comportamiento fuera del equilibrio mediante el concepto de producción de entropía. La relación de incertidumbre termodinámica revela que reducir el ruido en corrientes de canal único exige un mayor gasto de energía libre química por parte de la proteína operativa. Cada incremento en la precisión regulatoria conlleva un costo metabólico ineludible. Al examinar las fluctuaciones en registros eléctricos de patch clamp, los investigadores pueden deducir si una molécula altera el consumo de combustible o estabiliza estados intermedios.
Analogía Cotidiana: Considere una rueda hidráulica de madera montada sobre un canal de corriente rápida donde la caída del agua obliga a las paletas a girar en una sola dirección. La corriente cinética constante suministra energía mecánica continua, garantizando que los cangilones eleven agua hacia campos más altos con regularidad predecible. Si el estanque se seca y la rueda reposa en agua estancada, las brisas térmicas aleatorias podrían mover las paletas hacia adelante y hacia atrás, pero el dispositivo pierde toda capacidad de bombear agua cuesta arriba. En biofísica celular, el canal CFTR funciona exactamente como esta rueda hidráulica, utilizando el flujo químico de la hidrólisis de ATP para mantener un ciclo de compuerta unidireccional con alta precisión operativa. Sin ese aporte energético, el canal cae en fluctuaciones erráticas y reversibles que no logran sostener gradientes iónicos. Sin embargo, esta analogía mecánica cesa de ser precisa a escala molecular, porque las ruedas hidráulicas microscópicas operan bajo abrumadoras colisiones brownianas que intentan revertir continuamente las paletas.
Guía de Controles Interactivos: Antes de analizar las ecuaciones maestras de redes fuera del equilibrio, los lectores pueden explorar la termodinámica de compuerta mediante dos controles interactivos. El primer control, la Afinidad Química Impulsora, establece la diferencia de potencial químico suministrada por el combustible celular de ATP en relación con la energía térmica. El segundo control, la Tasa de Hidrólisis Catalítica del Sitio 2, dicta la rapidez con la que la proteína escinde la molécula de combustible para iniciar el cierre del canal. Ajustar estas dos entradas actualiza inmediatamente cuatro lecturas termodinámicas fundamentales. La lectura de Afinidad del Ciclo de Schnakenberg rastrea la fuerza impulsora química neta alrededor del ciclo de compuerta de cuatro estados. La Producción de Entropía en Estado Estacionario indica la velocidad a la que la energía metabólica se disipa en el entorno celular. La Cota Inferior de Precisión TUR calcula el límite mínimo de varianza para corrientes iónicas. Finalmente, el Estado del Balance Detallado confirma si el canal opera como un motor irreversible o ha colapsado hacia el equilibrio térmico.
Riesgo Biofísico y Clínico: La demostración matemática de la relación de incertidumbre termodinámica estableció un equilibrio fundamental entre la precisión cinética y la producción de entropía en todos los procesos markovianos de tiempo continuo. Para los biofísicos que estudian la fibrosis quística, este principio transforma la interpretación de los registros de electrofisiología de canal único. CFTR de tipo silvestre no opera como un poro pasivo en equilibrio; más bien, hidroliza MgATP en el Sitio 2 con una afinidad impulsora de aproximadamente 19.5 k_B T. Este impulso químico rompe el balance detallado microscópico y sostiene un flujo de probabilidad irreversible alrededor del ciclo de compuerta. Cuando mutaciones de compuerta como G551D desacoplan la unión de ATP de la apertura del poro, el canal pierde su regularidad direccional y falla al conducir corrientes fisiológicas de cloruro. Las implicaciones traslacionales de este marco fuera del equilibrio radican en el desarrollo racional de fármacos. Al analizar la varianza y la dispersión de Fano de los tiempos de residencia en patch-clamp, los investigadores pueden determinar si potenciadores de moléculas pequeñas como ivacaftor restauran el recambio catalítico normal o simplemente atrapan al canal en una conformación abierta artificial. Esta claridad cuantitativa evita conclusiones engañosas durante el cribado preclínico de fármacos.
Derivación Formal: La termodinámica estocástica de la compuerta de CFTR se formula sobre un proceso de saltos de Markov en tiempo continuo definido en cuatro estados conformacionales discretos: C_1 (cerrado apo), C_2 (cerrado con ATP unido), O_1 (abierto prehidrolítico) y O_2 (abierto poshidrolítico). Según el teorema de descomposición de ciclos de Jürgen Schnakenberg, la afinidad termodinámica A_C asociada a un ciclo dirigido C está determinada por la razón logarítmica de las tasas de transición directa e inversa: A_C = k_B T ln(Π_C w_{ij}^+ / w_{ij}^-). En un canal en equilibrio que obedece reversibilidad microscópica, el balance detallado exige que cada flujo de transición individual se anule idénticamente (J_{ij} = 0), forzando A_C = 0. Por el contrario, CFTR de tipo silvestre acopla la transición C_2 → O_1 → O_2 a la energía libre química de la hidrólisis de MgATP. En condiciones fisiológicas donde [ATP]/([ADP][Pi]) citosólico se mantiene lejos del equilibrio, la afinidad del ciclo equivale a la diferencia de potencial químico Δμ_ATP ≈ 19.5 k_B T (~12.0 kcal/mol), sosteniendo una corriente de probabilidad estacionaria no nula J_C alrededor del circuito. En 2015, Andre Barato y Udo Seifert demostraron la Relación de Incertidumbre Termodinámica (TUR), evidenciando que para cualquier corriente observable en estado estacionario J integrada en un intervalo de tiempo τ, la varianza relativa satisface Var(J) / ⟨J⟩² ≥ 2 k_B / (Σ · τ), donde Σ = J_C A_C / T es la tasa total de producción de entropía. En nuestro motor de simulación computacional dentro de src/scripts/living-proof-observatory.js, implementamos esta relación resolviendo la ecuación maestra algebraica Q^T π = 0 junto con la función de tasa de grandes desviaciones. El núcleo demuestra que reducir el factor de Fano de los tiempos de residencia en canal único exige incrementos cuadráticos en la disipación termodinámica. Cuando la mutación catalítica E1371Q anula la hidrólisis del Sitio 2 (k_hyd → 0), la disipación del ciclo colapsa a cero (Σ → 0), restaurando el balance detallado y bloqueando al canal en una trampa de equilibrio artificial.
Frontera Traslacional: Aplicar la termodinámica fuera del equilibrio a la farmacología de CFTR proporciona una claridad diagnóstica crítica para mutaciones de compuerta que conservan el tráfico normal a la membrana pero alteran la apertura del canal. En pacientes que portan la variante de Clase III G551D, la glicina se sustituye por aspartato en la interfaz del dímero NBD1–NBD2, impidiendo la dimerización dependiente de ATP y reduciendo la probabilidad de apertura en más del noventa y cinco por ciento. La terapia moduladora con ivacaftor actúa como un potenciador que elude el requerimiento natural de ATP, estimulando la apertura del canal incluso en ausencia de unión a nucleótidos. Sin embargo, el análisis termodinámico demuestra que la compuerta rescatada por ivacaftor exhibe una varianza del tiempo de residencia fundamentalmente distinta en comparación con los canales de tipo silvestre, reflejando un paisaje energético alterado que no restaura el recambio enzimático normal. Estas distinciones biofísicas definen los límites exactos del modelado estocástico. Si bien las cotas TUR permiten a los investigadores detectar si un nuevo potenciador actúa catalíticamente o mediante atrapamiento alostérico pasivo, no predicen desenlaces clínicos pulmonares globales ni reemplazan las pruebas farmacocinéticas empíricas. Además, los modelos de compuerta fuera del equilibrio son inherentemente inaplicables a genotipos sin sentido de Clase I o mutaciones de procesamiento graves donde ninguna proteína funcional llega a la membrana plasmática. En nuestro flujo de trabajo de ciencia abierta, cada parámetro termodinámico se coteja con registros de patch-clamp verificados mediante citas inteligentes de Scite y recibos DOI. Al respetar estos límites experimentales, nuestros modelos proporcionan conocimientos mecanicistas rigurosos sin ofrecer garantías clínicas no verificadas.
Horizonte terapéutico: Conecta los teoremas de la mecánica estadística fuera del equilibrio (descomposición de Schnakenberg y cota TUR de Barato–Seifert) con los registros experimentales de patch-clamp en CFTR. Investigación independiente, no asesoramiento médico.
Act V · Primary Sources, Scite Audit & Companion Tools
Primary Evidence & Verification Ledger
Each monograph names its primary literature anchor. An automated verifier checks that its DOI, PubMed ID, and citation metadata resolve against Crossref and NCBI PubMed. Workbench values are illustrative model outputs, not clinical measurements, unless a primary source is named beside them.
- Primary Literature Receipt: Barato & Seifert (2015), Physical Review Letters 114(15):158101 — Thermodynamic Uncertainty Relation for Biomolecular Processes (PMID 25933341). (DOI 10.1103/PhysRevLett.114.158101 · PubMed 25933341).
- Citation Context (Scite): Inspect how later literature cites 10.1103/PhysRevLett.114.158101 on Scite.ai ·
Crossref + Scite verified (data/folio-citations.cache.json · fetched 2026-11-08 · 1,293 citing · 31 supporting · 1,180 mentioning · 0 contrasting). - Open-Source Computational & Epistemic Engine: COMPAS & Python Optimal Transport (POT) Wasserstein Geodesics (DOI 10.1016/j.cad.2015.01.006 · Scite Report).
- Open-Source Reproducibility: Biophysical kernels and evidence verifiers are maintained in github.com/senseiewok/cf-research and github.com/senseiewok/cf-skills.
For builders: method, companion SKILL.md, Python verifier, and edition history (REV 1.2)
cf-stochastic-thermodynamics-tur-SKILL.md · Kernel src/scripts/living-proof-observatory.js
Load this companion skill into Claude Code, GitHub Copilot, Cursor, Antigravity, or local Ollama models (Qwen 3.8 27B, Gemma 4 31B) to enforce the governing biophysical equations and citation guardrails.
View raw SKILL.md source
---
name: cf-stochastic-thermodynamics-tur
description: AI-agnostic statistical physics skill for Schnakenberg Markov cycle affinity, non-equilibrium ATP hydrolysis in CFTR, and Barato-Seifert Thermodynamic Uncertainty Relation (TUR) bounds.
license: CC-BY-4.0
---
# L.02 — Schnakenberg Cycle Affinity, Entropy Production & TUR Bounds
## Governing Biophysical / Mathematical Law
A_C = k_B T ln(Π_C wᵢⱼ⁺ / w_{ij}^−) = Δμ_ATP · Var(J) / ⟨J⟩² ≥ 2 k_B / Σ_cycle
## Primary Evidence Receipt
- **Citation:** Barato & Seifert (2015), Physical Review Letters 114(15):158101 — Thermodynamic Uncertainty Relation for Biomolecular Processes (PMID 25933341).
- **DOI:** https://doi.org/10.1103/PhysRevLett.114.158101
- **PubMed PMID:** https://pubmed.ncbi.nlm.nih.gov/25933341/
- **Answer Tier:** T2 · Sourced Stochastic Thermodynamics Model
## Verification Rules
- Compute Schnakenberg cycle affinity as A_C = k_B * T * ln(product(w_forward / w_reverse)) = Delta mu_ATP.
- Enforce the Barato-Seifert TUR bound Var(J) / <J>^2 >= 2 * k_B / Sigma for any non-equilibrium gating current.
- Explain why E1371Q abolishes Site-2 hydrolysis (k_hyd -> 0), restoring detailed balance for cryo-EM trapping (PDB 6MSM, 8EIQ).
## Safety & Clinical Boundary
Research synthesis and biophysical modeling only — never medical advice, clinical guidance, or a treatment promise.
REV 1.2 (2026-11-04): Informatico magazine parity, Split Workbench (Plates A & B + 4-stage stepper), and live Crossref/PubMed citation gate verification.
Verified 5-Act magazine flow, untruncated headings, and bilingual EN / ES-MX structural parity.
Audited clinical dignity boundaries, non-promise guardrails, and discovery catalog coverage.
check_folio_citations.py
Verified DOI 10.1103/PhysRevLett.114.158101 against Crossref and NCBI PubMed metadata.
src/scripts/living-proof-observatory.js
Coupled governing relation A_C = k_B T ln(Π_C wᵢⱼ⁺ / w_{ij}^−) = Δμ_ATP to interactive Plates A & B.