Why structural divergence varies among residues in enzyme evolution: contributions of mutation, stability, and activity constraints

Dit onderzoek toont aan dat variaties in structurele divergentie tussen enzymresiduen worden bepaald door een familie-specifiek evenwicht tussen mutatie, stabiliteit en activiteit, waarbij het MSA-model deze selectiedrukken kan reconstrueren uit divergentieprofielen.

De kern

Waarom sommige delen van een enzym veranderen en andere niet: Een verhaal over mutaties, stabiliteit en werk

Stel je voor dat een enzym een complexe machine is, zoals een oude, goed onderhouden fiets. Elke schroef, elk tandwiel en elke ketting is een aminozuur (een 'residu'). In de loop van de evolutie ondergaan deze fietsen kleine reparaties en aanpassingen (mutaties). Maar als je naar een hele reeks van deze fietsen kijkt, zie je een interessant patroon: sommige onderdelen zijn na duizenden jaren nog steeds bijna hetzelfde, terwijl andere onderdelen er heel anders uitzien.

De vraag die Julian Echave en Mathilde Carpentier in hun onderzoek stellen, is: Waarom is dat zo? Waarom veranderen sommige onderdelen van de machine wel en andere niet?

Ze hebben een nieuw model ontwikkeld om dit te verklaren, gebaseerd op drie krachten die in het spel zijn: Mutatie, Stabiliteit en Activiteit. Laten we dit uitleggen met een paar creatieve vergelijkingen.

1. De Drie Krachten in het Spel

Stel je de evolutie van een enzym voor als een spelletje waarbij je blindelings onderdelen van je fiets vervangt.

• De Kracht van de Mutatie (Het Willekeurige Vervangen):
  Soms vervang je een onderdeel puur door toeval. Als je een losse bout vervangt door een andere bout, kan het zijn dat de fiets er net iets anders uitziet. Sommige onderdelen van de fiets (zoals de handvatten) zijn flexibel en kunnen veel verschillende vormen aan zonder dat de fiets kapot gaat. Andere onderdelen (zoals de as van het wiel) zijn stijf; als je die vervangt, stort de hele fiets in elkaar.

  • In het kort: Mutaties zorgen voor variatie, maar ze slaan harder toe op de flexibele delen van het eiwit dan op de stijve delen.

• De Kracht van de Stabiliteit (Het Niet-Kapot-Maken):
  Je wilt dat je fiets blijft staan. Als je een onderdeel vervangt dat ervoor zorgt dat het frame breekt, gooi je die fiets weg. In de natuur betekent dit: als een mutatie het eiwit instabiel maakt (het 'smelt' of valt uit elkaar), wordt die mutatie niet doorgegeven.

  • In het kort: De natuur selecteert tegen mutaties die het eiwit instabiel maken. Dit houdt bepaalde delen van het eiwit 'strak' in vorm.

• De Kracht van de Activiteit (Het Goed Werken):
  Een fiets moet ook kunnen rijden. Als je een onderdeel vervangt waardoor de ketting niet meer goed loopt, is de fiets nutteloos. Voor een enzym betekent dit: als een mutatie het actieve centrum (waar het werk gebeurt) verstoort, werkt het enzym niet meer goed. De natuur gooit deze mutaties ook weg.

  • In het kort: De natuur selecteert tegen mutaties die het werk van het enzym verstoren.

2. Het Grote Geheim: Het Evenwicht Verschilt per Familie

Het verrassende resultaat van dit onderzoek is dat elk enzym een eigen verhaal heeft.

In sommige families van enzymen is het vooral de stabiliteit die bepaalt hoe de onderdelen eruitzien. Het is alsof je een fiets hebt die in een storm wordt gebruikt; alles moet extreem stevig zijn, dus verandert er weinig.
In andere families is het vooral de activiteit die telt. Stel je een racefiets voor: het frame mag misschien wat buigen, maar de versnellingen en de ketting moeten perfect blijven om snel te blijven.
En in weer andere families is het vooral de willekeurige mutatie (de flexibiliteit van het materiaal) die bepaalt hoe het eruitziet, omdat de andere regels minder streng zijn.

De auteurs hebben 34 verschillende groepen (families) van enzymen onderzocht. Ze ontdekten dat:

• Mutaties altijd een rol spelen (het is de basis van verandering).
• Maar of stabiliteit of activiteit de baas is, verschilt enorm per groep. Soms is stabiliteit de belangrijkste, soms activiteit, en soms spelen ze allebei een kleine rol.

3. De "Vorm" van het Eiwit Vertelt het Verhaal

