Markov State Model for the forced unfolding of a small peptide

Dit artikel toont aan dat een dynamische grofkorrelige techniek gebaseerd op Markov-staatmodellering, waarbij donor-acceptorafstanden van helicale waterstofbruggen als collectieve variabelen worden gebruikt, de atomaire details van het mechanische ontvouwingproces van een klein peptide dat geen eenvoudig twee-toestands- of coöperatief mechanisme volgt, nauwkeurig kan reconstrueren.

De kern

Stel je een klein, opgerold veertje voor dat bestaat uit een kleine keten van bouwstenen (een peptide). Wetenschappers willen begrijpen hoe dit veertje ontrolt wanneer je eraan trekt, net als het rekken van een stukje taffy.

Meestal gebruiken wetenschappers om dit te bestuderen krachtige computers om de beweging van elk afzonderlijk atoom te simuleren. Maar er is een probleem: het echte leven verloopt langzaam, terwijl computersimulaties vaak extreem snel moeten verlopen om binnen een redelijke tijd te eindigen. Het is alsof je probeert een film van een slak die kruipt te bekijken door de video met 100x snelheid af te spelen; je mist dan alle subtiele details van hoe het zijn poten beweegt.

Om dit op te lossen, hebben de onderzoekers in dit artikel een "slimme shortcut"-methode ontwikkeld die een Markov-toestandsmodel wordt genoemd. Denk aan deze methode niet als een video in hoge snelheid, maar als een stroomschema van mogelijkheden. In plaats van elke kleine wiebel van elk atoom bij te houden, groepeert de methode de vormen van het peptide in "toestanden" (zoals "opgerold", "half-ontrold" of "volledig uitgerekt") en berekent de kansen om van de ene toestand naar de andere te springen.

Hier is hoe ze dit toepasten op hun specifieke puzzel:

1. De verkeerde kaart versus de juiste kaart
Bij eerdere experimenten met eenvoudigere veertjes konden wetenschappers gewoon de totale lengte van het veertje meten (afstand van einde tot einde) om te weten wat er gebeurde. Als het veertje langer werd, ontrolde het.
Echter, dit specifieke peptide is lastig. Het ontrolt niet zomaar in een simpele, rechte lijn. Het heeft een "tussenliggende" toestand waarbij de uiteinden open zijn, maar het midden nog opgerold blijft.

• De analogie: Stel je een rits voor. Als je alleen de totale lengte van het jasje meet, kun je niet zeggen of de rits halverwege open is of dat het jasje gewoon vreemd is gevouwen. De lengte alleen is een slechte kaart.
• De oplossing: De onderzoekers realiseerden zich dat ze naar de "ritsen" binnenin het veertje moesten kijken – de waterstofbruggen die de spiralen bij elkaar houden. Ze hielden de afstand tussen de specifieke delen van deze bindingen (donor-acceptor-afstanden) bij om een veel duidelijker beeld te krijgen.

2. Het bouwen van het stroomschema
Ze voerden duizenden computersimulaties uit om te zien hoe het peptide bewoog.

• Ze gebruikten een wiskundige truc (TICA genoemd) om de complexe data te vereenvoudigen, net als een chef die een saus reduceert om de essentiële smaak te krijgen.
• Ze ontdekten dat ze door naar de totale lengte plus drie specifieke patronen van de interne bindingen te kijken, een betrouwbaar stroomschema konden maken. Dit stroomschema voorspelde nauwkeurig hoe het peptide zich gedroeg, zelfs wanneer het vastzat in die lastige "tussenliggende" toestand.

3. Het trek-experiment
Ze simuleerden het uit elkaar trekken van het peptide op verschillende snelheden:

• Snel trekken: Alsof je een tapijt onder een tafel vandaan trekt. Het peptide knapt gewelddadig open en de gemeten krachten zijn enorm.
• Langzaam trekken: Alsof je zachtjes taffy rekkt. Het peptide heeft tijd om te ontspannen en zijn natuurlijke pad te vinden.
• Het resultaat: Hun "slimme shortcut" (het Markov-model) werkte perfect voor langzaam trekken. Het kon de zachte, realistische krachten voorspellen die onmogelijk te simuleren zijn met standaardmethoden, omdat die methoden te lang zouden duren om uit te voeren.

4. Wat ze ontdekten
De studie onthulde dat dit peptide niet in één keer uit elkaar valt.

• Het pad: Het begint meestal met het openen aan één uiteinde (de "N-terminus") en ontrolt vervolgens als een rits.
• De valstrik: Soms blijft het vastzitten in een tussenliggende toestand waarbij de uiteinden open zijn, maar het centrum nog steeds een strakke spiraal is. Dit verklaart waarom het proces complexer is dan een simpele "aan/uit"-schakelaar.

