Name: Modelowanie farmakoforów dla celów z obszernymi bibliotekami ligandów: studium przypadku SARS-CoV-2 Mpro (Wideo) | JoVE Journal
Uploaded: 2025-09-26T14:30:00.000+00:00
Duration: 5 min 50 s
Description: Poznaj modelowanie farmakoforów dla SARS-CoV-2 Mpro przy użyciu dużych bibliotek ligandów. Ta metoda obliczeniowa identyfikuje kluczowe cechy interakcji molekularnych, usprawniając wirtualne przesiewanie i racjonalny projekt leków.

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

Niniejsze badanie zostało częściowo sfinansowane przez PAPIIT UNAM IV200121 (M.A.V-V.), projekt secihti Cátedras CONACYT 639 (L.C-B. i M.A.V-V.), LANCAD-UNAM-DGTIC-386 (L.C-B.) oraz Zakład Badań i Rozwoju w Bioterapeutykach (UDIBI).
Autorzy wyrażają wdzięczność Eduardo Orozco za jego cenne wsparcie i wnikliwy wkład w prace nad tym manuskryptem.

1. Metoda 1
<ol>
	<li>Przygotowanie ligandów do generowania farmakoforów konsensusowych<ol>
		<li>Dopasuj wszystkie kompleksy białko-ligand za pomocą oprogramowania PyMOL11.</li>
		<li>Wyodrębnij każdy wyrównany konformer liganda i zapisz go jako osobny plik w formacie SDF. UWAGA: Inne formaty, takie jak MOL, MOL2 i PDB, mogą być również używane dla opisanego tutaj protokołu.</li>
	</ol>
	</li>
	<li>Generowanie plików JSON farmakoforów przy użyciu programu Pharmit12
	<ol>
		<li>Prześlij każdy plik liganda indywidualnie do Pharmit, korzystając z opcji Load Features (patrz Spis materiałów , aby uzyskać link do Pharmit).</li>
		<li>Użyj opcji Zapisz sesję , aby pobrać odpowiedni plik JSON farmakoforu.</li>
	</ol>
	</li>
	<li>Organizowanie plików JSON do użycia w ConPhar<ol>
		<li>Przechowuj wszystkie pobrane pliki JSON w jednym folderze. Pliki te zostaną przesłane do środowiska Google Colab w następnej metodzie.</li>
	</ol>
	</li>
