Obliczeniowe przewidywanie preferencji aminokwasów potencjalnie wielospecyficznych domen wiążących peptydy zaangażowanych w interakcje białko-białko

Podsumowanie

Opisujemy metodologię opartą na dywersyfikacji sekwencji w celu oszacowania preferencji aminokwasowych wielospecyficznych miejsc wiązania w interakcjach białko-białko (PPI). W ramach tej strategii tysiące potencjalnych ligandów peptydowych są generowane i badane in silico, pokonując w ten sposób niektóre ograniczenia dostępnych metod eksperymentalnych.

Streszczenie

Wiele interakcji białko-białko wiąże krótkie segmenty białka z domenami wiążącymi peptydy. Zazwyczaj takie oddziaływania wymagają rozpoznania motywów linearnych o zmiennej konserwacji. Połączenie wysoce konserwatywnych i bardziej zmiennych regionów w tych samych ligandach często przyczynia się do wielospecyficzności wiązania, co jest powszechną właściwością enzymów i białek sygnalizacji komórkowej. Charakterystyka preferencji aminokwasowych domen wiążących peptydy jest ważna dla projektowania mediatorów interakcji białko-białko (PPI). Metody obliczeniowe stanowią skuteczną alternatywę dla często kosztownych i kłopotliwych technik eksperymentalnych, umożliwiając projektowanie potencjalnych mediatorów, które można później zweryfikować w dalszych eksperymentach. W tym miejscu opisaliśmy metodologię wykorzystującą zastosowanie Pepspec w pakiecie do modelowania molekularnego Rosetta do przewidywania preferencji aminokwasowych domen wiążących peptydy. Metodologia ta jest przydatna, gdy struktura białka receptorowego i charakter liganda peptydowego są znane lub można je wywnioskować. Metodologia rozpoczyna się od dobrze scharakteryzowanej kotwicy z liganda, która jest przedłużana przez losowe dodawanie reszt aminokwasowych. Powinowactwo wiązania wytworzonych w ten sposób peptydów jest następnie oceniane za pomocą dokowania peptydów z elastycznym szkieletem w celu wybrania peptydów o najlepszych przewidywanych wynikach wiązania. Peptydy te są następnie wykorzystywane do obliczania preferencji aminokwasowych i opcjonalnie do obliczania macierzy pozycja-masa (PWM), którą można wykorzystać w dalszych badaniach. Aby zilustrować zastosowanie tej metodologii, wykorzystaliśmy interakcję między podjednostkami ludzkiego czynnika regulacyjnego interferonu 5 (IRF5), wcześniej znanego jako wielospecyficzny, ale globalnie kierowany przez krótki konserwatywny motyw zwany pLxIS. Oszacowane preferencje aminokwasowe były zgodne z wcześniejszą wiedzą na temat powierzchni wiązania IRF5. Pozycje zajmowane przez fosforylujące reszty seryny wykazywały wysoką częstotliwość asparaginianu i glutaminianu, prawdopodobnie dlatego, że ich ujemnie naładowane łańcuchy boczne są podobne do fosfoseryny.

Wprowadzenie

Oddziaływanie między dwoma białkami często obejmuje wiązanie krótkich odcinków aminokwasów z domenami wiążącymi peptydy, co przypomina interfejsy białko-peptyd. Białka receptorowe uczestniczące w takich oddziaływaniach białko-białko (PPI) często posiadają zdolność rozpoznawania określonego zestawu nakładających się, lecz zróżnicowanych sekwencji ligandów, co jest właściwością znaną jako wielospecyficzność^1,2. Rozpoznawanie wielospecyficzne jest cechą wielu białek komórkowych, ale jest szczególnie wyraźne w przypadku enzymów i białek sygnałowych komórki³. Białka oddziałujące z wielospecyficznymi miejscami wiązania często posiadają w swojej sekwencji kombinację regionów bardziej i mniej konserwatywnych^4,5,6. W takim scenariuszu bardziej konserwatywne motywy sekwencji biorą udział w rygorystycznych oddziaływaniach molekularnych. Z kolei sekwencje bardziej zmienne oddziałują z powierzchniami w miejscu wiązania receptora, które są w pewnym stopniu bardziej tolerancyjne. Zazwyczaj te mniej konserwatywne, lecz wciąż istotne funkcjonalnie segmenty są pętlami pozbawionymi zdefiniowanych wzorów struktury drugorzędowej lub przyjmują jeszcze bardziej dynamiczne konformacje, takie jak te typowe dla białek wewnętrznie nieuporządkowanych⁷.

