Accessing pore-blocker bound and open conformations of TMEM16A using PIP2-assisted adaptive sampling

Deze studie toont aan dat het gebruik van het lipide PIP2 in combinatie met adaptieve sampling de overgang van het TMEM16A-kanal van een gesloten naar een stabiele open conformatie mogelijk maakt, waardoor de bindingsplaatsen van pore-blokkers kunnen worden geïdentificeerd voor gerichte medicijnontwikkeling.

De kern

Titel: Hoe een kleine lipide-sleutel een vergrendelde deur opent voor medicijnen

Stel je voor dat je lichaam een enorme stad is, en de cellen zijn de huizen. Om deze huizen veilig te houden, hebben ze poorten: ionkanalen. Deze poorten laten bepaalde stoffen (zoals zouten) binnen of buiten, wat essentieel is voor spierbeweging, zenuwsignalen en meer.

Een van deze poorten heet TMEM16A. Als deze poort te vaak of op het verkeerde moment openstaat, kan dat leiden tot ziektes zoals astma, beroertes of kanker. Wetenschappers willen daarom medicijnen maken die deze poort kunnen sluiten of regelen. Maar er is een groot probleem: om een goede sleutel (een medicijn) te maken, moet je precies weten hoe de poort eruitziet als hij open is.

Helaas is het heel moeilijk om de "open" stand van deze poort vast te leggen met een camera (zoals een microscopie). Het is alsof je probeert een foto te maken van een deur die zo snel opengaat en dichtslaat, dat je camera te traag is.

De oplossing: Een slimme digitale simulatie

In dit onderzoek hebben de auteurs een digitale oplossing bedacht. Ze hebben een computermodel gemaakt van de poort en een slimme truc gebruikt om de deur te laten openen, zodat ze de "open" stand konden bestuderen.

Hier is hoe ze het deden, vertaald naar alledaagse taal:

1. De sleutel die de deur opent (PIP2)

De poort heeft twee sleutels nodig om open te gaan: Calcium (een mineraal) en PIP2. PIP2 is een klein vetmolecuul dat van nature in de celwand zit.

• De analogie: Stel je voor dat de poort een zware, roestige schuifdeur is. Calcium is de persoon die duwt, maar PIP2 is de olie die de scharnieren soepel maakt. Zonder die olie blijft de deur vastzitten.
• De onderzoekers gebruikten PIP2 in hun computerprogramma om de deur soepel te laten glijden van de gesloten naar de open stand.

2. De digitale dans (Adaptive Sampling)

Zelfs met de olie (PIP2) is het moeilijk om de deur volledig open te krijgen in een simulatie. De deur wil vaak terugsluiten.

• De analogie: Stel je voor dat je een danspartner hebt die steeds terugtrekt als je te dichtbij komt. De onderzoekers gebruikten een slim algoritme genaamd FAST. Dit is alsof je een danspartner hebt die je moedigt om steeds een stapje verder te gaan. Als de deur begint te openen, geeft het programma een beloning en laat het de simulatie doorgaan. Als de deur dichtgaat, start het opnieuw.
• Door dit duizenden keren te doen, konden ze eindelijk de poort volledig open zien staan.

3. De verrassing: De "tussenstand" is de echte doelwit

Toen ze de poort volledig open zagen, probeerden ze er medicijnen in te stoppen om te zien of ze de poort konden blokkeren.

• Het probleem: De poort was te open! De medicijnen vielen er zo doorheen of konden niet goed vasthouden. Het was alsof je probeert een deur te blokkeren met een steen, maar de deur staat zo wijd open dat de steen er gewoon doorheen valt.
• De ontdekking: Ze realiseerden zich dat de medicijnen waarschijnlijk niet op het moment van de maximale opening werken, maar in de tussenstand (de deuropening die half open is).
• De conclusie: Ze zagen dat medicijnen zoals 1PBC en A9C perfect passen in deze "half-open" stand. Ze vonden precies waar ze vasthaken (aan bepaalde aminozuren, die je kunt zien als de grepen van de deur).

4. Waarom werkt dit op de ene ziekte en niet op de andere?

De onderzoekers keken ook naar een andere poort, TMEM16B, die heel veel lijkt op TMEM16A. Een medicijn genaamd Ani9 blokkeert TMEM16A, maar werkt niet op TMEM16B.