</ol>
2. Metoda 2
<ol>
	<li>Konfigurowanie środowiska Google Colab<ol>
		<li>Uruchom nowy notatnik Google Colab: Otwórz Google Colab w przeglądarce internetowej i utwórz nowy notatnik oraz dostosuj ustawienia, aby korzystać z wcześniejszej wersji, wybierając opcję Środowisko uruchomieniowe → Zmień wersję środowiska uruchomieniowego → 2025.07.</li>
		<li>Zainstaluj narzędzia Conda i PyMOL. Kod wymagany do zainstalowania Conda i PyMOL w Google Colab jest dostępny tutaj w formie kursywy w celach informacyjnych, a pełny skrypt wykonywalny jest również dostępny w pliku pomocniczym 1. 
		# Zainstaluj CondaColab: Włącza obsługę środowiska Conda w Google Colab 
		z IPython.utils importuj operacje we/wy 
		importowanie tqdm.notebook 
		Importowanie systemu operacyjnego 
		Razem = 100 z tqdm.notebook.tqdm(total=total) jako pbar: za pomocą io.capture_output() w postaci, w jakiej została przechwycona: # Zainstaluj CondaColab !install -q condacolab Import Condacolab condacolab.install() pbar.aktualizacja(10) # Zaktualizuj ścieżkę Pythona, aby zlokalizować zainstalowane pakiety 
		Import sys sys.path.append('/usr/local/lib/python3.7/site-packages/') pbar.aktualizacja(20) # Zainstaluj pakiet PyMOL za pomocą mamba z kanału Schrödinger %shell mamba install -c schrodinger pymol-bundle --yes pbar.aktualizacja(90)</li>
		<li>Sprawdź pomyślne wykonanie: Uruchom komórkę, klikając ikonę odtwarzania lub naciskając Shift + Enter. Po pomyślnym wykonaniu nad komórką pojawi się poziomy zielony pasek (patrz Rysunek 1).</li>
	</ol>
	</li>
	<li>Zainstaluj pakiet ConPhar Python i zaimportuj wymagane moduły<ol>
		<li>Zainstaluj i zaimportuj ConPhar. Kod wymagany do zainstalowania pakietu ConPhar i zaimportowania niezbędnych modułów do Google Colab jest tutaj podany kursywą w celach informacyjnych, a pełny skrypt wykonywalny jest również dostępny w pliku pomocniczym 1. 
		# Zainstaluj pakiet do analizy farmakoforów ConPhar import pymol !install conphar z konfaru. Farmakofory importują parse_json_pharmacophore, show_pharmacophoric_descriptors, save_pharmacophore_to_pymol, save_pharmacophore_to_json, compute_concensus_pharmacophore Importowanie systemu operacyjnego Importowanie pand jako PD UWAGA: Instrukcja importu z conphar. Farmakofory są podzielone na wiele linii dla jasności, ale powinny być wprowadzone jako pojedyncza ciągła linia. Narzędzie ConPhar (https://github.com/AngelRuizMoreno/ConcensusPharmacophore) jest aktywnie rozwijane. Obecny protokół korzysta ze stabilnej wersji (0.1.2), która została zweryfikowana pod kątem opisanej procedury. Użytkownicy są zachęcani do korzystania z tej wersji, aby zapewnić odtwarzalność.</li>
		<li>Potwierdź pomyślną instalację: Uruchom komórkę, klikając ikonę odtwarzania lub naciskając Shift + Enter. Po pomyślnej instalacji i zaimportowaniu narzędzi ConPhar pojawi się komunikat potwierdzający (patrz Rysunek 2).</li>
	</ol>
	</li>
	<li>Wczytywanie indywidualnych modeli farmakoforów z plików JSON<ol>
		<li>Utwórz folder dla plików JSON farmakoforu. Kod wymagany do utworzenia folderu do przechowywania plików wejściowych farmakoforu w Google Colab jest tutaj podany kursywą w celach informacyjnych, a pełny skrypt wykonywalny jest również dostępny w pliku pomocniczym 1. 
		# Utwórz folder do przechowywania wejściowych plików JSON 
		os.makedirs("JSON_FOLDER", exist_ok=Prawda) UWAGA: To polecenie tworzy folder automatycznie, jeśli jeszcze nie istnieje.</li>
		<li>Prześlij pliki JSON do folderu: Kliknij ikonę folderu w lewym panelu Colab, otwórz nowo utworzony folder i kliknij prawym przyciskiem myszy, aby wybrać Prześlij. Dodaj wymagane pliki JSON (patrz rysunek 3). UWAGA: Upewnij się, że pliki są zgodne z oczekiwanym formatem wygenerowanym przez firmę Pharmit.</li>
	</ol>
	</li>
	<li>Analizowanie i konsolidowanie funkcji farmakoforycznych<ol>
		<li>Wyodrębnij cechy farmakoforyczne z przesłanych plików. Kod wymagany do przeanalizowania przesłanych plików JSON, wyodrębnienia cech farmakoforycznych i zapisania ich w jednej ramce danych w Google Colab jest dostępny tutaj w formie kursywy w celach informacyjnych, a pełny skrypt wykonywalny jest również dostępny w pliku pomocniczym 1. 
		p4_table=pd. Ramka danych() Dla pliku w OS.listdir('/content/JSON_FOLDER'): Jeśli '.json' w pliku: próbować: p4,lig,rec=parse_json_pharmacophore(f"/treść/JSON_FOLDER/{plik}") P4['ligand']=file.replace('.json','') p4_table=pd.concat([p4_table,p4],ignore_index=Prawda) z wyjątkiem wyjątku: przechodzić p4_table</li>
		<li>Uruchom komórkę, klikając ikonę odtwarzania lub naciskając Shift + Enter , aby wygenerować skonsolidowaną ramkę danych (patrz Rysunek 4). Powstała w ten sposób skonsolidowana ramka danych kompiluje wszystkie cechy farmakoforyczne wyodrębnione z poszczególnych ligandów w ujednoliconą tabelę, ułatwiając grupowanie i analizę statystyczną. UWAGA: Skrypt zawiera podstawową obsługę wyjątków w celu pominięcia zniekształconych plików JSON podczas przetwarzania, aby zapobiec przerwom w przepływie pracy. Biorąc pod uwagę dużą liczbę plików, wcześniejsze sprawdzenie każdego z nich jest niepraktyczne; Zamiast tego skrypt można zmodyfikować tak, aby wydrukował nazwę dowolnego pliku, którego nie można załadować, dzięki czemu użytkownik może go sprawdzić i poprawić indywidualnie.</li>
	</ol>
	</li>
	<li>Wygeneruj i zapisz konsensus Farmakofor<ol>
		<li>Wyświetlić wszystkie deskryptory farmakoforyczne. Kod wymagany do wizualizacji deskryptorów farmakoforycznych wyodrębnionych z plików wejściowych w Google Colab jest dostępny tutaj w formie kursywy w celach informacyjnych, a pełny skrypt wykonywalny jest również dostępny w pliku pomocniczym 1. 
		show_pharmacophoric_descriptors(p4_table)</li>
		<li>Uruchom komórkę, klikając ikonę odtwarzania lub naciskając Shift + Enter , aby zwizualizować grupowanie cech farmakoforycznych (patrz rysunek 5). Wyniki grupowane w klastry cech grupują podobne cechy farmakoforyczne w wielu ligandach na podstawie ich pozycji przestrzennej, umożliwiając identyfikację konserwatywnych wzorców interakcji.</li>
		<li>Zapisz model farmakoforu w formacie PyMOL. Kod wymagany do wygenerowania i zapisania modelu farmakoforu konsensusu w formacie zgodnym z PyMOL w Google Colab jest dostępny tutaj w formie kursywy w celach informacyjnych, a pełny skrypt wykonywalny jest również dostępny w pliku pomocniczym 1. 
		save_pharmacophore_to_pymol(p4_table, out_file='ConPhar_pymol.pse')</li>
		<li>Uruchom komórkę, aby utworzyć odpowiedni plik .pse (patrz Rysunek 6).</li>
		<li>Zapisz model farmakoforu w formacie JSON. Kod wymagany do wygenerowania i zapisania modelu farmakoforu konsensusu w formacie zgodnym z Pharmit w Google Colab jest dostępny tutaj w formie kursywy w celach informacyjnych, a pełny skrypt wykonywalny jest również dostępny w pliku pomocniczym 1. ​save_pharmacophore_to_json(p4_table,out_file= "ConPhar_pharmit.json")</li>
		<li>Uruchom komórkę, aby utworzyć odpowiedni plik .json (patrz Rysunek 6).</li>
		<li>Generuj dane wyjściowe z klastrami cech i dendrogramy. Kod wymagany do wygenerowania plików farmakoforowych konsensusu pogrupowanych według typu funkcji, w tym formatów kompatybilnych z PyMOL i Pharmit oraz wizualizacji dendrogramu w Google Colab jest dostępny tutaj w formie kursywy w celach informacyjnych, a pełny skrypt wykonywalny jest również dostępny w pliku pomocniczym 1. concensus,links=compute_concensus_pharmacophore (p4_table,save_data_per_descriptor=Prawda,out_folder='/treść')</li>
		<li>Uruchom komórkę, klikając ikonę odtwarzania lub naciskając Shift + Enter (patrz Rysunek 7). Wyniki wyświetlają wizualizacje dendrogramu, które w postaci drzewa reprezentują hierarchiczne relacje między zgrupowanymi cechami farmakoforycznymi, pomagając użytkownikom interpretować stopień podobieństwa i bliskość przestrzenną między klastrami. UWAGA: Ten krok zapisuje wiele plików wyjściowych i figur w określonym folderze.</li>
		<li>Eksportuj wyniki farmakoforu konsensusu do pliku CSV. Kod wymagany do zapisania ostatecznej tabeli farmakoforów konsensusu w formacie CSV w Google Colab jest tutaj dostępny w formie kursywy w celach informacyjnych, a pełny skrypt wykonywalny jest również dostępny w pliku pomocniczym 1. 
		concensus.to_csv('consensus_result.csv', index=Fałsz)</li>
		<li>Uruchom komórkę, klikając ikonę odtwarzania lub naciskając Shift + Enter , aby wyeksportować dane farmakoforu konsensusu do pliku CSV o nazwie consensus_result.csv w celu dalszej analizy (patrz Rysunek 8). UWAGA: Plik CSV ułatwia dalsze zastosowania, takie jak analiza statystyczna lub wizualizacja w arkuszu kalkulacyjnym.</li>
	</ol>
	</li>
	<li>Wykorzystanie konsensusu Farmakofor do wirtualnych badań przesiewowych<ol>
		<li>Eksportuj farmakofor konsensusu w formacie JSON. Upewnij się, że plik zawiera współrzędne przestrzenne dla wszystkich obiektów. Ten format jest bezpośrednio kompatybilny z internetowymi narzędziami do badań przesiewowych, takimi jak Pharmit12. Aby przeprowadzić wirtualne badanie przesiewowe, prześlij plik JSON farmakoforu na serwer Pharmit za pomocą opcji Load Features . Serwer umożliwia udoskonalanie modelu poprzez modyfikowanie, dodawanie lub usuwanie funkcji, a także sprawdzanie pod kątem dużych bibliotek złożonych, takich jak PubChem13lub ZINC14. UWAGA: Wszystkie skrypty wymagane do odtworzenia analiz opisanych w tym badaniu są dostępne w pliku pomocniczym 1 i można uzyskać do nich interaktywny dostęp za pośrednictwem https://colab.research.google.com/drive/1Zj8DjazV6mJM-we9iOpoGTZt8JfW4I7D?usp=sharing. Czytelnicy są zachęcani do stworzenia osobistej kopii notatnika Colab w celu uruchamiania i modyfikowania analiz bez zmiany oryginalnego skryptu.</li>
	</ol>
	</li>
