# WFS1 Wolframin — S443R Variant Card

**Molecular Atlas Pilot · RareResearch.AI · Generated by wolfram-variant-card skill**

*Serine → Arginine at position 443 inside wolframin's fourth transmembrane helix (TM4). ClinVar Pathogenic/Likely pathogenic. AlphaMissense 0.999 (near-maximum pathogenicity score), DynaMut2 ΔΔG -0.31 kcal/mol (destabilising). A pathogenic variant whose mechanism is charge-into-membrane disruption.*

---

## Identity

| Field | Value |
|---|---|
| **Variant** | S443R (p.Serine443Arginine) |
| **DNA change** | c.1329C>G |
| **Gene · Protein** | WFS1 · Wolframin (890 aa) |
| **UniProt** | O76024 · WFS1_HUMAN |
| **ClinVar accession** | VCV001315675 |
| **Amino acid change** | Serine (S) → Arginine (R) — a small polar hydroxyl-bearing residue replaced by a large, positively-charged residue with a long alkyl chain and a guanidinium group (the strongest hydrogen-bond donor in the amino acid alphabet). |

---

## Structural Context

| Field | Value |
|---|---|
| **AlphaFold model** | AF-O76024-F1, v6 |
| **pLDDT at residue 443** | **88.19** — well-folded |
| **Domain** | TM4 (427-447), helical transmembrane |
| **Position context** | TM4 (residues 427–447) · position 443 is bilayer-embedded near the lumenal end of the helix. The lipid environment penalizes charged side chains heavily. |
| **IDR flag** | No — pLDDT above 50 threshold |

**UniProt features at this position:**

  - Chain: 1-890 Wolframin
  - Transmembrane: 427-447 Helical
  - Natural variant: 443-443 in WFS1

> Position 443 sits inside TM4, one of wolframin's eleven transmembrane helices. The AlphaFold model places S443 within 5 Å of THR442 (2.5 Å) and TYR444 (2.5 Å), and into a packed environment with PHE439 (3.5 Å), SER446 (3.6 Å), PHE365 (3.7 Å, from TM3), and THR440 (3.9 Å). The wild-type serine's small polar hydroxyl fits well in this membrane-embedded context, possibly forming a hydrogen bond with SER446 or with the backbone carbonyl of a nearby residue.

Replacing serine with arginine here introduces three layered structural costs. First, the volume difference: arginine is one of the larger amino acids by side-chain volume, while serine is among the smaller. The local packing has to accommodate roughly four-fold more side-chain mass. Second, the charge: arginine's guanidinium group carries a positive charge that is thermodynamically penalized in the bilayer hydrophobic core. Third, the H-bonding capacity changes — the lost serine hydroxyl is replaced by the guanidinium's strong H-bond donor character, which would prefer to engage water or polar partners outside the membrane.

DynaMut2's |ΔΔG| of 0.31 kcal/mol underestimates the structural disruption. The variant probably forces local rearrangement: the arginine side chain extends toward the membrane-water interface (where its charge can be partially satisfied), pulling nearby residues out of their wild-type positions. The PHE365 contact from TM3 (3.7 Å) is particularly important — that's a helix-helix interaction, and disrupting it perturbs the relative geometry of TM3 and TM4.

AlphaMissense's score of 0.999 reflects the severity of this mechanism even with a small structural ΔΔG. The variant is pathogenic because it disrupts helix-helix packing in the membrane, not because it unfolds the protein.

---

## Computational Predictions

### AlphaMissense
| Field | Value |
|---|---|
| am_pathogenicity | **0.9985** |
| am_class | **LPath** |
| Interpretation | Likely pathogenic (threshold 0.564) |

### DynaMut2
| Field | Value |
|---|---|
| ΔΔG (kcal/mol) | **-0.31 (Destabilising)** |
| Job ID | 177990909658 |
| Result URL | https://biosig.lab.uq.edu.au/dynamut2/results_prediction/177990909658 |

---

## Clinical Evidence

| Field | Value |
|---|---|
| Classification | **Pathogenic/Likely pathogenic** |
| Review status | criteria provided, multiple submitters, no conflicts |
| Last evaluated | 2025/07/13 00:00 |
| Inheritance | Inheritance not specified in this ClinVar entry. Given the multiple submitters with consistent Pathogenic/Likely pathogenic classification, the variant likely contributes to the WFS1 spectrum across both AR and AD presentations. |
| WFS1 variant landscape | S443R is 1 of ~326 pathogenic-spectrum variants in WFS1 (out of 2,243 catalogued in ClinVar) |

- (no specific conditions catalogued for S443R — ClinVar Pathogenic/Likely pathogenic by review evidence)

---

## 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 only
Stable fold + functional site hit  →  CATEGORY 4 — site-specific docking
```

## Final Schema Categorization

**Category 3/4 — Most Druggable**

<strong>Category 3/4 — Most Druggable.</strong> |ΔΔG| = 0.31 kcal/mol — fold survives. AlphaMissense 0.999 (near-maximum) confirms severe functional consequence despite modest structural cost.<br/><br/>The mechanism is charge-into-membrane plus TM3-TM4 helix-helix interface disruption. S443R is one of several Atlas variants where a charged residue is introduced into a bilayer-embedded position (compare T641K, V536E). Across these cases, the therapeutic target is the helix-helix interface that the wild-type residue stabilized, not the helix itself.<br/><br/>The therapeutic strategy is site-specific: a small molecule that occupies the TM3-TM4 interface near the PHE365 contact would compensate for the disrupted packing. The atlas-derived structural framework makes this target geometry visible — pre-atlas, the small ΔΔG would have deprioritized this variant for structure-based drug design.

**Why this card matters.** S443R is the cleanest example in the Atlas of the "charge-into-membrane" mechanism class. The protein folds, the variant is pathogenic, and the mechanism is helix-helix interface disruption inside the bilayer. Drug discovery for this class of variant targets the TM-TM interface, not the individual TM helix. The Atlas surfaces this whole class by surfacing the cross-helix neighbor contacts in the PDB analysis.

---

## Files in this folder

- `AF-O76024-F1-model_v6.pdb` — AlphaFold structure
- `S443R_molstar_viewer.html` — interactive 3D viewer (auto-highlights position 443 with ball-and-stick + neighbors within 5Å)
- `S443R_variant_card.md` — this card (source of truth)
- `S443R_variant_card.html` — styled printable card
- `S443R_dynamut2_summary.html` — clean offline DynaMut2 result card
- `dynamut2_result.json` — structured result data
- `dynamut2_result_page.html` — local snapshot of the Biosig result page (asset URLs absolutized)
- `S443R_wildtype_interactions.pse` / `S443R_mutant_interactions.pse` — PyMOL sessions

---

*Generated by wolfram-variant-card skill · RareResearch.AI Molecular Atlas*
*Every assumption documented. Every score sourced.*
