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 wibracji 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), przy założeniu gęstości białka wynoszącej ~1.4 g/cm3,58 lipidów 1.0 g/cm3,59 oraz stosunku lipidów do białek 1/3 (w/w) w błonie purpurowej60.
Suchy film nawodniono nadmiarem 4 M NaCl, 100 mM NaPi o pH 7,4 (Rysunek 3). Poziom nawodnienia oraz efektywne stężenie białka w objętości badanej przez falę ewanescentną można wyznaczyć na podstawie współczynnika skalowania niezbędnego do cyfrowego usunięcia pasm absorpcyjnych wody. Optymalny współczynnik skalowania 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 badanej przez falę ewanescentną. Na podstawie poziomu nawodnienia można wywnioskować, że w tym konkretnym przypadku grubość filmu próbki zwiększyła się ok. 6-krotnie po rehydratacji, z ~1 do ~6 μm. Dla kąta padania 45°, typowego dla naszej i większości konfiguracji ATR, głębokość przenikania pola ewanescentnego (dp) dla nawodnionego filmu białkowego 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ę dwukrotnie powyżej dp 61, wnioskujemy, że ilość białka wykorzystanego w eksperymencie ATR mogłaby w zasadzie zostać zredukowana 5-krotnie bez żadnej istotnej utraty sygnału.
W przypadku zastosowania buforów o niższej sile jonowej film ulega większemu pęcznieniu, co prowadzi do zmniejszenia ilości białka w objętości badanej przez pole ewanescentne (Rycina 2). Dokładna zależność między pęcznieniem filmu a siłą jonową buforu zależy, między innymi, od rodzaju lipidów. Na przykład, podobne pęcznienie filmu, jakie uzyskano w tym przypadku dla bR w błonie purpurowej przy użyciu 4 M NaCl, zaobserwowano dla białka błonowego zrekonstytuowanego w polarnych lipidach E. coli przy zastosowaniu zaledwie 0,1 M NaCl 62. Jeśli próbka zajmuje mniej niż 5% badanej objętości, 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 w celu osiągnięcia stabilizacji (Rycina 4). W przypadku bR w błonie purpurowej proces ten jest monoeksponencjalny, z stałą czasową wynoszącą 12 min: uzyskanie stabilnego rehydratowanego filmu zajmuje nie więcej niż 30-60 min
Rysunek 5 przedstawia widmo absorpcji IR uwodnionego filmu ChR2 uzyskane metodą transmisji. W tym przypadku uwodnienie osiągnięto poprzez wystawienie suchego filmu na atmosferę o kontrolowanej wilgotności, zapewnioną 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 mogliśmy oszacować, korzystając z widm współczynnika ekstynkcji, przeskalowanych tak, aby pasowały do widma eksperymentalnego. W przypadku wody przyjęliśmy widmo współczynnika ekstynkcji z literatury63, a dla DM zmierzyliśmy go z roztworu o stężeniu 100 mg/ml (Rysunek 6). Wykorzystaliśmy również współczynniki ekstynkcji dla dwóch reprezentatywnych białek błonowych: mitochondrialnego przenośnika ADP/ATP64 oraz bR (Rysunek 7). Współczynnik skalowania wskazuje na powierzchniową gęstość 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łka/detergentu/wody obliczono na 1/60/2000, przyjmując odpowiednio masy cząsteczkowe 35 kDa, 480 Da i 18 Da.
Wykres 3D z typowego eksperymentu FT-IR w trybie step-scan z rozdzielczością czasową dla bR przedstawiono na Rysunku 10A; dane uzyskano metodą ATR podczas 200 min akwizycji. Widma można wyodrębnić dla konkretnych momentów czasowych, na przykład w czasie, gdy spodziewane jest osiągnięcie najwyższej populacji przez produkty pośrednie L, M i N fotocyklu bR (Rysunek 11). Ich cechy widmowe 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 jego zanik 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, będącego donorem protonów dla zasady Schiffa68. Zanik 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że zostać zmodyfikowany, aby uwzględnić zależność szumu od liczby fal oraz aby zniwelować fluktuacje/dryfty 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. (Abstraktowe) komponenty spektro-temporalne w U i V pojawiają się w kolejności malejącej istotności dla opisu danych eksperymentalnych w sensie najmniejszych kwadratów, co jest określone przez przypisaną im wartość osobliwą. Rycina 13A przedstawia wartości osobliwe jako funkcję numeru komponentu oraz reprodukuje pierwsze osiem kolumn/komponentów macierzy U (widma abstraktowe, Rycina 13B) oraz V (przebiegi czasowe abstraktowe, Rycina 13C). Sygnał jest skoncentrowany w pierwszych pięciu komponentach (co jest oczekiwane dla fotocyklu zawierającego pięć stanów pośrednich), natomiast komponenty powyżej są w dużej mierze zdominowane przez szum losowy i inne źródła błędów. Dane eksperymentalne zrekonstruowano, używając jedynie pierwszych pięciu kolumn U i V oraz pierwszych pięciu kolumn i wierszy macierzy S. Dane zrekonstruowane za pomocą SVD stanowią najlepszą aproksymację metodą najmniejszych kwadratów danych eksperymentalnych do macierzy o randze pięciu. Zrekonstruowane dane wykazują poprawioną jakość (Rycina 10B). Mianowicie, znacznie zredukowano szum, a także niektóre ledwo zauważalne fluktuacje i dryfty linii bazowej (częste w spektroskopii z rozdzielczością czasową typu step-scan 25). Przetwarzanie SVD jest szczególnie korzystne dla poprawy jakości eksperymentalnych przebiegów czasowych w zakresie milisekundowym (czerwone linie na Rycynie 12).
Dane czasowo-rozdzielcze ze skanowania skokowego (step-scan) dla fotocyklu ChR2 uzyskano, stosując opisany powyżej protokół dla eksperymentów transmisyjnych (Rysunek 14A), co wiązało się z około 120 h skumulowanych pomiarów. Dane czasowo-rozdzielcze zebrane metodą step-scan obejmują zakres od 6,25 µsec do 125 msec. Fotocykl ChR2 wymaga około 60 s do pełnego powrotu do stanu wyjściowego. Końcowa część fotocyklu może zostać objęta analizą za pomocą szybkiego skanowania FT-IR (rapid-scan), a zestawy danych step-scan i rapid-scan mogą zostać połączone, jak przedstawiono w innym opracowaniu48. Zmiany spektralne ChR2 są około 10-krotnie mniejsze od tych uzyskanych dla bR, co utrudnia pomiary. W szczególności, przy tak niskich zmianach absorpcji, wyraźne stają się oscylacje linii bazowej w zakresie milisekundowym. Wynikają one z niewielkich drgań ruchomego lustra w skali milisekundowej25. Oscylacje te, wraz z częścią szumu, mogą zostać w dużej mierze usunięte za pomocą SVD, co znacząco poprawia wygląd danych (Rysunek 14B).