</ol>
<img alt="Figure 1" class="xfigimg" src="/files/ftp_upload/68933/68933fig1.jpg"> 
Rysunek 1: Zrzut ekranu przedstawiający pierwszy blok kodu używany do skonfigurowania środowiska Google Colab przez zainstalowanie CondaColab. Poziomy niebieski pasek postępu pojawia się nad komórką podczas wykonywania i zmienia kolor na zielony po pomyślnym zakończeniu. <a href="https://www.jove.com/files/ftp_upload/68933/68933fig1large.jpg" target="_blank">Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.</a>
<img alt="Figure 2" class="xfigimg" src="/files/ftp_upload/68933/68933fig2.jpg"> 
Rysunek 2: Zrzut ekranu drugiego bloku kodu używanego do instalacji pakietu ConPhar w Google Colab. <a href="https://www.jove.com/files/ftp_upload/68933/68933fig2large.jpg" target="_blank">Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.</a>
<img alt="Figure 3" class="xfigimg" src="/files/ftp_upload/68933/68933fig3.jpg"> 
Rysunek 3: Zrzut ekranu przedstawiający pliki JSON farmakoforu pomyślnie przesłane do wyznaczonego folderu w Google Colab. <a href="https://www.jove.com/files/ftp_upload/68933/68933fig3large.jpg" target="_blank">Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.</a>
<img alt="Figure 4" class="xfigimg" src="/files/ftp_upload/68933/68933fig4.jpg"> 
Rysunek 4: Zrzut ekranu przedstawiający kod używany do wyodrębniania cech farmakoforycznych z plików JSON i konsolidowania ich w jedną ramkę danych. <a href="https://www.jove.com/files/ftp_upload/68933/68933fig4large.jpg" target="_blank">Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.</a>
<img alt="Figure 5" class="xfigimg" src="/files/ftp_upload/68933/68933fig5.jpg"> 
Rysunek 5: Wizualizacja grupowania cech farmakoforycznych wygenerowanych podczas budowy modelu farmakoforu konsensusu. <a href="https://www.jove.com/files/ftp_upload/68933/68933fig5large.jpg" target="_blank">Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.</a>
<img alt="Figure 6" class="xfigimg" src="/files/ftp_upload/68933/68933fig6.jpg"> 
Rysunek 6: Zrzut ekranu przedstawiający pomyślne wygenerowanie plików farmakoforowych konsensusu w formatach kompatybilnych z PyMOL i Pharmit. <a href="https://www.jove.com/files/ftp_upload/68933/68933fig6large.jpg" target="_blank">Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.</a>
<img alt="Figure 7" class="xfigimg" src="/files/ftp_upload/68933/68933fig7.jpg"> 
Rysunek 7: Zrzut ekranu przedstawiający generowanie plików farmakoforów konsensusu pogrupowanych według typu funkcji, w tym danych wyjściowych dla wizualizacji PyMOL, Pharmit i dendrogramu. <a href="https://www.jove.com/files/ftp_upload/68933/68933fig7large.jpg" target="_blank">Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.</a>
<img alt="Figure 8" class="xfigimg" src="/files/ftp_upload/68933/68933fig8.jpg"> 
Rysunek 8: Zrzut ekranu przedstawiający pomyślny eksport danych farmakoforu konsensusu do pliku consensus_result.csv. <a href="https://www.jove.com/files/ftp_upload/68933/68933fig8large.jpg" target="_blank">Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.</a>

<ol><li>Leach, A. R., Gillet, V. J., Lewis, R. A., Taylor, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Three-dimensional+pharmacophore+methods+in+drug+discovery%5BTitle%5D)+AND+%22J+Med+Chem%22%5BJournal%5D)">Three-dimensional pharmacophore methods in drug discovery.</a> J Med Chem. 53 (2), 539-558 (2010).</li>
<li>Giordano, D., Biancaniello, C., Argenio, M. A., Facchiano, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Drug+design+by+pharmacophore+and+virtual+screening+approach%5BTitle%5D)+AND+%22Pharmaceuticals+%28Basel%29%22%5BJournal%5D)">Drug design by pharmacophore and virtual screening approach.</a> Pharmaceuticals (Basel). 15 (5), 646(2022).</li>
<li>Schaller, D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Next+generation+3D+pharmacophore+modeling%5BTitle%5D)+AND+%22WIREs+Comput+Mol+Sci%22%5BJournal%5D)">Next generation 3D pharmacophore modeling.</a> WIREs Comput Mol Sci. 10 (4), e1468(2020).</li>
<li>Wolber, G., Langer, T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((LigandScout%3A+3-D+pharmacophores+derived+from+protein-bound+ligands+and+their+use+as+virtual+screening+filters%5BTitle%5D)+AND+%22J+Chem+Inf+Model%22%5BJournal%5D)">LigandScout: 3-D pharmacophores derived from protein-bound ligands and their use as virtual screening filters.</a> J Chem Inf Model. 45 (1), 160-169 (2005).</li>
<li>Ekins, S., Mestres, J., Testa, B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((In+silico+pharmacology+for+drug+discovery%3A+Methods+for+virtual+ligand+screening+and+profiling%5BTitle%5D)+AND+%22Br+J+Pharmacol%22%5BJournal%5D)">In silico pharmacology for drug discovery: Methods for virtual ligand screening and profiling.</a> Br J Pharmacol. 152 (1), 9-20 (2007).</li>
<li>Wolber, G., Dornhofer, A. A., Langer, T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Efficient+overlay+of+small+organic+molecules+using+3D+pharmacophores%5BTitle%5D)+AND+%22J+Comput+Aided+Mol%22%5BJournal%5D)">Efficient overlay of small organic molecules using 3D pharmacophores.</a> J Comput Aided Mol. 20 (12), 773-788 (2006).</li>
<li>Schuster, D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((3D+pharmacophores+as+tools+for+activity+profiling%5BTitle%5D)+AND+%22Drug+Discov+Today+Technol%22%5BJournal%5D)">3D pharmacophores as tools for activity profiling.</a> Drug Discov Today Technol. 7 (4), e205-e211 (2010).</li>
<li>Qing, X., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Pharmacophore+modeling%3A+advances%2C+limitations%2C+and+current+utility+in+drug+discovery%5BTitle%5D)+AND+%22J+Recept+Ligand+Channel+Res%22%5BJournal%5D)">Pharmacophore modeling: advances, limitations, and current utility in drug discovery.</a> J Recept Ligand Channel Res. 7, 81-92 (2014).</li>
<li>Ullrich, S., Nitsche, C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+SARS-CoV-2+main+protease+as+drug+target%5BTitle%5D)+AND+%22Bioorg+Med+Chem+Lett%22%5BJournal%5D)">The SARS-CoV-2 main protease as drug target.</a> Bioorg Med Chem Lett. 30 (17), 127377(2020).</li>
<li>Ruiz-Moreno, A. J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Consensus+pharmacophore+strategy+for+identifying+novel+SARS-CoV-2+Mpro+inhibitors+from+large+chemical+libraries%5BTitle%5D)+AND+%22J+Chem+Inf+Model%22%5BJournal%5D)">Consensus pharmacophore strategy for identifying novel SARS-CoV-2 Mpro inhibitors from large chemical libraries.</a> J Chem Inf Model. 64 (6), 1984-1995 (2024).</li>
<li>Schrödinger, L., DeLano, W. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aSchr%C3%B6dinger%2C+L+%22PyMOL%2C+2020%22">PyMOL, 2020</a>. , http://www.pymol.org/pymol (2020).</li>
<li>Sunseri, J., Koes, D. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Pharmit%3A+Interactive+exploration+of+chemical+space%5BTitle%5D)+AND+%22Nucleic+Acids+Res%22%5BJournal%5D)">Pharmit: Interactive exploration of chemical space.</a> Nucleic Acids Res. 44 (W1), W442-W448 (2016).</li>
<li>Kim, S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((PubChem+2025+update%5BTitle%5D)+AND+%22Nucleic+Acids+Res%22%5BJournal%5D)">PubChem 2025 update.</a> Nucleic Acids Res. 53 (D1), D1516-D1525 (2025).</li>
<li>Tingle, B. I., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((ZINC-22%3A+A+free+multi-billion-scale+database+of+tangible+compounds+for+ligand+discovery%5BTitle%5D)+AND+%22J+Chem+Inf+Model%22%5BJournal%5D)">ZINC-22: A free multi-billion-scale database of tangible compounds for ligand discovery.</a> J Chem Inf Model. 63 (4), 1166-1176 (2023).</li>
<li>Hazemann, J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Accelerating+the+hit-to-lead+optimization+of+a+SARS-CoV-2+Mpro+inhibitor+series+by+combining+high-throughput+medicinal+chemistry+and+computational+simulations%5BTitle%5D)+AND+%22J+Med+Chem%22%5BJournal%5D)">Accelerating the hit-to-lead optimization of a SARS-CoV-2 Mpro inhibitor series by combining high-throughput medicinal chemistry and computational simulations.</a> J Med Chem. 68 (8), 8269-8294 (2025).</li>
<li>Momin, Y., Beloshe, V. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Pharmacophore+modeling+in+drug+design%5BTitle%5D)+AND+%22Adv+Pharmacol%22%5BJournal%5D)">Pharmacophore modeling in drug design.</a> Adv Pharmacol. 103, 313-324 (2025).</li>
<li>Zdrazil, B., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+ChEMBL+database+in+2023%3A+A+drug+discovery+platform+spanning+multiple+bioactivity+data+types+and+time+periods%5BTitle%5D)+AND+%22Nucleic+Acids+Res%22%5BJournal%5D)">The ChEMBL database in 2023: A drug discovery platform spanning multiple bioactivity data types and time periods.</a> Nucleic Acids Res. 52 (D1), D1180-D1192 (2024).</li>
<li>Zhou, L., Griffith, R., Gaeta, B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Combining+spatial+and+chemical+information+for+clustering+pharmacophores%5BTitle%5D)+AND+%22BMC+Bioinformatics%22%5BJournal%5D)">Combining spatial and chemical information for clustering pharmacophores.</a> BMC Bioinformatics. 15 (Suppl 16), S5(2014).</li>
<li>Guner, O. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((History+and+evolution+of+the+pharmacophore+concept+in+computer-aided+drug+design%5BTitle%5D)+AND+%22Curr+Top+Med+Chem%22%5BJournal%5D)">History and evolution of the pharmacophore concept in computer-aided drug design.</a> Curr Top Med Chem. 2 (12), 1321-1332 (2002).</li>
<li>Wieder, M., Perricone, U., Seidel, T., Boresch, S., Langer, T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Comparing+pharmacophore+models+derived+from+crystal+structures+and+from+molecular+dynamics+simulations%5BTitle%5D)+AND+%22Monatsh+Chem%22%5BJournal%5D)">Comparing pharmacophore models derived from crystal structures and from molecular dynamics simulations.</a> Monatsh Chem. 147 (3), 553-563 (2016).</li>
<li>Santos, R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+comprehensive+map+of+molecular+drug+targets%5BTitle%5D)+AND+%22Nat+Rev+Drug+Discov%22%5BJournal%5D)">A comprehensive map of molecular drug targets.</a> Nat Rev Drug Discov. 16 (1), 19-34 (2017).</li>
</ol>