Identyfikacja potencjalnych ligandów peptydowych miejsc wiązania jest zazwyczaj pierwszym krokiem w projektowaniu mediatorów zdolnych do zakłócania odpowiadających im PPI⁸. Jednakże w przypadku ligandów wielospecyficznych miejsc wiązania często mało prawdopodobne jest znalezienie jednego, najczęściej występującego reszty aminokwasowej w większości pozycji sekwencji. Zamiast tego miejsca te mogą wykazywać szczególne preferencje dla konkretnej klasy aminokwasów, w zależności od ich właściwości chemicznych, np. aminokwasów kwasowych i naładowanych ujemnie, takich jak asparaginian lub glutaminian, objętościowych aminokwasów aromatycznych, takich jak fenyloalanina, lub bardziej hydrofobowych reszt, takich jak alifatyczne aminokwasy alanina, walina, leucyna lub izoleucyna³. Kilka metod eksperymentalnych może dostarczyć informacji na temat preferencji aminokwasowych miejsc wiązania białek, w tym ewolucja sterowana⁹, mutageneza z wielokodonowym skanowaniem¹⁰ oraz głębokie skanowanie mutacyjne¹¹. Wszystkie te metody opierają się na podejściu dywersyfikacji sekwencji, które polega na wprowadzaniu mutacji do oryginalnych ligandów i dalszej analizie ich wpływu na funkcję białka receptorowego (szeroki przegląd temu zagadnieniu poświęcono w pracy Bratulic i Badran¹²). Jednak metody te często wymagają analizy dużych bibliotek sekwencji, co czyni je bardziej uciążliwymi, kosztownymi i czasochłonnymi.

Metody obliczeniowe służące do wnioskowania o preferencjach aminokwasowych wielospecyficznych miejsc wiązania mają potencjał do obejścia ograniczeń metod laboratoryjnych. Wśród nich podejście polegające na dywersyfikacji sekwencji in silico ocenia wpływ energetyczny szerokiego zakresu zamian aminokwasowych w sekwencji ligandu w celu scharakteryzowania plastyczności strukturalnej PPI¹³. Metoda ta rozpoczyna się od struktury lub modelu ligandu peptydowego związanego z miejscem wiązania receptora, a następnie wprowadza mutacje do sekwencji ligandu. Następnie wykorzystuje się statystyczne i energetyczne funkcje oceny, aby przeanalizować wpływ tych mutacji na stabilność i powinowactwo wiązania. Zbiór najlepiej ocenionych sekwencji ligandów wynikający z fazy oceny może zostać następnie wykorzystany do obliczenia preferencji aminokwasowych. Strategia ta pozwala na efektywne przetwarzanie bardzo dużej liczby sekwencji ligandów. Dzięki temu może ona zapewnić pełniejsze i bardziej spójne wnioskowanie o preferencjach aminokwasowych w porównaniu z wynikami uzyskanymi z ograniczonej liczby sekwencji, które zazwyczaj można przetworzyć w podejściach laboratoryjnych.