Hoe weten ze dit? Ze hebben een slimme rekenmethode gebruikt (het MSA-model). Ze kijken naar de "vorm" van de veranderingen in de enzymen.

• Als een enzym erg flexibel is (zoals een zachte rubberen fiets), zie je dat de mutaties zich vooral daar manifesteren.
• Als een enzym heel streng wordt geselecteerd op zijn werk (activiteit), zie je dat de onderdelen die belangrijk zijn voor het werk, juist heel weinig veranderen, terwijl de rest wel mag veranderen.

Het mooie is: door naar de uiteindelijke vorm van het eiwit te kijken, kunnen de onderzoekers terugrekenen hoe streng de natuur op dat moment heeft geselecteerd. Het is alsof je naar de slijtage van een fiets kijkt en kunt zeggen: "Ah, deze fiets is veel gebruikt in de modder (stabiliteit was belangrijk)" of "Deze fiets is gebruikt op een racecircuit (activiteit was belangrijk)".

Conclusie: Een Unieke Identiteit voor Elk Enzym

Kortom, dit onderzoek laat zien dat er geen één groot recept is voor hoe eiwitten evolueren. Elk enzym heeft zijn eigen "persoonlijkheid" die wordt gevormd door de balans tussen:

1. Hoe makkelijk het materiaal verandert (mutatie).
2. Hoe bang we zijn dat het in elkaar stort (stabiliteit).
3. Hoe bang we zijn dat het niet meer werkt (activiteit).

De auteurs hebben een model gebouwd dat deze complexe balans kan voorspellen. Het helpt ons niet alleen te begrijpen waarom eiwitten eruitzien zoals ze eruitzien, maar ook om te zien hoe de omgeving (de selectiedruk) de bouw van het leven heeft beïnvloed. Het is een beetje alsof we de blauwdrukken van de evolutie kunnen lezen, niet alleen door naar de genen te kijken, maar door naar de vorm van de machines zelf te kijken.

Technische samenvatting

Probleemstelling

In de evolutie van enzymen varieert de structurele divergentie (de mate waarin de 3D-structuur verandert) aanzienlijk tussen verschillende residuen binnen hetzelfde eiwit. Sommige posities blijven bijna onveranderd, terwijl andere sterk variëren, wat leidt tot residu-afhankelijke profielen van structurele divergentie. Hoewel empirische studies hebben aangetoond dat deze patronen correleren met flexibiliteit en afstand tot het actieve centrum, ontbreekt er een mechanistisch model dat uitlegt waarom deze verschillen bestaan en welke evolutionaire krachten (mutatie, stabiliteit, activiteit) hieraan bijdragen. Bestaande diffusiemodellen nemen aan dat alle residuen volgens hetzelfde proces evolueren en kunnen deze variatie dus niet verklaren.

Methodologie

De auteurs bouwen voort op het Mutation-Stability-Activity (MSA) model, een mechanistisch mutatie-selectiemodel dat eerder werd ontwikkeld voor sequentie-evolutie, en passen dit toe op structurele evolutie. De kern van de methodologie bestaat uit de volgende componenten:

1. Het MSA-model:

   • Evolutie wordt gemodelleerd als een proces waarbij mutaties worden geïntroduceerd en gefixeerd of verloren gaan op basis van hun effecten op eiwitstabiliteit ($\Delta\Delta G$) en katalytische activiteit ($\Delta\Delta G^\ddagger$).
   • De fixatiekans wordt bepaald door selectieparameters $a_S$ (voor stabiliteit) en $a_A$ (voor activiteit). Mutaties die de stabiliteit of activiteit niet compromitteren, fixeren met een neutrale snelheid; schadelijke mutaties fixeren met een exponentieel afnemende kans.

2. Linearly Forced Elastic Network Model (LFENM):

   • Om de fysieke effecten van mutaties te berekenen, gebruiken de auteurs LFENM. Het eiwit wordt gemodelleerd als een netwerk van $C_\alpha$-atomen verbonden door harmonische veren.
   • Een mutatie wordt gesimuleerd door willekeurige verstoringen toe te voegen aan de lengte van de veren rondom het gemuteerde residu.
   • Hieruit worden drie grootheden berekend:
     • $\Delta r_0$: De structurele verplaatsing (RMSD) door relaxatie van het netwerk.
     • $\Delta\Delta G$: De verandering in vouwingsvrije energie (stabiliteit).
     • $\Delta\Delta G^\ddagger$: De verandering in activeringsenergie (activiteit), gebaseerd op de kosten om het actieve centrum terug te brengen naar de optimale geometrie.