Autorzy deklarują brak konfliktu interesów.

Modelowanie farmakoforów pomaga w racjonalnym odkrywaniu leków16. Szczegółowo opisany tutaj protokół został wykorzystany do wygenerowania konsensusu farmakoforu dla głównej proteazy SARS-CoV-2 (Mpro), zaczynając od dziesiątek inhibitorów niekowalencyjnych, ale można go zastosować do dowolnego celu biologicznego o znanych konformacjach związanych z ligandem, zwłaszcza tych z obszernymi zestawami danych ligandów. W porównaniu z indywidualnymi modelami farmakoforowymi lub modelami opartymi wyłącznie na ligandach, to podejście oparte na konsensusie zwiększa niezawodność i zmniejsza błąd systematyczny poprzez integrację wielu konformacji ligandów. Konsensusowy model farmakoforowy wygenerowany przez ten protokół może być stosowany do różnych bibliotek molekularnych, w tym do publicznych baz danych, takich jak ChEMBL17i ZINC14, kolekcji komercyjnych i zastrzeżonych wewnętrznych bibliotek związków. Ta wszechstronność umożliwia eksplorację szerokich i zróżnicowanych przestrzeni chemicznych w wirtualnych kampaniach przesiewowych, zwiększając w ten sposób prawdopodobieństwo identyfikacji nowych związków bioaktywnych. W ten sposób ten przepływ pracy wspiera identyfikację nowych związków bioaktywnych o pożądanych profilach interakcji.
Jednym z krytycznych etapów tego protokołu jest dokładne przygotowanie i wyrównanie kompleksów ligand-cel. Precyzyjne wstępne wyrównanie ligandów w ich bioaktywnych konformacjach jest niezbędne do zapewnienia znaczącej i niezawodnej integracji cech, ponieważ niewspółosiowość może osłabić moc predykcyjną modelu10. Aby osiągnąć optymalne wyrównanie, zalecamy wykonanie superpozycji strukturalnej przy użyciu konserwatywnych reszt miejsca wiązania białek lub atomów szkieletu. W przypadku zestawów danych z różnymi ligandami wyrównanie kluczowych punktów interakcji lub cech farmakoforycznych zwiększa spójność w całym zestawie. Ponadto zaleca się staranną ręczną kontrolę i korektę konformacji ligandów odstających, aby zapobiec artefaktom, które mogłyby zmniejszyć dokładność modelu. Ponadto otwarty charakter ConPhar ułatwia integrację z innymi przepływami pracy obliczeniowej, a także zaawansowane korekty parametrów klastrowania. Doświadczeni użytkownicy mogą modyfikować parametry grupowania cech farmakoforycznych, aby uzyskać równowagę między ogólnością a swoistością modelu. Modyfikacje te mogą być kierowane przez rozmiar zestawu danych ligandów, różnorodność chemiczną lub złożoność docelową. Na przykład, gdy mamy do czynienia z bardzo zróżnicowanymi ligandami, stopniowe grupowanie lub analiza podgrup mogą lepiej uchwycić istotne cechy18. 
Narzędzie ConPhar jest obecnie aktywnie rozwijane. Zapewnia to znaczną elastyczność dla użytkowników, którzy chcą dostosować narzędzie do konkretnych potrzeb badawczych, w tym integrację z uzupełniającymi potokami oprogramowania. Ponieważ jednak oprogramowanie może nadal ulegać zmianom, zachęca się użytkowników do monitorowania oficjalnych aktualizacji i zapoznawania się z dokumentacją specyficzną dla wersji. Zgłaszanie napotkanych problemów przyczyni się do jego ciągłego ulepszania i szerszej walidacji przez społeczność. Nasz protokół zbuduje modele z kompleksowymi wzorcami interakcji, a tym samym będą mogły być wykorzystane do wirtualnego badania przesiewowego dużych bibliotek. W naszych reprezentatywnych wynikach wykorzystano otrzymany farmakofor do identyfikacji dwóch nowych potencjalnych ligandów dla regionu katalitycznego Mpro. Ostatecznie ta metoda farmakoforu opartego na konsensusie usprawnia identyfikację trafień i przyspiesza optymalizację potencjalnych klientów.
Identyfikujemy jednak następujące ograniczenia protokołu: i) zależność od wysokiej jakości struktur związanych z ligandami, ponieważ dane o niskiej rozdzielczości mogą zmniejszyć dokładność modelu; ii) ryzyko przeoczenia unikalnych i ważnych interakcji kluczowych dla określonych ligandów, ponieważ farmakofor konsensusu kładzie nacisk na cechy wspólne; oraz (iii) na odtwarzalność mogą mieć wpływ różne definicje cech w różnych narzędziach programowych19. Ograniczenia te można jednak przezwyciężyć poprzez udoskonalenie struktury układu białko-ligand za pomocą symulacji dynamiki molekularnej20, utrzymanie klastrów w oparciu o wiedzę użytkownika oraz konsekwentne stosowanie oprogramowania do identyfikacji cech farmakoforycznych.
Chociaż konsensusowy model farmakoforowy w tym protokole nie został zweryfikowany przy użyciu klasycznych wskaźników binarnych, jego zdolność predykcyjna została wykazana w poprzednim badaniu10. W tej pracy wykorzystano niezależny zestaw testowy 78 chemicznie zróżnicowanych ligandów do walidacji uniwersalnego farmakoforu wygenerowanego przy użyciu tej samej metodologii. Zestaw walidacyjny zawierał ligandy o podobieństwie ≤0,5 Tanimoto, masach cząsteczkowych od 200 do 700 g/mol, wiązaniach obrotowych ≤17 i co najmniej trzech cechach farmakoforycznych. Model z powodzeniem zidentyfikował znane substancje aktywne, co potwierdza jego przydatność w wirtualnych przepływach pracy związanych z badaniami przesiewowymi oraz jego niezawodność w przypadku chemicznie zróżnicowanych kandydatów.
Przedstawiony tutaj protokół umożliwia budowę modeli farmakoforowych, które ujmują kompleksowe wzorce interakcji charakterystyczne dla celu. Jest to szczególnie cenne w przypadku celów z obszernymi danymi ligandowymi, takich jak proteazy wirusowe, kinazy i receptory jądrowe21. Wygenerowane modele nadają się do wirtualnych badań przesiewowych dużych bibliotek związków, ułatwiając identyfikację nowych kandydatów bioaktywnych. Podkreślając konserwatywne i funkcjonalnie istotne cechy, podejście oparte na farmakoforach opartych na konsensusie wspiera zarówno wykrywanie trafień na wczesnym etapie, jak i późniejszą optymalizację potencjalnych klientów.

