The Structural Carbide Anchor and the Hydride Spring: Resolving the FeMo-Cofactor Belt-Sulfur Paradox via Rigidity Matroids, Proton-Coupled Hemilability, and Hydride Reductive Elimination

## Abstract The catalytic mechanism of nitrogenase—the primary enzymatic system (alongside homologous vanadium and iron-only nitrogenases) capable of reducing ambient dinitrogen ($\mathrm{N}_2$) to ammonia ($\mathrm{NH}_3$) at room temperature and pressure—presents a longstanding biophysical paradox. High-resolution crystallographic structures (PDB 3U7Q, 4TKV/4TKU, 6FEA) reveal that the iron-molybdenum cofactor (FeMo-co, $[\mathrm{Mo}\text{-}7\mathrm{Fe}\text{-}9\mathrm{S}\text{-}\mathrm{C}]$) contains an interstitial central $\mu_6$-carbide ($\mathrm{C}^{4-}$) and undergoes reversible displacement of a belt sulfur ($\mathrm{S2B}$) into the secondary coordination sphere during ligand binding. However, conventional gas-phase broken-symmetry density functional theory (BS-DFT) calculations predict that uncoordinating a bare belt $\mu_2\text{-S}^{2-}$ sulfide requires an insurmountable activation barrier of $+38\text{ to }+45\,\mathrm{kcal/mol}$, yielding an apparent turnover rate of $k \sim 10^{-17}\,\mathrm{s}^{-1}$ (one turnover every 1.07 billion years)—a 19-order-of-magnitude failure compared to the physiological turnover frequency ($k_{\mathrm{cat}} \approx 1.0\text{--}2.0\,\mathrm{s}^{-1}$, macrogated by Fe-protein docking at $4\text{--}6\,\mathrm{s}^{-1}$). Concurrently, synthetic multimetallic inorganic mimics lacking the central carbide decompose or undergo irreversible rhomboidal collapse upon belt-ligand displacement. While Skubi and Holland (2018) qualitatively hypothesized that the interstitial carbide might prevent cluster collapse upon sulfur loss into an adjacent "holding chamber" and that hydrides assist displacement, this conceptual picture has lacked quantitative mechanical and energetic proof. Here, we formulate an integrated biophysical constraint analysis that resolves this paradox by rigorously decoupling structural cluster resilience from catalytic activation kinetics across five coupled physical pillars:1. **Infinitesimal Rigidity Matroid Invariant (Tier M1):** An unbraced $\mathrm{Fe}_6$ triangular prism possesses $V=6, E=9$, leaving $3V - 6 - E = 3$ internal shearing degrees of freedom under Cauchy-Laman rigidity theory. The central $\mu_6$-carbide contributes six radial struts, increasing edge rank to 14 and eliminating two rectangular face-shearing collapse modes ($\dim(\ker R) - 6 = 3 \to 1$), with the single remaining torsional breathing mode mechanically damped by the capping protein coordination sphere. This mathematically proves why the interstitial carbide acts as a structural backstop that arrests cluster collapse when $\mathrm{S2B}$ uncoordinates.2. **Proton-Coupled Hemilability (Tier M2):** In the Thorneley-Lowe catalytic accumulation, proton delivery selectively protonates the belt sulfur ($\mu_2\text{-S}^{2-} + \mathrm{H}^+ \to \mu_2\text{-SH}^-$), lengthening $\mathrm{Fe}\text{--}\mathrm{S}$ bonds from $2.20\,\text{\AA}$ to $2.58\,\text{\AA}$ and reducing the mechanical hinge strain from $\sim 40\,\mathrm{kcal/mol}$ to $+13.5\,\mathrm{kcal/mol}$.3. **The Hydride Reductive-Elimination Spring (Tier M2):** Spectroscopic ENDOR measurements demonstrate that the Janus $E_4(4\mathrm{H})$ intermediate stores two electrons and two protons as bridging hydrides $[\mathrm{Fe}\text{-}(\mu\text{-H})\text{-}\mathrm{Fe}]_2$. Obligate reductive elimination of $\mathrm{H}_2$ upon $\mathrm{N}_2$ binding provides standard-state exergonic driving potential ($\Delta G_{\mathrm{re}}^\circ = -20.5\,\mathrm{kcal/mol}$ in gas phase; $-16.2\,\mathrm{kcal/mol}$ in aqueous solution), offsetting the $+13.5\,\mathrm{kcal/mol}$ hinge penalty. While in vivo cellular and matrix coupling tunes the net steady-state equilibrium $E_4(4\mathrm{H}) + \mathrm{N}_2 \rightleftharpoons E_4(2\mathrm{N}_2\mathrm{H}) + \mathrm{H}_2$ near thermoneutrality ($\Delta G^\circ \approx 0 \pm 2\,\mathrm{kcal/mol}$, Hoffman & Seefeldt 2020) to prevent irreversible thermodynamic trapping, the exothermic $\mathrm{H}\text{--}\mathrm{H}$ bond formation synchronously compensates the hinge strain, satisfying the semi-empirical kinetic consistency boundary of $\Delta G^\ddagger \le 17.2\,\mathrm{kcal/mol}$ required for physiological turnover ($k_{\mathrm{obs}} \approx 1.5\,\mathrm{s}^{-1}$).4. **Active-Site Pocket Leash & Holding Chamber (Tier M1/M2):** As structurally captured in CO-inhibited Mo-nitrogenase (PDB 4TKV) and resolved $\sim 7\,\text{\AA}$ away in V-nitrogenase (PDB 6FEA), the secondary coordination sphere retains displaced sulfur within a specialized cage. The constricted local exit throat ($d_{\mathrm{throat}} \approx 2.30\,\text{\AA}$) sterically blocks both naked ($4.10\,\text{\AA}$) and hydrated ($6.80\,\text{\AA}$) hydrosulfide ($\mathrm{HS}^-$), while imposing a $+15.2\,\mathrm{kcal/mol}$ Born desolvation barrier (and $> +29.7\,\mathrm{kcal/mol}$ total escape barrier), preventing loss to bulk solvent and ensuring reversible re-ligation.5. **Electronic Delocalization and Non-Redox Mo Ballast (Tier M2):** Anderson-Hasegawa double exchange and $28.5\%$ carbide $2p\text{--}3d$ covalency buffer electron density, suppressing vacuum spin-localization traps to $\le 1.8\,\mathrm{kcal/mol}$. Molybdenum, assigned as formal $\mathrm{Mo(III)}$ spin-coupled to the iron core ($S=3/2$), does not undergo net redox cycling across $E_0 \to E_4$ ($S=1/2$), functioning alongside bidentate $(R)$-homocitrate as a non-redox torsional ballast that dampens the residual $N_{\mathrm{floppy}} = 1$ breathing mode to $\pm 6.12^\circ$. We formulate three testable experimental predictions using nuclear resonance vibrational spectroscopy (NRVS), site-directed secondary-sphere mutagenesis ($\alpha\text{-Arg359}$), and isotopic exchange to definitively differentiate this structural-kinetic model from competing static hypotheses.