Mapping protein neutral networks from predicted secondary structure

Dit onderzoek construeert een operationeel genotype-phenotype-kaart voor influenza hemagglutinine op basis van voorspelde secundaire structuren en onthult dat neutraliteit lokaal gestructureerd maar globaal beperkt is, waardoor de structurele evolvabiliteit aanzienlijk kleiner is dan bij RNA-systemen.

De kern

De Onzichtbare Landkaart van de Influenza-virus: Een Verhaal over Mutaties en Toevalligheid

Stel je voor dat het DNA van een virus (in dit geval de griep) een enorme bibliotheek is met miljarden boeken. Elk boek is een andere versie van het virus. De meeste boeken lijken op elkaar, maar sommige hebben een paar woorden veranderd. De vraag die de onderzoekers stellen is: Als we een woord in een boek veranderen, verandert het verhaal dan nog steeds, of blijft het verhaal hetzelfde?

In de wetenschap noemen we dit een Genotype-Phenotype-kaart.

• Genotype: De letters in het boek (de DNA-sequentie).
• Phenotype: Het verhaal dat eruit komt (de vorm en functie van het eiwit).

Deze studie kijkt specifiek naar het hemagglutinine (HA) eiwit van het griepvirus. Dit is het "sleutel" die het virus gebruikt om in onze cellen te komen. De onderzoekers wilden weten: hoe makkelijk kan dit eiwit veranderen zonder dat het zijn functie verliest?

Hier is wat ze ontdekten, vertaald naar alledaagse beelden:

1. Het is geen open veld, maar een labyrint

Bij RNA-virussen (zoals sommige andere virussen) is het alsof je door een open veld loopt. Als je een stap zet, kun je vaak nog steeds in dezelfde richting blijven lopen, zelfs als je ver weg bent van waar je begon. Er zijn enorme "paden" (neutrale netwerken) waar je heel lang kunt wandelen zonder dat je valt.

Bij eiwitten zoals HA is het echter meer als een labyrint met hoge muren.

• De onderzoekers ontdekten dat er weliswaar paden zijn, maar deze zijn vaak kort en eindig.
• Je kunt een paar stappen zetten zonder dat het eiwit "breekt", maar als je te ver gaat, stort de structuur in.
• De analogie: Stel je voor dat je een toren bouwt met LEGO-blokken. Bij RNA kun je de toren heel hoog maken en nog steeds nieuwe blokken toevoegen zonder dat hij omvalt. Bij HA (eiwitten) is de toren veel stugger. Als je één blokje verplaatst, moet je heel voorzichtig zijn, want de toren is al snel instabiel.

2. De "Ster"-vorm van de netwerken

De onderzoekers zagen dat de mogelijke paden eruitzien als sterren.

• Er is één centraal punt (een heel stabiel virus) waar veel paden naartoe leiden.
• Maar deze paden overlappen elkaar nauwelijks. Het is alsof je in een bos bent met veel kleine, afzonderlijke kluizenaars. Je kunt binnen één kluisje veel rondlopen, maar om naar een andere kluis te komen, moet je vaak een hele lange, moeilijke weg nemen of een heel specifiek "bruggetje" vinden.
• Dit betekent dat het virus niet zomaar overal naartoe kan evolueren. Het zit vast in zijn eigen kleine hoekje.

3. Niet alle plekken zijn even belangrijk

Niet alle onderdelen van het eiwit zijn even gevoelig voor veranderingen.

• De "Stam" (Stem): Het onderste deel van het eiwit is als een betonnen fundering. Hier kunnen er veel veranderingen plaatsvinden (mutaties) zonder dat het eiwit kapot gaat. Het is robuust.
• De "Lussen" en "Kleefpunten": De bovenkant en de delen die contact maken met het menselijk lichaam zijn als glazen bloemen. Hier mag je niets veranderen. Als je hier één letter in het DNA verandert, breekt het eiwit of werkt het niet meer.
• De les: Het virus kan op sommige plekken vrij spelen, maar op andere plekken moet het zich heel streng aan de regels houden.

4. Kleine stapjes, geen sprongen

Wat gebeurt er als het virus toch verandert?