Modelowanie farmakoforowe jest podstawową techniką w komputerowym wspomaganiu projektowania leków, umożliwiającą identyfikację istotnych cech molekularnych odpowiedzialnych za aktywność biologiczną 1,2. Farmakor określa przestrzenny układ cech, takich jak donory wiązań wodorowych, akceptory, pierścienie aromatyczne i regiony hydrofobowe, wymagane do interakcji molekularnej między ligandem a celem biologicznym3. Gdy dostępnych jest wiele kompleksów ligand-cel, zarówno z krystalografii eksperymentalnej, jak i modelowania molekularnego, wspólne wzorce interakcji można zintegrować z modelami farmakoforowymi konsensusu, zwiększając niezawodność modelu i zwiększając dokładność wirtualnych badań przesiewowych 4,5.
Pomimo użyteczności farmakoforów opartych na konsensusie, ich wytwarzanie pozostaje wyzwaniem technicznym, zwłaszcza gdy ligandy są zróżnicowane strukturalnie. Dostępnych jest kilka narzędzi programowych do generowania farmakoforów z poszczególnych ligandów, ale brakuje ustandaryzowanych procedur integracji wielu zestawów funkcji w spójny model 6,7. Co więcej, niewiele platform oferuje kompleksowe przepływy pracy, które obsługują grupowanie, filtrowanie i eksportowanie funkcji w formatach kompatybilnych z narzędziami do przesiewania i wizualizacji8. Ograniczenia te ograniczyły szerokie zastosowanie farmakoforów konsensusu w wysiłkach związanych z odkrywaniem leków na dużą skalę.
Aby wypełnić tę lukę, ConPhar został opracowany jako nowatorskie narzędzie typu open source, zaprojektowane specjalnie do systematycznej ekstrakcji, grupowania i modelowania konsensusu cech farmakoforycznych z obszernych zestawów wstępnie wyrównanych kompleksów ligand-cel. W przeciwieństwie do istniejącego oprogramowania, ConPhar oferuje elastyczne dostrajanie parametrów, automatyczną integrację funkcji i kompatybilność z wieloma formatami wyjściowymi, ułatwiając generowanie solidnych modeli konsensusu odpowiednich dla wirtualnych potoków przesiewowych. W ten sposób narzędzie to pozwala przezwyciężyć wcześniejsze wąskie gardła w obsłudze dużych i zróżnicowanych chemicznie bibliotek ligandów, zwiększając odtwarzalność i skalowalność w przepływach pracy związanych z modelowaniem farmakoforów.
Aby zademonstrować to podejście, przedstawiamy powtarzalny protokół konstruowania konsensusowych modeli farmakoforowych na podstawie zbiorów kompleksów ligand-cel. Przepływ pracy integruje narzędzia typu open source do ekstrakcji cech farmakoforów, grupowania, wizualizacji i dalszych zastosowań. W ramach studium przypadku zastosowaliśmy protokół do głównej proteazy SARS-CoV-2 (Mpro), krytycznego celu terapeutycznego z obszernymi danymi strukturalnymi9. Opracowaliśmy zestaw danych 100 ligandów niekowalencyjnych skrystalizowanych wspólnie z Mpro (wpisy PDB na dzień 27 maja 2025 r.), z wyłączeniem form apo i redundantnych kompleksów. Poszczególne cechy farmakoforyczne zostały wyodrębnione i połączone w model konsensusu za pomocą ConPhar, narzędzia zaprojektowanego specjalnie do grupowania cech10. Ten przypadek ilustruje zdolność protokołu do ujawniania konserwatywnych wzorców interakcji i wspierania racjonalnego badania przesiewowego ultradużych bibliotek.

