RNA Riboswitches · NGS Bridge Amplification · Microtubule Curvature (Domains 1–3)
Three companion preprints applying the TOGT/GTCT contact-geometric operator framework (Principia Orthogona series, DOI 10.5281/zenodo.19117399) to three biological domains. Each preprint contains the problem statement, operator mapping, falsifiable quantitative predictions, experimental protocols, and 6-panel domain-specific figures.
Domain 1: Aminoglycoside resistance riboswitches (AAC/AAD leader RNA; κ* modulation; antagonist OFF-lock predictions).
Domain 2: NGS bridge amplification on Illumina flow cells (κNGS derivable from DNA persistence length and primer spacing; Whitney A₁ bifurcation at cluster quality boundary).
Domain 3: Microtubule curvature dynamics in volatile anesthesia and neuroprotection (classical κ shift mechanism distinguishable from Orch-OR; epothilone B LORR delay prediction).
Also includes the dm³ Python simulation (DOP853, rtol=1e-10) reproducing all figures. Lean 4 source: github.com/TOTOGT/AXLE.
This preprint series applies the Topographical Orthogenetic Theory (TOGT) and Generative Time Circuit Theorem (GTCT) — a formally verified contact-geometric framework — to three biological systems sharing a common operator structure. The framework is part of the Principia Orthogona series, a multi-volume mathematical research programme by G6 LLC.
All three domains are mapped to the same non-commutative operator chain on a contact 3-manifold (dm³):
where C = Compression, K = Curvature intensification, F = Fold (Whitney A₁ singularity at κ*), U = Unfold/Stabilization. The fold operator F is the critical event: the Whitney A₁ singularity at the curvature threshold κ* produces the characteristic bifurcation observed in all three biological systems.
The mathematical substrate is a three-dimensional contact manifold with:
The canonical invariant triple is \((T^*, \mu_{\max}, \tau) = (2\pi, -2, 2)\) with stability radius[Ch 10] \(\varepsilon_0 = 1/3\).
Aminoglycoside resistance riboswitches in the 5' leader of aac/aad resistance genes. Operator mapping of the OFF→ON conformational switch. Three falsifiable predictions for antagonist design.
NGS bridge amplification on Illumina flow cells. κNGS derivable from DNA persistence length and primer spacing without free parameters. Whitney A₁ bifurcation at cluster quality boundary.
Microtubule curvature dynamics in volatile anesthesia, neuroprotection, and tauopathies. Classical κ shift mechanism distinguishable from Orch-OR. Epothilone B LORR delay prediction.
The central mathematical result unifying all three domains:
Any two dm³ systems with matching invariants \((\mu_{\max}, \omega, \kappa^*)\) are categorically equivalent under explicit contact morphisms. That is, if
then there exists an explicit contact morphism \(\Phi: M_1 \to M_2\) such that \(\Phi^* \alpha_2 = f \cdot \alpha_1\) for some smooth positive function \(f\), and the operator chains intertwine: \(\Phi \circ G_1 = G_2 \circ \Phi\).
Consequence: predictions derived in one domain transfer to any other domain with matching invariants, subject to explicit morphism scaling.
| Component | Lean File | Status |
|---|---|---|
| Core operator chain (C, K, F, U) | AXLE_v5_1.lean | 0 sorry |
| Coherence Bridge (Theorem 5.4) | AXLE_v5_1.lean | Verified |
| Whitney A₁ singularity conditions | AutophagyDm3_v2.lean | 18 theorems |
| Gronwall radius ε₀ = 1/3 | Multiple | Proved |
| Domain axiom: Mather step | AXLE Issue #14 | Open |
| Domain axiom: Poincaré–Bendixson | AXLE Issue #14 | Open |
| File | Description | Size |
|---|---|---|
dm3_simulation.py | Python DOP853 reference integrator (rtol=1e-10), generates all 4 figures | 14.7 kB |
fig_domain1_riboswitch.pdf | 6-panel domain figure — Riboswitch | 70.6 kB |
fig_domain2_ngs.pdf | 6-panel domain figure — NGS | 63.1 kB |
fig_domain3_mt.pdf | 6-panel domain figure — Microtubule | 69.5 kB |
togt_domain1_riboswitch.pdf | Full paper — Domain 1 | 349.4 kB |
togt_domain2_ngs.pdf | Full paper — Domain 2 | 331.0 kB |
togt_domain3_mt.pdf | Full paper — Domain 3 | 347.8 kB |
| Role | DOI / URL |
|---|---|
| Series root / concept DOI | 10.5281/zenodo.19117399 |
| Principia Orthogona Vol. I | 10.5281/zenodo.20298665 |
| Principia Orthogona Vol. II | 10.5281/zenodo.20755436 |
| GCM paper (dm³ toy model) | 10.5281/zenodo.19379385 |
| GTCT / Ring 5 | 10.5281/zenodo.20239928 |
| AXLE formal verification hub | github.com/TOTOGT/AXLE |
TOGT/GTCT — Domain 1 of 3
Aminoglycoside-sensing riboswitches in the 5' leader of aac/aad resistance genes undergo an inducible conformational switch — OFF (SD2 sequestered) to ON (SD2 exposed) — that activates resistance enzyme expression only in the presence of drug. No crystal structure of the full ON state exists; the switching kinetics have not been measured directly; and the universality of the switch across integron contexts is not established.
We apply the TOGT/GTCT contact-geometric operator framework to map the riboswitch conformational cycle to the operator chain \(G = U \circ F \circ K \circ C\) on a dm³ contact manifold. The fold operator \(F\) (Whitney A₁ singularity) maps to the stem competition bifurcation; the curvature threshold \(\kappa^*\) maps to the free-energy barrier between OFF and ON stems.
We derive three falsifiable quantitative predictions: (P1) the switching free-energy barrier is linear in \(\kappa^*_{\text{ribo}}\), measurable by SHAPE/DMS probing; (P2) an effective antagonist raises \(\kappa^*\) at the fold point — steric bulk at the switching junction suffices, full pocket occupancy is not necessary; (P3) the dose-response curve (OFF-lock probability vs. antagonist concentration) follows a dm³ sigmoid with Hill coefficient determined by \(\mu_{\max} = -2\) and \(\kappa^*_{\text{ribo}}\).
Riboswitches de aminoglicosídeos na região 5' dos genes de resistência aac/aad realizam uma troca conformacional induzível — OFF (SD2 sequestrado) para ON (SD2 exposto) — que ativa a expressão da enzima de resistência somente na presença do fármaco. Aplicamos o arcabouço TOGT/GTCT de geometria de contato formalmente verificada — cadeia de operadores \(G = U\circ F\circ K\circ C\) em uma variedade dm³ com forma de contato \(\alpha = dz - r^2 d\theta\) — ao ciclo conformacional do riboswitch.
Aminoglycosides (gentamicin, tobramycin, amikacin) are broad-spectrum antibiotics targeting the 16S rRNA decoding site. Resistance is primarily mediated by acetyltransferases (AAC) and adenylyltransferases (AAD) encoded in class 1 integrons. Jia et al. [2013] established that these genes carry a structured RNA element in their 5' leader that acts as an inducible riboswitch: the anti-Shine–Dalgarno helix (anti-SD2) sequesters the ribosome binding site in the drug-free state (OFF), and drug binding triggers a conformational shift that exposes SD2 (ON), enabling resistance enzyme expression only under selective pressure.
This inducibility minimises the metabolic cost of resistance in drug-free environments — a key factor in integron persistence in clinical settings across P. aeruginosa, E. coli, K. pneumoniae, and Campylobacter spp.
Four critical gaps remain open:
(1) ON-state structure — No crystal structure of the full ON-state conformation has been reported; the aptamer binding geometry is inferred from probing data.
(2) Switching kinetics — The rate of the aptamer→expression-platform conformational change has not been directly measured.
(3) Cross-integron universality — Whether \(\kappa^*_{\text{ribo}}\) is conserved across integron contexts and host organisms is unknown.
(4) Controversy — Roth and Breaker [2013] argued that integron attI1 site association, not RNA switching, explains resistance gene expression patterns. If correct, the riboswitch interpretation requires revision.
A small-molecule resistance blocker must: (a) bind the riboswitch with sufficient affinity; (b) lock the OFF state under physiological drug concentrations; (c) show selectivity over endogenous RNA; and (d) be cell-permeable and non-toxic. Ribocil [Howe et al., 2015] established proof-of-concept for the FMN riboswitch. No equivalent has been published for the aminoglycoside riboswitch.
The dm³ contact manifold \(M = \mathbb{R}^2_{>0} \times \mathbb{R}\) is equipped with contact form:
The operator chain acts on trajectories in \(M\):
The canonical invariant triple is \((T^*, \mu_{\max}, \tau) = (2\pi, -2, 2)\) with:
The fold operator \(F\) is characterized by the Whitney A₁ singularity at \(\kappa^*\). The potential \(V(q) = q^3 - 3q\) at \(q = 1\) satisfies:
giving \(\mu_{\text{canonical}} = -V''(1)/2 = -3\), rescaled to \(\mu_{\max} = -2\) in the dm³ toy model.
| Operator | Biological Realisation |
|---|---|
| C (Compression) | Aminoglycoside drug binds aptamer domain; RNA begins to deform |
| K (Curvature intensification) | Stem competition increases — anti-SD2 vs. aptamer-stabilised helix compete |
| F (Fold at κ*) | Whitney A₁ bifurcation at stem competition threshold; SD2 sequestration breaks |
| U (Unfold/Stabilise) | Expression platform stabilises in ON conformation; ribosome accesses SD2 |
The curvature threshold \(\kappa^*_{\text{ribo}}\) is identified with the free-energy barrier between OFF and ON conformations:
where \(\ell_{\text{helix}}\) is the characteristic helix length at the switching junction. This is a dimensionless curvature invariant and is directly measurable from SHAPE/DMS probing experiments.
The switching free-energy barrier \(\Delta G^{\ddagger}_{\text{OFF→ON}}\) is linear in \(\kappa^*_{\text{ribo}}\):
where \(\lambda_{\text{ribo}}\) is the dm³ Lipschitz constant for the riboswitch domain and \(c_0\) is the baseline barrier.
Experimental test: SHAPE/DMS probing of aac/aad leader RNA at varying magnesium concentrations (which tune \(\kappa^*_{\text{ribo}}\)) should yield a linear relationship between the SHAPE reactivity at the switching junction and the free-energy barrier estimated by smFRET.
Falsification criterion: A statistically significant (p < 0.05) nonlinear dependence (e.g., quadratic term dominant) would falsify P1 within the dm³ framework.
An effective antagonist raises \(\kappa^*\) at the fold point. Steric bulk at the switching junction suffices — full pocket occupancy is not necessary.
The antagonist efficacy \(E\) satisfies:
where \(\Delta\kappa^* = \kappa^*_{\text{blocked}} - \kappa^*_{\text{ribo}} > 0\) is the increase in curvature threshold induced by the antagonist.
Experimental test: Screen steric blockers (e.g., PNA clamps, intercalators) at the switching junction using an in-vitro translation (IVT) assay reporting SD2 accessibility. Efficacy should correlate with \(\Delta\kappa^*\) estimated from SHAPE, not with binding affinity to the aptamer pocket.
Falsification criterion: If a high-affinity aptamer binder consistently outperforms a steric junction blocker with lower aptamer affinity, P2 is falsified.
The dose-response curve (OFF-lock probability vs. antagonist concentration) follows a dm³ sigmoid with Hill coefficient determined (not fitted) by \(\mu_{\max} = -2\):
The Hill coefficient \(n \approx 3.64\) is derived from the dm³ invariants, not fitted to data. This is the strongest prediction of the framework — a sharp, numerical, parameter-free claim.
Experimental test: Measure the dose-response of an OFF-lock compound using an IVT or cell-based luminescence assay at ≥ 8 concentrations spanning at least 3 decades. Fit a Hill equation and compare the fitted \(n\) to the dm³ prediction.
Falsification criterion: A fitted Hill coefficient \(n < 2.5\) or \(n > 5\) would falsify P3 at the dm³ level (the prediction tolerates ±30% for domain-specific corrections).
The core mathematical claims are verified in Lean 4 / Mathlib4 via the AXLE repository. The domain-specific
axioms (Mather's theorem for the Whitney fold, and Poincaré–Bendixson for limit cycle existence) remain as
open obligations, declared with trivial placeholders, not hidden sorry.
| # | Theorem | Statement | Status |
|---|---|---|---|
| 1 | contactCoeff_neg | c(ρ) = −2ρ < 0 for ρ > 0 | Proved |
| 2 | V_critical_at_one | V′(1) = 0 | Proved |
| 3 | V_factored | V(q)+2 = (q−1)²(q+2) | Proved |
| 4 | mu_dm3_neg | μ_max = −2 < 0 (transverse attraction) | Proved |
| 5 | gronwall_radius | ε₀ = 1/3 | Proved |
| 6 | whitneyFold_from_kinase_data | Mather's theorem + biological data | Open (trivial) |
| 7 | limitCycle_exists_ribo | Poincaré–Bendixson or Lyapunov construction | Open (trivial) |
Source: fig_domain1_riboswitch.pdf (Zenodo 10.5281/zenodo.20559510) · Generated by dm3_simulation.py (DOP853, rtol=1e-10)
The dm³ framework provides a principled geometric account of the riboswitch switch that is both mathematically rigorous and experimentally falsifiable. The key strength is that all three predictions involve specific numerical values or functional forms derivable from first principles — not fitted parameters.
The Hill coefficient prediction (n ≈ 3.64) is particularly diagnostic: it differs from the empirically common n = 2 (cooperative dimer) and n = 4 (cooperative tetramer), and its origin in \(\mu_{\max} = -2\) and the dm³ geometry is traceable. A confirmation would strongly support the contact-geometric framework; a falsification would constrain its domain of applicability.
The antagonist mechanism prediction (P2) is clinically significant: it implies that junction-targeted steric blockers — potentially simpler to design than aptamer-pocket binders — may be sufficient for OFF-lock activity. This opens a new structural class for antibiotic adjuvant design.
We have applied the TOGT/GTCT contact-geometric operator framework to the aminoglycoside riboswitch conformational cycle and derived three falsifiable quantitative predictions: (P1) free-energy barrier linearity in κ*_ribo, (P2) steric-bulk sufficiency for antagonist efficacy, and (P3) the dm³ sigmoid with Hill coefficient n ≈ 3.64. Each prediction comes with an explicit experimental protocol and falsification criterion. The mathematical claims are machine-verified in Lean 4 with zero sorry in the core chain. Two domain-specific open obligations are clearly marked.
We present the Topographical Orthogenetic Theory (TOGT) and Generative Time Circuit Theorem (GTCT) — a formally verified contact-geometric framework — and apply it to three biological systems sharing a common operator structure: (i) aminoglycoside resistance riboswitches, (ii) NGS bridge amplification on Illumina flow cells, and (iii) microtubule curvature dynamics in volatile anesthesia and neuroprotection.
The mathematical substrate is a three-dimensional contact manifold (dm³) with contact form \(\alpha = dz - r^2 d\theta\) and non-commutative operator chain \(G = U \circ F \circ K \circ C\), where the fold operator \(F\) is a Whitney A₁ singularity at the curvature threshold \(\kappa^*\). The Coherence Bridge Theorem (Theorem 5.4), machine-verified in Lean 4 / Mathlib4 with zero sorry obligations in the core chain, proves that any two dm³ systems with matching invariants \((\mu_{\max}, \omega, \kappa^*)\) are categorically equivalent under explicit contact morphisms.
Each biological domain is mapped to the operator chain; cross-domain transfer predictions are stated as falsifiable hypotheses with explicit experimental protocols and falsification criteria. The paper includes dm³ phase portrait and spiral return simulations (DOP853, rtol=1e-10), tikz operator diagrams, and domain-specific figures. The framework is distinguished from quantum-mechanical models of consciousness (Orch-OR): TOGT operates at the level of classical contact geometry. Two domain axioms remain explicitly open.
A amplificação em ponte (bridge amplification) em flow cells da Illumina — tecnologia reconhecida pelo Prêmio Princesa das Astúrias 2026 — é um processo topográfico de superfície cujo desempenho depende de dois limiares geométricos: uma densidade mínima de primers abaixo da qual a formação de clusters falha, e uma densidade máxima acima da qual o carregamento policlonal domina. Nenhum dos dois limiares foi derivado a partir de primeiros princípios geométricos.
Aplicamos o arcabouço TOGT/GTCT de geometria de contato formalmente verificada — cadeia de operadores \(G = U\circ F\circ K\circ C\) em uma variedade dm³ com forma de contato \(\alpha = dz-r^2d\theta\) — à amplificação em ponte. O operador de dobramento F (singularidade Whitney A₁) mapeia-se ao ponto de bifurcação de ramificação de clusters; o limiar de curvatura \(\kappa^*\) mapeia-se à densidade crítica de primers.
The dm³ system is defined on the contact 3-manifold \(M = \mathbb{R}^2_{>0} \times \mathbb{R}\) with coordinates \((r, \theta, z)\) and contact form:
The dm³ toy model equations (exact, Principia Orthogona Vol. II §4.3):
Parameters: \((\mu_{\max}, \omega, \beta) = (-2, 1, 1)\). The limit cycle is \(\Gamma = \{r=1\}\) with period \(T^* = 2\pi\). The canonical invariant triple is \((T^*, \mu_{\max}, \tau) = (2\pi, -2, 2)\).
The non-commutative operator chain \(G = U \circ F \circ K \circ C\) acts on trajectories in \(M\):
| Operator | Name | Mathematical Role | Geometric Action |
|---|---|---|---|
| C | Compression | Contractive, injective map | Reduces r toward the fold threshold |
| K | Curvature | Curvature intensification | Increases curvature; drives system toward κ* |
| F | Fold | Whitney A₁ singularity at κ* | Non-injective; produces bifurcation at fold |
| U | Unfold | Φ-decrease + stable branch selection | Stabilises onto the limit cycle Γ |
Non-commutativity theorem: \(F \circ K \neq K \circ F\) — operator order is physically meaningful. Swapping C and K produces different biological outcomes (as demonstrated empirically in zeolite catalysis, Domain 0 of the Principia Orthogona series).
The global attractor of the full dm³ system is the resonant orbit \(\Gamma_{12}\). Any trajectory starting in the basin \(\{r^* < r < r_{\text{outer}}\}\) converges to \(\Gamma\) with exponential rate \(|\mu_{\max}| = 2\).
\(|\kappa| \uparrow \kappa^* \iff \mu_{\max} < 0 \iff \tau = \sqrt{c/\kappa_{\text{noise}}} \in (0, \infty)\). The curvature threshold \(\kappa^*\) and the embodiment threshold \(\tau\) are two parameterizations of the same event.
The four dm³ bifurcations (Contact Hopf, Saddle-node, Neimark–Sacker, Slow-fast crossover) correspond bijectively to the Whitney A₁–A₃ singularity types.
The stationary SDE measure concentrates on \(\Gamma\) for noise amplitude below the embodiment threshold \(\tau = 2\) and spreads above \(\tau\).
Any two dm³ systems with matching invariants \((\mu_{\max}, \omega, \kappa^*)\) are categorically equivalent under explicit contact morphisms. Cross-domain transfer predictions follow as corollaries.
NGS bridge amplification on Illumina flow cells is a surface-based DNA amplification technique in which primer-coated flow cell surfaces support the formation of DNA clusters through sequential extension of bridged single-stranded templates. The technology was recognised by the Premio Princesa de Asturias 2026.
Two critical geometric thresholds govern cluster formation:
Lower threshold: Minimum primer density below which cluster formation fails — insufficient template bridging probability.
Upper threshold: Maximum primer density above which polyclonal loading dominates — multiple templates per cluster region degrade sequencing quality.
Neither threshold has been derived from first geometric principles.
| Operator | NGS Biological Realisation |
|---|---|
| C | DNA denaturation — template collapses to surface as single strand |
| K | Template bending — ssDNA curves toward surface-anchored primer; persistence length effect |
| F | Whitney A₁ bifurcation at κ*_NGS — cluster branching threshold; primer density bifurcation point |
| U | Cluster stabilisation — amplified cluster reaches stable monoclonal configuration |
The curvature threshold for NGS bridge amplification is derivable from first principles:
where \(L_p\) is the DNA persistence length (≈ 50 nm for double-stranded DNA, ≈ 1 nm for single-stranded DNA) and \(d_{\text{primer}}\) is the primer surface spacing.
\(\kappa^*_{\text{NGS}}\) is derivable from DNA persistence length and primer spacing without free parameters:
Experimental test: Vary flow cell primer density systematically using diluted primer mixtures. Measure cluster formation efficiency vs. density. The onset of failure should occur at \(d_{\text{primer}} \approx 1/\sqrt{\kappa^*_{\text{NGS}} \cdot L_p}\).
Falsification: If the experimentally measured critical density deviates from the prediction by more than 50% across three independent flow cell preparations, P1 is falsified.
The cluster quality curve (quality score vs. primer density) exhibits a Whitney A₁ bifurcation structure with asymmetry ratio ≥ 2:
Experimental test: Scan primer density from 0.1× to 10× optimal across ≥ 8 points. Fit a cusp catastrophe or Whitney A₁ profile to the quality vs. density curve.
Falsification: A symmetric quality curve (ratio < 1.2) or a monotone relationship (no bifurcation) would falsify P2.
The functional form of the cluster quality vs. primer density curve is identical (up to contact morphism scaling) to the riboswitch OFF-lock dose-response curve (Domain 1):
where \(\Phi_{\text{NGS→Ribo}}\) is the explicit contact morphism guaranteed by the Coherence Bridge Theorem (Theorem 5.4).
Experimental test: After fitting the riboswitch dose-response (Domain 1), use the dm³ morphism to predict the NGS quality curve shape without additional free parameters.
Falsification: A statistically significant mismatch (χ² test, p < 0.01) between the morphism-scaled riboswitch curve and the NGS quality curve would falsify P3.
The three domains share the same invariant triple \((\mu_{\max}, \omega, \kappa^*)\) under appropriate scaling. This is the content of the Coherence Bridge Theorem:
| Invariant | Domain 1 (Riboswitch) | Domain 2 (NGS) | Domain 3 (Microtubule) |
|---|---|---|---|
| \(\mu_{\max}\) | −2 (canonical) | −2 (canonical) | −2 (canonical) |
| \(\omega\) | ≈ 1 (aptamer oscillation) | ≈ 1 (cluster oscillation) | ≈ 1 (MT dynamics) |
| \(\kappa^*\) | κ*_ribo (from ΔG‡) | κ*_NGS (from L_p, d_primer) | κ*_MT (from tubulin mechanics) |
| Contact form | α = dz − r²dθ | α = dz − r²dθ | α = dz − r²dθ |
| Morphism | Identity (base) | Φ_NGS→Ribo | Φ_MT→Ribo |
All figures are generated by dm3_simulation.py using the DOP853 integrator (Hairer–Nørsett–Wanner)
with rtol=1e-10, atol=1e-12. The simulation integrates 66 orbits from the dm³ toy model equations,
verifying:
• \(\mu_{\max} = -2\) confirmed numerically (transverse Lyapunov exponent)
• Inner basin boundary \(r^* \approx 0.776\) (Gronwall asymmetry correction, AXLE Issue #13)
• Stability hierarchy: \(\varepsilon_0 = 1/3 < 2/3 < r^* \approx 0.776 < \kappa^* \approx 0.882 < 1\)
• Spiral return confirmed: \(x_0 \to G^{64}(x_0) \to G^{64}(x_{64}) = x_0'\) with \(x_0' \neq x_0\)
# dm3_simulation.py (excerpt)
import numpy as np
from scipy.integrate import ode
def dm3_rhs(t, y):
r, theta, z = y
dr = r*(1 - r**2) + 2*(r-1)*np.exp(-z)
dtheta = 1.0
dz = r**2 - 2*(r-1)**2*np.exp(-z)
return [dr, dtheta, dz]
# DOP853 integrator (Hairer-Norsett-Wanner)
solver = ode(dm3_rhs).set_integrator('dop853',
rtol=1e-10, atol=1e-12, nsteps=50000)
| File | Description |
|---|---|
togt_domain1_riboswitch.pdf | Full paper — Domain 1 (350 kB) |
togt_domain2_ngs.pdf | Full paper — Domain 2 (331 kB) |
togt_domain2_ngs_impa.pdf | IMPA/SBM submission version — Domain 2 (209 kB) |
togt_domain3_mt.pdf | Full paper — Domain 3 (349 kB) |
The TOGT/GTCT framework is explicitly distinguished from quantum-mechanical models of consciousness (Orch-OR, Penrose–Hameroff): TOGT operates at the level of classical contact geometry. The microtubule curvature predictions (Domain 3) involve macroscopic curvature fields measurable by electron microscopy — no quantum coherence claims are made. This distinction is important for the domain of applicability of the framework and for the experimental falsifiability of the predictions.
Two domain axioms remain explicitly open (declared with trivial, not hidden sorry):
Axiom A.NGS: The mTORC1-analogous kinase governing the Whitney fold in the NGS domain satisfies Mather's transversality condition at the cluster bifurcation point. (Blocker: Mather's theorem requires specific smoothness conditions on the primer-DNA interaction potential.)
Axiom A.NGS.2: The NGS bridge amplification flow satisfies the Poincaré–Bendixson conditions for limit cycle existence in the effective 2D system. (Blocker: requires Mathlib semiflow definition and Poincaré–Bendixson theorem.)
TOGT/GTCT — Domain 3 of 3
We apply the TOGT/GTCT contact-geometric framework to model microtubule dynamic instability as an operator sequence \(G = U \circ F \circ K \circ C\) on a dm³ contact manifold.
We derive three falsifiable classical predictions:
(P1) Anesthetic binding produces a measurable downward shift \(\Delta\kappa_{\text{MT}}\) linear in the minimum alveolar concentration (MAC).
(P2) Epothilone B LORR (Loss of Righting Reflex) delay is proportional to κ elevation with proportionality constant derivable from the dm³ Lipschitz bound.
(P3) The dose-response functional form is identical (up to contact morphism scaling) to the riboswitch OFF-lock curve (Domain 1).
All predictions are classically falsifiable without quantum measurement. Extends to tauopathies and TBI. Part of the Principia Orthogona series. XII Bienal da SBM, Natal, 2026.
Microtubules (MTs) are dynamic cytoskeletal polymers composed of α/β-tubulin dimers assembled into 13-protofilament hollow cylinders. Their dynamic instability — stochastic switching between growth (polymerization) and shrinkage (catastrophe) — is critical for mitotic spindle formation, axonal transport, and neuronal connectivity.
MT dynamics are sensitive to small-molecule perturbations at multiple levels:
Volatile anesthetics (isoflurane, sevoflurane, halothane) have been proposed to modulate consciousness through effects on MT curvature dynamics, independently of known membrane and ion-channel targets.
Neuroprotective agents (epothilone B, paclitaxel/taxol) stabilize MTs against catastrophe, with therapeutic applications in TBI, neurodegeneration, and axonal regeneration.
Tauopathies (Alzheimer's disease, frontotemporal dementia, CTE) involve hyperphosphorylated tau that detaches from MTs, leading to catastrophe and axonal disintegration.
The Penrose–Hameroff Orchestrated Objective Reduction (Orch-OR) hypothesis proposes that quantum coherence in MT tubulin dimers underlies consciousness. The TOGT/GTCT framework takes an explicit classical position: MT curvature dynamics sufficient to account for anesthetic modulation can be described within classical contact geometry without invoking quantum coherence.
This is not a claim that quantum effects are absent — it is a claim that the classical mechanism is sufficient for the falsifiable predictions in this paper, and that these predictions are experimentally distinguishable from Orch-OR predictions.
(1) The molecular target of volatile anesthetics on MTs is not established — proposed targets include the hydrophobic pocket on β-tubulin, the GTP-binding site, and lateral contacts between protofilaments.
(2) The quantitative relationship between anesthetic concentration and MT curvature changes has not been measured directly at single-filament resolution.
(3) The mechanism by which tau hyperphosphorylation disrupts the κ threshold is not understood geometrically.
| Operator | MT Biological Realisation |
|---|---|
| C | GTP-tubulin addition — compression of the protofilament tip; lattice strain accumulates |
| K | Curvature intensification at the MT tip — flared protofilament ends; GTP hydrolysis front approaches |
| F | Whitney A₁ fold at κ*_MT — catastrophe point; transition from growth to shrinkage |
| U | Rescue and stabilisation — GTP cap replenishment, return to growth phase (±anesthetic or epothilone B modulation) |
The curvature threshold \(\kappa^*_{\text{MT}}\) is identified with the catastrophe threshold:
where \(R_{\text{critical}}\) is the critical radius of MT-tip curvature at which the catastrophe transition occurs. This value corresponds to a protofilament tip angle of approximately 25° — consistent with cryo-EM measurements of curled protofilaments at catastrophe ends.
Anesthetic binding produces a measurable downward shift in the MT curvature threshold, linear in the minimum alveolar concentration (MAC):
where \(\lambda_{\text{anes}}\) is the dm³ Lipschitz constant for the anesthetic-MT interaction, derivable from the dm³ framework as:
where \(E_{\text{MT}}\) is the MT bending modulus (≈ 26 pN·μm²).
Experimental test: Single-filament MT dynamics assay (TIRF or cryo-EM) in the presence of isoflurane or sevoflurane at 0.5, 1.0, 1.5, and 2.0 MAC equivalents. Measure catastrophe frequency vs. anesthetic concentration. The slope should equal \(\lambda_{\text{anes}}\).
Falsification criterion: A nonlinear relationship (significant quadratic term, p < 0.05) or a slope inconsistent with the predicted \(\lambda_{\text{anes}}\) (beyond ±40%) falsifies P1.
Distinguishing from Orch-OR: The Orch-OR prediction requires quantum coherence times correlated with anesthetic potency. P1 only requires classical curvature measurements — no quantum measurement is needed to falsify it.
Epothilone B (EpoB) raises the MT curvature threshold \(\kappa^*_{\text{MT}}\), delaying the Loss of Righting Reflex (LORR) under volatile anesthetic challenge. The LORR delay time is proportional to \(\Delta\kappa^*_{\text{EpoB}}\) with proportionality constant from the dm³ Lipschitz bound:
where \(T_{\text{MT}} \approx 2\pi / \omega_{\text{MT}}\) is the characteristic MT oscillation period (on the order of minutes for in vivo dynamics).
Experimental test: Rodent LORR assay (Drosophila or mouse) with pre-administration of EpoB at doses below the anti-mitotic threshold (to avoid confounds from MT stabilization in dividing cells). Compare LORR onset time under standard isoflurane dose with and without EpoB. The delay should scale linearly with EpoB dose.
Falsification criterion: No statistically significant (p < 0.05) EpoB dose-dependent LORR delay at sub-anti-mitotic doses falsifies P2.
The dose-response functional form for MT catastrophe probability vs. anesthetic concentration is identical (up to contact morphism scaling) to the riboswitch OFF-lock dose-response curve (Domain 1):
This follows directly from the Coherence Bridge Theorem (Theorem 5.4): both systems have \(\mu_{\max} = -2\) and share the dm³ Whitney A₁ singularity structure.
In particular, the Hill coefficient for the MT dose-response satisfies the same dm³ relation:
Falsification criterion: A fitted Hill coefficient outside the range [1.5, 6.0] for the MT catastrophe dose-response, or a statistically significant mismatch (χ² test, p < 0.01) with the morphism-scaled riboswitch curve, would falsify P3.
In Alzheimer's disease and other tauopathies, hyperphosphorylated tau detaches from MTs, reducing the effective curvature threshold. Within the dm³ framework, tau hyperphosphorylation shifts \(\kappa^*_{\text{MT}}\) downward:
where \(\Delta\kappa_{\text{phos}}\) is the per-phosphorylation-site reduction in the curvature threshold, and \([\text{tau-P}]\) is the local concentration of hyperphosphorylated tau. When \(\kappa^*_{\text{tau-P}} < \kappa^*_{\text{thermal}}\) (the thermal noise threshold), spontaneous catastrophe occurs — modelling the progressive axonal degeneration of tauopathy.
TBI produces acute mechanical distortion of axonal MTs. Within the dm³ framework, this corresponds to a sudden large perturbation of the trajectory in phase space, pushing the system outside the Gronwall basin \(\varepsilon_0 = 1/3\). If the perturbation is outside the outer basin but within \(r^* \approx 0.776\), the system returns to the limit cycle (recoverable injury). If outside \(r^*\), the trajectory diverges (severe injury).
This predicts a sharp threshold for TBI recovery based on the magnitude of the mechanical perturbation — a prediction testable in in vitro MT stretching experiments.
Source: fig_domain3_mt.pdf (Zenodo 10.5281/zenodo.20559510) · Generated by dm3_simulation.py (DOP853, rtol=1e-10)
Core dm³ theorems are verified in Lean 4 via the AXLE repository (0 sorry in the core chain). Domain-specific open obligations:
| Obligation | Blocker | Status |
|---|---|---|
| MT fold from tubulin mechanics | Mather's theorem applied to MT bending potential | Open (trivial) |
| MT limit cycle existence | Poincaré–Bendixson for MT dynamic instability 2D projection | Open (trivial) |
| Tau-κ* coupling formula | Requires mTOR-analogous signalling cascade analysis | Planned |
The contact-geometric framework makes three qualitatively distinct predictions for microtubule dynamics under anesthetic and stabilizing-agent perturbations. The key feature is that all three predictions derive from a single geometric parameter (\(\kappa^*_{\text{MT}}\)) and two framework invariants (\(\mu_{\max} = -2\) and \(\varepsilon_0 = 1/3\)).
The EpoB LORR prediction (P2) is particularly important because it bridges in vitro MT biochemistry to a well-established in vivo phenotype. If confirmed, it would provide strong evidence that the classical curvature mechanism (not quantum coherence) accounts for MT-mediated anesthetic modulation of consciousness.
The extension to tauopathies points toward a unified geometric account of MT-related neurodegeneration that encompasses both chronic (tauopathy) and acute (TBI) pathological modes. The framework predicts a sharp boundary — not a gradient — between recoverable and irreversible MT network disruption, which could have implications for TBI treatment windows.
We have applied the TOGT/GTCT contact-geometric operator framework to microtubule curvature dynamics and derived three classically falsifiable predictions: (P1) anesthetic-induced κ shift linear in MAC, (P2) EpoB LORR delay with dm³ Lipschitz proportionality, and (P3) cross-domain dose-response morphism with the riboswitch domain. The framework is operationally distinguishable from Orch-OR — no quantum measurement is required for any of the three tests. Extensions to tauopathies and TBI are outlined. Lean 4 verification is complete for the core chain; two domain-specific axioms remain open.