$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Rysunek 1 pokazuje widmo absorpcyjne suchej warstwy bR osadzonej na powierzchni diamentowego elementu wewnętrznego odbicia używanego do spektroskopii ATR. Charakterystyczne pasma od drgań wiązania peptydowego (amid A, amid I i amid II) są wyraźnie rozróżnialne. Przybliżoną grubość suchej warstwy można oszacować na ~1 μm, biorąc pod uwagę ilość dodanego białka (18 μg) i powierzchnię ATR (~0,2 cm2) oraz przyjmując gęstość białka jako ~1,4 g/cm 3,58 i lipidów jako 1,0 g/cm 3,59 i stosunek lipidów do białka 1/3 (w/w) w fioletowej błonie60.
Sucha warstwa została ponownie nawodniona nadmiarem 4 M NaCl, 100 mM NaPi przy pH 7,4 (Rysunek 3). Poziom nawodnienia i efektywne stężenie białka w objętości badanej przez falę ulotną można wywnioskować na podstawie współczynnika skalującego potrzebnego do cyfrowego usunięcia pasm absorpcyjnych z wody. Optymalny współczynnik skalowania wynosił 0,87, co oznacza, że w tym przypadku bufor zajmuje 87%, a próbka 13% objętości przy powierzchni. Biorąc pod uwagę gęstość białek i lipidów oraz stosunek lipidów do białek w błonach purpurowych (patrz wyżej), możemy wywnioskować efektywne stężenie białka 125 mg/ml w objętości badanej przez falę ulotną. Z poziomu nawodnienia możemy wywnioskować, że w tym konkretnym przypadku grubość warstwy próbki zwiększyła się ~6 razy po ponownym uwodnieniu, z ~1 do ~6 μm. Dla kąta padania 45°, typowego dla naszych i większości układów ATR, głębokość penetracji pola ulotnego (dp) dla uwodnionego filmu białkowego waha się od 0,3 do 0,6 μm w przedziale 1,800-850 cm-1 61. Ponieważ pole ulotne jest prawie niewrażliwe na próbkę znajdującą się dwa razy powyżej d p 61, wnioskujemy, że ilość białka użytego w eksperymencie ATR może być w zasadzie zmniejszona 5-krotnie bez znaczącej utraty sygnału.
Podczas korzystania z o niższej sile jonowej, film rozszerza się bardziej, a ilość białka w objętości badanej przez ulotne pole jest zmniejszona (Rysunek 2). Dokładna zależność między pęcznieniem filmu a siłą jonową buforu będzie zależeć między innymi od charakteru lipidów. Na przykład, podobne pęcznienie filmu, jak uzyskane tutaj dla bR w fioletowej błonie przy użyciu 4 M NaCl, uzyskano dla białka błonowego odtworzonego w polarnych lipidach E. coli przy użyciu zaledwie 0,1 M NaCl 62. Jeśli próbka zajmuje mniej niż 5% badanej objętości, należy rozważyć ponowne wykonanie etapu hydratacji przy użyciu buforu o wyższej sile jonowej. Obrzęk filmu po uwodnieniu wymaga pewnego czasu, aby osiągnąć stabilizację (ryc. 4). W przypadku bR w fioletowej membranie proces jest jednowykładniczy, ze stałą czasową 12 min: nie potrzeba więcej niż 30-60 minut, aby uzyskać stabilną ponownie uwodnioną warstwę
Rysunek 5 pokazuje widmo absorpcyjne IR uwodnionej warstwy ChR2 uzyskanej przez transmisję. W tym przypadku nawodnienie osiągnięto poprzez wystawienie suchej warstwy na działanie atmosfery o kontrolowanej wilgotności zapewnianej przez mieszaninę glicerolu i wody. Widmo można rozłożyć na wkład wody, detergentu i białka. Mogliśmy oszacować ilość każdego z nich za pomocą widm współczynnika ekstynkcji, przeskalowanych tak, aby pasowały do spektrum eksperymentalnego. Dla wody wzięliśmy spektrum współczynnika ekstynkcji z literatury63, a dla DM zmierzyliśmy je z roztworu 100 mg/ml (ryc. 6). Użyliśmy również współczynników ekstynkcji dla dwóch reprezentatywnych białek błonowych: mitochondrialnego nośnika ADP / ATP64 i bR (Figura 7). Współczynnik skalowania wskazuje gęstość powierzchniową wody i masy DM wynoszącą odpowiednio 260 μg/cm2 i 200 μg/cm2 w folii. Pozostała absorpcja (ryc. 5, czerwona linia) pochodzi głównie z białka, którego gęstość powierzchniową szacuje się na 250 μg/cm2 przy użyciu współczynnika ekstynkcji amidu II z dwóch różnych białek błonowych (patrz ryc. 6). Stosunek molowy białka/detergentu/wody obliczono na 1/60/2000 przy masie cząsteczkowej odpowiednio 35 kDa, 480 Da i 18 Da.
Wykres 3D z typowego eksperymentu FT-IR z rozdzielczością czasową na bR jest przedstawiony na rysunku 10A, uzyskanym przez ATR i obejmującym 200 minut zbierania danych. Widma można wyodrębnić w określonych momentach, na przykład, gdy oczekuje się, że produkty pośrednie L, M i N fotocyklu bR osiągną swoją największą populację (ryc. 11). Ich cechy spektralne zostały obszernie opisane13,41,65 i nie będą tutaj dalej omawiane. Rysunek 12 przedstawia ślady czasowe dla wybranych liczb falowych. Mianowicie, wzrost kinetyki przy 1,762 cm-1 (t1/2 ~ 60 μs) informuje o dynamice protonacji Asp85 z siatkówkowej zasady Schiffa66, a jej rozpad do zera wskazuje na jej deprotonację po odzyskaniu stanu67. Ujemny wzrost kinetyki na poziomie 1740 cm-1 (t1/2 ~ 1 ms) informuje o deprotonacji łańcucha bocznego Asp96, dawcy protonów dla zasady Schiffa68. Rozpad intensywności do zera informuje o jego retonacji z cytoplazmy67. Dynamikę zmian konformacyjnych białek można badać za pomocą zmian absorpcji w regionie amidu I i amidu II, który osiąga maksymalną zmianę po ~3 ms w temperaturze pokojowej67,69.
Zastosowanie SVD do problemów spektroskopowych zostało sprawdzone przed55,56. Krótko mówiąc, SVD rozkłada dane eksperymentalne na czynniki pierwsze, ułożone w macierz A, jako: A = U S VT. Ta faktoryzacja może być modyfikowana w celu uwzględnienia zależności szumu od liczby falowej i ukarania fluktuacji/dryftów w linii bazowej57. Kolumny U i V zawierają widma ortonormalne i wektory śladów czasowych odpowiednio dla każdego składnika SVD. S jest macierzą diagonalną zawierającą tzw. wartości osobliwe. (Abstrakcyjne) składowe spektrologiczno-czasowe w U i V wydają się mieć malejące znaczenie dla opisu danych eksperymentalnych w sensie najmniejszych kwadratów, określonych ilościowo przez przypisaną im wartość osobliwą. Rysunek 13A przedstawia wartości osobliwe w funkcji liczby składnika i odtwarza pierwsze osiem kolumn/składowych U (widma abstrakcyjne, rysunek 13B) i V (abstrakcyjne ślady czasowe, rysunek 13C). Sygnał jest skoncentrowany w pierwszych pięciu składowych (zgodnie z oczekiwaniami dla fotocyklu zawierającego pięć stanów pośrednich), przy czym składowe powyżej są w dużej mierze zdominowane przez losowe szumy i inne źródła błędów. Dane eksperymentalne zostały zrekonstruowane przy użyciu tylko pierwszych pięciu kolumn U i V oraz pierwszych pięciu kolumn i wierszy S. Dane zrekonstruowane przez SVD reprezentują najlepsze przybliżenie metodą najmniejszych kwadratów danych eksperymentalnych do macierzy rangi piątej. Zrekonstruowane dane wykazują poprawę jakości (Rysunek 10B). Mianowicie, szum jest w znacznym stopniu zredukowany, podobnie jak pewne ledwo zauważalne wahania i dryfty w linii podstawowej (powszechne w czasowo-rozdzielczej spektroskopii krokowej 25). Przetwarzanie SVD jest szczególnie korzystne dla poprawy jakości eksperymentalnych śladów czasowych w zakresie milisekund (czerwone linie na rysunku 12).
Dane skanowania krokowego z rozdzielczością czasową dla fotocyklu ChR2 uzyskano przy użyciu wyżej opisanego protokołu eksperymentów transmisyjnych (Rysunek 14A), co w przybliżeniu obejmowało 120 godzin skumulowanych pomiarów. Dane czasowo-rozdzielcze zebrane za pomocą skanowania krokowego rozciągają się od 6,25 μs do 125 ms. Fotocykl ChR2 wymaga około 60 s do pełnego wyzdrowienia. Najnowsza część fotocyklu może być pokryta za pomocą szybkiego skanowania FT-IR, a zestawy danych zarówno step-scan, jak i rapid-scan są połączone, jak pokazano w innym miejscu48. Zmiany widmowe ChR2 są ~10-krotnie mniejsze niż te uzyskiwane z bR, co sprawia, że pomiary są trudniejsze. Szczególnie oscylacje linii podstawowej w zakresie milisekund stają się wyraźnie widoczne przy tych niskich zmianach absorpcji. Są one spowodowane niewielkimi oscylacjami w ruchomym zwierciadle w milisekundowej skali czasowej25. Oscylacje, wraz z częścią szumu, mogą być w znacznym stopniu usunięte przez SVD, co znacznie poprawia wygląd danych (Rysunek 14B).