3. Data en Analyse:

   • De studie omvat 34 families van functioneel geconserveerde, monomeer, enkel-domein enzymen uit de M-CSA database.
   • Voor elke familie wordt een waargenomen structureel divergentieprofiel gegenereerd door de $C_\alpha$ RMSD te meten tussen homologue structuren en een referentie-enzym.
   • Bayesiaanse inferentie wordt gebruikt om de parameters $a_S$ en $a_A$ te schatten door de voorspelde profielen van het MSA-model aan te passen aan de waargenomen profielen.
   • Geneste modelvergelijkingen: Er worden vier modellen vergeleken:
     • $M_0$: Uniforme divergentie (null-model).
     • $MM$: Alleen mutatie-effecten ($a_S = a_A = 0$).
     • $MS$: Mutatie + stabiliteitselectie ($a_A = 0$).
     • $MSA$: Volledig model (mutatie + stabiliteit + activiteit).
   • Decompositie: De voorspelde profielen worden opgesplitst in drie componenten ($\phi_M, \phi_S, \phi_A$) om de relatieve bijdrage van elke beperking te kwantificeren via de standaarddeviatie over de residuen.

Belangrijkste Bijdragen

• Extensie van het MSA-model: Voor het eerst wordt het MSA-model succesvol toegepast om structurele divergentieprofielen te voorspellen, niet alleen sequentie-divergentie.
• Mechanistische verklaring: Het biedt een fysisch onderbouwd mechanisme voor de variatie in structurele divergentie, waarbij de interactie tussen de mechanische respons van het eiwit en evolutionaire selectie wordt gekwantificeerd.
• Kwantificering van selectiedruk: Het model stelt de auteurs in staat om familiespecifieke selectiestrengths ($a_S$ en $a_A$) af te leiden uit structurele data, wat een brug slaat tussen biophysica en evolutionaire biologie.

Resultaten

1. Modelnauwkeurigheid: Het MSA-model (met slechts twee vrije parameters) voorspelt de waargenomen structurele divergentieprofielen even nauwkeurig als een complexer empirisch model (M12). De gemiddelde correlatie is 0,66.
2. Bijdrage van beperkingen: Alle drie de beperkingen (mutatie, stabiliteit, activiteit) dragen bij aan de structurele divergentieprofielen:
   • Mutatie: Draagt altijd substantieel bij (gemiddeld ~47% van het profiel). Dit komt door de heterogeniteit in de flexibiliteit van het eiwit.
   • Stabiliteit en Activiteit: Hun bijdrage varieert sterk per familie (van verwaarloosbaar tot dominant). Gemiddeld dragen ze respectievelijk ~33% en ~20% bij.
3. Variatie tussen families: Er is geen universele hiërarchie van beperkingen. In sommige families is stabiliteit de dominante factor, in andere activiteit.
   • De bijdrage van mutatie wordt bepaald door de architectuur van het eiwit (specifiek de heterogeniteit in flexibiliteit).
   • De bijdragen van stabiliteit en activiteit worden voornamelijk bepaald door de sterkte van de selectie ($a_S$ en $a_A$), niet door de verdeling van de mutatie-effecten zelf.
4. Correlaties:
   • De mutatiebijdrage correleert sterk met de heterogeniteit in flexibiliteit ($\rho = 0,99$).
   • De stabiliteitsbijdrage correleert sterk met de parameter $a_S$ ($\rho = 0,93$).
   • De activiteitsbijdrage correleert sterk met de parameter $a_A$ ($\rho = 0,94$).

Significantie en Conclusie

Dit onderzoek toont aan dat structurele divergentieprofielen niet alleen informatie bevatten over de architectuur van enzymen, maar ook coderen voor het selectieregime waaronder ze evolueren.

• Oplossing voor een open probleem: Het model suggereert dat de variatie in de snelheid van structurele divergentie tussen verschillende enzymfamilies kan worden verklaard door verschillen in de selectiestrengths ($a_S$ en $a_A$), een probleem dat eerder onopgelost bleef.
• Biologische connectie: De geschatte parameters $a_S$ en $a_A$ kunnen worden gekoppeld aan biologische variabelen zoals expressieniveau (voor stabiliteit) en metabole flux of functionele optimalisatie (voor activiteit).
• Toekomstperspectief: Het MSA-model fungeert als een krachtig "nulmodel" om afwijkingen te detecteren die wijzen op extra beperkingen (zoals cofactor-binding of allosterie) die nog niet in het model zijn opgenomen.

Samenvattend biedt dit werk een mechanistisch inzicht in hoe de fysieke eigenschappen van eiwitten en evolutionaire selectiedruk samenwerken om de structurele diversiteit van enzymen te vormen.