🌍 The Dynamic Conformation Adsorption Model 💾

This repository contains the core predictive framework and Python model simulation for the Dynamic Conformation Adsorption Model (DCAM). DCAM resolves traditional linear decay contradictions in soil bioinformatics by modeling extracellular DNA (exDNA) adsorption kinetics via a first-principles physical chemistry approach.

📋 Executive Summary

Traditional soil bioinformatics models often present a structural paradox: they treat the accumulation and decay of extracellular (naked) DNA on mineral surfaces as a static linear process. This creates a mathematical contradiction where genetic information is predicted to decay uniformly, failing to account for the planet-wide genetic homogeneity observed in highly isolated, arid soil micro-environments.

DCAM resolves this paradox by replacing arbitrary decay constants with a first-principles physical chemistry framework. By coupling ambient ionic strength ($I$) to the geometric footprint ($\sigma$) of the DNA polymer, DCAM demonstrates that drying soils actively compress molecular data and throttle chemical degradation via steric hindrance, transforming dry mineral matrices into protective informational vaults.

🔺 The Master Kinetic Equation

The net change in active, readable genetic information bound to clay mineral surfaces over time ($t$) is governed by the following core differential equation:

$$\frac{dD_{surf}}{dt} = \alpha I_{in} \left(\frac{C_{ion}}{K_{bridge} + C_{ion}}\right) \theta_{clay} - \left[ \gamma_{chem} + \left( \frac{V_{max} \cdot B_{live}}{K_m + D_{surf}} \right) \right] D_{surf}$$

Variable and Parameter Definitions

• $D_{surf}$: Density of active, bound genetic information ($\text{base pairs / }\mu\text{m}^2$).
• $I_{in}$: Influx rate of extracellular DNA via cellular lysis ($\text{base pairs / day}$).
• $C_{ion}$: Concentration of divalent cations ($Ca^{2+}, Mg^{2+}$) driving ion-bridging.
• $\alpha$: Scaling coefficient representing DNA-clay surface affinity ($\left[ \frac{1}{\mu\text{m}^2} \right]$).
• $B_{live}$: Biomass density of active, surface-surfing microbes ($\text{cells / }\mu\text{m}^2$).

⚙️ Core Mechanistic Upgrades (The New Logic)

1. The Dynamic Space Governor ($\theta_{clay}$)

Rather than treating clay adsorption sites like a static grid, DCAM models the remaining capacity ($\theta_{clay}$) as an environmental function of molecular shape:

$$\theta_{clay} = 1 - \frac{D_{surf} \cdot \sigma(I)}{M_{max}}$$

Where $M_{max}$ is the maximum electrical saturation threshold of the clay matrix, and $\sigma(I)$ is the dynamic footprint area of the DNA molecule governed by the ionic strength ($I$) of the evaporating soil solution:

$$\sigma(I) = \sigma_{min} + (\sigma_{max} - \sigma_{min})e^{-\kappa I}$$

As soil moisture decreases, salts concentrate, causing $I$ to spike. This forces the DNA polymer to coil tightly into a compressed conformation, driving $\sigma \to \sigma_{min}$. Consequently, $\theta_{clay}$ increases, automatically opening up physical storage space on the mineral surface during droughts.

2. Chemically Derived, Shielded Decay ($\gamma_{chem}$)

To eliminate arbitrary information loss, the chemical degradation rate constant ($\gamma_{chem}$) is derived directly from collision theory and the Arrhenius equation, scaled by a Probability of Exposure Factor ($P_{exposed}$):

$$\gamma_{chem} = A \cdot e^{-\frac{E_a}{R \cdot T}} \cdot \left( \frac{\sigma}{\sigma_{max}} \right)^\beta$$

• $A$: Pre-exponential frequency factor ($\text{day}^{-1}$).
• $E_a$: Activation energy for phosphodiester bond hydrolysis ($\sim 100,000 \text{ J/mol}$).
• $R, T$: Universal gas constant ($8.314 \text{ J/mol}\cdot\text{K}$) and Temperature ($298.15 \text{ K}$).
• $\beta$: Steric hindrance exponent ($\beta = 3.0$).

When the DNA molecule coils tightly ($\sigma \to \sigma_{min}$), the inner core of the macromolecule is physically protected from chemical collisions. The term $\left(\frac{\sigma}{\sigma_{max}}\right)^\beta$ collapses toward zero, effectively freezing chemical decay during dry periods.