<table><tbody><tr><td>Współpraca z Google</td><td>Google LLC</td><td>https://colab.google</td><td>Platforma internetowa do uruchamiania notatników w języku Python; używany tutaj do wykonywania skryptów ConPhar i PyMOL do modelowania farmakoforów konsensusu.</td></tr><tr><td>Przemysł farmaceutyczny</td><td>Laboratorium Koesa na Uniwersytecie w Pittsburghu</td><td>https://pharmit.csb.pitt.edu/search.html</td><td>Platforma internetowa do interaktywnych wirtualnych badań przesiewowych farmakoforów i eksploracji przestrzeni chemicznej, umożliwiająca przesyłanie plików ligandów i generowanie funkcji farmakoforu</td></tr><tr><td>PyMOL</td><td>Schr&ouml; dinger, Inc.</td><td>https://www.pymol.org</td><td>Oprogramowanie do wizualizacji molekularnej i wyrównywania strukturalnego używane do białek&ndash; Superpozycja i analiza zespolonego liganda</td></tr><tr><td>Obsługa skryptu przepływu pracy Pharmacophore File 1</td><td>Google LLC</td><td>https://colab.research.google.com/drive/1Zj8DjazV6mJM-we9iOpoGTZt8JfW4I7D?usp=sharing</td><td>Pełny skrypt wykonywalny do instalacji Conda, PyMOL i uruchamiania przepływu pracy farmakoforu w Google Colab. Dostarczany również jako plik pomocniczy 1.</td></tr></tbody></table>

modelowanie farmakoforów dla celów z obszernymi bibliotekami ligandów: studium przypadku sars-cov-2 mpro

Pharmacophore Modeling for Targets with Extensive Ligand Libraries: A Case Study on SARS-CoV-2 Mpro

Farkofor określa przestrzenny układ cech molekularnych wymaganych do optymalnych interakcji między związkiem a jego celem biologicznym. Modele te można uzyskać, analizując interakcje międzycząsteczkowe między celem a zestawem znanych ligandów w ich konformacjach wiązania. Farmakofor oparty na konsensusie integruje wspólne cechy wielu ligandów, zmniejszając błąd systematyczny modelu i zwiększając moc predykcyjną. Jednak wygenerowanie solidnego farmakoforu opartego na konsensusie z dużego i zróżnicowanego chemicznie zestawu ligandów stanowi wyzwanie techniczne.
W tym miejscu przedstawiamy protokół budowy farmakoforów konsensusu przy użyciu ConPhar, narzędzia informatycznego typu open source zaprojektowanego do identyfikacji i grupowania cech farmakoforycznych w wielu kompleksach związanych z ligandami. Protokół obejmuje generowanie, udoskonalanie i stosowanie modeli do wirtualnego badania przesiewowego ultradużych bibliotek molekularnych. Jako studium przypadku zastosowaliśmy metodę do głównej proteazy SARS-CoV-2 (Mpro), używając stu niekowalencyjnych inhibitorów skrystalizowanych z celem. Uzyskany w ten sposób model farmakoforowy pozwolił uchwycić kluczowe cechy interakcji w regionie katalitycznym Mpro i umożliwił identyfikację nowych potencjalnych ligandów.
Strategia ta ma szerokie zastosowanie do każdego celu biologicznego, dla którego dostępne są konformacje związane z ligandem. Jest to szczególnie cenne dla celów z obszernymi zestawami danych ligandów i wspiera racjonalne odkrywanie leków poprzez usprawnienie identyfikacji nowych kandydatów o pożądanych profilach interakcji.

Sto kompleksów Mpro skrystalizowanych z różnymi niekowalencyjnymi inhibitorami zostało wyrównanych, jak pokazano na rysunku 9A. Każdy ligand został następnie wyodrębniony jako pojedynczy plik (rysunek 9B), a następnie przesłany na serwer Pharmit. Sesja została zapisana w celu wygenerowania odpowiedniego pliku JSON (patrz Rysunek 9C).
<img alt="Figure 9" class="xfigimg" src="/files/ftp_upload/68933/68933fig9.jpg"> 
Rysunek 9: Przygotowanie ligandów do modelowania farmakoforów w drodze konsensusu. (A) Strukturalne wyrównanie trzech reprezentatywnych kompleksów Mpro skrystalizowanych z niekowalencyjnymi inhibitorami. (B) Superpozycja 100 ligandów uwzględnionych w tym badaniu. (C) Przesłanie każdego liganda na serwer Pharmit za pomocą opcji "Załaduj funkcje", a następnie wygenerowanie pliku JSON za pomocą opcji "Zapisz sesję". Obie opcje są oznaczone czerwonymi strzałkami. <a href="https://www.jove.com/files/ftp_upload/68933/68933fig9large.jpg" target="_blank">Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.</a>
Kompletny zestaw plików JSON wykorzystano do wygenerowania modelu farmakoforu konsensusu, który zawierał 1450 cech farmakoforycznych pogrupowanych w 110 klastrów: 23 aromatyczne (Aro), 30 akceptorów wiązań wodorowych (HBA), 16 donorów wiązań wodorowych (HBD), 36 hydrofobowych (Hyd) i 5 ujemnie naładowanych (Anion) klastrów (Rysunek 10A, Tabela uzupełniająca 1). Zidentyfikowaliśmy większe klastry dla każdej funkcji w celu ich włączenia do modelu konsensusu. Spośród klastrów Aro wybraliśmy te, które zawierają ponad 20 elementów. Z grup HBA, HBD i Hyd zachowaliśmy klastry z co najmniej 50 elementami. Żaden z klastrów anionów nie miał więcej niż 4 członków, więc zostały one wyłączone z modelu konsensusu (rysunek 10B). Te progi liczbowe (&gt;20 członków dla Aro; ≥50 dla HBA, HBD i Hyd) zostały empirycznie zdefiniowane na podstawie rozkładu rozmiarów klastrów, aby nadać priorytet najbardziej konserwatywnym i gęsto zaludnionym wzorcom interakcji.
<img alt="Figure 10" class="xfigimg" src="/files/ftp_upload/68933/68933fig10.jpg"> 
Rycina 10: Grupowanie cech farmakoforycznych ze 100 kompleksów Mpro-ligand. (A) Rozkład 1450 cech farmakoforycznych: aromatyczny (Aro) w kolorze fioletowym, akceptor wiązań wodorowych (HBA) w kolorze pomarańczowym, donor wiązań wodorowych (HBD) w kolorze białym, hydrofobowy (Hyd) w kolorze zielonym i ujemnie naładowany (Anion) w kolorze czerwonym. (B) Klastry zawierające ≥10 członków są pokazane dla Aro, HBA, HBD i Hyd; wyświetlane są wszystkie klastry anionów. Największe klastry - zdefiniowane jako te, które mają ≥20 elementów dla Aro i ≥50 członków dla HBA, HBD i Hyd oraz ≥4 członków dla anionów - są wyróżnione kropkowanymi kulami. Liczba członków w każdym z największych klastrów jest wyraźnie wskazana. <a href="https://www.jove.com/files/ftp_upload/68933/68933fig10large.jpg" target="_blank">Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.</a>
Tak więc nasz konsensusowy model farmakoforowy obejmował 11 cech: 3 cechy aromatyczne (Aro), Aro 1, 2 i 3, pochodzące z klastrów mających odpowiednio 20, 24 i 45 członków; 4 akceptory wiązań wodorowych (HBA), HBA 1, 2, 3 i 4, z klastrów liczących odpowiednio 51, 55, 77 i 81 członków; 2 donory wiązań wodorowych (HBD) z klastrów liczących odpowiednio 51 i 71 członków; i 2 cechy hydrofobowe (Hyd) z gromad liczących odpowiednio 50 i 52 członków (patrz rysunek 11A).
Model farmakoforowy został wykorzystany do przeszukania bazy danych PubChem przy użyciu strategii opisanej wcześniej10, ale nie znaleziono żadnych dopasowań. Aby zwiększyć elastyczność wyszukiwania, usunięto cechę aromatyczną Aro 1 – pochodzącą z najmniejszego i najmniej reprezentatywnego klastra. Ta korekta doprowadziła do zidentyfikowania dwóch trafień, w tym konformera związku z CID PubChem 101267741 i 10285538 (patrz rysunek 11B). Strukturę chemiczną 2D 101267741 złożonego przedstawiono na rysunku 11C. Co ciekawe, ten zidentyfikowany związek dobrze mieści się w kieszeni na wiązania Mpro. Analiza porównawcza z kokrystalizowanym ligandem 38a15 (PDB ID: 9HAJ) pokazuje, że konformer pasujący do farmakoforu 101267741 jest zakopany jeszcze głębiej niż 38a w podkieszeniach S1 i S2 Mpro. Jednak 38a zajmuje szerszy obszar kieszeni, głównie ze względu na boczny łańcuch etylokarboksyamidowy, który rozciąga się w kierunku kieszeni S1′. Niemniej jednak, kluczowe oddziaływania międzycząsteczkowe 38a obejmują inne ugrupowania, a nie ten łańcuch boczny. Natomiast 101267741 złożony tworzy 11 oddziaływań międzycząsteczkowych, w tym siedem wiązań wodorowych i cztery kontakty hydrofobowe. W porównaniu z trzema wiązaniami wodorowymi i dwoma oddziaływaniami hydrofobowymi utworzonymi przez 38a, tryb wiązania związku 101 wydaje się silniejszy (ryc. 11D).
<img alt="Figure 11" class="xfigimg" src="/files/ftp_upload/68933/68933fig11.jpg"> 
Rysunek 11: Wyszukiwanie farmakoforów w bazie danych PubChem. (A) Kompletny model farmakoforu konsensusu zawiera 11 cech: 3 aromatyczne (fioletowy), 4 akceptory wiązań wodorowych (pomarańczowy), 2 donory wiązań wodorowych (biały) i 2 hydrofobowe (zielony). (B) Zidentyfikowany konformer związku PubChem 101267741 zgodny ze zredukowanym modelem farmakoforu po usunięciu cechy Aro 1 oraz (C) Zidentyfikowano dwuwymiarową strukturę chemiczną związku. (D) Analiza porównawcza sposobów wiązania związku 38a (magenta; od PDB ID: 9HAJ) i związku PubChem CID 101267741 (zielony) do kieszeni katalitycznej SARS-CoV-2 Mpro. Odpowiednie oddziaływania międzycząsteczkowe dla każdego związku są pokazane po prawej stronie. Wiązania wodorowe i oddziaływania hydrofobowe są pokazane jako linie przerywane odpowiednio w kolorze niebieskim i szarym. <a href="https://www.jove.com/files/ftp_upload/68933/68933fig11large.jpg" target="_blank">Kliknij tutaj, aby zobaczyć większą wersję tego rysunku.</a>
Plik pomocniczy 1: Skrypt przepływu pracy farmakoforu. <a href="https://www.jove.com/files/ftp_upload/68933/Supporting%20File%201%20Pharmacophore%20Workflow%20Script.txt" target="_blank">Kliknij tutaj, aby pobrać ten plik.</a>
Tabela 1: Grupowanie Podsumowanie cech farmakoforycznych. Liczba klastrów, pręty na grupę, współrzędne środka ciężkości i promienie klastra <a href="https://www.jove.com/files/ftp_upload/68933/Table_S1.xlsx" target="_blank">Kliknij tutaj, aby pobrać ten plik.</a>

