$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Rysunek 1 przedstawia widmo absorpcyjne suchego filmu bR osadzonego na powierzchni diamentowego elementu wewnętrznego odbicia stosowanego w spektroskopii ATR. Wyraźnie widoczne są charakterystyczne pasma drgań wiązania peptydowego (amid A, amid I i amid II). Przybliżoną grubość suchego filmu można oszacować na ~1 μm, biorąc pod uwagę ilość dodanego białka (18 μg) oraz powierzchnię ATR (~0.2 cm2), a także przyjmując gęstość białka na poziomie ~1.4 g/cm3,58, lipidów 1.0 g/cm3,59 oraz stosunek lipidów do białek 1/3 (w/w) w błonie purpurowej60.
Suchą błonę rehydratowano nadmiarem 4 M NaCl, 100 mM NaPi o pH 7.4 (Rysunek 3). Poziom hydratacji oraz efektywne stężenie białka w objętości sondowanej przez falę ewanescentną można wyprowadzić z czynnika skalującego niezbędnego do cyfrowego usunięcia pasm absorpcyjnych wody. Optymalny czynnik skalujący wyniósł 0.87, co oznacza, że w tym przypadku bufor zajmuje 87%, a próbka 13% objętości w pobliżu powierzchni. Biorąc pod uwagę gęstość białek i lipidów oraz stosunek lipidów do białek w purpurowych błonach (patrz wyżej), możemy wyznaczyć efektywne stężenie białka wynoszące 125 mg/ml w objętości sondowanej przez falę ewanescentną. Na podstawie poziomu hydratacji można wywnioskować, że w tym konkretnym przypadku grubość błony próbki zwiększyła się ~6 razy podczas rehydratacji, z ~1 do ~6 μm. Dla kąta padania 45°, typowego dla naszej i większości konfiguracji ATR, głębokość penetracji pola ewanescentnego (dp) dla hydratowanej błony białkowej w zakresie 1,800-850 cm-1 waha się od 0.3 do 0.6 μm 61. Ponieważ pole ewanescentne jest niemal niewrażliwe na próbkę znajdującą się na wysokości dwukrotności dp 61, wnioskujemy, że ilość białka użyta w eksperymencie ATR mogłaby w zasadzie zostać zredukowana 5x bez żadnej znaczącej utraty sygnału.
Przy użyciu buforów o niższej sile jonowej film bardziej pęcznieje, a ilość białka w objętości sondowanej przez pole ewanescentne ulega zmniejszeniu (Rysunek 2). Dokładna zależność między pęcznieniem filmu a siłą jonową buforu będzie zależeć, między innymi, od rodzaju lipidów. Na przykład podobne pęcznienie filmu, jakie uzyskano tutaj dla bR w błonie purpurowej przy użyciu 4 M NaCl, osiągnięto dla białka błonowego zrekonstytuowanego w polarnych lipidach E. coli przy użyciu zaledwie 0,1 M NaCl 62. Jeśli próbka zajmuje mniej niż 5% objętości sondowanej, należy rozważyć powtórzenie etapu hydratacji z użyciem buforu o wyższej sile jonowej. Pęcznienie filmu po hydratacji wymaga pewnego czasu do osiągnięcia stabilizacji (Rysunek 4). Dla bR w błonie purpurowej proces ten jest monoeksponencjalny, ze stałą czasową wynoszącą 12 min: uzyskanie stabilnego rehydratowanego filmu zajmuje nie więcej niż 30-60 min
Rysunek 5 przedstawia widmo absorpcyjne IR uwodnionego filmu ChR2 uzyskane metodą transmisji. Uwodnienie osiągnięto poprzez wystawienie suchego filmu na działanie atmosfery o kontrolowanej wilgotności, zapewnionej przez mieszaninę glicerolu i wody. Widmo można rozłożyć na wkłady pochodzące z wody, detergentu i białka. Ilość każdego z tych składników oszacowano przy użyciu widm współczynnika ekstynkcji, przeskalowanych tak, aby pasowały do widma eksperymentalnego. Dla wody przyjęto widmo współczynnika ekstynkcji z literatury63, a dla DM zmierzono je z roztworu o stężeniu 100 mg/ml (Rysunek 6). Wykorzystano również współczynniki ekstynkcji dla dwóch reprezentatywnych białek błonowych: mitochondrialnego transportera ADP/ATP64 oraz bR (Rysunek 7). Współczynnik skalowania wskazuje na gęstość powierzchniową masy wody i DM w filmie wynoszącą odpowiednio 260 μg/cm2 i 200 μg/cm2. Pozostała absorpcja (Rysunek 5, czerwona linia) pochodzi głównie z białka, którego gęstość powierzchniową oszacowano na 250 μg/cm2, wykorzystując współczynnik ekstynkcji pasma amidu II dla dwóch różnych białek błonowych (patrz Rysunek 6). Stosunek molowy białko/detergent/woda obliczono na 1/60/2000, przyjmując odpowiednio masy cząsteczkowe 35 kDa, 480 Da i 18 Da.
Wykres 3D z typowego eksperymentu czasowo-rozdzielczego FT-IR z metodą step-scan na bR przedstawiono na Rysunku 10A; dane uzyskano metodą ATR w czasie akwizycji wynoszącym 200 min. Widma można wyodrębnić dla konkretnych momentów czasowych, na przykład gdy produkty pośrednie L, M i N fotocyklu bR osiągają najwyższą populację (Rysunek 11). Ich cechy spektralne zostały szeroko opisane13,41,65 i nie będą tutaj dalej omawiane. Rysunek 12 przedstawia przebiegi czasowe dla wybranych liczb falowych. Mianowicie, wzrost kinetyki przy 1,762 cm-1 (t1/2 ~60 μsec) informuje o dynamice protonowania Asp85 z zasady Schiffa retinalu66, a jej spadek do zera wskazuje na deprotonowanie podczas powrotu do stanu podstawowego67. Ujemny wzrost kinetyki przy 1,740 cm-1 (t1/2 ~1 msec) informuje o deprotonowaniu łańcucha bocznego Asp96, który jest donorem protonu dla zasady Schiffa68. Spadek intensywności do zera informuje o jego reprotonowaniu z cytoplazmy67. Dynamikę zmian konformacyjnych białka można badać poprzez zmiany absorpcji w obszarze amidu I i amidu II, które osiągają maksymalną zmianę przy ~3 msec w temperaturze pokojowej67,69.
Zastosowanie SVD w problemach spektroskopowych było już wcześniej omawiane55,56. W skrócie, SVD rozkłada dane eksperymentalne, uporządkowane w macierzy A, jako: A = U S VT. Rozkład ten można zmodyfikować, aby uwzględnić zależność szumu od liczby fal oraz aby zminimalizować wpływ fluktuacji/dryftów linii bazowej57. Kolumny macierzy U i V zawierają odpowiednio ortonormalne widma oraz wektory przebiegów czasowych dla każdego komponentu SVD. S jest macierzą diagonalną zawierającą tzw. wartości osobliwe. (Abstrakcyjne) komponenty spektro-temporalne w macierzach U i V pojawiają się w kolejności malejącej istotności dla opisu danych eksperymentalnych w sensie najmniejszych kwadratów, co jest określane przez przypisaną im wartość osobliwą. Rycina 13A przedstawia wartości osobliwe jako funkcję numeru komponentu oraz reprodukuje pierwsze osiem kolumn/komponentów macierzy U (abstrakcyjne widma, Rycina 13B) i V (abstrakcyjne przebiegi czasowe, Rycina 13C). Sygnał jest skoncentrowany w pierwszych pięciu komponentach (zgodnie z oczekiwaniami dla fotocyklu zawierającego pięć stanów pośrednich), podczas gdy komponenty powyżej są w dużej mierze zdominowane przez szum losowy i inne źródła błędów. Dane eksperymentalne zostały zrekonstruowane przy użyciu tylko pierwszych pięciu kolumn macierzy U i V oraz pierwszych pięciu kolumn i wierszy macierzy S. Dane zrekonstruowane metodą SVD stanowią najlepszą aproksymację danych eksperymentalnych w sensie najmniejszych kwadratów do macierzy o randzie pięć. Zrekonstruowane dane wykazują poprawioną jakość (Rycina 10B). A konkretnie, znacząco zredukowano szum, a także niektóre ledwo zauważalne fluktuacje i dryfty linii bazowej (częste w spektroskopii step-scan z rozdzielczością czasową 25). Przetwarzanie SVD jest szczególnie korzystne dla poprawy jakości eksperymentalnych przebiegów czasowych w zakresie milisekund (czerwone linie na Rycini 12).
Dane z pomiarów step-scan z rozdzielczością czasową dla fotocyklu ChR2 uzyskano zgodnie z opisanym powyżej protokołem dla eksperymentów transmisyjnych (Rysunek 14A), co obejmowało łącznie około 120 h pomiarów. Zakres czasowy danych zebranych metodą step-scan wynosi od 6,25 µsec do 125 msec. Fotocykl ChR2 wymaga około 60 s do pełnego powrotu do stanu wyjściowego. Końcową fazę fotocyklu można objąć pomiarami rapid-scan FT-IR, a następnie połączyć zbiory danych step-scan i rapid-scan, zgodnie z opisem przedstawionym w innym opracowaniu48. Zmiany spektralne ChR2 są około 10-krotnie mniejsze niż te uzyskiwane dla bR, co czyni pomiary bardziej wymagającymi. W szczególności przy tak niskich zmianach absorbancji wyraźnie widoczne stają się oscylacje linii bazowej w zakresie milisekundowym. Wynikają one z niewielkich drgań ruchomego lustra w skali czasowej milisekundowej25. Oscylacje te, wraz z częścią szumu, można w dużej mierze usunąć za pomocą SVD, co znacząco poprawia wygląd danych (Rysunek 14B).