Aplikacja Pepspec z pakietu modelowania molekularnego Rosetta¹⁴ jest narzędziem, które przeprowadza dywersyfikację sekwencji jako kluczowy krok w trybie projektowania peptydów. Aplikacja ta wymaga struktury lub modelu białka receptorowego z przyłączonym peptydem o długości nawet jednego aminokwasu, który służy jako kotwica dla kolejnych etapów. Sekwencja przyłączonego peptydu jest następnie przedłużana (jeśli to konieczne) i dywersyfikowana w celu wygenerowania licznej grupy domniemanych ligandów peptydowych. Powinowactwo wiązania tych peptydów jest następnie oceniane za pomocą dokowania peptydów z elastycznym szkieletem, aby wybrać te z najlepszymi przewidywanymi wynikami wiązania. Chociaż głównym wynikiem działania tej aplikacji są najlepsze kandydaty peptydowe wybrane na koniec fazy projektowania, znacznie większy zbiór peptydów zaakceptowanych w tej fazie może być również wykorzystany do obliczenia preferencji aminokwasowych w docelowym miejscu wiązania. Preferencje aminokwasowe oblicza się jako częstotliwość występowania każdego aminokwasu w danej pozycji sekwencji ligandu, przedstawiając je w formie macierzy wag pozycyjnych (PWM) lub bardziej wizualnego logo sekwencyjnego.

W niniejszym artykule opisujemy protokół służący do szacowania preferencji aminokwasowych powierzchni wiążącej białka receptorowego uczestniczącego w PPI. Protokół koncentruje się na PPI, w których wiadomo, że liniowy fragment białka-liganda wiąże się z białkiem receptorowym, zatem scenariusz ten można modelować jako interfejs białko-peptyd. W takim przypadku konserwowane motywy liganda zazwyczaj oddziałują z określonymi kieszeniami w miejscu wiązania receptora, choć cały segment liganda zaangażowany w PPI może zawierać obszary mniej konserwowane. Schemat blokowy podsumowujący główne etapy protokołu przedstawiono na Rysunku 1. Protokół rozpoczyna się od struktury 3D kompleksu białko-białko, a następnie redukuje białko-ligand do potencjalnie najlepiej oddziałującego segmentu, pozostawiając białko receptorowe w niezmienionej formie. Najlepiej oddziałujący segment jest wyznaczany przy użyciu serwera BUDE Alanine Scan¹⁵, który przeprowadza obliczeniowe skanowanie mutagenezą alaninową w celu zidentyfikowania reszt typu hot-spot pomiędzy dwiema oddziałującymi białkami. W tym podejściu reszty liganda są pojedynczo zastępowane alaniną, a oszacowana zmiana wolnej energii lub stabilności kompleksu (ΔΔG) jest następnie wykorzystywana do wnioskowania o znaczeniu odpowiadającej reszty dla docelowego PPI. Po wyznaczeniu najlepiej oddziałującego segmentu, jego kompleks z białkiem receptorowym jest wykorzystywany jako struktura bazowa przesyłana do Pepspec w celu przeprowadzenia dywersyfikacji sekwencji.

Rycina 1: Przegląd głównych etapów protokołu zaproponowanego w niniejszej pracy. Numery odpowiadają numerom kroków w sekcji protokołu. Ryciny zostały przygotowane z wykorzystaniem kompleksu białko-białko jako przykładu opisanego w tekście. W kompleksie tym łańcuch białkowy uznany za receptor jest pokazany na różowo, natomiast łańcuch uznany za ligand jest pokazany na jasnoniebiesko, z wyróżnionym na czerwono przewidywanym najlepiej oddziałującym segmentem. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Jednym z ograniczeń sugerowanego protokołu jest wymóg posiadania rozdzielonej struktury interfejsu białko-peptyd. Protokół może alternatywnie rozpocząć się od modelu docelowego interfejsu białko-peptyd, choć konkretne kroki modelowania nie zostały tutaj opisane. Ponadto, choć protokół można przeprowadzić na komputerze osobistym z dowolnym systemem operacyjnym, do kroków obejmujących aplikacje Rosetta wymagane jest środowisko Linux. Ze względu na dużą liczbę iteracji zazwyczaj wykonywanych przez program Pepspec, do etapu dywersyfikacji sekwencji wysoce zaleca się użycie klastra komputerowego.