• De onderzoekers zagen dat het virus alleen kleine, voorzichtige stapjes kan maken.
• Als het eiwit verandert, blijft het resultaat er vaak heel veel op lijken. Het is alsof je een foto van een gezicht hebt en je verandert de kleur van de ogen een beetje. Het is nog steeds hetzelfde gezicht.
• Grote, radicale veranderingen (zoals een nieuw type griep dat er totaal anders uitziet) zijn extreem zeldzaam en moeilijk te bereiken via kleine mutaties. Het virus zit vast in een "bubbel" van kleine variaties.

Waarom is dit belangrijk?

Deze studie legt uit waarom het griepvirus zo lastig te verslaan is, maar ook waarom het niet zomaar een heel nieuw monster wordt.

• Het virus kan zich aanpassen (door kleine mutaties op de veilige plekken) om ons immuunsysteem te ontlopen.
• Maar het kan niet zomaar alles veranderen. De fysieke bouw van het eiwit (de "architectuur") beperkt waar het naartoe kan.

Samenvattend:
Het evolueren van het griepvirus is niet als het vrij rondzwemmen in een oceaan (zoals bij sommige andere virussen). Het is meer als het lopen door een gigantisch, complex kasteel met smalle gangen. Je kunt veel doen binnen één kamer, maar om naar een andere vleugel te komen, moet je door specifieke deuren. De muren van het kasteel (de fysieke structuur van het eiwit) bepalen waar je naartoe kunt, niet alleen de keuzes die het virus maakt.

Dit helpt wetenschappers beter te begrijpen hoe het virus zich ontwikkelt en waarom sommige veranderingen wel mogelijk zijn en andere niet.

Technische samenvatting

Probleemstelling

