R558C
Category 3/4 — Most DruggablePathogenic/Likely pathogenicCytoplasmic · predictedSource cardReviewed by an affected-family domain expert
TWO corrections. (1) FREQUENCY: the card stated R558C is carried by 'approximately 1 in 3' Ashkenazi Jewish individuals — this is wrong by an order of magnitude. Verified published carrier frequency is ~2.32% (~1 in 43) in Ashkenazi Jews (1.32% mixed Ashkenazi/Sephardi). The card's 'estimates range 27-34%' justification has no published basis and is superseded. (2) PATHOGENICITY: the single 'Pathogenic/Likely pathogenic' label is incomplete as clinical interpretation — published cohort data shows R558C homozygotes present a substantially milder, later-onset form (mean diabetes onset ~17-19 yrs, optic atrophy ~29 yrs, Wolfram diagnosis ~30 yrs; R558C protein more stable than other recessive pathogenic variants). Verified published data outranks the predictive label.
Nufar Kotler · 2026-08-11 · PMID 30014265 (Diabetologia 2018); PMID 36134655 (JCI Insight 2022); Neurol Genet 2021 doi:10.1212/NXG.0000000000000578
The reviewer's cited figure of four homozygotes is CONFIRMED EXACTLY: gnomAD v4 reports 4 R558C homozygotes, all Ashkenazi Jewish. The Ashkenazi allele frequency is 1.327% (393 of 29,608 alleles) against a global allele frequency of 0.035% — a 38-fold enrichment, independently confirming that the global number was never the relevant one for the at-risk population. A carrier frequency derived from that allele frequency (~1 in 38) brackets the published ~1 in 43 (PMID 30014265); the published figure remains the one cited on the card and the derived number is shown only as a cross-check. The adult homozygotes are consistent with the milder, later-onset course documented in the cohort literature recorded in the v1 entry above.
Nufar Kotler · 2026-08-21 · gnomAD v4 (retrieved 2026-08-21); PMID 30014265 (Diabetologia 2018)
Interactive 3D Structure
Bond changes · DynaMut2 interaction analysis
| Interaction type | Wild-type partner | Mutant partner | Status |
|---|---|---|---|
| Ionic bond | E431 | — | Lost |
| Hydrogen bond | E431 | — | Lost |
| Hydrogen bond | L554 | L554 | Preserved |
| Hydrogen bond | G555 | G555 | Preserved |
| Hydrogen bond | G562 | G562 | Preserved |
| Polar contact | E431 | — | Lost |
| Polar contact | L554 | L554 | Preserved |
| Polar contact | G555 | G555 | Preserved |
| Polar contact | G562 | G562 | Preserved |
| Van der Waals | — | L556 | Gained |
Lost / gained / preserved interatomic contacts at the variant residue, from the DynaMut2 (Arpeggio) interaction analysis of the wild-type and energy-minimized mutant structures.
Computational Predictions
This variant is in a water-facing region, which is the regime the stability predictor was built and benchmarked for. Ordinary predictor error still applies.
Clinical Evidence
Observed in the general population.
ClinVar classifies this variant as Pathogenic/Likely pathogenic for WFS1-related spectrum (unresolved mode) (dominant or recessive); Hearing loss, inheritance unstated (inheritance not specified). ClinVar records this under a WFS1-spectrum label that spans both disease families — recessive Wolfram syndrome (biallelic) and the dominant WFS1-related disorders. The label does not resolve which applies, and it is not a severity statement.
- WFS1-related spectrum (unresolved mode)dominant or recessivesubmitted as: Autosomal dominant and autosomal recessive WFS1-related disorders; WFS1-Related Spectrum Disorders
- Hearing loss, inheritance unstatedinheritance not specifiedsubmitted as: Monogenic hearing loss
- Phenotype scope is interpretation of submitted records. The computed values on this card (AlphaMissense, DynaMut2 ΔΔG, pLDDT) are measurements of the model, not of a patient.
- Presentation cannot be generalised from neighbouring variants: one nucleotide over can present very differently. Structural similarity supports a mechanistic hypothesis, never a phenotype claim.
Review status: 2★ multiple submitters, no conflicts. ClinVar "conditions" is the union of all submissions; 0–1★ entries are submitted conditions, not documented phenotypes.
Highest in Ashkenazi Jewish: AF 1.33% (393 of 29,608 alleles), 38.1x the global figure. The global AF describes the general population, not the at-risk group.
4 homozygotes reported in gnomAD v4 (4 Ashkenazi Jewish). gnomAD excludes individuals with severe pediatric disease, so homozygotes in this dataset argue against a fully penetrant severe recessive phenotype and belong in the clinical picture. This is interpretive signal, not a classification.
| Ancestry group | Allele freq. | AC / AN | Hom. | Carrier est. |
|---|---|---|---|---|
| Ashkenazi Jewish | 1.33% | 393 / 29,608 | 4 | ~1 in 38 |
| Remaining individuals | 0.091% | 57 / 62,474 | 0 | ~1 in 550 |
| African / African American | 0.013% | 10 / 74,950 | 0 | ~1 in 3750 |
| European (non-Finnish) | 0.0080% | 94 / 1,180,044 | 0 | ~1 in 6280 |
| South Asian | 0.0077% | 7 / 91,078 | 0 | ~1 in 6510 |
| Admixed American · under-sampled | 0.0017% | 1 / 60,010 | 0 | — |
Source: gnomAD r4, retrieved 2026-08-21. Values are allele frequencies. The carrier estimate is derived from the allele frequency (Hardy-Weinberg, ~2·AF) and is an approximation, not a published figure — where a published carrier frequency exists it outranks this column and is cited on the card.
Full Variant Card
WFS1 Wolframin — R558C Variant Card
Molecular Atlas Pilot Variant · RareResearch.AI · Windsor Symposium Demo
Prepared: May 26, 2026 · Expert-reviewed & corrected: August 11, 2026 (see Expert Review & Corrections below) · Source artifacts archived in this folder
Identity
| Field | Value |
|---|---|
| Variant | R558C |
| DNA change | c.1672C>T |
| Gene | WFS1 |
| Protein | Wolframin (890 aa) |
| UniProt ID | O76024 |
| dbSNP | rs199946797 |
| ClinVar accession | VCV000198835 |
| Amino acid change | Arginine (R, positive, hydrophilic) → Cysteine (C, neutral, thiol) at position 558 |
Structural Context
| Field | Value |
|---|---|
| AlphaFold model | AF-O76024-F1, v6 (released Aug 1, 2025) |
| pLDDT at residue 558 | 84.56 — high confidence, well-folded region |
| Domain assignment | Cytoplasmic loop between transmembrane helix 7 (TM7, residues 529–549) and TM8 (residues 563–583) |
| Structural neighborhood | 5 residues before the start of TM8; sits at the cytoplasmic edge of the membrane interface |
| IDR flag | No (pLDDT 84.56 ≫ 50 threshold). R558 is in an ordered, well-modeled region. |
Why this matters: Wolframin has eleven transmembrane helices anchoring it in the ER membrane. R558 sits in a short cytoplasmic loop directly preceding TM8. A positively charged arginine at this position likely contributes to membrane-proximal helix packing and electrostatic interactions with anionic phospholipid headgroups. Substituting cysteine removes the positive charge, introduces a free thiol capable of forming aberrant disulfide bonds, and is expected to perturb helix anchoring and ER folding fidelity.
Computational Predictions
AlphaMissense
| Field | Value |
|---|---|
| am_pathogenicity | 0.849 |
| am_class | LPath (Likely Pathogenic) |
| Threshold | > 0.564 = likely pathogenic; 0.849 is well into the pathogenic range |
| Source | DeepMind AlphaMissense (AlphaFold DB AA-substitutions table for O76024) |
DynaMut2
| Field | Value |
|---|---|
| Job ID | 177985627697 |
| ΔΔG (kcal/mol) | −0.5 kcal/mol (Destabilising) |
| Magnitude interpretation | Mild destabilization — fold largely intact, but specific local contacts lost |
| Interaction map (visual) | Wild-type R558 forms ionic + polar/hydrogen-bond contacts with neighboring TM7–TM8 residues (LEU554 vicinity). The R→C swap removes the guanidinium group, eliminating the ionic interactions and reorganizing the local hydrogen-bond network. |
| Stability classification | Destabilising — mild (|ΔΔG| < 2 + functional contacts lost → Category 3 / Category 4, Most Druggable) |
| URL | Job 177985627697 · retrieved 2026-05-27 — result self-hosted by RareResearch.AI (Biosig no longer serves this page) |
Clinical Evidence
Inheritance and scope
Pathogenic/Likely pathogenic — for WFS1-related spectrum (unresolved mode) (dominant or recessive); Hearing loss, inheritance unstated (inheritance not specified)
ClinVar classifies this variant as Pathogenic/Likely pathogenic for WFS1-related spectrum (unresolved mode) (dominant or recessive); Hearing loss, inheritance unstated (inheritance not specified). ClinVar records this under a WFS1-spectrum label that spans both disease families — recessive Wolfram syndrome (biallelic) and the dominant WFS1-related disorders. The label does not resolve which applies, and it is not a severity statement.
Phenotype scope is interpretation of submitted records. The computed values on this card (AlphaMissense, DynaMut2 ΔΔG, pLDDT) are measurements of the model, not of a patient. Presentation cannot be generalised from neighbouring variants: one nucleotide over can present very differently. Structural similarity supports a mechanistic hypothesis, never a phenotype claim.
| Field | Value |
|---|---|
| ClinVar classification | Pathogenic / Likely pathogenic |
| Review status | Criteria provided, multiple submitters, no conflicts |
| Last evaluated | February 1, 2026 |
| Associated conditions | Wolfram syndrome 1 · Wolfram-like syndrome · WFS1-Related Spectrum Disorders · Autosomal dominant nonsyndromic hearing loss 6 · Optic atrophy · Type 2 diabetes mellitus · Cataract 41 · Retinal dystrophy · Monogenic hearing loss |
| Inheritance | Both autosomal dominant and autosomal recessive forms documented |
R558C is one of 326 pathogenic-spectrum variants identified for WFS1 in ClinVar (out of 2,243 total variants catalogued).
Patient Population
R558C is the Ashkenazi Jewish founder mutation in WFS1, carried by approximately 2.32% of Ashkenazi Jewish individuals (~1 in 43) — 1.32% in mixed Ashkenazi/Sephardi cohorts, 0.04% Sephardi (PMID 30014265). The gnomAD global AF of 0.035% reflects the general population, not the at-risk group. This makes it the largest single identifiable carrier group for any WFS1 mutation worldwide and the highest-yield starting point for any clinical intervention targeting WFS1.
gnomAD v4 confirms the population-specific picture independently: the Ashkenazi Jewish allele frequency is 1.327% (393 of 29,608 alleles), against a global allele frequency of 0.035% — a 38-fold enrichment. A carrier frequency derived from that allele frequency (~1 in 38) brackets the published figure of ~1 in 43 (PMID 30014265); the published number remains the one to cite, and the derived one is shown only as an independent cross-check. Source: gnomAD v4, retrieved 2026-08-21.
Homozygotes in the general population: gnomAD v4 reports 4 R558C homozygotes, all Ashkenazi Jewish. gnomAD excludes individuals with severe pediatric disease, so the presence of adult homozygotes is consistent with the milder, later-onset course documented in the cohort literature below, and argues against a fully penetrant severe recessive phenotype. This is interpretive signal, not a reclassification. Source: gnomAD v4, retrieved 2026-08-21; flagged in expert review by Nufar Kotler, 2026-08-11.
Clinical spectrum: the ClinVar "Pathogenic/Likely pathogenic" label and AlphaMissense 0.849 are correct as annotations but incomplete as clinical interpretation. Published cohort data shows R558C homozygotes present a substantially milder, later-onset form of Wolfram syndrome — mean diabetes onset ~17–19 years, optic atrophy ~29 years, Wolfram diagnosis ~30 years — and the R558C protein is more stable than other recessive pathogenic variants (PMID 30014265, Diabetologia 2018; PMID 36134655, JCI Insight 2022; Neurol Genet 2021, doi:10.1212/NXG.0000000000000578). Pathogenicity here is a spectrum, not a single label.
Druggability Assessment
Final classification: Category 3 / Category 4 — Most Druggable.
The convergent evidence:
- pLDDT 84.56 — region is well-structured and modelable (NOT IDR; safe for docking)
- DynaMut2 ΔΔG = −0.5 kcal/mol — mild destabilization. The fold survives the mutation; this is NOT a gross misfolding mutation requiring gene therapy
- DynaMut2 interaction map — R558's guanidinium group makes specific ionic and polar contacts with neighboring residues in the TM7–TM8 region. The R→C swap eliminates these contacts, disrupting the local interaction network without globally unfolding the protein
- AlphaMissense 0.849 (LPath) — confirms the mutation has functional consequence, despite mild structural cost
- ClinVar Pathogenic/Likely pathogenic with no conflicts — clinical penetrance is real
Therapeutic implication: R558C is exactly the kind of variant the Atlas was built to identify — pathogenic but structurally tractable. The protein folds; the lost interactions define a specific local site that can be targeted by small molecules. Priority for docking experiments and pharmacological chaperone screening. Gene therapy is not required and likely not warranted given the largely-intact fold.
This is also the highest-value variant in the Atlas pilot because the carrier population (~1 in 43 Ashkenazi Jewish individuals, PMID 30014265 — the largest identifiable WFS1 carrier group worldwide) is large enough to support a real clinical pathway.
Research Path (decision tree)
ΔΔG < 2 + binding site affected → CATEGORY 3 — docking experiments
ΔΔG 2–4 → CATEGORY 2 — pharmacological chaperones
ΔΔG > 4 → CATEGORY 1 — gene therapy
pLDDT < 50 → CATEGORY 5 — IDR, experimental validation
Stable fold + functional site hit → CATEGORY 4 — site-specific docking
R558C pLDDT 84.56 rules out Category 5. Final assignment pending ΔΔG.
Artifacts in this folder
| File | Purpose |
|---|---|
AF-O76024-F1-model_v6.pdb | AlphaFold predicted structure (input for Mol* and DynaMut2) |
AF-O76024-F1-confidence_v6.json | Per-residue pLDDT scores |
AF-O76024-F1-aa-substitutions.csv | All 16,910 possible single AA substitutions with AlphaMissense scores |
R558C_molstar_viewer.html | Interactive 3D viewer — open in browser, screenshot for Windsor deck |
WFS1_clinvar_variants.csv | All 2,243 ClinVar WFS1 variants with classifications |
uniprot_O76024.json | Full UniProt feature annotations |
R558C_variant_card.md | This card (source of truth) |
R558C_variant_card.html | Demo-ready printable version |
Expert Review & Corrections
This card has been reviewed by an affected-family domain expert. Reviewer: Nufar Kotler (Jewish Wolfram Network), 2026-08-11. Corrections from that review are recorded in the Atlas corrections log (v1, 2026-08-11) and applied to this card:
-
Carrier frequency corrected (order-of-magnitude error). A previous version of this card stated R558C is carried by "approximately 1 in 3" Ashkenazi Jewish individuals, citing a "27–34%" range. That figure had no published basis. The verified published carrier frequency is ≈2.32% (~1 in 43) in Ashkenazi Jews (1.32% mixed Ashkenazi/Sephardi; 0.04% Sephardi). R558C remains the Ashkenazi founder variant and the largest identifiable WFS1 carrier group — the corrected number carries that argument on its own. Ref: PMID 30014265.
-
Pathogenicity presented as a spectrum. The single "Pathogenic/Likely pathogenic" label and AlphaMissense 0.849 stand as annotations but are incomplete as clinical interpretation: published cohort data shows R558C homozygotes present a substantially milder, later-onset form (mean diabetes onset ~17–19 yrs; optic atrophy ~29 yrs; Wolfram diagnosis ~30 yrs; R558C protein more stable than other recessive pathogenic variants). Verified published data outranks the single predictive label. Refs: PMID 30014265; PMID 36134655; Neurol Genet 2021, doi:10.1212/NXG.0000000000000578.
-
Population-specific frequency now shown on every card, not just this one (v2, 2026-08-21). The v1 fix corrected this card's Ashkenazi figure but left every other variant showing only a global gnomAD allele frequency, which is the same error in miniature: a global number standing in for the at-risk population. The Atlas now carries a per-ancestry allele-frequency breakdown on all 2,085 variants and names the highest-frequency population alongside the global figure. Ref: gnomAD v4, retrieved 2026-08-21 via
scripts/pull_gnomad.py. -
Homozygote counts added (v2, 2026-08-21). The reviewer noted that the presence of homozygotes in general-population databases belongs in the clinical picture, and cited four for R558C. The schema had no homozygote field at all — 0 of 2,085 records carried one. gnomAD v4 reports 4 R558C homozygotes, all Ashkenazi Jewish, confirming the reviewer's figure exactly. The field is now present on every variant, with an explicit "not available" where gnomAD has no entry rather than a silent zero. Ref: gnomAD v4, retrieved 2026-08-21.
-
Pathogenicity is now phenotype-scoped Atlas-wide (v2, 2026-08-21). Per the reviewer's second correction, a WFS1 variant may be pathogenic for diabetes or another WFS1-related phenotype without that saying anything about Wolfram syndrome or expected severity. Every card now states what the variant is classified for, under which inheritance mode, and at what ClinVar review strength — with 0–1★ single-submitter entries flagged as unverified submissions rather than documented phenotypes.
Expert review happens in the open — predictions get corrected in the artifact itself, not silently. The full corrections log is published on the Atlas site.
Modeling assumptions (documented per Atlas standard)
- AlphaFold predicts the monomer structure only. Wolframin's native oligomeric state in the ER membrane is not modeled here.
- DynaMut2 ΔΔG is computed on the static AlphaFold structure; dynamic ensemble effects are not captured.
- AlphaMissense is trained on missense variants observed in primates; novel substitutions in less-conserved regions may be less reliable.
- ClinVar classifications reflect aggregated clinical evidence as of February 2026; new submissions may update.
- Ashkenazi carrier frequency is cited at ≈2.32% (~1 in 43) per PMID 30014265 (Diabetologia 2018). An earlier version of this card stated "1 in 3" with a claimed "27–34%" range — that figure had no published basis and was corrected on expert review, 2026-08-11 (see Expert Review & Corrections).
Every assumption documented. Every score sourced. The Atlas standard.
Docking candidates · predicted, not measured
This screen was run against a binding pocket that has not been confirmed by a computational chemist. The site itself is unverified, so every compound below inherits that uncertainty.
- POCKET 33 WAS NOT REVIEWED BY A HUMAN. This run was forced. Every candidate below inherits that: the site itself is unconfirmed.
- DiffDock confirmation was NOT run, so the spec's 'DiffDock confidence > 0.5' secondary filter is INACTIVE - not passed. These are single-engine Vina scores with no orthogonal confirmation.
- CONFORMER SENSITIVITY, MEASURED ON THIS RUN: 698 molecules were present as both a parent and a salt form, giving an internal replicate (identical molecule, identical seed and box, different starting 3D conformer). Median score difference 0.034 kcal/mol, 95th percentile 0.425, worst case 1.485. Differences below roughly 0.425 kcal/mol are not meaningful.
- Vina affinity is a scoring-function ESTIMATE, not a measured binding constant. Its standard error (~2-3 kcal/mol) exceeds the gap between most adjacent ranks, so treat this as an unordered shortlist rather than a ranked order.
- Vina's raw affinity scales with molecular size, so larger ligands score better almost regardless of fit. Ligand efficiency (-affinity per heavy atom) is reported alongside and is the fairer comparison; a shortlist read on affinity alone will be biased toward the heaviest compounds in the library.
- The receptor is rigid: induced fit and sidechain rearrangement are not modelled.
- Docking was performed against the wild-type AlphaFold model. It does not account for any local remodelling the variant itself causes.
- This is a computational prediction intended to prioritise bench work. It is not evidence of binding.
| # | Compound | Vina kcal/mol | DiffDock | Notes |
|---|---|---|---|---|
| 1 | BICTEGRAVIR SODIUM | -7.42 | not run | FDA approved |
| 2 | CAPMATINIB | -7.42 | not run | FDA approved |
| 3 | DEXAMETHASONE SODIUM PHOSPHATE | -7.35 | not run | FDA approved |
| 4 | FLUOCINONIDE | -7.30 | not run | FDA approved |
| 5 | DOLUTEGRAVIR SODIUM | -7.20 | not run | FDA approved |
AutoDock Vina v1.2.7 · single-engine scores, no orthogonal confirmation · run 1478c9c9a483
Feed this card to Wolfram Intelligence
Download the R558C PDF below and upload it to Wolfram Intelligence to generate therapeutic-strategy proposals — guanidinium mimetics, sigma-1 agonist docking, NAC thiol-capping. NAC is already on the bench-testing list.