W artykule przedstawiono protokół konstruowania modelu farmakoforu konsensusu poprzez integrację cech molekularnych z wielu ligandów. Metoda ta ma zastosowanie w pracach nad odkrywaniem leków ukierunkowanych na dowolny cel biologiczny o znanych konformacjach związanych z ligandem, umożliwiając identyfikację kluczowych cech interakcji dla wirtualnych badań przesiewowych i racjonalnego projektowania leków.

Modelowanie farmakoforów dla celów z obszernymi bibliotekami ligandów: studium przypadku SARS-CoV-2 Mpro

A pharmacophore defines the spatial arrangement of molecular features required for optimal interactions between a compound and its biological target. ...

Nasze badania koncentrują się na modelowaniu farmakoforów z dużymi bibliotekami ligandów. W tym badaniu skupiamy się na SARS-CoV-2 Mpro, aby zidentyfikować wzorce strukturalne, które kierują odkrywaniem leków. Najnowsze postępy w modelowaniu farmakoforów wykorzystują większe i bardziej zróżnicowane zestawy danych, które zwiększają dokładność identyfikacji kluczowych czynników molekularnych do odkrywania leków.Symulacje dynamiki molekularnej i uczenie maszynowe są obecnie wykorzystywane do rozwoju modelowania farmakoforów i poprawy wyników odkrywania leków. Kluczowym wyzwaniem jest zarządzanie dużymi i różnorodnymi zestawami danych ligandów, przy jednoczesnej poprawie dokładności ekstrakcji w przyszłości i obniżeniu kosztów obliczeniowych. Wykazaliśmy, że drogie biblioteki ligandów poprawiają przyszłą dokładność farmakoforów, umożliwiając bardziej wiarygodne prognozy i lepszą identyfikację potencjalnych leków.Aby rozpocząć, uruchom oprogramowanie PyMOL. Dopasuj wszystkie kompleksy ligandów białkowych w oprogramowaniu. Dla każdego wyrównanego kompleksu wyodrębnij konformer liganda i zapisz go jako osobny plik w formacie SDF.Prześlij każdy plik liganda indywidualnie, aby go farmić, korzystając z opcji Załaduj funkcje. Po załadowaniu użyj opcji Zapisz sesję, aby pobrać plik farmakoforowy w formacie JSON. Przechowuj wszystkie pobrane pliki JSON farmakoforu razem w jednym folderze do późniejszego wykorzystania w środowisku Google Colab.Uruchom Google Colab w przeglądarce internetowej i utwórz nowy notatnik. Zainstaluj Conda i PyMOL, a następnie uruchom komórkę, klikając ikonę odtwarzania lub naciskając Shift i Enter, aby zweryfikować pomyślne wykonanie. Nad komórką pojawi się zielony pasek, jeśli zostanie wykonany poprawnie.Teraz zainstaluj pakiet CANFAR Python i zaimportuj wymagane moduły. Uruchom komórkę, klikając ikonę odtwarzania, aby potwierdzić pomyślną instalację. Po pomyślnym wykonaniu pojawi się komunikat z potwierdzeniem.Następnie utwórz nowy folder do przechowywania plików JSON farmakoforu. Otwórz folder, kliknij w nim prawym przyciskiem myszy i wybierz prześlij, aby przesłać pliki JSON. Aby przeanalizować i skonsolidować cechy farmakoforu, wyodrębnij cechy ze wszystkich przekazanych plików.Kliknij ikonę odtwarzania, aby uruchomić komórkę i wygenerować skonsolidowaną ramkę danych zawierającą cechy z każdego pojedynczego ligandu. Wyświetl wszystkie deskryptory farmakoforów, a następnie uruchom komórkę, aby wyświetlić wzorce grupowania. Teraz zapisz model farmakoforu w formacie PyMOL.Uruchom komórkę, aby wygenerować plik pse, a następnie zapisz model jako plik JSON, uruchamiając odpowiednią komórkę kodu. Generuj dane wyjściowe klastrowane funkcji i dendogramy przy użyciu dostarczonego kodu. Kliknij przycisk odtwarzania i wyeksportuj dane farmakoforu konsensusu do pliku CSV, wykonując polecenie eksportu.Wykorzystaj wygenerowany model farmakoforu konsensusu do wirtualnych badań przesiewowych. 100 głównych kompleksów proteazy lub Mpro skrystalizowanych z różnymi niekowalencyjnymi inhibitorami zostało ustawionych w celu wizualizacji różnorodności ligandów w kieszeni wiążącej. Każdy ligand został indywidualnie wyekstrahowany z wyrównanych kompleksów w celu wygenerowania oddzielnych plików do wirtualnych badań przesiewowych.Wyekstrahowane ligandy zostały przesłane na serwer farmaceutyczny i wykorzystane do wygenerowania ustrukturyzowanych plików JSON do modelowania farmakoforów. Pliki JSON dostarczyły 1450 cech farmakoforowych, pogrupowanych w 110 klastrów, w tym 23 aromatyczne, 30 akceptorów wiązań wodorowych, 16 donorów wiązań wodorowych, 36 klastrów hydrofobowych i 5 klastrów anionowych. Aby zbudować model konsensusu, zachowano tylko najgęściej zaludnione klastry, podczas gdy klastry anionowe zostały wykluczone z powodu niewystarczającej liczby członków.Ostateczny model farmakoforu konsensusu składał się z 11 cech, w tym trzech aromatycznych, dla akceptorów wiązań wodorowych, dwóch donorów wiązań wodorowych i dwóch cech hydrofobowych. Po usunięciu najmniej reprezentatywnej cechy aromatycznej, aro pierwszej, zidentyfikowano dwa związki kandydujące. Dwuwymiarowa struktura chemiczna związku 101 267 741 została zobrazowana w celu potwierdzenia jego zgodności z modelem farmakoforowym.Związek 101 267 741 wpasowywał się głęboko w kieszeń wiążącą Mpro, zajmując kieszenie podrzędne S1 i S2 skuteczniej niż ligand referencyjny 38a, który rozprzestrzeniał się w kierunku kieszeni S1. Związek 101 267 741 utworzył 11 oddziaływań międzycząsteczkowych, siedem wiązań wodorowych i cztery kontakty hydrofobowe.