Evolutionaire uitkomsten worden beperkt door het feit dat slechts een subset van fenotypes zowel functioneel als toegankelijk is. Hoewel genotype-fenotype (GP) kaarten voor RNA uitgebreid zijn gekarakteriseerd, blijft dit voor eiwitten grotendeels onontgonnen vanwege de complexiteit van eiwitvouwing. In tegenstelling tot RNA, waar baseparing relatief eenvoudig is, omvatten eiwitvouwingen langeafstandsinteracties tussen chemisch diverse residuen. Dit kan leiden tot gefragmenteerde neutrale netwerken (NC's) en beperkte mutatie-toegankelijkheid. De vraag is of de eigenschappen die bij RNA zijn waargenomen (zoals grote, percolerende neutrale netwerken die toegang geven tot verre fenotypes) ook gelden voor echte eiwitten.

Methodologie

De auteurs hebben een operationele GP-kaart geconstrueerd voor het influenza hemagglutinine (HA) eiwit, waarbij voorspelde secundaire structuur dient als een grofkorrelig fenotype.

1. Data: Er werden 19.289 volledige HA-sequenties (lengte L=566) uit openbare databases (NCBI, GISAID) gehaald. Deze werden vertaald naar secundaire structuurstrings (Helix, Strand, Coil) met behulp van Porter5.
2. Definitie van de GP-kaart: De kaart $f: G \to \Phi$ koppelt elke aminozuursequentie aan een secundaire structuur. Een neutrale component (NC) is een maximaal verbonden subset van sequenties die naar hetzelfde fenotype leiden via enkelvoudige mutaties.
3. Schatting van NC-grootte en Robuustheid:
   • Omdat een volledige enumeratie van de ruimte ($20^{566}$) onmogelijk is, gebruikten ze statistische schatters.
   • Site-scanning: Deterministisch scannen van sequenties om te zien hoeveel alternatieve aminozuren op elke positie neutraal zijn.
   • Exhaustieve lokale enumeratie: Voor een steekproef van 20 genotypen per stam werden alle 19 mogelijke enkelvoudige mutaties per positie getest om lokale tolerantie en netwerktopologie te bepalen.
   • Schattingen: De grootte van de NC ($S_{NC}$) en de robuustheid ($r_{NC}$) werden berekend op basis van de "site versatility" (het effectieve aantal toegestane aminozuren per positie).
4. Netwerkanalyse: De lokale topologie van de NC's werd gereconstrueerd als grafen (knooppunten = sequenties, randen = enkelvoudige mutaties) om connectiviteitspatronen te analyseren.

Belangrijkste Bijdragen

• Eerste operationele GP-kaart voor een groot eiwit: De studie biedt een empirisch raamwerk om de structuur van neutrale netwerken in eiwitten te analyseren, vergelijkbaar met wat er voor RNA is gedaan.
• Kwantificering van beperkingen: Het toont aan dat eiwitneutraliteit fundamenteel anders is dan RNA-neutraliteit, voornamelijk door de aard van de vouwing en de grootte van de genotype-ruimte.
• Locale vs. Globale structuur: Het onderscheidt tussen lokale, dichte clusters van neutraliteit en de globale, beperkte toegankelijkheid van nieuwe structuren.

Resultaten

1. Sterke Fenotypische Bias:

   • De geschatte grootte van NC's varieert over meer dan 200 ordes van grootte ($10^{176}$ tot $10^{444}$).
   • Een klein aantal secundaire structuren domineert de sequentieruimte (Zipf-achtige verdeling), terwijl de meeste structuren slechts kleine gebieden innemen.

2. Robuustheid en Schaling:

   • Er is een positieve correlatie tussen NC-grootte en mutatie-robuustheid, maar deze relatie is veel zwakker dan bij RNA.
   • De helling van de relatie ($\beta \approx 0.001$) is veel ondieper dan bij RNA ($\beta \approx 0.1$). Dit komt doordat HA-netwerken een verwaarloosbaar klein deel van de enorme genotype-ruimte innemen; ze zijn "boundary-dominated" (grens-gedomineerd) en daardoor minder robuust.

3. Lokale Topologie:

   • NC's vertonen een ster-achtige, modulaire structuur. Ze bestaan uit dichte, positie-gerichte clusters (cliques) die verbonden zijn door een klein aantal goed verbonden "seed"-genotypen.
   • Er is weinig overlap tussen deze clusters, wat suggereert dat de connectiviteit heterogeen en zwak gepereerd is.
   • De netwerken zijn lokaal dicht, maar globaal gefragmenteerd.

4. Evolueerbaarheid en Structuurtransities:

   • Toegankelijke nieuwe fenotypes (via enkelvoudige mutaties) zijn sterk redundant en structuur-geproximeerd.
   • De meeste niet-neutrale mutaties leiden tot structuren die zeer vergelijkbaar zijn met het oorspronkelijke fenotype (kleine Hamming-afstand).
   • "Nieuwheid" is voornamelijk incrementeel; grote structurele sprongen zijn zeldzaam en moeilijk toegankelijk.

5. Invloed van Structurele Context:

   • Robuustheid is niet uniform verdeeld. Helix-structuren in de stam-domeinen zijn robuuster dan lussen of $\beta$-platen.
   • Functioneel kritieke gebieden (zoals fusiepeptiden en antigenische sites) tonen een lagere robuustheid.

Betekenis en Conclusie

De studie concludeert dat neutraliteit in eiwitten lokaal gestructureerd maar globaal beperkt is. In tegenstelling tot RNA, waar grote neutrale netwerken vaak door de hele sequentieruimte percoleren en verre fenotypes toegankelijk maken, zijn eiwitnetwerken in HA gefragmenteerd en beperkt door lokale positietolerantie.

Dit impliceert dat:

• De evolvabiliteit van eiwitten beperkt is tot incrementele, redundante veranderingen.
• Historische contingentie (de volgorde van mutaties) cruciaal is, omdat vroege mutaties de toegankelijke paden in de sequentieruimte onherroepelijk kunnen veranderen.
• De intrinsieke biofysische beperkingen van eiwitvouwing de evolutionaire trajecten sterker beperken dan eerder werd aangenomen, wat verklaart waarom structurele innovatie in HA zeldzaam is ondanks sterke immuunselectie.

Deze bevindingen verfijnen de neutrale netwerktheorie en benadrukken dat de mechanismen van neutraliteit fundamenteel verschillen tussen RNA en eiwitten.