💻 Validated Simulation Results

The model was verified using a 44-day weather cycle: a 30-day progressive drying phase followed by an immediate flash flood on Day 31.

Day	Soil State	Ionic ($I$)	Footprint ($\sigma$)	Memory Stored ($D$)	Notes
1	Drying	0.13	1.1397	1508.09
10	Drying	0.40	0.4436	1678.70
20	Drying	0.70	0.2544	1910.85
30	Drying	1.00	0.2121	2147.96	Peak Vault Capacity
31	Flooded	0.05	1.6018	2114.23	Rapid Release
35	Flooded	0.05	1.6018	1988.26
44	Flooded	0.05	1.6018	1751.00	Biological Drain

Analysis of the Output

• Days 1–30 (Drying Vault): As the soil dries out, the molecular footprint ($\sigma$) shrinks drastically from 1.14 down to 0.21. Because the DNA folds tightly, it shields itself from hydrolysis. Stored memory steadily and safely climbs from 1508.09 to a peak of 2147.96 base pairs, functioning as an archival safe.
• Day 31 (The Flood Flush): When moisture returns, the sudden dilution causes the ionic strength to drop to 0.05, forcing the DNA to rapidly uncoil ($\sigma = 1.60$). The molecule loses its steric protection, and the waking microbial population ($B_{live}$) immediately begins downloading and consuming the exposed genetic archive, draining the surface storage down to 1751.00 by Day 44.

🎯 Definitive Empirical Predictions

DCAM offers three highly distinct, testable hypotheses for field researchers:

1. The Flash-Drought Spike: Environments subjected to rapid desiccation will exhibit higher concentrations of intact extracellular DNA than environments subjected to gradual cooling, due to the rapid locking of the conformational vault.
2. Mineral Sinks: Fine-grained, high-surface-area clays (e.g., Smectite, Bentonite) will behave as non-linear information buffers, retaining ancestral genetic sequences across climate eras, whereas coarse sandy soils will exhibit low informational thresholds and rapid data wash-out.
3. Post-Rain Transformation Burst: Horizontal Gene Transfer (HGT) frequencies among soil bacteria will peak in a tight non-linear window immediately following a rewetting event, driven by the structural expansion and sudden exposure of the stored DNA vault.

⚠️ Key Assumptions & Model Limitations

To maintain a parsimonious and mathematically traceable framework, DCAM operates under the following boundary conditions:

1. Divalent Cation Homogeneity: The ion-bridging mechanism assumes a soil solution dominated by divalent cations ($Ca^{2+}, Mg^{2+}$). It does not currently account for high-salinity competitive adsorption from monovalent ions ($Na^{+}$), which may destabilize the bridge under extreme sodic conditions.
2. Static Surface Area ($M_{max}$): The total electrical surface area capacity is treated as a fixed constant based on soil mineralogy class. Real-world 2:1 expanding lattice clays (e.g., Smectite/Bentonite) exhibit physical swelling and shrinking cycles that may dynamically alter available binding coordinates.
3. Homogeneous Micro-environments: The simulation assumes idealized, uniform contact between the DNA polymer and flat aluminum-silicate mineral plates, omitting macroscopic soil structural aggregates or organic matter competitive binding.

📢 Future Research & Open Call for Contributions

DCAM provides a baseline predictive architecture. We welcome extensions and contributions from soil scientists, biophysicists, and computational ecologists to address the following vectors:

• Sodic Soil Mechanics: Introducing a competitive ion-exchange coefficient to account for monovalent vs. divalent ratios in high-salt environments.
• Dynamic Mineralogy: Coupling $M_{max}$ to a matrix hydration variable to simulate the physical swelling and shrinking of expansive clays.
• Empirical Field Validation: Integrating field data from environmental DNA (eDNA) profiles—specifically tracking historical relic DNA vs. active microbial biomass fluctuations across stark wetting/drying cycles.

📄 Citation

Reed, Jonathan ƒ(n). (2026). The Soil Memory Overwrite Paradox: Modeling Extracellular DNA Adsorption Kinetics via the Dynamic Conformation Adsorption Model (1.0). Zenodo. https://doi.org/10.5281/zenodo.20289730

---

⚖️ Licensed under CC BY 4.0 International