RareResearch.AI
← Back to atlas

R558C

Category 3/4 — Most DruggablePathogenic/Likely pathogenicCytoplasmic · predictedSource card
ArginineCysteine at position 558 · Connecting loop · WFS1 (Wolframin)
Expert Review · corrections applied

Reviewed 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)

Full corrections log →

Interactive 3D Structure

Wild-type reference
Wild-type R558 — ionic bond to E431
Fullscreen ↗
DynaMut2 mutant · R558C
Mutant C558 — ionic bond to E431 lost (3 contacts lost)
Fullscreen ↗

Bond changes · DynaMut2 interaction analysis

3 lost1 gained6 preserved
Interaction typeWild-type partnerMutant partnerStatus
Ionic bondE431Lost
Hydrogen bondE431Lost
Hydrogen bondL554L554Preserved
Hydrogen bondG555G555Preserved
Hydrogen bondG562G562Preserved
Polar contactE431Lost
Polar contactL554L554Preserved
Polar contactG555G555Preserved
Polar contactG562G562Preserved
Van der WaalsL556Gained

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

DynaMut2 ΔΔG
-0.50kcal/mol
Destabilising — mild
AlphaMissense
0.849
LPath
AlphaFold pLDDT
85
model confidence
Schema
Cat 3/4
Category 3/4 — Most Druggable
ΔΔG confidence: Standard confidence · Cytoplasmic

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

ClinVar classificationPathogenic/Likely pathogenic
Review statuscriteria provided, multiple submitters, no conflicts
Associated conditionsAutosomal dominant and autosomal recessive WFS1-related disorders; WFS1-Related Spectrum Disorders; Monogenic hearing loss
Population frequency (gnomAD v4)Low frequency · AF 0.035%
cDNA changec.1672C>T
ClinVar accessionVCV000198835
Last evaluated2026/02/01 00:00

Observed in the general population.

Classified for2★ documented assertion

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.

Population frequency
Global (all ancestries)
AF 0.035% · 562 / 1,612,648 alleles
Homozygotes
4
Highest-frequency population
Ashkenazi Jewish · AF 1.33%

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 groupAllele freq.AC / ANHom.Carrier est.
Ashkenazi Jewish1.33%393 / 29,6084~1 in 38
Remaining individuals0.091%57 / 62,4740~1 in 550
African / African American0.013%10 / 74,9500~1 in 3750
European (non-Finnish)0.0080%94 / 1,180,0440~1 in 6280
South Asian0.0077%7 / 91,0780~1 in 6510
Admixed American · under-sampled0.0017%1 / 60,0100

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

FieldValue
VariantR558C
DNA changec.1672C>T
GeneWFS1
ProteinWolframin (890 aa)
UniProt IDO76024
dbSNPrs199946797
ClinVar accessionVCV000198835
Amino acid changeArginine (R, positive, hydrophilic) → Cysteine (C, neutral, thiol) at position 558

Structural Context

FieldValue
AlphaFold modelAF-O76024-F1, v6 (released Aug 1, 2025)
pLDDT at residue 55884.56 — high confidence, well-folded region
Domain assignmentCytoplasmic loop between transmembrane helix 7 (TM7, residues 529–549) and TM8 (residues 563–583)
Structural neighborhood5 residues before the start of TM8; sits at the cytoplasmic edge of the membrane interface
IDR flagNo (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

FieldValue
am_pathogenicity0.849
am_classLPath (Likely Pathogenic)
Threshold> 0.564 = likely pathogenic; 0.849 is well into the pathogenic range
SourceDeepMind AlphaMissense (AlphaFold DB AA-substitutions table for O76024)

DynaMut2

FieldValue
Job ID177985627697
ΔΔG (kcal/mol)−0.5 kcal/mol (Destabilising)
Magnitude interpretationMild 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 classificationDestabilising — mild (|ΔΔG| < 2 + functional contacts lost → Category 3 / Category 4, Most Druggable)
URLJob 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.

FieldValue
ClinVar classificationPathogenic / Likely pathogenic
Review statusCriteria provided, multiple submitters, no conflicts
Last evaluatedFebruary 1, 2026
Associated conditionsWolfram 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
InheritanceBoth 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

FilePurpose
AF-O76024-F1-model_v6.pdbAlphaFold predicted structure (input for Mol* and DynaMut2)
AF-O76024-F1-confidence_v6.jsonPer-residue pLDDT scores
AF-O76024-F1-aa-substitutions.csvAll 16,910 possible single AA substitutions with AlphaMissense scores
R558C_molstar_viewer.htmlInteractive 3D viewer — open in browser, screenshot for Windsor deck
WFS1_clinvar_variants.csvAll 2,243 ClinVar WFS1 variants with classifications
uniprot_O76024.jsonFull UniProt feature annotations
R558C_variant_card.mdThis card (source of truth)
R558C_variant_card.htmlDemo-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:

  1. 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.

  2. 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.

  3. 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.

  4. 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.

  5. 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)

  1. AlphaFold predicts the monomer structure only. Wolframin's native oligomeric state in the ER membrane is not modeled here.
  2. DynaMut2 ΔΔG is computed on the static AlphaFold structure; dynamic ensemble effects are not captured.
  3. AlphaMissense is trained on missense variants observed in primates; novel substitutions in less-conserved regions may be less reliable.
  4. ClinVar classifications reflect aggregated clinical evidence as of February 2026; new submissions may update.
  5. 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

pocket 33· 1,566 compounds screened
Pocket not reviewed

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.

  1. POCKET 33 WAS NOT REVIEWED BY A HUMAN. This run was forced. Every candidate below inherits that: the site itself is unconfirmed.
  2. 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.
  3. 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.
  4. 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.
  5. 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.
  6. The receptor is rigid: induced fit and sidechain rearrangement are not modelled.
  7. Docking was performed against the wild-type AlphaFold model. It does not account for any local remodelling the variant itself causes.
  8. This is a computational prediction intended to prioritise bench work. It is not evidence of binding.
#CompoundVina kcal/molDiffDockNotes
1BICTEGRAVIR SODIUM-7.42not runFDA approved
2CAPMATINIB-7.42not runFDA approved
3DEXAMETHASONE SODIUM PHOSPHATE-7.35not runFDA approved
4FLUOCINONIDE-7.30not runFDA approved
5DOLUTEGRAVIR SODIUM-7.20not runFDA approved
Full assay input document (PDF) →

AutoDock Vina v1.2.7 · single-engine scores, no orthogonal confirmation · run 1478c9c9a483

Therapeutic Strategy Handoff · prediction

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.

Download R558C PDF card ↓Strategies are AI-generated predictions, not validated therapeutics.