Samenvattend
Het artikel laat zien dat je voor complexe, wiebelende moleculen niet alleen de totale lengte kunt meten om ze te begrijpen. Je moet kijken naar de interne verbindingen. Door een "stroomschema"-benadering te gebruiken die zich richt op deze interne verbindingen, creëerden de onderzoekers een methode die langzame, realistische trek-experimenten op een computer kan simuleren. Hierdoor kunnen ze de gedetailleerde stappen zien van hoe een molecuul ontrolt, wat eerder te langzaam was om te bekijken met standaard computersimulaties.

Technische samenvatting

Technische Samenvatting: Markov-toestandsmodel voor de geforceerde ontvouwing van een klein peptide

Probleemstelling
Enkelmolecuul-krachtspectroscopie-experimenten en hun computationele tegenhangers, krachtsonderzochte moleculaire dynamica (FPMD)-simulaties, zijn essentieel voor het bestuderen van de mechanische ontvouwing van biomoleculen. Er bestaat echter een aanzienlijke beperking: de trekkingssnelheden in standaard FPMD-simulaties zijn vaak vele ordes van grootte hoger dan die welke in experimenten toegankelijk zijn. Hoewel grofkorrelige technieken (coarse-graining) frequent worden ingezet om deze tijdschaalkloof te overbruggen, vertrouwen velen op structurele vereenvoudigingen of impliciete oplosmiddelenmodellen die atomaire details van conformatieovergangen verdoezelen.

Een specifieke uitdaging doet zich voor wanneer het bestudeerde systeem niet ontvouwt via een eenvoudig tweestapsmechanisme, of wanneer de eind-tot-eind-afstand ($r_{ee}$) niet voldoet als voldoende ordeparameter. In dergelijke complexe scenario's kan het aannemen van een enkele reactiecoördinaat leiden tot niet-Markoviaanse dynamica, waardoor standaard extrapolatiemethoden ongeldig worden. Eerder werk van de auteurs heeft dynamische grofkorreliging succesvol toegepast op een eenvoudig tweestaps calixareen-systeem, maar de toepasbaarheid van deze methode op complexere systemen met intermediaire toestanden en ontoereikende ordeparameters bleef ongetest.

