Experimentally accessible measurement of irreversibility in stochastic systems by categorizing single-molecule displacements

Technische Samenvatting: Experimenteel Toegankelijke Meting van Irreversibiliteit in Stochastische Systemen

Probleemstelling
Het kwantificeren van de irreversibiliteit en dissipatie van niet-evenwichtsprocessen is fundamenteel voor het begrijpen van het gedrag, de capaciteiten en de efficiëntie van stochastische systemen, met name in biologische contexten zoals eiwitvouwing. Hoewel de stochastische thermodynamica een rigoureus theoretisch kader biedt (bijv. fluctuatietheorema's, Kullback-Leibler-divergenties), is het toepassen van deze concepten op experimentele gegevens uitdagend. Bestaande methoden vereisen vaak:

1. Volledige kennis van de krachten die op het systeem inwerken en de diffusiecoëfficiënt.
2. Toegang tot volledige trajectgegevens, wat vaak onmogelijk is vanwege de eindige temporele resolutie.
3. Specifieke condities (bijv. stationariteit) of complexe gegevensverwerking (bijv. hyperparameterinstellingen voor machine learning).
4. Zowel voorwaartse als achterwaartse protocollen voor bepaalde schatters (bijv. de relatie van Crooks).

Bijgevolg is er behoefte aan een model-vrije, experimenteel toegankelijke methode om irreversibiliteit te meten in systemen met beperkte informatie, zoals bij single-molecule tracking waarbij krachten en diffusiecoëfficiënten onbekend of moeilijk te schatten zijn.

Methodologie
De auteurs introduceren een techniek om irreversibiliteit in overdempte Langevin-systemen te kwantificeren op basis van de categorisering van single-molecule verplaatsingen over eindige tijdsintervallen ($\Delta t$).

1. Categorisering van Verplaatsing: In plaats van volledige trajecten te analyseren, observeert de methode de initiële positie ($x'$) en de finale positie ($x$) van een deeltje na een tijd $\Delta t$, waarbij de verplaatsing wordt gedefinieerd als $\ell = x - x'$. Deze verplaatsingen worden gecategoriseerd in klassen ($C$) op basis van hun ruimtelijke locatie ten opzichte van gedefinieerde grenzen (bijv. het kruisen van een grens $k$ van links naar rechts, van rechts naar links, of het binnenblijven in de bulk).
2. Irreversibiliteitsmaat ($\Sigma_{\Delta t}$): De auteurs definiëren een maat voor irreversibiliteit als de Kullback-Leibler (KL) divergentie tussen de gezamenlijke waarschijnlijkheidsverdelingen van voorwaartse verplaatsingen en hun tijd-omgekeerde tegenhangers over deze klassen:
   $$\Sigma_{\Delta t}(C_n) \equiv k_B \sum_{C \in C_n} \int d\ell \, q_{\Delta t}(\ell, C) \ln \frac{q_{\Delta t}(\ell, C)}{q_{\Delta t}(-\ell, \tilde{C})}$$
   waarbij $q_{\Delta t}(\ell, C)$ de waarschijnlijkheidsdichtheid is van een verplaatsing $\ell$ behorend tot klasse $C$, en $\tilde{C}$ de tijd-omgekeerde klasse is.
3. Theoretische Grenzen:
   • Stationaire Systemen: Voor stationaire dynamica vormt $\Sigma_{\Delta t}$ een rigoureuze ondergrens voor de gemiddelde entropieproductie $\langle \Delta S_{tot} \rangle$. Deze grens wordt nauwer naarmate het aantal verplaatsingsklassen toeneemt (een hiërarchie van grenzen).
   • Tijdsafhankelijke Systemen: Voor niet-stationaire protocollen biedt de maat een ondergrens in de limiet van een verdwijnende $\Delta t$.
   • Lokale Dissipatie: Door zich te richten op verplaatsingen die een specifieke grens kruisen, kan de methode de lokale dissipatiesnelheid schatten.
4. Experimentele Implementatie: De methode wordt toegepast op single-molecule magnetische tweepers experimenten op het eiwit taline R3 (wild-type en mutant). Het systeem wordt onderworpen aan cyclische lineaire krachtrampen (ontvouwing en hervouwing). De irreversibiliteit wordt berekend door verplaatsingen in vier klassen in te delen ten opzichte van een grens tussen de gevouwen en ontvouwen toestanden.

Kernbijdragen en Resultaten

• Model-Vrije Schatting: De methode vereist geen voorafgaande kennis van de deterministische krachten, de diffusiecoëfficiënt of het onderliggende energielandschap. Het steunt uitsluitend op de statistieken van geobserveerde verplaatsingen.
• Validatie op Gesimuleerde Systemen:
  • Bij een gedreven deeltje op een ring met een vierkante potentiaal slaagde de methode erin de entropieproductie te vangen. De auteurs toonden aan dat het verhogen van het aantal verplaatsingsklassen (het fijner partitioneren van de ruimte) de schatting aanzienlijk verbeterde, vooral voor systemen nabij evenwicht waar de Thermodynamic Uncertainty Relation (TUR) alleen een losse grens bood.
  • De methode toonde aan accuraat het ruimtelijke profiel van lokale dissipatie in kaart te brengen.
• Toepassing op Eiwitvouwing (Talin R3):
  • Wild-type vs. Mutant: De methode maakte onderscheid tussen het wild-type taline R3 en een mechanisch stabielere mutant (R3IVVI). Het identificeerde correct dat de mutant een hogere totale irreversibiliteit vertoont.
  • Asymmetrie in Vouwingsdynamica: Een belangrijke bevinding was de detectie van asymmetrie in de dynamica van de mutant: ontvouwing bleek systematisch irreversibeler dan hervouwing. Dit werd niet waargenomen bij het wild-type.
  • Landschapsinferentie: Door de experimentele resultaten te vergelijken met simulaties van asymmetrische dubbel-well potentialen, koppelden de auteurs de geobserveerde asymmetrie in irreversibiliteit aan het onderliggende vrije-energielandschap. Het gedrag van de mutant suggereerde een asymmetrisch landschap waarbij de ontvouwen well breder en verder van de barrière verwijderd is dan de gevouwen well.
• Robuustheid: De methode bleek robuust tegen meetruis (gesimuleerde Gaussische ruis) en beperkte datasets, waarbij het in staat bleef systemen te rangschikken op basis van dissipatie en het identificeren van het tijdstip van irreversibele gebeurtenissen (pieken in de irreversibiliteitsgraad die overeenkomen met transitiegebeurtenissen).

Betekenis en Claims
Het artikel beweert een "rigoureus instrument te bieden voor het rationaliseren van de functie van mechanisch actieve eiwitten" door een kwantitatief kader te bieden voor niet-evenwichtsthermodynamica in biomoleculaire transities.

• Experimentele Toegankelijkheid: De primaire betekenis ligt in de haalbaarheid van de methode voor experimentele systemen waarbij volledige trajectreconstructie of kennis van krachten niet beschikbaar is. Het omzeilt de noodzaak om de exponentieel grote ruimte van trajecten te bemonsteren door zich te concentreren op grofmazige initiële en finale toestanden.
• Gevoeligheid voor Landschapskenmerken: De auteurs demonstreren dat het kwantificeren van irreversibiliteit subtiele eigenschappen van het energielandschap kan onthullen (bijv. asymmetrieën in barrièreafstanden) die niet gemakkelijk toegankelijk zijn via standaard evenwichtsmetingen of de oppervlakte van hysteresislussen alleen.
• Vergelijking met Bestaande Methoden: De auteurs positioneren hun methode als superieur aan de Thermodynamic Uncertainty Relation (TUR) voor ruimtelijk heterogene systemen, en robuuster/interpreteerbaarder dan machine learning-benaderingen of methoden die zowel voorwaartse als achterwaartse protocollen vereisen (zoals de relatie van Crooks of specifieke binning-methoden). Ze merken op dat hoewel hun maat een ondergrens is, het de relatieve trends en verschillen tussen systemen (bijv. wild-type versus mutant, verschillende trekknelheden) met hoge getrouwheid vastlegt.

Het artikel concludeert dat deze benadering onderzoekers in staat stelt om te beoordelen hoe ver een systeem zich van het evenwicht bevindt en om trade-offs tussen dissipatie en prestaties te karakteriseren, waarbij specifiek de bruikbaarheid wordt benadrukt voor het ontleden van de thermodynamische signaturen van eiwitvouwingsmechanismen.