Zastosowanie sugerowanego protokołu zilustrowano na przykładzie oceny preferencji aminokwasowych powierzchni wiążącej białka IRF5, należącego do ludzkiej rodziny czynników regulacyjnych interferonu (IRF). Wybrano to białko jako przykład, ponieważ podczas jego aktywacji dwie podjednostki wiążą się, tworząc dimer o dobrze scharakteryzowanej strukturze¹⁶. W dimerach IRF wiązanie można modelować jako interfejs białko-peptyd, w którym jedna podjednostka stanowi powierzchnię wiążącą, a druga oddziałuje poprzez region zawierający krótki konserwatywny motyw zwany pLxIS^17,18. Ponadto wiązanie z podjednostkami IRF jest wielospecyficzne, w związku z czym mogą one tworzyć homodimery, heterodimery oraz kompleksy z innymi białkami komórkowymi, znanymi jako koaktywatory¹⁸.

Protokół

1. Wstępne przygotowanie interfejsu białko-peptyd

1. Pobieranie struktury kompleksu białko-białko
   1. Przejdź do strony głównej Protein Data Bank (PDB) (https://www.rcsb.org/) i wpisz w głównym polu wyszukiwania identyfikator PDB dla struktury kompleksu białko-białko (Rysunek 2A). Identyfikator PDB dla struktury dimeru IRF5, wykorzystany jako przykład w tej pracy, to 3DSH¹⁹.
   2. Na stronie głównej wybranej struktury kliknij Download Files (Rysunek 2B), a następnie Biological Assembly 1 (PDB - gz) (Rysunek 2C).
      UWAGA: W bazie danych PDB struktury wielu kompleksów białkowych utworzonych z identycznych monomerów są reprezentowane jako zgrupowania biologiczne (biological assemblies), w których w pliku PDB zapisana jest struktura tylko jednego monomeru (jednostka asymetryczna). Strukturę multimeru, w tym przypadku dimeru IRF5, należy pobrać jako zgrupowanie biologiczne zawierające dwie kopie jednostki asymetrycznej. Aby ułatwić kolejne kroki tego protokołu, dwa monomery są najpierw rozdzielane, a następnie przypisuje się im różne identyfikatory łańcuchów.
   3. Otwórz pobraną strukturę w programie UCSF Chimera²⁰ i kliknij Tools > Structure Editing > Change Chain IDs. W tym przykładzie oba łańcuchy w zgrupowaniu biologicznym nazywają się A. Zmień nazwę drugiego łańcucha (oznaczonego jako #0.2) na B i kliknij OK.
   4. Kliknij Favorites > Model Panel, a następnie wybierz model zawierający oba łańcuchy. Kliknij przycisk Group/Ungroup, aby rozdzielić każdy łańcuch do osobnego modelu. Następnie zaznacz oba modele i kliknij przycisk Copy/Combine. Wprowadź nową nazwę dla połączonego modelu, zaznacz opcję Close Source Models i kliknij OK.
   5. Kliknij Select > Chain i potwierdź, że każdy łańcuch w dimerze jest teraz zidentyfikowany inną literą, mianowicie A i B.
   6. Użyj funkcji File > Save PDB, aby zapisać zmodyfikowaną strukturę w nowym pliku PDB, który będzie wykorzystany w kolejnych krokach protokołu (tutaj zastosowano nazwę IRF5_dimer.pdb).

Rycina 2: Strona Protein Data Bank (PDB) dla struktury wykorzystanej jako reprezentatywny przykład w niniejszej pracy. (A) Pole wyszukiwania do wprowadzenia kodu dostępu PDB docelowej struktury. (B) Menu umożliwiające pobranie struktury w kilku formatach. (C) Opcje pobierania zespołów biologicznych w przypadku, gdy struktura została zapisana jako jednostka asymetryczna (szczegóły w kroku 1.1.2). Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

2. Identyfikacja docelowego segmentu w białku ligandu
   1. Przejdź do serwera BUDE Alanine Scan (https://pragmaticproteindesign.bio.ed.ac.uk/balas/). Kliknij przycisk Choose File w sekcji Structure Upload i wybierz plik PDB zapisany w kroku 1.1.6.
   2. Na następnej stronie sprawdź, czy struktura została poprawnie wczytana (Rycina 3A), a następnie wprowadź nazwę zadania w serwerze (Rycina 3B).
   3. Wyznacz łańcuchy z pliku PDB, które będą traktowane jako receptor (A) i ligand (B) (Rycina 3C). Następnie kliknij przycisk Start Scan, aby przesłać zadanie.
   4. Po zakończeniu zadania kliknij Show Results, aby otworzyć stronę z wynikami (Rycina 4).
      UWAGA: Na stronie wyników reszty struktury ligandu są pokolorowane zgodnie z szacowaną zmianą energii swobodnej (ΔΔG); reszty o wyższych wartościach są zaznaczone na czerwono.
   5. Z listy reszt wybierz odcinek aminokwasów, który według przewidywań będzie najlepiej oddziaływał z docelową powierzchnią wiązania. Upewnij się, że reszty te grupują najwyższe wartości różnicy energii swobodnej (ΔΔG). W tym przykładzie wybrano segment między resztami Leu424 a Ser436 (zaznaczony czerwoną ramką na prawym panelu Ryciny 4).
3. Przygotowanie interfejsu białko-peptyd do dywersyfikacji sekwencji
   1. Otwórz w programie Chimera plik PDB zapisany w kroku 1.1.6 i upewnij się, że w strukturze podjednostek docelowych nie brakuje atomów ani wiązań.
   2. Usuń wszystkie małe cząsteczki, jony i rozpuszczalniki, które współkrystalizowały z oryginalną strukturą. Aby to zrobić, kliknij Select > Residues, a następnie zaznacz wszystkie cząsteczki inne niż standardowe aminokwasy. Następnie kliknij Actions > Atoms/Bonds i Delete.
   3. Przytnij łańcuch ligandu do najlepiej oddziałującego segmentu wybranego w kroku 1.2.5. Aby to zrobić, kliknij Favorites oraz Sequence , a następnie kliknij łańcuch uznany za ligand (B). W panelu Sequence przeciągnij myszą, aby zaznaczyć wszystkie reszty z wyjątkiem tych między pozycjami 424 a 436. Aby usunąć te reszty, kliknij Actions > Atoms/Bonds i Delete.
   4. Użyj opcji File > Save PDB, aby zapisać zmodyfikowaną strukturę w osobnym pliku PDB, który będzie wykorzystany w kolejnych krokach protokołu (tutaj użyto nazwy IRF5_interface.pdb).

Rysunek 3: Wybór receptora i liganda w serwerze BUDE Alanine Scan. (A) Graficzna reprezentacja kompleksu białko-białko. (B) Pole tekstowe do wprowadzenia nazwy zadania w serwerze. (C) Panel do interaktywnego wyboru łańcuchów, które zostaną uznane za receptor i ligand (szczegóły w kroku 1.2). Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rycina 4: Strona wyników serwera BUDE Alanine Scan. Potencjalny najlepiej oddziałujący segment w sekwencji ligandu zaznaczono czerwoną ramką. W lewym panelu reszta o wyższym przewidywanym udziale energetycznym (Leu433) jest wyróżniona na zielono. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

2. Dywersyfikacja sekwencji

UWAGA: W poniższych krokach rosetta_main odnosi się do głównego katalogu instalacyjnego Rosetta, który zazwyczaj znajduje się w /opt/rosetta_src__bundle/main/, gdzie oznacza zainstalowaną wersję Rosetta. Zakłada się również, że aplikacje Rosetta są dostępne w całym systemie; jeśli tak nie jest, należy podać pełną ścieżkę do plików wykonywalnych. W przypadku kompilacji ze źródeł, pliki wykonywalne te znajdują się w katalogu /rosetta_main/source/bin/.

1. Wstępna optymalizacja łańcuchów bocznych aminokwasów
   1. Skopiuj edytowaną strukturę zapisaną w kroku 1.3.4 do lokalizacji w systemie Linux dostępnej dla aplikacji Rosetta.
   2. Użyj aplikacji FixBB z pakietu Rosetta, aby przeprowadzić repakowanie (repack) wszystkich łańcuchów bocznych aminokwasów struktury bazowej przed dywersyfikacją sekwencji. W tej operacji orientacja wszystkich łańcuchów bocznych aminokwasów jest optymalizowana w celu zminimalizowania energii i poprawy stabilności kompleksu. Aby to zrobić, wykonaj następujące polecenie:
      
      UWAGA: To polecenie generuje plik PDB nazwany zgodnie z oryginalną strukturą z dodatkowym sufiksem numerycznym (IRF5_interface_0001.pdb w tym przykładzie).
   3. Aby ułatwić kolejny krok protokołu, zmień nazwę repakowanego pliku PDB, dodając sufiks _repack, korzystając z następującego polecenia:
      mv IRF5_interface_0001.pdb IRF5_repack.pdb
2. Dywersyfikacja sekwencji
   1. Uruchom Pepspec w trybie projektowania (design mode), aby przeprowadzić właściwy etap dywersyfikacji sekwencji, używając następującego polecenia:
      
      Poniżej znajdują się opcje ogólne:
      • -s wskazuje plik wejściowy (repakowany plik PDB wygenerowany w kroku 2.1.3).
      • -o wskazuje prefiks dla nazw plików wyjściowych.
      • - database wskazuje ścieżkę do głównej bazy danych Rosetta 3.
      • -ex1, -ex2 oraz extrachi_cutoff są opcjami biblioteki rotamerów (szczegóły znajdują się w dokumentacji Pepspec).
      • -overwrite instruuje aplikację, aby nadpisała ewentualne istniejące wyniki wygenerowane w poprzednich iteracjach.
      Poniżej znajdują się opcje odnoszące się bezpośrednio do dywersyfikacji sekwencji:
      • -pepspec:pep_chain wskazuje łańcuchy PDB uznane za ligand ('b' w tym przykładzie).
      • -pepspec:native_pep_anchor wskazuje resztę aminokwasową użytą jako kotwicę (w tym przykładzie resztę Leu w pozycji 10 peptydu ligandu).
      • -pepspec:n_peptides określa liczbę struktur peptydowych do wygenerowania.
      • -pepspec:no_prepack_prot instruuje aplikację, aby pominęła repakowanie w wejściowej strukturze bazowej (ponieważ zostało to wykonane wcześniej w kroku 2.1).
        UWAGA: Głównym wynikiem działania Pepspec jest katalog zawierający pliki PDB dla peptydów powstałych w fazie projektowania, nazwane przy użyciu prefiksu wyjściowego z sufiksem .pdbs (IRF5.pdbs w przykładzie). Dodatkowo Pepspec generuje wszystkie zaakceptowane sekwencje peptydowe przetestowane w kroku dywersyfikacji sekwencji oraz odpowiadające im wyniki energii Rosetta w pliku tekstowym z rozdzielaczem tabulacyjnym, nazwanym zgodnie z prefiksem wyjściowym z sufiksem .spec (IRF5.spec w przykładzie). Ponieważ protokół opisany w niniejszej pracy ma na celu oszacowanie preferencji aminokwasowych, a nie właściwe projektowanie peptydów, kolejne kroki wykorzystują plik IRF5.spec zamiast struktur PDB w katalogu .pdbs .

3. Ocena preferencji aminokwasowych

1. Obliczanie PWM
   1. Aby wygenerować PWM, należy użyć skryptu gen_pepspec_pwm.py dołączonego do pakietu Rosetta. Aby uruchomić ten skrypt, należy użyć następującego polecenia:
      
      gdzie:
      • IRF5.spec to plik wyjściowy Pepspec wygenerowany w kroku 2.2.
      • -1 wskazuje, że w sekwencji nie ma dodatkowych reszt N-końcowych, w związku z czym pozycje w PWM są indeksowane od 1.
      • 0.2 instruuje skrypt, aby uwzględnił tylko 20% najlepiej ocenianych peptydów z wyniku Pepspec (wartość domyślna to 0.1, co odpowiada 10%)
      • interface_score instruuje skrypt, aby uszeregował peptydy na podstawie wyniku interfejsu (interface score), który jest jedną z różnych punktacji Rosetta zawartych w pliku wyjściowym Pepspec.
        UWAGA: Skrypt ten generuje dwa pliki wyjściowe: jeden dla obliczonego PWM (z rozszerzeniem .pwm) oraz drugi dla sekwencji podzbioru peptydów użytych do obliczenia PWM (z rozszerzeniem .seq). Nazwy tych plików zawierają również wynik oraz ułamek peptydów użytych do rankingu. W tym przykładzie pliki te nazywają się odpowiednio IRF5_interface_score_0.2.pwm oraz IRF5_interface_score_0.2.seq.
2. Generowanie logo sekwencji
   1. Przejdź do serwera WebLogo (https://weblogo.berkeley.edu/logo.cgi)²¹ i kliknij przycisk Choose File obok opcji Upload Sequence Data. Prześlij plik z sekwencjami peptydów wygenerowany w kroku 3.1.1 (w tym przykładzie IRF5_interface_score_0.2.seq).
   2. Wybierz żądany format i rozmiar logo zgodnie z długością danych wejściowych. W przykładzie zastosowano format PDF i rozmiar 15 cm x 12 cm. Kliknij Create Logo.

Wyniki

W niniejszym artykule opisano protokół służący do przewidywania preferencji aminokwasowych powierzchni wiążącej IRF5, białka należącego do rodziny czynników transkrypcyjnych znanych jako ludzkie czynniki regulacyjne interferonu. Białka te są regulatorami wrodzonych i nabytych odpowiedzi immunologicznych i uczestniczą w różnicowaniu oraz aktywacji wielu komórek odpornościowych. Podjednostki IRF posiadają wysoce plastyczne i wielospecyficzne powierzchnie wiążące, dzięki czemu są zdolne do tworzenia homodimerów, heterodimer...

Dyskusja

W niniejszym artykule opisano protokół służący do oszacowania preferencji aminokwasowych potencjalnie wielospecyficznych miejsc wiązania w oparciu o dywersyfikację sekwencji in silico. Opracowano niewiele narzędzi obliczeniowych do oszacowania preferencji aminokwasowych granic faz białko-peptyd 14,25,26. Narzędzia te mają charakter predykcyjny, ale różnią się algorytmami obliczeniowymi używanymi do wykonywania ich przewidywań i ...

Materiały

Lista materiałów użytych w tym artykule

Nazwa	Firma	Numer katalogowy	Komentarze
BUDE Alanine Scan Server	Uniwersytet w Edynburgu	https://pragmaticproteindesign.bio.ed.ac.uk/balas/	doi: 10.1021/acschembio.9b00560
Rosetta Modeling Software	Rosetta Commons	https://www.rosettacommons.org/software	doi: 10.1002/prot.22851
UCSF Chimera	Uniwersytet Kalifornijski San Francisco	https://www.cgl.ucsf.edu/chimera/	doi: 10.1002/ JCC.20084