Methodologie
De auteurs onderzochten de mechanische ontvouwing van een $\beta$-alanine-octameer in methanol, een systeem dat bekend staat om zijn complexe, meerstaps ontvouwingstraject met een metastabiel intermediair stadium. Het onderzoek maakte gebruik van een combinatie van atomaire FPMD-simulaties en een dynamische grofkorreligingsbenadering gebaseerd op Markov-toestandsmodellen (MSM's).

1. Systeem en simulatieopzet:

   • Model: Een $\beta$-alanine-octameer dat een stabiele $3_{14}$-helix vormt, gestabiliseerd door zes waterstofbruggen (H-bruggen).
   • FPMD-simulaties: Uitgevoerd met GROMACS 2020.1 en het GROMOS 53a6-krachtenveld bij $T=240$ K. De simulaties maakten gebruik van een krachtstroomprotocol met een constante trekkingssnelheid en een harmonische veer ($K=1$ N/m).
   • Evenwichtssimulaties: Uitgevoerd bij $T=298$ K om thermische ontvouwingstrajecten te analyseren.

2. Collectieve variabelen (CV's) en dimensiereductie:

   • In het besef dat $r_{ee}$ alleen ontoereikend is, selecteerden de auteurs de donor-acceptor-afstanden van de zes helicale H-bruggen als primaire collectieve variabelen.
   • Time-lagged Independent Component Analysis (TICA) werd toegepast om de dimensie van de dynamica te reduceren. De auteurs identificeerden dat de drie meest relevante onafhankelijke componenten (IC1, IC2, IC3), gecombineerd met $r_{ee}$, meer dan 95% van de cumulatieve kinetische variantie vingen en de Markoviaanse aard van de dynamica herstelden.

3. Constructie van het Markov-toestandsmodel:

   • De gereduceerde viervoudige configuratieruimte werd gediskretiseerd in 500 toestanden met behulp van $k$-means clustering.
   • Voor thermische ontvouwing werd een overgangsmatrix geconstrueerd en grofkorrelig gemaakt met behulp van Perron Cluster Cluster Analysis (PCCA+) om 7 metastabiele toestanden te identificeren.
   • Voor krachtstroom-simulaties werd een multi-ensemble-benadering gebruikt. De bemonsteringscoördinaat ($r_s$) werd gediskretiseerd in 21 ensembles. Umbrella Sampling (US)-berekeningen werden uitgevoerd voor elk ensemble om de Potentiaal van Gemiddelde Kracht (PMF) en overgangssnelheden te berekenen.
   • De Transition Matrix Reweighting Analysis Method (TRAM) werd ingezet om overgangssnelheidsmatrices ($W^{(m)}$) voor elk ensemble te schatten, rekening houdend met de tijdsafhankelijkheid van het trekprotocol.

4. Reconstructie van de krachtstroom:

   • De tijdsontwikkeling van toestandspopulaties werd berekend door het systeem te propageren door de reeks ensembles die door de trekkingssnelheid zijn gedefinieerd.
   • De gemiddelde kracht werd afgeleid uit de tijdsafhankelijke populaties en de gemiddelde eind-tot-eind-afstand, waardoor reconstructie van Kracht-Extensiecurven (FEC's) mogelijk werd bij trekkingssnelheden die veel trager zijn dan die haalbaar in directe FPMD.

Belangrijkste resultaten

• Validatie van de ordeparameter: Het onderzoek bevestigde dat het gebruik van alleen $r_{ee}$ leidt tot niet-Markoviaanse dynamica. Echter, de opname van de drie uit TICA afgeleide onafhankelijke componenten (IC1, IC2, IC3) naast $r_{ee}$ maakte de kinetiek succesvol Markoviaans, waarmee de multidimensionale benadering werd gevalideerd.
• Ontvouwingstrajecten:
  • Thermische ontvouwing ($T=298$ K): Het peptide ontvouwt voornamelijk beginnend bij het N-terminus op een "rits-achtige" manier. Het proces vertoont heterogeniteit, inclusief misgevouwen structuren, en kan niet worden gekarakteriseerd als een eenvoudig coöperatief proces.
  • Mechanische ontvouwing ($T=240$ K): Onder externe kracht verloopt het systeem via een metastabiel intermediair stadium waarbij de innerste H-brug intact blijft terwijl de buitenste lussen open zijn. Dit intermediair is duidelijker aanwezig in de PMF verkregen uit Umbrella Sampling dan in standaard MD-simulaties.
• Kracht-Extensiecurves:
  • De MSM-benadering slaagde erin de kwalitatieve kenmerken van FEC's (specifiek het optreden van twee krachtpieken) te reproduceren die in FPMD-simulaties werden waargenomen voor trage trekkingssnelheden.
  • Snelheidsafhankelijkheid: Voor zeer trage trekkingssnelheden ($v \le 10^{-4}$ m/s) bereikt het systeem een quasi-evenwichtstoestand waarbij de piek van het intermediaire stadium afvlakt, consistent met evenwichtsgedrag waarbij ontvouwing begint bij het N-terminus zonder significante accumulatie van intermediairen.
  • Beperkingen: Bij hoge trekkingssnelheden onderschat de MSM-benadering de gemiddelde kracht niet, maar overschat deze deze. Deze discrepantie ontstaat omdat het wisselen tussen ensembles sneller plaatsvindt dan de interne overgangssnelheden tussen toestanden, en de aanname van Markoviaanse kinetiek bezwijkt als de tijd voor het wisselen van ensembles korter is dan de MSM-lag-tijd.

Betekenis en claims
Het artikel claimt succesvol een dynamische grofkorreligingstechniek, die eerder was gevalideerd op eenvoudige tweestapsystemen, uit te breiden naar een complex peptidesysteem dat een meerstaps ontvouwingstraject vertoont. De primaire bijdrage is het aantonen dat door het selecteren van geschikte collectieve variabelen (H-brugafstanden) en het uitvoeren van dimensiereductie (TICA), men een geldig Markov-toestandsmodel kan construeren, zelfs wanneer de eind-tot-eind-afstand een ontoereikende ordeparameter is.

De auteurs stellen dat deze methode toelaat om observabelen te berekenen in het regime van trage trekking, wat computationeel ontoegankelijk is voor directe FPMD-simulaties. Het onderzoek benadrukt dat hoewel de methode robuust is voor matige tot kleine trekkingssnelheden, deze beperkingen ondervindt bij zeer hoge snelheden door het bezwijken van de Markoviaanse aanname ten opzichte van de snelheid van protocolwisseling. Uiteindelijk biedt het werk een raamwerk voor het nabootsen van mechanische ontvouwingsprocessen in systemen waar eenvoudige ordeparameters falen, en overbrugt het de kloof tussen simulatietijdschalen en experimentele krachtspectroscopie.