Rysunek 1. Widmo absorpcyjne bR w fioletowej membranie wysuszonej na wierzchu diamentowej powierzchni akcesorium ATR. Wskazane są pasma od drgań wiązania peptydowego (amid I, II i A).

Rysunek 2. Widmo absorpcyjne warstwy bR po ponownym uwodnieniu przy użyciu o różnej sile jonowej (rozcieńczenia 4 M NaCl, 100 mM Na2HPO4 / NaH2PO4 przy pH 7,4). Zwrócić uwagę na zmniejszenie pasma amidu II, tj. zwiększone pęcznienie filmu, przy malejącej sile jonowej buforu. (Wstawić) Względna ilość białka w pobliżu powierzchni, określona ilościowo przez intensywność absorbancji na poziomie 1,541 cm-1 (maksimum amidu II) po odjęciu wkładu absorpcyjnego buforu.

Rysunek 3. Widmo absorpcyjne warstwy bR ponownie uwodnionej buforem objętościowym (4,3 M siła jonowa) i po odjęciu wkładu buforu. Współczynnik odejmowania 0,87 został wybrany, aby usunąć silną absorpcję wody w zakresie 3 700-3 000 cm-1 i uzyskać płaską linię podstawową w zakresie 2 700-1 800 cm-1.

Rysunek 4. Widmo absorpcyjne bR po uwodnieniu buforem o sile jonowej 4,3 M (absorpcja buforu została odjęta dla jasności, jak na rysunku 3). Wkładka pokazuje ewolucję pęcznienia filmu po ponownym uwodnieniu, a następnie absorbancję białka amidu II na poziomie 1,541 cm-1. Dopasowanie do pojedynczego wykładnika wskazuje stałą czasową stabilizacji pęcznienia filmu wynoszącą 12 minut.

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

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