Rycina 1. Widmo absorpcyjne bR w błonie purpurowej wysuszonej na powierzchni diamentowej akcesorium ATR. Oznaczono pasma drgań wiązania peptydowego (amid I, II i A).

Rycina 2. Widmo absorpcyjne filmu bR po rehydratacji przy użyciu buforów o różnych mocnościach jonowych (rozcieńczenia 4 M NaCl, 100 mM Na)2HPO4/NaH2PO4 przy pH 7,4). Zauważenie zmniejszenia intensywności pasma amidowego II, t.j., zwiększone pęcznienie błony wraz ze spadkiem siły 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 pasma amidowego II) po odjęciu wkładu absorpcyjnego buforu.

Rycina 3. Widmo absorpcyjne błony bR rehydratowanej w buforze 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 3 700-3 000 cm-1 oraz uzyskania płaskiej linii bazowej w zakresie 2 700-1 800 cm-1.

Rysunek 4. Widmo absorpcyjne bR po rehydratacji w buforze o sile jonowej 4,3 M (dla przejrzystości odjęto absorpcję bufora, analogicznie jak na Rysunku 3). Wstawka przedstawia ewolucję pęcznienia filmu po rehydratacji, a następnie absorbancję amidu II białka przy 1 541 cm-1. Dopasowanie do funkcji pojedynczej 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 wymierania wody (przerywana linia pomarańczowa) i DM (przerywana linia cyjanowa) zostało przeskalowane i odjęte (linia czerwona).

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