• De analogie: Stel je voor dat TMEM16A en TMEM16B twee bijna identieke huizen zijn. Het medicijn Ani9 is een sleutel die alleen in het slot van het eerste huis past.
• Door hun simulatie zagen ze dat het slot van het tweede huis (TMEM16B) net iets anders is gevormd. Het slot is te smal of heeft een ander profiel, waardoor de sleutel er niet in past. Dit verklaart waarom het medicijn veilig is voor het ene orgaan en niet voor het andere.

Waarom is dit belangrijk?

Dit onderzoek is een doorbraak voor drie redenen:

1. Het vinden van onzichtbare deuren: Ze hebben een manier gevonden om de "open" stand van een poort te zien die wetenschappers nog nooit met een echte camera hebben kunnen vastleggen.
2. Beter medicijnontwerp: Nu weten ze precies hoe de poort eruitziet als hij half open is (waar de medicijnen werken). Dit helpt farmaceuten om medicijnen te bouwen die perfect in dat slot passen, net als een op maat gemaakte sleutel.
3. Minder bijwerkingen: Door te begrijpen waarom een medicijn op de ene poort werkt en op de andere niet, kunnen ze medicijnen maken die alleen de ziekte bestrijden en geen schade aanrichten aan gezonde delen van het lichaam.

Kortom:
De onderzoekers hebben een digitale "olie" (PIP2) gebruikt om een vastzittende poort te laten openen in een computer. Ze ontdekten dat de poort in een tussentijdse stand het beste werkt als doelwit voor medicijnen. Dit opent de deur voor nieuwe, veiligere behandelingen voor ziektes zoals astma en beroertes.

Technische samenvatting

Probleemstelling

Ionenkanalen zijn cruciale doelwitten voor de geneesmiddelenontwikkeling, maar er is een significant obstakel: de database van ionenkanaalstructuren bevat voornamelijk thermodynamisch stabiele conformaties die experimenteel zijn opgelost (bijv. gesloten of gedeeltelijk geopende toestanden). Veel "druggable" (drug-bespeelbare) toestanden, zoals de volledig open conformatie of specifieke intermediaire toestanden die nodig zijn voor blokkering, zijn experimenteel moeilijk te vangen of niet stabiel onder standaard simulatieomstandigheden. Dit beperkt de mogelijkheid tot structureel gebaseerd drugdesign, vooral voor kanalen zoals TMEM16A (een calcium-geactiveerd chloridekanaal), waarvan overactivatie geassocieerd is met ziekten zoals pulmonale hypertensie en ischemische beroerte. Bestaande modellen van de open toestand zijn vaak onvolledig of vertonen afwijkingen in geleidbaarheid in vergelijking met experimentele data.

Methodologie

De auteurs hebben een geavanceerde computationele aanpak ontwikkeld die drie kerncomponenten combineert:

1. PIP2-gestuurde Adaptive Sampling (FAST):

   • De studie maakt gebruik van Fluctuation Amplification of Specific Traits (FAST), een adaptieve sampling-techniek.
   • Als katalysator voor de conformatieovergang werd PIP2 (phosphatidylinositol-4,5-bisphosphate) gebruikt, een endogeen regulatormolecuul dat bekend staat als een allosterische modulator van TMEM16A.
   • Eerst werden Coarse-Grained (CG) MD-simulaties uitgevoerd om het PIP2-bindingsgebied te identificeren (tussen TM3, TM4 en TM5).
   • Vervolgens werd het systeem omgezet naar een All-Atom (AA) representatie. De FAST-algoritme werd toegepast om de afstand tussen de residuen die de porie bekleden te maximaliseren, waardoor het systeem wordt gedwongen om de gesloten toestand te verlaten en de open toestand te bereiken.

2. Markov State Modelling (MSM):

   • De gegenereerde trajecten werden geanalyseerd met MSM om de vrije-energielandschappen (Free Energy Landscapes - FEL) te construeren.
   • Hierdoor konden drie macrotoestanden worden geïdentificeerd: gesloten (C), intermediair (I) en open (O).

3. Versterkte Sampling voor Drugbinding (AWH):

   • Om de bindingsplaatsen van verschillende blokkers te bepalen, werd Accelerated Weighted Histogram (AWH) gebruikt.
   • Dit werd toegepast op de intermediaire en open toestanden om de vrije energie van de binding van vier verschillende geneesmiddelen te berekenen: 1PBC, A9C, niclosamide en Ani9.