Rysunek 1. Widmo absorpcyjne bR w purpurowej membranie osuszonej na diamentowej powierzchni akcesorium ATR. Oznaczono pasma drgań wiązania peptydowego (amid I, II i A).

Rycina 2. Widmo absorpcyjne błony bR po rehydratacji przy użyciu buforów o różnych mocnościach jonowych (rozcieńczenia 4 M NaCl, 100 mM Na2HPO4/NaH2PO4 przy pH 7.4). Zwróć uwagę na spadek pasma amidu II, t.j. zwiększone pęcznienie błony, wraz ze spadkiem mocności jonowej buforu. (Wstawka) Względna ilość białka w pobliżu powierzchni, określona na podstawie intensywności absorbancji przy 1,541 cm-1 (maksimum amidu II) po odjęciu wkładu absorpcyjnego buforu.

Rycina 3. Widmo absorpcyjne filmu bR rehydratyzowanego buforem objętościowym (siła jonowa 4,3 M) po odjęciu wkładu buforu. Współczynnik odejmowania, 0,87, wybrano w celu usunięcia silnej absorpcji wody w zakresie 3700-3000 cm-1 oraz uzyskania płaskiej linii bazowej w zakresie 2700-1800 cm-1.

Rycina 4. Widmo absorpcyjne bR po rehydratacji buforem o sile jonowej 4,3 M (dla przejrzystości odjęto absorpcję buforu, analogicznie jak na Rycini 3). Wstawka przedstawia przebieg pęcznienia filmu po rehydratacji, a następnie absorbancję amidu II białka przy 1 541 cm-1. Dopasowanie do pojedynczej funkcji wykładniczej wskazuje, że stała czasowa stabilizacji pęcznienia filmu wynosi 12 min.

