Zrozumienie reprezentacji świata molekularnego jest kluczowe dla zdobycia wiedzy eksperckiej w naukach biomolekularnych1, ponieważ interpretacja takich obrazów jest niezbędna do zrozumienia funkcji biologicznych2. Pierwsze zetknięcie ucznia z makrocząsteczkami zazwyczaj następuje poprzez dwuwymiarowe ilustracje w podręcznikach przedstawiające błony komórkowe, organelle, makrocząsteczki itp., jednak rzeczywistość biologiczna polega na tym, że są to struktury trójwymiarowe, a zrozumienie ich właściwości wymaga sposobów wizualizacji i wyciągania wniosków z modeli 3D.
W związku z tym szczególną uwagę poświęcono rozwojowi kompetencji w zakresie wizualizacji biomolekularnej w zaawansowanych kursach nauk o życiu na poziomie molekularnym, a w licznych artykułach wskazywano na znaczenie i trudności związane z nauczaniem oraz ocenianiem umiejętności wizualizacji1,3,4,5,6,7,8,9. Odpowiedzią na te publikacje był wzrost liczby interwencji dydaktycznych, zazwyczaj prowadzonych w ciągu jednego semestru w jednej instytucji, w których programy wizualizacji molekularnej i modele są wykorzystywane do wyjaśniania trudnych koncepcji2,10,11,12,13,14,15. Ponadto badacze starali się scharakteryzować sposób, w jaki studenci wykorzystują programy i/lub modele wizualizacji biomolekularnej do zgłębiania konkretnych zagadnień16,17,18,19. Nasza grupa, BioMolViz, opisała ramy koncepcyjne (Framework), które dzielą nadrzędne motywy kompetencji wizualnych na cele i zamierzenia kształcenia, aby wspomagać tego typu interwencje20,21; prowadzimy również warsztaty szkolące kadrę naukową w zakresie stosowania tych ram w procesie projektowania wstecznego (backward design) ocen służących do pomiaru umiejętności wizualizacji22.
Kluczową umiejętnością we wszystkich tych badaniach jest zdolność do manipulowania strukturami makrocząsteczek przy użyciu programów do wizualizacji biomolekularnej. Narzędzia te zostały opracowane niezależnie na różnych platformach, dlatego ich obsługa i zastosowanie mogą być dość specyficzne. Wymaga to instrukcji dostosowanych do konkretnego programu, a zidentyfikowanie oprogramowania, z którym użytkownik czuje się komfortowo, jest istotne dla ułatwienia jego dalszego wdrażania.
Poza podstawowymi operacjami manipulowania strukturami w 3D (obracaniem, zaznaczaniem i modyfikowaniem modelu), głównym celem jest modelowanie miejsca aktywnego białka. Proces ten pozwala uczącemu się rozwinąć zrozumienie trzech nadrzędnych zagadnień opisanych w BioMolViz Framework: oddziaływań molekularnych, ligandów / modyfikacji oraz zależności między strukturą a funkcją20,21.
Cztery popularne programy do wizualizacji biomolekularnej to: Jmol/JSmol23, iCn3D24, PyMOL25 oraz UCSF Chimera26,27. Osobowom rozpoczynającym pracę z programem Chimera zalecamy korzystanie z UCSF ChimeraX, czyli nowej generacji programu do wizualizacji molekularnej Chimera, która jest obecnie wspieraną wersją tego oprogramowania.
W niniejszym protokole demonstrujemy, jak wykorzystać każdy z tych czterech programów do modelowania miejsca aktywnego ludzkiej glukokinazy w kompleksie z analogiem substratu (PDB ID: 3FGU) oraz do prezentowania pomiarów ilustrujących specyficzne oddziaływania wiążące28. Model ten przedstawia kompleks katalityczny enzymu. Aby uchwycić miejsce aktywne w stanie przedkatalitycznym, do miejsca aktywnego glukokinazy przyłączono niehydrolizowalny analog ATP. Ten ester kwasu fosfoaminofosfonowego z adenylanem (ANP) zawiera wiązanie fosfor-azot zamiast zwykłego wiązania fosfor-tlen w tej pozycji. Miejsce aktywne zawiera również glukozę (oznaczoną w modelu jako BCG) i magnez (oznaczony jako MG). Dodatkowo w strukturze znajduje się jon potasu (K), będący wynikiem zastosowania chlorku potasu w rozpuszczalniku do krystalizacji. Jon ten nie jest krytyczny dla funkcji biologicznej i znajduje się poza miejscem aktywnym.