4. Validatie:

   • De geleidbaarheid van de gegenereerde open toestand werd getest door een elektrisch veld toe te passen en ionenpermeatie te meten, vergeleken met experimentele data.
   • De specificiteit van de drug Ani9 voor TMEM16A versus TMEM16B werd onderzocht door middel van sequentie-uitlijning en simulaties van beide subtypes.

Belangrijkste Resultaten

• Generatie van een stabiele open toestand: Met PIP2-assisted FAST sampling slaagden de auteurs erin om een thermodynamisch stabiele open toestand van TMEM16A te genereren. Deze toestand vertoont drie cruciale structurele kenmerken:
  1. Een vergrote buitenste poort (outer gate).
  2. Een knik (kink) in de TM6-helix bij residu G654.
  3. Een breuk in de $\alpha$-helixstructuur in TM4 bij residu I551 (hoewel een volledige overgang naar een $\pi$-helix niet werd waargenomen).
• Validatie van geleidbaarheid: De berekende stroom-doorlaatbaarheid (conductance) van deze nieuwe open toestand stemde nauwkeurig overeen met experimentele waarden bij potentiaalverschillen tot +300 mV, wat aantoont dat het een fysiologisch relevante, volledig geleidende toestand is.
• De rol van PIP2 en Ca2+: Simulaties zonder PIP2 of Ca2+ toonden aan dat de open toestand dan in een gebied met hoge vrije energie ligt en dus onstabiel is. PIP2 verlaagt de energiebarrière en stabiliseert de open toestand, terwijl Ca2+ essentieel is voor het initiëren van de opening.
• Drugbinding en "Induced Fit":
  • 1PBC en A9C: Deze blokkers binden niet stabiel in de volledig open toestand (ze dissocieerden snel). Ze binden echter stabiel in de intermediaire toestand. Dit suggereert dat de drug-bound conformatie een "induced-fit" mechanisme vertegenwoordigt die niet identiek is aan de volledig open, geleidende toestand.
  • Niclosamide: De studie weerlegt eerdere hypothese dat niclosamide bindt op dezelfde plek als in TMEM16F. In TMEM16A bindt het in de intermediaire toestand aan de buitenste poort, met specifieke interacties met residuen zoals R535 en E623. Dit verklaart het remmende effect op TMEM16A.
  • Ani9 (Subtype-specificiteit): Ani9 is een specifieke blokker voor TMEM16A, maar niet voor TMEM16B. De studie toont aan dat Ani9 twee mogelijke bindingsplekken heeft: de buitenste poort (in de open toestand) en de binnenste poort (in de intermediaire toestand). De specificiteit wordt bepaald door een waterstofbrug met residu S592 (aanwezig in TMEM16A, afwezig in TMEM16B) en door de ruimtelijke beperking van de binnenste poort in TMEM16B, die te smal is voor de drug.

Bijdragen en Significantie

1. Nieuwe Computationele Paradigma: De studie introduceert een krachtige methode waarbij endogene liganden (zoals PIP2) worden gebruikt als "computational tools" om adaptieve sampling te sturen naar experimenteel onbereikbare conformaties. Dit overbrugt de kloof tussen bestaande PDB-structuren en de dynamische realiteit van ionenkanalen.
2. Structuur-specifiek Drugdesign: Het onderzoek benadrukt dat de "meest open" toestand van een kanaal niet altijd de "druggable" toestand is. Veel blokkers (zoals 1PBC en A9C) vereisen een intermediaire conformatie voor binding. Dit is cruciaal voor het ontwikkelen van nieuwe geneesmiddelen die specifiek gericht zijn op ziekte-gerelateerde conformaties.
3. Mechanistisch Inzicht: De studie biedt gedetailleerd inzicht in de moleculaire mechanismen van subtype-specificiteit (TMEM16A vs. TMEM16B) en de verschillende werkingsmechanismen van bestaande drugs (activatie vs. inhibitie van niclosamide).
4. Toekomstperspectief: De methode biedt een blauwdruk voor het ontdekken van nieuwe therapeutische doelwitten voor ionenkanalen die betrokken zijn bij ernstige ziekten, door gebruik te maken van AI-gestructureerde modellen (zoals AlphaFold) gecombineerd met fysisch gebaseerde versterkte sampling.

Kortom, deze studie demonstreert hoe geavanceerde moleculaire dynamica, gekoppeld aan biologische kennis van lipiden, kan worden gebruikt om de dynamische landschappen van ionenkanalen volledig te ontrafelen en zo de basis te leggen voor de volgende generatie structureel gebaseerde geneesmiddelen.