Rysunek 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 względnej energii 1 000 impulsów laserowych oraz dopasowanie do rozkładu gaussowskiego z odchyleniem standardowym 0,05.

Rycina 9. Indukowane światłem zmiany absorbancji bR przy trzech reprezentatywnych liczbach falowych w równych odstępach czasu (linie zielona, czarna i niebieska) oraz po uśrednianiu quasi-logarytmicznym do ~20 punktów na dekadę (linie czerwone). Zwróć uwagę na redukcję 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, zarejestrowana metodą ATR przy pH 7,4 (4 M NaCl, 100 mM NaPi). A) Surowe dane. B) Dane zrekonstruowane przy użyciu pięciu komponentów SVD. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 11. Widma różnicowej absorpcji FT-IR fotocyklu bR w trzech wybranych momentach czasu, w których największe wzbogacenie wykazują produkty pośrednie L (12,5 µsec), M (300 µsec) i N (6 msec).

Rysunek 12. Transfer protonów i dynamika szkieletu białkowego w fotocyklu bR rozdzielone za pomocą czasowo rozdzielonej różnicowej spektroskopii FT-IR z przesuwaniem krokowym (step-scan). Zmiany absorbancji przy 1 762 cm-1 informują o dynamice protonowania/deprotonowania Asp85, a przy 1 741 cm-1 o dynamice deprotonowania/reprotonowania Asp96 (wliczając zmiany wiązań wodorowych przed 300 μsec). Przebiegi czasowe przy 1 670 i 1 555 cm-1 informują o zmianach w wibracjach amidowych I i II, które obie są wrażliwe na konformację szkieletu peptydowego. Czarne linie odpowiadają surowym danym, a czerwone danym przetworzonym metodą SVD. Proszę kliknąć tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 13. Dekompozycja według wartości osobliwych (SVD) eksperymentalnych danych FT-IR z rozdzielczością czasową, uzyskanych metodą skanowania skokowego dla fotocyklu bR (patrz Rysunek 10A). SVD wykonano po zważeniu danych eksperymentalnych, uwzględniając zależność odchylenia standardowego szumu od liczby falowej (Rysunek 15); połączone z 1Proszę podać tekst źródłowy, który należy przetłumaczyć. 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 wykazują spadek wykładniczy, zgodnie z oczekiwaniami 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 liniami czerwonymi, a pozostałe liniami czarnymi. Aby wyświetlić powiększoną wersję tej figury, kliknij tutaj.

Rysunek 14. Trójwymiarowa reprezentacja zmian absorbancji IR indukowanych światłem dla ChR2, zarejestrowana metodą skanowania skokowego w trybie transmisji. Dane ze skanowania skokowego 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ć powiększoną wersję tego rysunku.

Rysunek 15. Szacowany poziom szumu w FT-IR. Widmo różnicowe absorpcji przy rozdzielczości czasowej 6,25 μsec i rozdzielczości widmowej 8 cm-1 dla eksperymentu obejmującego 1 koaddycję/pozycję lustra (500 fotoreakcji) lub ~200 koaddycji/pozycję lustra (105 fotoreakcji). Wartości przedstawiono dla konfiguracji z tłumionym całkowitym odbiciem (ATR) oraz transmisji.