Rysunek 7. Przepuszczalność filtra optycznego stosowanego w pomiarach FT-IR metodą skanowania i skanowania w czasie.

Rysunek 8. Wydajność impulsów laserowych dostarczanych przez drugą harmoniczną (532 nm) lasera Nd:YAG. Histogram zmienności energii względnej 1000 impulsów laserowych i dopasowanie do rozkładu Gaussa z odchyleniem standardowym 0,05.

Rysunek 9. Wywołane światłem zmiany absorbancji bR przy trzech reprezentatywnych liczbach falowych w jednakowych odstępach czasowych (linie zielone, czarne i niebieskie) i po quasi-logarytmicznym uśrednieniu do ~20 punktów/dekadę (czerwone linie). Zwróć uwagę na redukcję szumów po uśrednieniu logarytmicznym, ujawniającym oscylacje śladów czasu w zakresie milisekund.

Rysunek 10. 3D przedstawienie zmian absorbancji bR wywołanych światłem zarejestrowanych przez ATR przy pH 7,4 (4 M NaCl, 100 mM NaPi). A) Dane surowe. ZA)Dane zrekonstruowane za pomocą pięciu komponentów SVD. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 11. Widma absorpcji różnic FT-IR fotocyklu bR w trzech wybranych momentach, w których produkty pośrednie L (12,5 μs), M (300 μs) i N (6 ms) są najbardziej wzbogacone.

