| Orientation | PCA frame canonicalization → barrel-axis / hydrophobic-belt SVD fit (orientor.py, barrel.py) | Rotation-invariant: 5× random-rotation extracellular-set Jaccard ≥ 0.95 (all 11 antigens = 1.00). Barrel quality: ring-CV (roundness), lipid−pore ΔKD gate ≥ 0.5, aromatic Trp/Tyr girdle fraction. | 3,5 |
| Orientation | Surface-antigen anchor-relative outward axis (orientor.py) | Exposed at RSA ≥ 0.20; distal-axis fraction ≥ 0.5 = outward face. | 2 |
| Membrane | Asymmetric Gram-negative outer-membrane bilayer model (membrane.py) | Hydrophobic core 20–35 Å; interfacial band 8 Å/leaflet; LPS headgroup shield 4 Å; periplasmic buffer 1.5 Å. | 1,4 |
| Membrane | Kyte–Doolittle hydropathy + side-chain radial facing (membrane.py) | KD ∈ [−4.5, +4.5]; lipid-facing radial dot ≥ +0.20, pore/water-facing ≤ −0.20. Side (EC/periplasmic) labelled only when combined confidence ≥ 70%. | 5 |
| Accessibility | Shrake–Rupley SASA, 1.4 Å probe (sasa.py) | RSA = SASA ÷ Tien-2013 maximum allowed ASA; antibody-accessible gate RSA ≥ 0.20. | 2,16 |
| B-cell mining | Per-residue linear propensity (mine_membrane_epitopes.py) | P = 0.40·hydrophilicity + 0.20·flexibility + 0.40·RSA, each ∈ [0,1]; hydrophilicity = (4.5 − KD)/9, flexibility = Karplus–Schulz normalised. | 5,7,2 |
| B-cell mining | Sliding window + accessibility gate (mine_membrane_epitopes.py) | 15-mer windows; β-barrel gate = antibody-accessible fraction ≥ 0.50 (LPS-shielded counts 0.5×); non-barrel surface gate = raw RSA ≥ 0.25; predicted Sec signal peptide excluded before mining. | 6 |
| Topology | Extracellular-side call (topology.py) | Loop-architecture: extracellular loop ≥ 8 residues, decisive median+long-loop margin ≥ 0.15; tiebreakers = positive-inside (Arg/Lys/His density) and N/C-terminus periplasmic prior. | 6 |
| T-cell | MHC-I / MHC-II prescreen | NetMHCpan-4.1 / NetMHCIIpan-4.0 percentile rank (lower = stronger); HLA-agnostic and not membrane-gated (peptides are processed intracellularly). | 8 |
| Conservation | Strain-ortholog conservancy (analysis/conservation.py) | Mean sequence identity % and conservancy % across strain orthologs; epitope retained on the mean-identity gate. | — |
| Construct | Adjuvant / helper / linker architecture (construct/) | β-defensin (TLR1/2 agonist) · HP91 (innate) · PADRE (universal CD4 helper) · EAAAK rigid / AAY / GPGPG / KK linkers; every epitope keeps an exact residue link to its parent antigen. | 9,11,12,13 |
| Validation | Physicochemistry + fold + innate docking + display (validation_map.py) | ProtParam-equivalent: instability index (stable < 40), GRAVY, aliphatic index, theoretical pI. Fold confidence: ESMFold/AlphaFold3 pLDDT & pTM. Innate docking: ClusPro TLR1/TLR2 (PDB 2Z81) center score + interface residues. Surface display: protruding-residue clustering. | 14,15,17 |
| Dynamics | Elastic-network flexibility — GNM / ANM (dynamics/enm.py) | Per-residue mean-square fluctuation (RMSF proxy) from inverse-eigenvalue-weighted modes; GNM cutoff 7.5 Å, ANM 15 Å; mode collectivity κ; GNM B-factor recovery (Pearson r) as a validation check. Structure-derived, not MD. | 18,19,20 |
| Dynamics | Receptor–vaccine interface (dynamics/interface.py) | Inter-body Cα contacts (8 Å), interface residues, buried SASA (Shrake–Rupley), and an ANM flexibility overlay (is the interface a rigid binding-competent scaffold). | 16,19 |
| Dynamics | TITO transfer-operator surrogate (dynamics/timescales.py) | Brownian-mode RELATIVE implied timescales (τ_i/τ_j = λ_j/λ_i, exact & calibration-free) + a parametric N-ns horizon. The absolute ns axis, macrostate transition matrix, and bound-state occupancy are NOT produced — they require MD + a trained operator (see the 200 ns MD→MSM scaffold). | 21,19 |