01
Introduction et Motivation Scientifique
Le paradigme quaternaire et ses limites
Depuis Watson et Crick (1953), l'ADN est décrit comme une séquence de quatre bases azotées — adénine (A), thymine (T), cytosine (C), guanine (G) — appariées selon des règles de complémentarité strictes. Ce paradigme quaternaire a conduit à des avancées extraordinaires : décryptage du code génétique, séquençage haut débit, édition génomique. Il constitue la fondation de toute la bioinformatique moderne.
Cependant, ce modèle traite chaque base comme un symbole discret équivalent. Il fait abstraction de la richesse physique des liaisons hydrogène qui stabilisent la double hélice. Ces liaisons ne sont pas des entités binaires on/off — elles sont dynamiques, caractérisées par des longueurs variables, des asymétries de densité électronique, et des couplages mutuels mesurables par spectroscopie et cristallographie.
Question centrale
Les liaisons hydrogène des paires de bases de l'ADN constituent-elles un espace d'états ternaire discret {−1, 0, +1} physiquement mesurable, contrôlable et exploitable comme support d'information indépendant de l'identité de la base ?
Qudits-36 propose de répondre à cette question par une formalisation mathématique rigoureuse. L'hypothèse n'est pas que l'ADN est un ordinateur quantique — c'est une proposition plus modeste et plus précise : les états de résonance des liaisons hydrogène constituent une couche informationnelle supplémentaire, encodable en logique ternaire équilibrée, dont la richesse dépasse de plusieurs ordres de grandeur le modèle quaternaire classique.
EVK
EVK
Espace Vectoriel du Codon
Modélise chaque codon comme vecteur dans la somme directe d'espaces ternaires T₃ⁿ, un par base du triplet.
ERK
ERK
Espace de Résonance du Codon
Formalise les états discrets {−1, 0, +1} des liaisons H comme configurations thermodynamiquement accessibles à 300K.
MRK
MRK
Matrice de Résonance du Codon
Encode les couplages inter-liaisons comme opérateur linéaire — matrice bloc symétrique semi-définie positive.
02
Fondements Physiques et Biologiques
Structure de la double hélice · Liaisons hydrogène Watson-Crick
La double hélice d'ADN est stabilisée par deux familles d'interactions non-covalentes : les interactions de stacking π-π entre bases empilées le long de l'axe hélicoïdal, et les liaisons hydrogène entre bases complémentaires sur un même plan perpendiculaire à cet axe. Ce sont ces dernières qui constituent le substrat de Qudits-36.
2.1 Géométrie des paires de bases Watson-Crick
Les deux familles de paires de bases se distinguent par le nombre de leurs liaisons hydrogène, déterminé par la géométrie de reconnaissance moléculaire :
2.2 Discrétisation ternaire — justification physique
La liaison hydrogène X─H···Y est un système dynamique dont la longueur r(H···Y) fluctue en fonction de la température, du pH, et des interactions avec des protéines avoisinantes. Les données cristallographiques de la PDB (résolution ≤ 1.5 Å) montrent des distributions approximativement gaussiennes centrées autour de 2.0 Å avec un écart-type σ ≈ 0.15–0.20 Å.
La discrétisation ternaire de Qudits-36 propose de regrouper ces états continus en trois classes discrètes selon leur écart à la longueur d'équilibre :
| État ternaire |
Longueur H···Y |
Signification physique |
Énergie ΔE |
−1 Comprimée |
~1.80 Å |
Compression — répulsion de Pauli + terme harmonique |
+0.800 kcal/mol |
0 Équilibre |
~2.00 Å |
État de repos — énergie de référence |
0.000 kcal/mol |
+1 Étirée |
~2.20 Å |
Extension — terme harmonique seul |
+0.400 kcal/mol |
Fondement énergétique · CHARMM36
Les énergies sont dérivées du potentiel harmonique du champ de forces CHARMM36 (Best et al., JCTC 2012) avec kHB = 20 kcal/mol/Ų et Δd = ±0.20 Å. La dissymétrie entre les états −1 et +1 (ΔΔE = 0.4 kcal/mol) est physiquement justifiée : la compression crée une répulsion de Pauli stérique absente en extension. Cette valeur reste à calibrer par calcul DFT explicite sur les paires AT et CG isolées. [FLAG_CALIBRATE]
03
EVK — Espace Vectoriel du Codon
Somme directe des espaces ternaires · T₃² ⊕ T₃³
3.1 Alphabet ternaire équilibré
L'alphabet ternaire équilibré T₃ = {−1, 0, +1} est muni d'une loi d'addition définie par :
Éq. (1)
a ⊕₃ b = ((a + b + 1) mod 3) − 1, ∀ a, b ∈ T₃
Addition ternaire équilibrée. Cette opération fait de T₃ un groupe abélien cyclique d'ordre 3.
La table de vérité complète confirme la clôture de T₃ sous ⊕₃ :
| ⊕₃ | −1 | 0 | +1 |
|---|
| −1 | +1 | −1 | 0 |
| 0 | −1 | 0 | +1 |
| +1 | 0 | +1 | −1 |
3.2 Espaces des paires de bases
À chaque paire de bases on associe un espace vectoriel ternaire dont la dimension est égale au nombre de ses liaisons hydrogène :
Éq. (2)
VAT = T₃² ⟹ |VAT| = 3² = 9 états vecteur : (h₁, h₂) ∈ T₃²
VCG = T₃³ ⟹ |VCG| = 3³ = 27 états vecteur : (h₁, h₂, h₃) ∈ T₃³
Espaces ternaires des paires AT et CG. Chaque coordonnée hᵢ ∈ {−1, 0, +1} représente l'état de la i-ème liaison hydrogène.
3.3 Espace du codon — Somme directe
L'espace d'un codon B₁B₂B₃ est la somme directe des espaces de ses trois bases. Cette construction est algébriquement rigoureuse — la somme directe de deux espaces vectoriels V ⊕ W est l'espace produit cartésien V × W muni des opérations composante par composante.
Éq. (3)
Vcodon = VB₁ ⊕ VB₂ ⊕ VB₃
dim(Vcodon) = n₁ + n₂ + n₃
|Vcodon| = |VB₁| × |VB₂| × |VB₃| = 3n₁+n₂+n₃
Vecteur codon : concaténation des vecteurs des trois bases. L'espace résultant est un produit de groupes cycliques Z₃.
| Composition | Espace | Dimension | États | ×ADN classique |
|---|
| AT+AT+AT | T₃² ⊕ T₃² ⊕ T₃² | ℝ⁶ | 729 | ×11 |
| AT+AT+CG | T₃² ⊕ T₃² ⊕ T₃³ | ℝ⁷ | 2 187 | ×34 |
| AT+CG+CG | T₃² ⊕ T₃³ ⊕ T₃³ | ℝ⁸ | 6 561 | ×103 |
| CG+CG+CG | T₃³ ⊕ T₃³ ⊕ T₃³ | ℝ⁹ | 19 683 | ×307 |
04
ERK — Espace de Résonance du Codon
États thermodynamiquement accessibles · CHARMM36 + MP2/DFT
4.1 États quantiques de repos — données MP2/DFT
L'assignation des états de repos de chaque liaison H repose sur l'asymétrie de transfert de proton δ, calculée en MP2/DFT pour les paires Watson-Crick en géométrie optimisée. Cette asymétrie mesure la localisation préférentielle du proton dans le double puits de potentiel de la liaison.
Règle d'assignation · Asymétrie δ
Soit δ = (affinité protonique donneur − affinité protonique accepteur) / (somme). Alors :
δ < −0.10 → état −1 (proton vers le donneur, liaison courte)
−0.10 ≤ δ ≤ +0.10 → état 0 (quasi-symétrique, équilibre)
δ > +0.10 → état +1 (proton vers l'accepteur, liaison longue)
| Liaison | Atomes | δ (MP2) | État | ρBCP (e/a₀³) | Référence |
|---|
| AT · H₁ | N6(A)─H···O4(T) | −0.08 | 0 | 0.031 | Hobza & Šponer 1999 |
| AT · H₂ | N1(A)···H─N3(T) | −0.15 | −1 | 0.038 | Hobza & Šponer 1999 |
| CG · H₁ | O6(G)···H─N4(C) | +0.12 | +1 | 0.042 | Šponer et al. 2018 |
| CG · H₂ | N1(G)─H···N3(C) | −0.20 | −1 | 0.045 | Šponer et al. 2018 |
| CG · H₃ | N2(G)─H···O2(C) | −0.06 | 0 | 0.028 | Šponer et al. 2018 |
4.2 Calcul de l'énergie de résonance
Pour un vecteur d'état v = (t₁, t₂, ..., tₙ) ∈ T₃ⁿ, l'énergie de résonance du codon est la somme des contributions individuelles de chaque liaison H :
Éq. (4)
EERK(v) = Σᵢ ΔE(tᵢ)
avec ΔE(−1) = ½ · kHB · (Δd)² + EPauli = +0.800 kcal/mol
ΔE( 0) = 0.000 kcal/mol (référence)
ΔE(+1) = ½ · kHB · (Δd)² = +0.400 kcal/mol
kHB = 20 kcal/mol/Ų [CHARMM36, Best et al. 2012]
Δd = 0.20 Å [déviation par état ternaire]
EPauli = 0.40 kcal/mol [FLAG_CALIBRATE — DFT requis]
Énergie de résonance du codon. La dissymétrie ΔΔE = 0.40 kcal/mol entre états −1 et +1 reflète la répulsion de Pauli lors de la compression de la liaison H.
4.3 Population de Boltzmann à 300K
Les probabilités d'occupation à l'équilibre thermique (300K, kBT = 0.596 kcal/mol) sont données par la distribution de Boltzmann :
Éq. (5)
P(t) = exp(−ΔE(t) / kBT) / Z
Z = Σs ∈ {−1,0,+1} exp(−ΔE(s) / kBT)
P(−1) ≈ 0.147 P(0) ≈ 0.564 P(+1) ≈ 0.288
Les trois états sont thermiquement accessibles à 300K. L'état d'équilibre (0) est le plus probable, mais les états −1 et +1 sont peuplés à hauteur de ~15% et ~29% respectivement.
05
MRK — Matrice de Résonance du Codon
Opérateur linéaire des couplages · Structure en blocs · Garanties mathématiques
La MRK est la contribution originale centrale de Qudits-36. Aucun outil bioinformatique existant ne calcule cet opérateur. Elle encode les dépendances entre les états de liaisons H adjacentes dans le codon — une information absente du modèle quaternaire classique.
5.1 Structure en blocs
Pour un codon de bases B₁B₂B₃ avec n₁, n₂, n₃ liaisons H respectivement, la matrice MRK de dimension n×n (n = n₁+n₂+n₃) est construite par blocs :
Éq. (6)
B₁ B₂ B₃
┌──────────┬──────────┬──────────┐
B₁ │ M₁₁ │ M₁₂ │ M₁₃ │
├──────────┼──────────┼──────────┤
B₂ │ M₂₁ │ M₂₂ │ M₂₃ │
├──────────┼──────────┼──────────┤
B₃ │ M₃₁ │ M₃₂ │ M₃₃ │
└──────────┴──────────┴──────────┘
Cyan : blocs diagonaux (intra-paire)
Violet : blocs hors-diagonaux (inter-bases)
Blocs diagonaux (intra-paire) :
Éq. (7)
Mii[a][b] = tₐ² · δab + κ(paire_i) · tₐtᵦ · (1−δab)
κAT = 0.30 [coopérativité intra-AT, Šponer et al. 2018]
κCG = 0.50 [coopérativité intra-CG, Šponer et al. 2018]
Couplage intra-paire. Le terme diagonal encode l'auto-énergie de chaque liaison H. Le terme hors-diagonal encode la coopérativité entre liaisons de la même paire.
Blocs hors-diagonaux (inter-bases) :
Éq. (8)
Mij[a][b] = λ · tᵢₐ · tⱼᵦ (i ≠ j)
Mji[b][a] = Mij[a][b] (symétrie)
λ = 0.12 [estimé depuis J-couplings NMR, Wijmenga & van Buuren 1998]
[FLAG_CALIBRATE — calibration NMR/DFT stacking requis]
Couplage inter-bases, médié par les interactions de stacking π-π. La symétrie est imposée explicitement pour garantir M = Mᵀ.
5.2 Garanties mathématiques
Une matrice de couplage physiquement valide doit satisfaire trois propriétés :
Théorème — Propriétés de la MRK
La matrice MRK construite selon les équations (7) et (8) satisfait :
P1. Symétrie : M = Mᵀ (les couplages sont réciproques)
P2. Semi-définie positivité : construite avec régularisation ε·I pour garantir toutes les valeurs propres ≥ 0 (l'énergie de couplage est non-négative)
P3. Rang variable : le rang peut atteindre n, selon la configuration des états. Pour CG+CG+CG à l'état maximal, rang = 6 sur 9.
Régularisation de Tikhonov :
Éq. (9)
Mreg = M + ε · In
ε = max(0, −λmin) + 10⁻⁸
où λmin = plus petite valeur propre de M avant régularisation
La régularisation annule les valeurs propres négatives résiduelles dues aux erreurs numériques, garantissant la semi-définie positivité sans modifier la structure des couplages significatifs.
5.3 Invariants et interprétation physique
| Invariant | Formule | Signification physique |
|---|
| Trace |
tr(M) = Σᵢ M[i][i] |
Énergie d'auto-couplage totale — somme des auto-énergies des liaisons H |
| Déterminant |
det(M) |
Indépendance des états. det = 0 si deux liaisons H ont des états parfaitement corrélés → redondance informationnelle |
| Frobenius |
‖M‖_F = √(Σᵢⱼ M[i][j]²) |
Intensité globale des couplages. Sensible aux configurations extrêmes. Utilisable comme signature de séquence. |
| Spectre |
λ₁ ≥ λ₂ ≥ ... ≥ λₙ |
Modes propres de couplage. Chaque valeur propre correspond à un mode de corrélation indépendant entre liaisons H. |
| Anisotropie |
λ_max / λ_min |
Ratio d'anisotropie informationnelle. Élevé : couplages très directionnels. Faible : couplages isotropes. |
| Signature |
(n⁺, n⁰, n⁻) |
Nombre de valeurs propres positives, nulles, négatives. Caractérise la géométrie de l'espace de couplage. |
Un résultat concret sur l'insuline humaine (19 premiers codons) : rang moyen = 5.16 sur 9 possible, anisotropie moyenne = 3.67, toutes les matrices sont symétriques et semi-définies positives. La norme de Frobenius moyenne est de 2.54. Ces valeurs constituent la première signature MRK d'une protéine humaine dans le cadre Qudits-36.
06
Modes d'Encodage
Référence quantique · Expérimental · Deux niveaux épistémiques
Une distinction fondamentale structure l'implémentation de Qudits-36 : les états des liaisons H peuvent être assignés théoriquement depuis la chimie quantique (mode de référence quantique) ou mesurés expérimentalement (mode expérimental). Ces deux modes ont des statuts épistémiques différents et ne doivent jamais être confondus dans un rapport scientifique.
| Critère | Mode référence quantique | Mode expérimental |
|---|
| Source des états |
Asymétrie δ, MP2/DFT, calculs ab initio publiés |
Mesure directe : spectroscopie FTIR, cristallographie X ≤ 1.0 Å |
| Statut |
Prédiction théorique |
Donnée expérimentale |
| Usage |
Simulation, exploration, comparaison, première approximation |
Validation de l'hypothèse, analyse de données réelles |
| Dépendance de la séquence |
Déterministe depuis l'identité de la base |
Indépendant — les états peuvent différer de la prédiction théorique |
| Validation croisée |
Comparable à ρBCP (AIM) |
Référence pour affiner les calculs théoriques |
Test falsifiable de Qudits-36
Soient deux protéines P₁ et P₂ de séquences ADN identiques, produites dans des conditions différant uniquement par les états des liaisons H (mesurés par FTIR). Si les matrices MRK correspondantes ont des rangs différents, et si les comportements biologiques de P₁ et P₂ diffèrent de façon statistiquement significative, l'hypothèse Qudits-36 est partiellement confirmée. Si les comportements sont identiques malgré des MRK différentes, l'hypothèse est réfutée.
07
Applications Théoriques
Cryptographie · Structure 3D · Fonctions biologiques
7.1 Cryptographie — MICA-OTP
Le schéma de chiffrement MICA-OTP est un masque jetable (One-Time Pad) opérant dans T₃ⁿ. Il hérite de la sécurité inconditionnelle au sens de Shannon (1949) : si la clé K est uniformément distribuée dans T₃ⁿ et utilisée une seule fois, aucune information sur le message M n'est extraite du chiffré C, quelle que soit la puissance de calcul de l'adversaire.
Éq. (10)
E(M, K) = M ⊕₃ K = (m₁⊕₃k₁, ..., mₙ⊕₃kₙ)
D(C, K) = C ⊖₃ K = C ⊕₃ (−K)
I(M ; C) = 0 ↔ H(M|C) = H(M) [Shannon 1949]
Chiffrement MICA-OTP. L'espace de clés 3ⁿ est 1.58× plus grand que l'espace binaire 2ⁿ pour la même longueur n.
7.2 Structure 3D — Hamiltonien Qudits
La MRK peut être intégrée comme terme supplémentaire dans un Hamiltonien de structure protéique. Pour chaque codon i avec matrice MRK M_i :
Éq. (11)
Etotal = ERamachandran + Eliaison + EMRK
EMRK(i) = α · (tr(M_i) − nᵢ) + β · ‖M_i‖_F
α = 0.10 kcal/mol [FLAG_CALIBRATE]
β = 0.05 kcal/mol [FLAG_CALIBRATE]
Hamiltonien Qudits v1.0. Le terme E_MRK est la première intégration de la matrice de résonance dans un champ de forces de structure protéique. Les constantes α et β seront calibrées sur des structures PDB connues.
La perturbation ERK sur les angles de Ramachandran traduit l'état dominant de résonance en déformation locale de la chaîne principale : état −1 (tendu) pousse vers la région hélice α (φ ≈ −65°, ψ ≈ −40°) ; état +1 (étiré) pousse vers le feuillet β (φ ≈ −120°, ψ ≈ +120°). Cette perturbation est de l'ordre de 5–15° selon la configuration, ce qui reste dans les tolérances de Ramachandran mais différencie statistiquement deux populations de structures.
7.3 Stockage haute densité
La capacité d'encodage théorique du système Qudits-36 en mode CG+CG+CG est :
Éq. (12)
Capacité Qudits-36 = log₂(19 683) ≈ 14.3 bits / codon
Capacité ADN classique ≈ log₂(64) = 6 bits / codon
Facteur d'enrichissement = 14.3 / 6 ≈ ×2.38 bits par codon
Comparaison de la capacité d'encodage théorique. Pour référence, Church et al. (Science 2012) ont démontré un stockage de 5.5 pétabits/mm³ dans l'ADN classique. Qudits-36 ouvre théoriquement un facteur ×2.38 supplémentaire.
08
Limites et Hypothèses Non Vérifiées
Honnêteté intellectuelle · Programme expérimental
La rigueur scientifique exige d'énoncer explicitement les limites du modèle. Les points suivants constituent des hypothèses de travail qui requièrent validation expérimentale.
| Hypothèse | Statut | Test proposé |
|---|
| 3 clusters discrets de longueurs H |
Non vérifiée |
Analyse statistique de structures PDB (résolution ≤ 1.5 Å) — test de trimodalité Hartigan + GMM |
| États thermodynamiquement stables |
Non vérifiée |
Simulations MD 100 ns en solution, analyse des distributions à l'équilibre |
| Couplages κ et λ dans la MRK |
Estimés |
DFT explicite sur paires AT/CG isolées et en contexte de stacking |
| Répulsion de Pauli E_Pauli = 0.40 kcal/mol |
Estimée |
Calcul MP2/aug-cc-pVTZ sur H₂O···H₂O et N─H···N à courte distance |
| MRK influence la structure 3D |
Non vérifiée |
Comparaison de structures PDB pour protéines de même séquence, conditions différentes |
| MRK corrèle avec la fonction biologique |
Non vérifiée |
Analyse statistique sur UniProt — corrélation signature MRK / annotation fonctionnelle GO |
Note épistémologique
Qudits-36 est une hypothèse théorique formalisée mathématiquement. Sa cohérence algébrique est établie. Sa validité physique est entièrement ouverte et constitue le programme de recherche central de kyriosMICA pour les cinq prochaines années. Chaque paramètre marqué [FLAG_CALIBRATE] dans le code source identifie précisément un point de calibration requis avant toute affirmation quantitative sur des données biologiques réelles.
09
Références
Bibliographie scientifique · Sources des paramètres physiques
[1]
Watson, J.D. & Crick, F.H.C.
A Structure for Deoxyribose Nucleic Acid.
Nature, 171, 737–738 (1953).
[2]
Shannon, C.E.
Communication Theory of Secrecy Systems.
Bell System Technical Journal, 28(4), 656–715 (1949).
MICA-OTP
[3]
Hobza, P. & Šponer, J.
Structure, Energetics, and Dynamics of the Nucleic Acid Base Pairs.
Chemical Reviews, 99(11), 3247–3276 (1999).
États QM · AT
[4]
Šponer, J., Bussi, G., Krepl, M., et al.
RNA Structural Dynamics As Captured by Molecular Simulations.
Chemical Reviews, 118(8), 4177–4338 (2018).
États QM · CG · κ_CG
[5]
Best, R.B., Zhu, X., Shim, J., et al.
Optimization of the Additive CHARMM All-Atom Protein Force Field.
Journal of Chemical Theory and Computation, 8(9), 3257–3273 (2012).
k_HB = 20 kcal/mol/Ų
[6]
Bader, R.F.W.
Atoms in Molecules: A Quantum Theory.
Oxford University Press (1990).
AIM · ρ_BCP
[7]
Koch, U. & Popelier, P.L.A.
Characterization of C-H-O Hydrogen Bonds on the Basis of the Charge Density.
Journal of Physical Chemistry A, 99(24), 9747–9754 (1995).
Critères AIM · ρ_BCP
[8]
Wijmenga, S.S. & van Buuren, B.N.M.
The Use of NMR Methods for Conformational Studies of Nucleic Acids.
Progress in NMR Spectroscopy, 32(4), 287–387 (1998).
λ inter-bases [FLAG_CALIBRATE]
[9]
Church, G.M., Gao, Y. & Kosuri, S.
Next-Generation Digital Information Storage in DNA.
Science, 337, 1628 (2012).
Stockage ADN
[10]
Adleman, L.M.
Molecular Computation of Solutions to Combinatorial Problems.
Science, 266(5187), 1021–1024 (1994).
[11]
Muthukrishnan, A. & Stroud, C.R.
Multivalued Logic Gates for Quantum Computation.
Physical Review A, 62(5), 052309 (2000).
Qudits quantiques
[12]
Berman, H.M., Westbrook, J., Feng, Z., et al.
The Protein Data Bank.
Nucleic Acids Research, 28(1), 235–242 (2000).
PDB · Validation cristallographique