Rysunek 5. Widmo absorpcyjne uwodnionej warstwy ChR2 zmierzone metodą transmisji (linia niebieska). Widmo współczynnika ekstynkcji wody (przerywana linia pomarańczowa) i DM (przerywana linia cyjanowa) zostało przeskalowane i odjęte (linia czerwona).

Rysunek 6. Odtworzone widma współczynnika ekstynkcji masowej w 25 °C dla wody w stanie ciekłym (http://www.ualberta.ca/~jbertie/JBDownload.HTM#Spectra) oraz dla uwodnionego filmu nośnika ADP/ATP (AAC)64. Widmo współczynnika ekstynkcji masowej zmierzono w roztworze dla DM (100 mg/ml) oraz w uwodnionym filmie dla bR.

Rycina 7. Transmitancja filtra optycznego zastosowanego w pomiarach FT-IR w trybie scan-scan z rozdzielczością czasową.

Rycina 8. Charakterystyka impulsów laserowych generowanych przez drugą harmoniczną (532 nm) lasera Nd:YAG. Histogram zmienności energii względnej 1 000 impulsów laserowych oraz dopasowanie do rozkładu Gaussa z odchyleniem standardowym 0,05.

Rycina 9. Zmiany absorbancji bR indukowane światłem przy trzech reprezentatywnych liczbach falowych w jednolitych odstępach czasu (linie zielona, czarna i niebieska) oraz po uśrednianiu kwazi-logarytmicznym do ~20 punktów/dekady (linie czerwone). Zauważalna jest redukcja szumów po uśrednianiu logarytmicznym, co ujawnia oscylacje przebiegów czasowych w zakresie milisekund.

Rysunek 10. Trójwymiarowa reprezentacja zmian absorbancji indukowanych światłem dla bR zarejestrowanych za pomocą ATR przy pH 7.4 (4 M NaCl, 100 mM NaPi). A) Dane surowe. B) Dane zrekonstruowane z pięciu komponentów SVD. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 11. Spektra różnicowego pochłaniania FT-IR fotocyklu bR w trzech wybranych punktach czasowych, w których najwyższe stężenie wykazują produkty pośrednie L (12.5 µsec), M (300 µsec) oraz N (6 msec).

