Stosując opisane powyżej metody, wygenerowaliśmy rozpuszczalny TCR (Tabela 1) oraz cząsteczki pHLA, aby przeprowadzić szczegółowe analizy molekularne rozpoznawania gp100280-288 przez limfocyty T CD8+. W celu wytworzenia nierozpuszczalnych ciał inkluzyjnych (IB) dla każdego oddzielnego łańcucha zarówno TCR (łańcuchy α i β), jak i pHLA (łańcuch α i β2m), zastosowano zmodyfikowany system ekspresyjny E. coli. Metoda ta ma zaletę niskiego kosztu i łatwości wdrożenia oraz pozwala na uzyskanie wysokich wydajności białka (100-500 mg/L hodowli). Ponadto nierozpuszczalne białka są bardzo stabilne przy przechowywaniu w temperaturze -80 °C. Następnie wykorzystano sprawdzoną technikę refoldingowania i oczyszczania w celu uzyskania funkcjonalnych, jednorodnych i rozpuszczalnych białek. Metoda ta jest przydatna do otrzymywania białek do eksperymentów biofizycznych, strukturalnych i komórkowych, a także odczynników, które mogą być wykorzystane w diagnostyce lub terapii.
W tym badaniu wykorzystano te białka do przeprowadzenia eksperymentów z mutagenezą skanującą alaniną w obrębie szkieletu peptydu, a powinowactwo wiązania TCR oceniono za pomocą eksperymentów z wykorzystaniem rezonansu plazmonowego powierzchni (SPR) (Tabela 2). Test ten pozwolił wskazać, które reszty w peptydzie były najważniejsze dla wiązania TCR. Wysokorozdzielcze analizy powinowactwa wiązania z wykorzystaniem tej techniki są niezwykle użyteczne w określaniu mechanizmów biologicznych kontrolujących oddziaływania białko-białko, a także w analizie powinowactwa wiązania cząsteczek terapeutycznych.
Następnie poddaliśmy krystalizacji rozpuszczalny TCR specyficzny dla czerniaka (PMEL17 TCR) w kompleksie z modyfikowanym pHLA pochodzącym z guza (A2-YLE-9V), aby zbadać sposób wiązania w rozdzielczości atomowej (Ryciny 1 i 2 oraz Tabela 3). Eksperymenty te zapewniły bezpośrednią wizualizację interfejsu wiązania między dwiema cząsteczkami, dostarczając kluczowych informacji o zasadach sterujących tą interakcją. Przeprowadziliśmy ponadto analizę termodynamiczną oddziaływania z wykorzystaniem metod SPR oraz ITC, co pozwoliło wykazać wkłady energetyczne umożliwiające wiązanie (Rycina 3). Analizy te zostały dodatkowo wsparte wysokorozdzielczym opisem powierzchni kontaktu między obiema białkami (Rycina 4 i Tabela 4).
Następnie rozwiązaliśmy struktury niepowiązanych cząsteczek pHLA, prezentujących zmutowane formy peptydów, co wykazało, że przełącznik molekularny może wyjaśniać, dlaczego niektóre mutacje znosiły wiązanie z TCR (Ryciny 5).
Podsumowując, techniki te dostarczyły nowych danych demonstrujących mechanizm wyjaśniający sposób, w jaki limfocyty T rozpoznają antygen pochodzący z czerniaka, który stanowi istotny cel dla terapii przeciwnowotworowych. W szerszym kontekście techniki te mogą być wykorzystane do badania niemal każdej interakcji receptor-ligand, odsłaniając nowe mechanizmy biologiczne, które mogłyby stać się celem innowacyjnych postępów terapeutycznych.