Rysunek 12. Transfer protonów i dynamika szkieletu białkowego w fotocyklu bR rozwiązana za pomocą czasowo-rozdzielczej spektroskopii różnicowej FT-IR ze skanowaniem krokowym. Zmiany absorbancji przy 1,762 cm-1 informują o dynamice protonacji/deprotonacji Asp85, a przy 1,741 cm-1 o dynamice deprotonacji/retonacji Asp96 (w tym zmiany wiązania H przed 300 μs). Ślady czasowe na wysokości 1670 i 1555 cm-1 wskazują na zmiany w drganiach amidu I i II, które są wrażliwe na konformację szkieletu peptydowego. Ślady wsteczne odpowiadają danym pierwotnym, a czerwone – danym poddanym obróbce SVD. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 13. Dekompozycja wartości osobliwych (SVD) eksperymentalnych danych FT-IR ze skanowaniem krokowym w czasie fotocyklu bR (patrz rysunek 10A). SVD wykonano po zważeniu danych doświadczalnych z uwzględnieniem zależności odchylenia standardowego szumu od liczby falowej (ryc. 15); w połączeniu zpochodną 1. w celu zmniejszenia wagi statystycznej wahań linii bazowej, jak opisano wcześniej57. ZA)Wykres względnych wartości osobliwych pierwszych 50 składników (czarne kółka). Pierwsze pięć składowych jest przypisanych do składowych sygnałowych (czerwone kółka). Reszta komponentów zanika wykładniczo, zgodnie z oczekiwaniami dla komponentów szumu (patrz przerywana szara linia). ZA)Pierwsze osiem widm abstrakcyjnych (od U1 do U8). C)Pierwsze osiem abstrakcyjnych śladów czasowych (odV 1 do V8). Abstrakcyjne widma i ślady czasowe przypisane do sygnału są przedstawione za pomocą czerwonych linii, a w przeciwnym razie za pomocą czarnych linii. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 14. 3D przedstawiające zmiany absorbancji w podczerwieni wywołane światłem dla ChR2 zarejestrowane metodą skanowania krokowego w trybie transmisji. Dane skanowania krokowego rozciągają się do 125 ms, obejmując tylko część fotocyklu. ZA)Surowe dane. ZA)Dane zrekonstruowane za pomocą pięciu komponentów SVD. Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.

Rysunek 15. Szacowany poziom szumów w FT-IR. Widmo absorpcji różnic przy rozdzielczości czasowej 6,25 μs i rozdzielczości spektralnej 8 cm-1 dla eksperymentu obejmującego 1 pozycję koaddycji/lustra (500 fotoreakcji) lub ~200 pozycji koaddycji/lustra (105 fotoreakcji). Podane są wartości dla tłumienia całkowitego odbicia (ATR) i konfiguracji skrzyni biegów.