Rysunek 12. Transfer protonów i dynamika szkieletu białkowego w fotocyklu bR rozstrzygnięta za pomocą spektroskopii różnicowej FT-IR w trybie step-scan z rozdzielczością czasową. Zmiany absorbancji przy 1,762 cm-1 informują o dynamice protonowania/deprotonowania Asp85, a przy 1,741 cm-1 o dynamice deprotonowania/reprotonowania Asp96 (w tym o zmianach wiązań wodorowych przed 300 μsec). Przebiegi czasowe przy 1,670 i 1,555 cm-1 informują o zmianach w drganiach amidowych I i II, które są wrażliwe na konformację szkieletu peptydowego. Linie czarne odpowiadają surowym danym, a czerwone danym przetworzonym metodą SVD. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rysunek 13. Rozkład według wartości osobliwych (SVD) eksperymentalnych danych FT-IR z rozdzielczością czasową i skanowaniem skokowym dla fotocyklu bR (patrz Rysunek 10A). SVD wykonano po ważeniu danych eksperymentalnych z uwzględnieniem zależności odchylenia standardowego szumu od liczby falowej (Rysunek 15); w połączeniu z 1szą pochodną w celu zmniejszenia wagi statystycznej fluktuacji linii bazowej, zgodnie z wcześniejszym opisem57. A) Wykres względnych wartości osobliwych dla pierwszych 50 komponentów (czarne kółka). Pierwsze pięć komponentów przypisano do komponentów sygnału (czerwone kółka). Pozostałe komponenty zanikają wykładniczo, co jest oczekiwane dla komponentów szumu (patrz przerywana szara linia). B) Pierwsze osiem widm abstrakcyjnych (U1 do U8). C) Pierwsze osiem abstrakcyjnych przebiegów czasowych (V1 do V8). Widma abstrakcyjne i przebiegi czasowe przypisane do sygnału przedstawiono czerwonymi liniami, a pozostałe czarnymi liniami. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rycina 14. Trójwymiarowa reprezentacja zmian absorbancji IR indukowanych światłem dla ChR2, zarejestrowanych metodą step-scan w trybie transmisji. Dane step-scan obejmują czas do 125 msec, co pokrywa jedynie część fotocyklu. A) Dane surowe. B) Dane zrekonstruowane przy użyciu pięciu komponentów SVD. Kliknij tutaj, aby wyświetlić większą wersję tej ryciny.

Rycina 15. Szacowany poziom szumu w spektrometrze FT-IR. Widmo różnicowe absorpcji przy rozdzielczości czasowej 6,25 μsec i spektralnej 8 cm-1 dla eksperymentu obejmującego 1 koaddycję/pozycję lustra (500 fotoreakcji) lub ~200 koaddycji/pozycji lustra (105 fotoreakcji). Wartości przedstawiono dla układu osłabionego całkowitego odbicia (ATR) oraz układu transmisyjnego.