Rycina 1: Struktury ATP/ANP. Struktura trifosforanu adenozyny (ATP) w porównaniu z estrem adenylanowym kwasu fosfoaminofosfonowego (ANP). Kliknij tutaj, aby wyświetlić większą wersję tej ryciny.
Protokół demonstruje wybór związanych ligandów kompleksu analogu substratu oraz identyfikację reszt centrum aktywnego znajdujących się w odległości 5 Å od związanego kompleksu, co pozwala wyodrębnić aminokwasy i cząsteczki wody zdolne do tworzenia istotnych oddziaływań molekularnych, w tym oddziaływań hydrofobowych i sił van der Waalsa.
Wyświetlanie jest początkowo skonfigurowane tak, aby większość białka była przedstawiona w formie schematycznej (cartoon), natomiast aminokwasy centrum aktywnego w formie pałeczkowej (stick), co pozwala na ukazanie odpowiednich atomów białka i wyróżnienie oddziaływań molekularnych. Po kroku 3 protokołu dla każdego programu zastosowano te reprezentacje, a widok białka jest podobny we wszystkich programach (Rycina 2). Na koniec protokołu schemat białka zostaje ukryty, aby uprościć widok i skupić się na centrum aktywnym.

Rycina 2: Porównanie struktur w różnych programach. Porównanie struktury 3FGU w każdym programie po etapie Dostosowania Reprezentacji (krok 2 lub 3 każdego protokołu). Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.
Kolorowanie CPK stosuje się do aminokwasów centrum aktywnego oraz związanych ligandów29,30. Ten schemat kolorystyczny pozwala odróżnić atomy różnych pierwiastków chemicznych w modelach cząsteczkowych przedstawionych w reprezentacjach liniowej, prętowej, kulkowo-prętowej oraz wypełniającej przestrzeń. W schemacie kolorowania CPK wodór jest biały, azot niebieski, tlen czerwony, siarka żółta, a fosfor pomarańczowy. Tradycyjnie do węgla używa się koloru czarnego, choć w nowoczesnych zastosowaniach kolorowanie węgla może być zróżnicowane.
Atomy wodoru nie są widoczne w strukturach krystalicznych, mimo że każdy z tych programów potrafi przewidzieć ich położenie. Dodanie atomów wodoru do dużej struktury makromolekularnej może ograniczyć widoczność, dlatego nie są one wyświetlane w tym protokole. W związku z tym wiązania wodorowe w tych strukturach będą przedstawiane poprzez pomiar odległości między centrami dwóch heteroatomów (np. tlenu do tlenu, tlenu do azotu).
Przegląd programów
Graficzne interfejsy użytkownika (GUI) do pobrania: PyMOL (wersja 2.4.1), ChimeraX (wersja 1.2.5) oraz Jmol (wersja 1.8.0_301) to narzędzia do modelowania molekularnego oparte na GUI. Te trzy interfejsy posiadają wiersze poleceń do wprowadzania kodu; wiele z tych samych funkcji jest dostępnych za pomocą menu i przycisków w GUI. Wspólną cechą wiersza poleceń w tych programach jest możliwość wczytywania i ponownego wykonywania poprzednich poleceń za pomocą klawiszy strzałek w górę i w dół na klawiaturze.
Interfejsy graficzne oparte na przeglądarce internetowej: iCn3D (I-see-in-3D) to przeglądarka oparta na WebGL służąca do interaktywnego wyświetlania trójwymiarowych struktur makromolekularnych i związków chemicznych w sieci, bez konieczności instalacji osobnej aplikacji. Nie korzysta ona z wiersza poleceń, choć pełna wersja internetowa posiada edytowalny dziennik poleceń. JSmol to wersja Jmol w języku JavaScript lub HTML5 przeznaczona do użytku na stronie internetowej lub w oknie przeglądarki, której działanie jest bardzo zbliżone do Jmol. JSmol może być wykorzystywany do tworzenia tutoriali online, w tym animacji.
Proteopedia31,32, FirstGlance in Jmol33 oraz interfejs internetowy JSmol (JUDE) w Milwaukee School of Engineering Center for BioMolecular Modeling są przykładami takich środowisk projektowych online opartych na Jmol34. Wiki Proteopedia to narzędzie dydaktyczne, które umożliwia użytkownikowi modelowanie struktury makrocząsteczki i tworzenie stron zawierających te modele w obrębie witryny35. Narzędzie do tworzenia scen w Proteopedia, zbudowane przy użyciu JSmol, integruje GUI z dodatkowymi funkcjami niedostępnymi w GUI Jmol.
Programy Jmol i iCn3D opierają się na języku programowania Java; JSmol wykorzystuje język Java lub HTML5, natomiast PyMOL i ChimeraX opierają się na języku programowania Python. Każdy z tych programów ładuje pliki z bazy danych białek, które można pobrać z RCSB Protein Data Bank, posługując się 4-cyfrowym alfanumerycznym identyfikatorem PDB ID36,37. Najczęstszymi typami plików są pliki Protein Data Bank (PDB) z rozszerzeniem .pdb oraz pliki Crystallographic Information File (CIF lub mmCIF) z rozszerzeniem .cif. Format CIF zastąpił PDB jako domyślny typ pliku w Protein Data Bank, jednak oba formaty plików funkcjonują w tych programach. Przy użyciu plików CIF w porównaniu do plików PDB mogą wystąpić niewielkie różnice w sposobie wyświetlania sekwencji/struktury; pliki te działają jednak podobnie, a różnice te nie będą tutaj szczegółowo omawiane. Molecular Modeling Database (MMDB), produkt National Center for Biotechnology Information (NCBI), jest podzbiorem struktur PDB, do których przypisano informacje kategoryczne (np. cechy biologiczne, konserwatywne domeny białkowe)38. Program iCn3D, produkt NCBI, umożliwia ładowanie plików PDB zawierających dane MMDB.
Aby wyświetlić model, użytkownik może pobrać odpowiedni plik z dedykowanej strony Protein Data Bank dla danej struktury (np. https://www.rcsb.org/structure/3FGU), a następnie użyć rozwijanego menu Plik programu, aby otworzyć strukturę. Wszystkie programy umożliwiają również ładowanie pliku struktury bezpośrednio przez interfejs; metoda ta została szczegółowo opisana w protokołach.
Interfejsy graficzne programów ChimeraX, Jmol i PyMOL zawierają jedno lub więcej okien konsoli, których rozmiar można zmienić, przeciągając narożnik. Programy iCn3D i JSmol są w całości osadzone w przeglądarce internetowej. Podczas korzystania z iCn3D użytkownik może potrzebować przewinięcia zawartości okien wyskakujących, aby wyświetlić wszystkie elementy menu, w zależności od rozmiaru i rozdzielczości ekranu.
Szczegółowe protokoły opisane w niniejszym dokumencie przedstawiają prostą metodę wizualizacji centrum aktywnego enzymu przy użyciu każdego z programów. Należy zauważyć, że w każdym programie istnieje wiele sposobów na wykonanie poszczególnych kroków. Na przykład w programie ChimeraX to samo zadanie można zrealizować za pomocą menu rozwijanych, paska narzędzi u góry lub wiersza poleceń. Użytkownicy zainteresowani szczegółowym poznaniem konkretnego programu są zachęcani do zapoznania się z dostępnymi w sieci samouczkami, podręcznikami i stronami Wiki dla tych programów39,40,41,42,43,44,45,46.
Istniejące instrukcje i samouczki do tych programów przedstawiają elementy zawarte w niniejszym protokole jako odrębne zadania. Aby wyświetlić centrum aktywne, użytkownik musi samodzielnie zsyntetyzować wymagane operacje na podstawie różnych instrukcji i samouczków. Niniejsza publikacja uzupełnia dostępne samouczki, prezentując liniowy protokół modelowania oznakowanego centrum aktywnego z oddziaływaniami molekularnymi, dostarczając użytkownikowi logiki modelowania centrum aktywnego, którą można zastosować do innych modeli i programów.

Rycina 3: Interfejs graficzny (GUI) programu ChimeraX. Interfejs graficzny programu ChimeraX z opisanymi menu rozwijanymi, paskiem narzędzi, przeglądarką struktur oraz linią komend. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 4: Interfejs graficzny (GUI) iCn3D. Interfejs GUI iCn3D z oznaczonymi menu rozwijanymi, paskiem narzędzi, przeglądarą struktury, dziennikiem poleceń, wyskakującym oknem zestawów zaznaczeń oraz wyskakującymi menu sekwencji i adnotacji. Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.

Rycina 5: Interfejs graficzny Jmol. Interfejs graficzny Jmol z oznaczonymi menu rozwijanymi, paskiem narzędzi, przeglądarą struktur, menu kontekstowym oraz konsolą/linią komend. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 6: Interfejs graficzny PyMOL. Interfejs graficzny PyMOL z oznakowanymi menu rozwijanymi, przeglądarą struktur, panelem nazw/obiektów, menu sterowania myszą oraz wierszem poleceń. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.