Rycina 1: Analiza wykresu gęstości. Lewa kolumna przedstawia mapy pominięć (omit maps), w których model został udoskonalony przy pominięciu peptydu. Gęstość różnicowa jest wyznaczona dla poziomu 3,0 sigma, kontury dodatnie są zaznaczone kolorem zielonym, a ujemne kolorem czerwonym. Prawa kolumna przedstawia mapę obserwowaną przy poziomie 1,0 sigma (wyświetloną jako szara siatka wokół reprezentacji proteinowych łańcuchów w formie prętów) po kolejnym udoskonaleniu z zastosowaniem automatycznych więzów symetrii niekrystalograficznej za pomocą programu REFMAC5.A) Model PMEL17 TCR-A2-YLE-9V z pętlami CDR3 TCR zaznaczonymi kolorem niebieskim (α łańcuch) i pomarańczowy (β łańcuch) oraz peptyd na zielono. (BModel A2-YLE z peptydem zaznaczonym na ciemnozielono.CModel A2-YLE-3A z peptydem zaznaczonym na pomarańczowo (w przypadku A2-YLE-3A w jednostce asymetrycznej znajdowały się 2 cząsteczki, lecz były one praktycznie identyczne pod względem map omit i map gęstości, dlatego przedstawiono tutaj tylko kopię 1).DModel A2-YLE-5A z peptydem zaznaczonym kolorem różowym. Przedrukowano za zgodą źródła31. Kliknij tutaj, aby wyświetlić powiększoną wersję tej figury.

Rycina 2: Przegląd TCR PMEL17 w kompleksie z A2-YLE-9V. (A) Schematyczny model kompleksu TCR PMEL17-A2-YLE-9V. TCR zaznaczono kolorem czarnym; pętle CDR TCR są przedstawione w kolorach (czerwony, CDR1α; ciemnozielony, CDR2α; niebieski, CDR3α; żółty, CDR1β; aqua, CDR2β; pomarańczowy, CDR3β); natomiast HLA-A*0201 przedstawiono w kolorze szarym. Peptyd YLE-9V przedstawiono jako zielone pręty. (B) Reprezentacje powierzchniowe i prętowe reszt pętli CDR TCR PMEL17 (oznaczone kolorami jak w A), które kontaktują się z powierzchnią A2-YLE (A2, kolor szary; YLE-9V, zielone pręty). Czarna linia ukośna wskazuje kąt przecięcia TCR względem długiej osi peptydu YLEPGPVTV (46,15°). (C) Ślad kontaktu TCR PMEL17 na powierzchni A2-YLE-9V (A2, kolor szary); kolor fioletowy i zielony (powierzchnia i pręty) wskazują odpowiednio reszty HLA-A*0201 i YLE kontaktujące się z TCR gp100. Wartość odcięcia 3,4 Å dla wiązań wodorowych i 4 Å dla kontaktów van der Waalsa. Przedrukowano za zgodą z odniesienia 31. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 3: Analiza termodynamiczna oddziaływania PMEL17 TCR-A2-YLE. (A) Równowagowe odpowiedzi wiązania PMEL17 TCR z A2-YLE w temperaturach 5, 12, 18, 25 i 37 °C dla dziewięciu do dziesięciu rozcieńczeń seryjnych TCR. Surowe i dopasowane dane SPR (przy założeniu wiązania Langmura 1:1) przedstawiono we wstawkach do każdej krzywej; posłużyły one do obliczenia wartości Kon i Koff przy użyciu algorytmu globalnego dopasowania (BIAevaluation 3.1). Tabela przedstawia stałe wiązania równowagowego (KD (E)) i wiązania kinetycznego (KD (K) = Koff/Kon) dla każdej temperatury. Wartości stałej wiązania równowagowego (KD, µM) obliczono za pomocą nieliniowego dopasowania (y = (P1x)/(P2+x)). (B) Parametry termodynamiczne obliczono zgodnie z równaniem Gibbsa-Helmholtza (ΔG° = ΔH° − TΔS°). Wolne energie wiązania, ΔG° (ΔG° = -RTlnKD), naniesiono na wykres w funkcji temperatury (K), stosując regresję nieliniową do dopasowania równania trójparametrycznego (y = ΔH°+ΔCp°*(x-298)-x*ΔS°-x*ΔCp°*ln(x/298)). Entalpia (ΔH°) i entropia (TΔS°) w temperaturze 298 K (25 °C) podano w kcal/mol i obliczono poprzez regresję nieliniową temperatury (K) w funkcji energii swobodnej (ΔG°). (C) Pomiary izotermicznej kalorymetrii miareczkowej (ITC) dla oddziaływania PMEL17 TCR-A2-YLE. Entalpia (ΔH°) i entropia (TΔS°) w temperaturze 298 K (25 °C) podano w kcal/mol. Przedrukowano za zgodą z referencji31. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 4: Pętle CDR PMEL17 koncentrują się na resztach peptydowych Pro4, Val7 i Thr8. (A) Schematyczne przedstawienie kontaktów między peptydem YLE-9V a resztami pętli CDR PMEL17 (kodowanie kolorami jak na Rysunku 2A). Liczby na dole panelu wskazują całkowitą liczbę kontaktów między TCR a peptydem. (B) Kontakty między TCR PMEL17 a peptydem YLE-9V (zielone pręty) pokazujące kontakty van der Waalsa (czarne przerywane linie) oraz wiązania wodorowe (czerwone przerywane linie) tworzone przez pętle CDR3α (niebieska), CDR1β (żółta), CDR2β (akwamarynowa) i CDR3β (pomarańczowa) TCR. W dolnym panelu znajduje się powiększony widok kontaktów odpowiednio między YLE Pro4, Val7 i Thr8 a resztami pętli CDR TCR (kolory prętów zgodnie z Rysunkiem 1A). Wartość odcięcia wynosi 3.4 Å dla wiązań wodorowych i 4 Å dla kontaktów van der Waalsa. Przedrukowano za zgodą z źródła31. Prosimy kliknąć tutaj, aby zobaczyć powiększoną wersję tego rysunku.

Rycina 5: Porównanie konformacyjne peptydów YLE, YLE-3A i A2-YLE-5A prezentowanych przez HLA-A*0201. (A) Dopasowanie peptydów YLE (ciemnozielone pręty) i YLE-3A (pomarańczowe pręty) poprzez nałożenie helisy α1 HLA-A*0201 (szary model wstążkowy). Obramowane reszty wskazują mutację Glu3 na alaninę. Wstawki pokazują, w jaki sposób substytucja Glu3Ala powoduje przesunięcie pozycji (czarna strzałka) sąsiedniej reszty Pro4 w strukturze A2-YLE-3A w porównaniu do struktury A2-YLE. (B) Dopasowanie peptydów YLE (ciemnozielone pręty) i YLE-5A (różowe pręty) poprzez nałożenie helisy α1 HLA-A*0201 (szary model wstążkowy). Obramowane reszty wskazują mutację glicyny 5 na alaninę. Przedrukowano za zgodą z referencji31. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
| TCR | CDR1α | CDR2α | CDR3α | CDR1β | CDR2β | CDR3β |
| PMEL17 | DSAIYN | IQSSQRE | CAVLSSGGSNYKLTFG | SGHTA | FQGTGA | CASSFIGGTDTQYFG |
| gp100 | TSINN | IRSNERE | CATDGDTPLVFG | LNHDA | SQIVND | CASSIGGPYEQYFG |
| MPD | KALYS | LLKGGEQ | CGTETNTGNQFYFG | SGHDY | FNNNVP | CASSLGRYNEQFFG |
| 296 | DSASNY | IRSNVGE | CAASTSGGTSYGKLTFG | MNHEY | SMNVEV | CASSLGSSYEQYFG |
Tabela 1: Dopasowanie regionów TCR CDR3 dla PMEL17, gp100, MPD oraz 296 TCR specyficznych dla gp100. Przedrukowano za zgodą autorów źródła31.
| Sekwencja peptydowa | Peptyd | TCR PMEL17 TRAV21 TRBV7-3 Stała powinowactwa KD | TCR gp100 TRAV17 TRBV19 Stała powinowactwa KD |
| YLEPGCTVTA | YLE | 7.6 ±2 μM | 26.5 ±2.3 μM |
| YLEPGPVTV | YLE-9V | 6.3 ±1.2 μM | 21.9 ±2.4 μM |
| ALEPGPVTA | YLE-1A | 15.9 ±4.1 μM | 60.6 ±5.4 μM |
| YLAPGPVTA | YLE-3A | Brak wiązania | Brak wiązania |
| YLEAGPVTA | YLE-4A | 19.7 ±1.3 μM | 144.1 ±7.8 μM |
| YLEPAPVTA | YLE-5A | >1 mM | >1 mM |
| YLEPGAVTA (pole brzuszne nakrywki) | YLE-6A | 11.4 ±2.7 μM | 954.9 ±97.8 μM |
| YLEPGPATA | YLE-7A | 31.1 ±4 μM | 102.0 ±9.2 μM |
| YLEPGPVAA | YLE-8A | 38.1 ±7.4 μM | 121.0 ±7.5 μM |
Tabela 2: Analiza powinowactwa (KD) TCR PMEL17 oraz TCR gp100 do wariantów peptydów gp100280-288. Przedrukowano za zgodą z referencji31.
| Parametry | PMEL17 TCR-A2-YLE-9V | A2-YLE | A2-YLE-3A | A2-YLE-5A |
| Kod PDB | 5EU6 | 5EU3 | 5EU4 | 5EU5 |
| Statystyki zbioru danych | | | | |
| Grupa przestrzenna | P1 | P1 21 1 | P1 | P1 21 1 |
| Parametry komórki elementarnej (Å) | a= 45.52, b= 54.41, c= 112.12, a=85.0°, b=81.6°, g=72.6° | a= 52.81, b= 80.37, c= 56.06, b=112.8° | a= 56.08, b= 57.63, c= 79.93, a=90.0°, b=89.8°, g=63.8° | a= 56.33, b= 79.64, c= 57.74, b=116.2° |
| Źródło promieniowania | DIAMOND I03 | DIAMOND I03 | DIAMOND I02 | DIAMOND I02 |
| Długość fali (Å) | 0.9763 | 0.9999 | 0.9763 | 0.9763 |
| Zmierzony zakres rozdzielczości (Å) | 51.87 – 2.02 | 45.25 – 1.97 | 43.39 – 2.12 | 43.42 – 1.54 |
| Zewnętrzna powłoka rozdzielczości (Å) | 2.07 - 2.02 | 2.02 – 1.97 | 2.18 – 2.12 | 1.58 - 154 |
| Obserwowane refleksy | 128,191 (8,955) | 99,442 (7,056) | 99,386 (7,463) | 244,577 (17,745) |
| Unikalne refleksy | 64,983 (4,785) | 30,103 (2,249) | 49,667 (3,636) | 67,308 (4,962) |
| Kompletność (%) | 97.7 (96.7) | 98.5 (99.3) | 97.4 (96.7) | 99.6 (99.9) |
| Multiplikacja | 2.0 (1.9) | 3.3 (3.1) | 2.0 (2.1) | 3.6 (3.6) |
| I/Sigma(I) | 5.5 (1.9) | 7.2 (1.9) | 6.7 (2.3) | 13 (2.3) |
| Rpim (%) | 5.7 (39.8) | 8.8 (44.7) | 8.7 (41.6) | 4.5 (35.4) |
| Rmerge (%) | 7.8 (39.6) | 9.8 (50.2) | 8.7 (41.6) | 5.0 (53.2) |
| Statystyki udoskonalenia | | | | |
| Rozdzielczość (Å) | 2.02 | 1.97 | 2.12 | 1.54 |
| Liczba wykorzystanych refleksów | 61688 | 28557 | 47153 | 63875 |
| Liczba refleksów w zbiorze Rfree | 3294 | 1526 | 2514 | 3406 |
| Rcryst (bez odcięcia) (%) | 18.1 | 19.7 | 17.2 | 17.0 |
| Rfree | 22.2 | 25.5 | 21.1 | 20.1 |
| Średnio kwadratowe odchylenie od geometrii idealnej | | | | |
| Długości wiązań (Å) | 0.018 (0.019)* | 0.019 (0.019)* | 0.021 (0.019)* | 0.018 (0.019)* |
| Kąty wiązań (°) | 1.964 (1.939)* | 1.961 (1.926)* | 2.067 (1.927)* | 1.914 (1.936)* |
| Całkowity błąd współrzędnych (Å) | 0.122 | 0.153 | 0.147 | .055 |
| Statystyki Ramachandrana | | | | |
| Obszary najbardziej preferowane | 791 (96%) | 371 (98%) | 749 (99%) | 384 (98%) |
| Obszary dopuszczalne | 32 (4%) | 6 (2%) | 10 (1%) | 5 (1%) |
| Wartości odstające | 2 (0%) | 3 (1%) | 1 (0%) | 2 (0%) |
Tabela 3: Statystyki redukcji i udoskonalenia danych (wymiana molekularna). Przedrukowano za zgodą z referencji31. Wartości w nawiasach dotyczą powłoki o najwyższej rozdzielczości.
| pozostałość HLA/peptydu | pozostałość TCR | Liczba vdW (≤4)Å) | Liczba wiązań wodorowych (≤3,4Å) |
| Gly62 | αGly98 | 3 | |
| αSer99 | 1 | |
| Arg65 | αSer99 | 2 | |
| Arg65 O | αAsn100 Nδ2 | 2 | 1 |
| Arg65 NH1 | βAsp58 Oδ2 | | 1 |
| βSer59 | 8 | |
| Lys66 | αGly98 | 1 | |
| αSer99 | 4 | |
| αAsn100 | 4 | |
| Ala69 | αAsn100 | 2 | |
| βAla56 | 2 | |
| Gln72 Nε2 | βGln51 O | 3 | 1 |
| βGly54 | 7 | |
| βAla55 | 1 | |
| Thr73 | βGln51 | 1 | |
| Val76 | βGln51 | 3 | |
| βGly52 | 2 | |
| Lys146 | βPhe97 | 3 | |
| βIle98 | 3 | |
| Ala150 | βIle98 | 1 | |
| βAsp102 | 3 | |
| Val152 | βIle98 | 1 | |
| Glu154 | αTyr32 | 1 | |
| Gln155 N | αTyr32 OH | 4 | 1 |
| Gln155 Oε1 | βThr101 N | 10 | 1 |
| Tyr1OH | αGly97 O | 1 | 1 |
| αGly98 | 1 | |
| αSer96 | 1 | |
| Glu3 | αTyr101 | 1 | |
| Pro4 | αSer96 | 1 | |
| αSer99 | 1 | |
| αAsn100 | 4 | |
| Pro4 O | αTyr101N | 14 | 1 |
| Gly5 | αTyr101 | 3 | |
| βGly100 | 2 | |
| Val7 | βIle98 | 7 | |
| βGly99 | 2 | |
| βGly100 | 2 | |
| Thr8 | βThr31 | 5 | |
| βGln51 | 1 | |
| βPhe97 | 1 | |
| Thr8 N | βIle98 O | 6 | 1 |
Tabela 4: Tabela kontaktów PMEL17 TCR-A2-YLE-9V. Przedrukowano za zgodą z referencji31.