pharmacophore-modeling-for-targets-with-extensive-ligand-libraries

Research

JoVE Journal

Biochemistry

Biochemia

To badanie miało na celu przedstawienie strategii identyfikacji interakcji lek-peptyd. Strategia obejmuje biopanowanie krótkich peptydów rozpoznających leki w oparciu o bioczujnik mikrowagi kwarcowej (QCM), a następnie analizę bioinformatyczną w celu ilościowej oceny uzyskanych informacji w celu rozpoznania leku i adnotacji miejsc wiązania leku na białkach.

Oparta na biosensorach strategia wysokoprzepustowej analizy biopanningowej i bioinformatycznej do globalnej walidacji interakcji lek-białko

Ten artykuł opisuje cały proces XChem do badania fragmentów na bazie kryształów, począwszy od złożenia wniosku o dostęp i wszystkich kolejnych kroków do rozpowszechniania danych.

Osiągnięcie efektywnego przesiewania fragmentów w zakładzie XChem w Diamond Light Source

Kluczową umiejętnością w modelowaniu biomolekularnym jest wyświetlanie i opisywanie aktywnych miejsc w białkach. Technika ta została zademonstrowana przy użyciu czterech popularnych darmowych programów do wizualizacji makromolekularnej: iCn3D, Jmol, PyMOL i UCSF ChimeraX.

Modelowanie miejsca aktywnego enzymu za pomocą darmowego oprogramowania do wizualizacji molekularnej

Ten protokół wprowadza dostępne narzędzia do modelowania ligandów małocząsteczkowych na mapach krioEM makromolekuł.

Modelowanie ligandów na mapach pochodzących z kriomikroskopii elektronowej

Eksperyment wykorzystany tutaj pokazuje metodę molekularnego dokowania w połączeniu z testem przesunięcia termicznego komórkowego w celu przewidywania i walidacji interakcji między małymi cząsteczkami a celami białkowymi.

Przewidywanie docelowego białka i walidacja związku małocząsteczkowego

Przedstawiamy tutaj protokół strukturalnie ukierunkowanego projektowania interakcji wiązania białko-ligand między Sanglifehriną A a Cyklofiliną A/Ganoderiol-F, niezbędny w odkrywaniu nowych leków. Stabilność zespołu ligand-receptor została oceniona za pomocą symulacji dynamiki molekularnej.

Struktura-kierowana projektowanie i opracowywanie nowych inhibitorów cyklofiliny A oraz pochodnych ganoderiol-F: podejście In-Silico

Metody obliczeniowe obiecują przyspieszyć odkrywanie leków, jednak często pomijają dynamiczną naturę struktur białkowych. W tym miejscu omawiamy analizę dokowania opartą na zespołach, aby pośrednio uwzględnić elastyczność białek, potencjalnie poprawiając dokładność i niezawodność wysiłków związanych z odkrywaniem leków.

Wykorzystanie elastyczności i dynamiki struktury białka docelowego w obliczeniowym odkrywaniu leków przy użyciu analizy dokowania opartej na zespołach

W tym badaniu wykorzystano strategie in-silico, aby zidentyfikować enumerowaną etrawirynę jako obiecujący środek terapeutyczny dla HIV. Nasze odkrycia dotyczące interakcji molekularnych i dynamiki wspierają racjonalne projektowanie nowych NNRTI jako możliwych alternatyw leczenia HIV.

Ilościowa zależność struktura-aktywność, przewidywanie aktywności i dynamika molekularna nienukleotydowych inhibitorów odwrotnej transkryptazy

Ten dwujęzyczny protokół zapewnia obliczeniowy proces odkrywania leków, oceniający interakcje białko-ligandy kinaz Polopodobnych 1 do 3 (PLK1–3) oraz właściwości absorpcji, dystrybucji, metabolizmu, wydalania, toksyczności i stabilności (ADMET-S) naturalnych cząsteczek pochodzących z bazy danych.

Dwujęzyczny proces obliczeniowy do identyfikacji potencjalnych inhibitorów PLK1 w amerykańskim języku migowym i angielskim

Ten protokół wykorzystuje zarówno koekspresję podjednostek, jak i mieszanie podjednostek postlizy w celu dokładniejszego zbadania złożenia rekombinowanego proteasomu.

Badanie składania proteasomu za pomocą rekombinowanych proteasomów archeonów i niedenaturujących PAGE: Argumenty przemawiające za podejściem łączonym

Poznaj modelowanie farmakoforów dla SARS-CoV-2 Mpro przy użyciu dużych bibliotek ligandów. Ta metoda obliczeniowa identyfikuje kluczowe cechy interakcji molekularnych, usprawniając wirtualne przesiewanie i racjonalny projekt leków.

Modelowanie farmakoforów dla celów z obszernymi bibliotekami ligandów: studium przypadku SARS-CoV-2 Mpro

Przeglądaj więcej filmów