Name: Kuratorstwo obliczeniowych bibliotek chemicznych zademonstrowanych za pomocą alfa-aminokwasów (Wideo) | JoVE Journal
Uploaded: 2022-04-13T13:30:00.000+00:00
Duration: 8 min 21 s
Description: Poznaj metody kuratorowania obliczeniowych bibliotek chemicznych z wykorzystaniem oprogramowania open-source. Protokół umożliwia generowanie struktur alfa-aminokwasów i ich przesiewanie in silico.

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

MAY dziękuje za sfinansowanie przez Fundację Carl-Zeiss. Wszystkie rysunki zostały wygenerowane przy użyciu programu Microsoft Excel.

UWAGA: Podsumowanie protokołu znajduje się na Rysunku 1, natomiast szczegółowe informacje o wykorzystanym oprogramowaniu zawiera Tabela materiałów.
<img alt="Schemat generowania struktury molekularnej, przedstawiający kroki: filtrowanie substruktur, zastępowanie pseudoatomów." class="xfigimg" src="/files/ftp_upload/63632/63632fig01.jpg" data-altupdatedat="1788282455000" data-alt="Rysunek 1"> 
Rysunek 1: Schemat podsumowujący protokół. <a href="https://www.jove.com/files/ftp_upload/63632/63632fig01large.jpg" target="_blank">Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.</a>
1. Pobieranie oprogramowania i plików
UWAGA: Wszystkie programy są bezpłatne do użytku indywidualnego i mogą być uruchamiane na komputerze osobistym.
<ol>
	<li>Utwórz nowy katalog dla tego projektu. Umieść w nim pliki i pliki wykonywalne, aby zapewnić do nich łatwy dostęp.</li>
	<li>Pobierz i zainstaluj niezbędne pakiety oprogramowania.
	<ol>
		<li>Pobierz najnowszą wersję MAYGEN jako plik .jar. 
		UWAGA: MAYGEN jest dostępny bezpłatnie jako plik .jar pod adresem https://github.com/MehmetAzizYirik/MAYGEN/releases</li>
		<li>Pobierz i zainstaluj oprogramowanie do zarządzania pakietami Conda oraz zestaw narzędzi chemoinformatycznych RDKit18. 
		UWAGA: RDKit będzie filtrować struktury molekularne generowane przez MAYGEN i najlepiej działa w środowisku Conda. Instrukcje dotyczące pobierania platformy Conda znajdują się na stronie https://conda.io/projects/conda/en/latest/user-guide/install/index.html. Instrukcje dotyczące instalacji RDKit i konfiguracji środowiska można znaleźć pod adresem https://www.rdkit.org/docs/Install.html.
		<ol>
			<li>Zainstaluj RDKit w głównym środowisku Conda zamiast w osobnym środowisku RDKit za pomocą wiersza poleceń Anaconda prompt. W systemach Windows wyszukaj „Anaconda prompt” i kliknij wynikający z tego skrót, aby go uruchomić. W systemach MacOS i Linux komunikuj się z Condą poprzez terminal bez uruchamiania dodatkowych programów. Następnie wpisz poniższą komendę i naciśnij Enter, aby ją wykonać, a na wszelkie pytania pojawiające się podczas instalacji odpowiedz twierdząco: 
			conda install -c rdkit rdkit. 
			Choć istnieje wiele bezpłatnych programów do obliczania deskryptorów, w tym przykładzie zastosowano PaDEL-Descriptor19, bezpłatny i szybki kalkulator deskryptorów molekularnych i odcisków palców (fingerprints).</li>
		</ol>
		</li>
		<li>Pobierz i zapisz plik .jar w folderze projektu. 
		UWAGA: PaDEL-Descriptor można pobrać bezpłatnie ze strony http://www.yapcwsoft.com/dd/padeldescriptor/.</li>
	</ol>
	</li>
	<li>Pobierz notatniki Jupyter oraz pliki tekstowe z wzorcami podstruktur z Plików uzupełniających 1-5. 
	UWAGA: Notatniki Jupyter można również pobrać z następującej strony GitHub: https://github.com/cmayerb1/AA-structure-manip.</li>
</ol>
2. Generowanie struktury przy użyciu MAYGEN
<ol>
	<li>W wierszu poleceń przejdź do katalogu zawierającego plik wykonywalny MAYGEN .jar.</li>
	<li>Dla każdego interesującego wzoru chemicznego uruchom MAYGEN za pomocą następującego polecenia: 
	java -jar [nazwa pliku MAYGEN .jar] -f [wzór chemiczny] -v -o [folder na wyniki MAYGEN] -m -sdf. 
	UWAGA: Spowoduje to zapisanie pliku .sdf w wyznaczonym folderze, nazwanego zgodnie z użytym wzorem.
	<ol>
		<li>Jeśli zamiast wzoru dyskretnego użyto wzoru rozmytego, zastąp flagę -f flagą -fuzzy i zamknij przedziały pierwiastków w nawiasach kwadratowych (np. użyj C[5-7]H[12-15], aby zapewnić, że wszystkie wygenerowane struktury posiadają od 5 do 7 atomów węgla i od 12 do 15 atomów wodoru).</li>
	</ol>
	</li>
</ol>
3. Filtrowanie związków z niepożądanymi podstrukturami
<ol>
	<li>Otwórz wiersz poleceń Anaconda (patrz krok 1.2.2.1) i przejdź do folderu zawierającego notatniki Jupyter pobrane z Supplemental File 1.</li>
	<li>Otwórz notatnik Jupyter do filtrowania podstruktur za pomocą następującego polecenia: 
	jupyter notebook [nazwa pliku notatnika]</li>
	<li>W wyznaczonej komórce na początku notatnika wprowadź jako ciągi znaków (w cudzysłowie) pełną ścieżkę do pliku wejściowego .sdf (wygenerowanego przez MAYGEN), pełną ścieżkę do pożądanego pliku wyjściowego .sdf oraz ścieżkę do pliku „badlist”. Przykład pliku badlist znajduje się w Supplemental File 2.
	<ol>
		<li>Jeśli niektóre podstruktury w filtrowanej bibliotece (tzw. goodlist) mają zostać zachowane, utwórz plik .txt z wzorami SMARTS20 dla tych podstruktur (goodlist) i wpisz ścieżkę do pliku goodlist w wyznaczonej linii na początku notatnika. Przykład pliku goodlist znajduje się w Supplemental File 3.</li>
	</ol>
	</li>
	<li>Zrestartuj jądro (kernel) notatnika i uruchom wszystkie komórki (z menu na górze wybierz Kernel, Restart &amp; Run All), aby uzyskać plik .sdf o żądanej nazwie w określonym folderze wyjściowym.</li>
	<li>Powtórz dwa poprzednie kroki dla każdego pliku struktur wygenerowanego przez MAYGEN w kroku 2.</li>
</ol>
4. (Opcjonalnie) Dodatkowe modyfikacje struktury
UWAGA: Czynności te zostały wykonane w tym przykładzie, lecz mogą nie być wymagane przy kurateli innych bibliotek.
<ol>
	<li>Zastępowanie pseudoatomami. 
	UWAGA: W tym kontekście pseudoatom to unikalny atom używany do reprezentacji większej podstruktury wspólnej dla wszystkich generowanych struktur, co skraca czas generowania w programie MAYGEN. Przykład zastępowania pseudoatomami znajduje się w Pliku uzupełniającym 4.
	<ol>
		<li>Otwórz wiersz poleceń Anaconda (patrz krok 1.2.2.1) i przejdź do folderu zawierającego notatniki Jupyter.</li>
		<li>Otwórz notatnik Jupyter do zastępowania pseudoatomami: 
		jupyter notebook [nazwa pliku notatnika]</li>
		<li>W wyznaczonej komórce na początku notatnika wprowadź pełną ścieżkę do pliku wejściowego .sdf oraz pełną ścieżkę do pożądanego pliku wyjściowego .sdf w formie ciągów znaków (w cudzysłowie).</li>
		<li>Zrestartuj jądro (kernel) notatnika i uruchom wszystkie komórki, aby uzyskać plik .sdf o żądanej nazwie w określonym folderze wyjściowym.</li>
	</ol>
	</li>
	<li>Blokowanie końców N- i C- aminokwasów 
	UWAGA: Ta procedura jest specyficzna dla alfa-aminokwasów i polega na dodawaniu bloków molekularnych do końców N- i C- szkieletów alfa-aminokwasów. Przykład blokowania aminokwasów znajduje się w Pliku uzupełniającym 5.
	<ol>
		<li>Otwórz wiersz poleceń Anaconda (patrz krok 1.2.2.1) i przejdź do folderu zawierającego notatniki Jupyter.</li>
		<li>Otwórz notatnik Jupyter do blokowania aminokwasów: 
		jupyter notebook [nazwa pliku notatnika]</li>
		<li>W wyznaczonej komórce na początku notatnika wprowadź pełną ścieżkę do pliku wejściowego .sdf oraz pełną ścieżkę do pożądanego pliku wyjściowego .sdf w formie ciągów znaków (w cudzysłowie).</li>
		<li>Zrestartuj jądro (kernel) notatnika i uruchom wszystkie komórki, aby uzyskać plik .sdf o żądanej nazwie w określonym folderze wyjściowym.</li>
	</ol>
	</li>
</ol>
5. Generowanie deskryptorów
<ol>
	<li>Przed generowaniem deskryptorów umieść wszystkie pliki .sdf, dla których mają zostać obliczone deskryptory, w jednym folderze. 
	UWAGA: Jeśli nie zostało to jeszcze zrobione, nadaj tym plikom opisowe nazwy, aby ułatwić filtrowanie po wygenerowaniu deskryptorów.</li>
	<li>Otwórz wiersz poleceń i przejdź do folderu zawierającego plik .jar programu PaDEL-Descriptor.</li>
	<li>Uruchom PaDEL-Descriptor dla zebranych plików .sdf, używając następującej komendy: 
	java -jar PaDEL-Descriptor.jar -dir [ścieżka do katalogu z plikami .sdf] -file [ścieżka do pliku .csv dla wyników] -2d -retainorder -usefilenameasmolname 
	UWAGA: Plik z wynikami będzie zawierał nazwę cząsteczki w pierwszej kolumnie oraz każdy deskryptor w kolejnych kolumnach.</li>
	<li>Wyeksportuj te dane do dowolnego programu arkusza kalkulacyjnego w celu dalszej analizy.</li>
</ol>

<ol><li>Ruddigkeit, L., van Deursen, R., Blum, L. C., Reymond, J. -L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Enumeration+of+166+billion+organic+small+molecules+in+the+chemical+universe+database+GDB-17%5BTitle%5D)+AND+%22Journal+of+Chemical+Information+and+Modeling%22%5BJournal%5D)">Enumeration of 166 billion organic small molecules in the chemical universe database GDB-17.</a> Journal of Chemical Information and Modeling. 52 (11), 2864-2875 (2012).</li>
<li>Buchanan, B. G., Feigenbaum, E. A. Dendral and Meta-Dendral: their applications dimension. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aBuchanan%2C+BG+%22Readings+in+Artificial+Intelligence%22">Readings in Artificial Intelligence</a>. Webber, B. L., Nilsson, N. J. , Morgan Kaufmann. 313-322 (1981).</li>
<li>Gugisch, R., et al. MOLGEN 5.0, A Molecular Structure Generator. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aGugisch%2C+R+%22Advances+in+Mathematical+Chemistry+and+Applications%22">Advances in Mathematical Chemistry and Applications</a>. Basak, S. C., Restrepo, G., Villaveces, J. L. , Bentham Science Publishers. 113-138 (2015).</li>
<li>Yirik, M. A., Steinbeck, C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Chemical+graph+generators%5BTitle%5D)+AND+%22PLOS+Computational+Biology%22%5BJournal%5D)">Chemical graph generators.</a> PLOS Computational Biology. 17 (1), 1008504(2021).</li>
<li>Jaghoori, M. M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((PMG%3A+multi-core+metabolite+identification%5BTitle%5D)+AND+%22Electronic+Notes+in+Theoretical+Computer+Science%22%5BJournal%5D)">PMG: multi-core metabolite identification.</a> Electronic Notes in Theoretical Computer Science. 299, 53-60 (2013).</li>
<li>Yirik, M. A., Sorokina, M., Steinbeck, C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((MAYGEN%3A+an+open-source+chemical+structure+generator+for+constitutional+isomers+based+on+the+orderly+generation+principle%5BTitle%5D)+AND+%22Journal+of+Cheminformatics%22%5BJournal%5D)">MAYGEN: an open-source chemical structure generator for constitutional isomers based on the orderly generation principle.</a> Journal of Cheminformatics. 13 (1), 48(2021).</li>
<li>Sims, C. C. Computational methods in the study of permutation groups. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aSims%2C+CC+%22Computational+Problems+in+Abstract+Algebra%22">Computational Problems in Abstract Algebra</a>. Leech, J. , Pergamon. 169-183 (1970).</li>
<li>Mat, W. -K., Xue, H., Wong, J. T. -F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+genomics+of+LUCA%5BTitle%5D)+AND+%22Frontiers+in+Bioscience%22%5BJournal%5D)">The genomics of LUCA.</a> Frontiers in Bioscience. 13, 5605-5613 (2008).</li>
<li>Fournier, G. P., Alm, E. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Ancestral+reconstruction+of+a+pre-LUCA+aminoacyl-tRNA+synthetase+ancestor+supports+the+late+addition+of+Trp+to+the+genetic+code%5BTitle%5D)+AND+%22Journal+of+Molecular+Evolution%22%5BJournal%5D)">Ancestral reconstruction of a pre-LUCA aminoacyl-tRNA synthetase ancestor supports the late addition of Trp to the genetic code.</a> Journal of Molecular Evolution. 80 (3-4), 171-185 (2015).</li>
<li>Higgs, P. G., Pudritz, R. E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+Thermodynamic+basis+for+prebiotic+amino+acid+synthesis+and+the+nature+of+the+first+genetic+code%5BTitle%5D)+AND+%22Astrobiology%22%5BJournal%5D)">A Thermodynamic basis for prebiotic amino acid synthesis and the nature of the first genetic code.</a> Astrobiology. 9 (5), 483-490 (2009).</li>
<li>Trifonov, E. N. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Consensus+temporal+order+of+amino+acids+and+evolution+of+the+triplet+code%5BTitle%5D)+AND+%22Gene%22%5BJournal%5D)">Consensus temporal order of amino acids and evolution of the triplet code.</a> Gene. 261 (1), 139-151 (2000).</li>
<li>Cleaves, H. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+origin+of+the+biologically+coded+amino+acids%5BTitle%5D)+AND+%22Journal+of+Theoretical+Biology%22%5BJournal%5D)">The origin of the biologically coded amino acids.</a> Journal of Theoretical Biology. 263 (4), 490-498 (2010).</li>
<li>Moosmann, B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Redox+biochemistry+of+the+genetic+code%5BTitle%5D)+AND+%22Trends+in+Biochemical+Sciences%22%5BJournal%5D)">Redox biochemistry of the genetic code.</a> Trends in Biochemical Sciences. 46 (2), 83-86 (2021).</li>
<li>Simkus, D. N., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Methodologies+for+analyzing+soluble+organic+compounds+in+extraterrestrial+samples%3A+amino+acids%2C+amines%2C+monocarboxylic+acids%2C+aldehydes%2C+and+ketones%5BTitle%5D)+AND+%22Life%22%5BJournal%5D)">Methodologies for analyzing soluble organic compounds in extraterrestrial samples: amino acids, amines, monocarboxylic acids, aldehydes, and ketones.</a> Life. 9 (2), 47(2019).</li>
<li>Criado-Reyes, J., Bizzarri, B. M., García-Ruiz, J. M., Saladino, R., Di Mauro, E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+role+of+borosilicate+glass+in+Miller-Urey+experiment%5BTitle%5D)+AND+%22Scientific+Reports%22%5BJournal%5D)">The role of borosilicate glass in Miller-Urey experiment.</a> Scientific Reports. 11 (1), 21009(2021).</li>
<li>Parker, E. T., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Primordial+synthesis+of+amines+and+amino+acids+in+a+1958+Miller+H2S-rich+spark+discharge+experiment%5BTitle%5D)+AND+%22Proceedings+of+the+National+Academy+of+Sciences+of+the+United+States+of+America%22%5BJournal%5D)">Primordial synthesis of amines and amino acids in a 1958 Miller H2S-rich spark discharge experiment.</a> Proceedings of the National Academy of Sciences of the United States of America. 108 (14), 5526-5531 (2011).</li>
<li>Bada, J. L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((New+insights+into+prebiotic+chemistry+from+Stanley+Miller%27s+spark+discharge+experiments%5BTitle%5D)+AND+%22Chemical+Society+Reviews%22%5BJournal%5D)">New insights into prebiotic chemistry from Stanley Miller's spark discharge experiments.</a> Chemical Society Reviews. 42 (5), 2186-2196 (2013).</li>
<li>RDKit: Open-source cheminformatics. , Available from: 
 <a target="_blank" href="http://www.rdkit.org">http://www.rdkit.org</a> (2021).</li>
<li>Yap, C. W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((PaDEL-descriptor%3A+An+open+source+software+to+calculate+molecular+descriptors+and+fingerprints%5BTitle%5D)+AND+%22Journal+of+Computational+Chemistry%22%5BJournal%5D)">PaDEL-descriptor: An open source software to calculate molecular descriptors and fingerprints.</a> Journal of Computational Chemistry. 32 (7), 1466-1474 (2011).</li>
<li>SMARTS - A language for describing molecular patterns. Daylight Chemical Information Systems, Inc. , Available from: 
 <a target="_blank" href="https://www.daylight.com/html/doc/theory/theory.smarts.html">https://www.daylight.com/html/doc/theory/theory.smarts.html</a> (2019).</li>
<li>Meringer, M., Cleaves, H. J., Freeland, S. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Beyond+terrestrial+biology%3A+charting+the+chemical+universe+of+%CE%B1-amino+acid+structures%5BTitle%5D)+AND+%22Journal+of+Chemical+Information+and+Modeling%22%5BJournal%5D)">Beyond terrestrial biology: charting the chemical universe of α-amino acid structures.</a> Journal of Chemical Information and Modeling. 53 (11), 2851-2862 (2013).</li>
<li>Zherebker, A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Speciation+of+organosulfur+compounds+in+carbonaceous+chondrites%5BTitle%5D)+AND+%22Scientific+Reports%22%5BJournal%5D)">Speciation of organosulfur compounds in carbonaceous chondrites.</a> Scientific Reports. 11 (1), 7410(2021).</li>
<li>Tanford, C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+hydrophobic+effect+and+the+organization+of+living+matter%5BTitle%5D)+AND+%22Science%22%5BJournal%5D)">The hydrophobic effect and the organization of living matter.</a> Science. 200 (4345), 1012-1018 (1978).</li>
<li>Grantham, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Amino+acid+difference+formula+to+help+explain+protein+evolution%5BTitle%5D)+AND+%22Science%22%5BJournal%5D)">Amino acid difference formula to help explain protein evolution.</a> Science. 185 (4154), 862-864 (1974).</li>
<li>Cordero, B., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Covalent+radii+revisited%5BTitle%5D)+AND+%22Dalton+Transactions%22%5BJournal%5D)">Covalent radii revisited.</a> Dalton Transactions. (21), 2832-2838 (2008).</li>
<li>Cleaves, H. J., Butch, C., Burger, P. B., Goodwin, J., Meringer, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((One+among+millions%3A+the+chemical+space+of+nucleic+acid-like+molecules%5BTitle%5D)+AND+%22Journal+of+Chemical+Information+and+Modeling%22%5BJournal%5D)">One among millions: the chemical space of nucleic acid-like molecules.</a> Journal of Chemical Information and Modeling. 59 (10), 4266-4277 (2019).</li>
</ol>

Autorzy nie mają do ujawnienia żadnych konfliktów interesów.

Jedną z cech "wczesnych" aminokwasów jest brak siarki. Wspomniane wcześniej metaanalizy generalnie uważają, że kodowane aminokwasy zawierające siarkę (Cys i Met) były stosunkowo późnymi dodatkami do kodu genetycznego, a wnioski te potwierdza brak aminokwasów zawierających siarkę w meteorytach i eksperymentach z lampami zapłonowymi. Jednak związki siarkoorganiczne są łatwo wykrywalne w kometach i meteorytach22, a ponowna analiza eksperymentów z lampą iskrową z użyciem gazuH2Swykazała aminokwasy i inne związki organiczne zawierające siarkę16. Rozważając alternatywny alfabet aminokwasów, warto zapoznać się z tym wzbogaconym w siarkę.
W powyższym protokole generowanie struktury i filtrowanie podstruktury są uważane za kroki krytyczne; W zależności od składu gotowej biblioteki struktur, badacz może być zmuszony do wykonania tylko tych dwóch kroków. Instrukcje i oprogramowanie do dodatkowych działań (zastępowanie pseudoatomów i dodawanie podstruktur (w tym przypadku capping aminokwasów)) są dołączone w celu uzyskania bardziej odpowiednich obliczeń deskryptorowych (capping zapewnia, że na obliczenia XLogP ma wpływ łańcuch boczny, a nie grupy aminowe lub karboksylowe w szkielecie) oraz szybszego generowania struktury za pomocą pseudoatomu, co omówiono bardziej szczegółowo poniżej. Dodatkowo wykonywane są tutaj obliczenia deskryptorowe, które w łatwy sposób pozwalają na wizualizację różnorodności generowanych struktur i porównanie efektów wzbogacania siarką w gotowych bibliotekach.
Podczas gdy PaDEL-Descriptor może obliczyć tysiące właściwości molekularnych, objętość cząsteczkowa (obliczona objętość van der Waalsa) i współczynnik podziału (jako XLogP) zostały tutaj użyte z dwóch różnych powodów. Po pierwsze, te dwa deskryptory mierzą właściwości molekularne (odpowiednio rozmiar i hydrofobowość), które są znane większości chemików i biologów. Po drugie, w przypadku aminokwasów te dwie właściwości są znaczące. Przez dziesięciolecia wiadomo było, że wielkość aminokwasów i hydrofobowość wpływają na termodynamikę fałdowania białek23. Te dwie właściwości pomagają wyjaśnić częstotliwości podstawiania aminokwasów, które były integralną częścią zrozumienia ewolucji białek24.
Powyższy przykład pokazuje, że w dwóch badanych deskryptorach (objętość cząsteczkowa i hydrofobowość) podstawienie dwuwartościowej siarki na węgiel i dwa wodory nie powoduje istotnych zmian. Niewielki, nieistotny wzrost średniej objętości cząsteczkowej spowodowany podstawieniem siarki (ryc. 3) można przypisać większemu promieniowi kowalencyjnemu siarki (~103 pm) w porównaniu z węglem sp3 (~75 pm) lub sp2 (~73 pm)25. Podobnie, podstawienie siarki ma minimalny wpływ na średnią XLogP (ryc. 4). Największy efekt zaobserwowano między bibliotekami VAIL i VAIL_S, prawdopodobnie ze względu na połączenie biblioteki VAIL szczególnie hydrofobowej (łańcuchy boczne to tylko węglowodory) i grup sulfhydrylowych, które byłyby znacznie bardziej kwaśne niż grupy metylowe, które miałyby zastąpić. Minimalny efekt podstawiania siarki jest widoczny na rysunku 2, gdzie biblioteki z podstawieniem siarki zajmują tę samą przestrzeń chemiczną, co analogiczne biblioteki bez podstawiania siarki.
Zmniejszenie liczby struktur (Rysunek 5A) i czasu potrzebnego do wygenerowania tych struktur (Rysunek 5B) przy użyciu pseudoatomu nie jest zaskakujące. Użycie pseudoatomu zmniejsza liczbę ciężkich atomów, które muszą zostać włączone do wykresu chemicznego, zmniejszając liczbę węzłów grafu i powodując wykładnicze skrócenie czasu generowania i liczby struktur. W tym przypadku wybór fosforu trójwartościowego jako pseudoatomu wynika z podstawowej biochemii (brak potranslacyjnego dodawania grup fosforanowych, żadne genetycznie kodowane aminokwasy nie zawierają fosforu) oraz wartościowości atomu, który by go zastąpił (fosfor trójwartościowy można łatwo zastąpić czterowartościowym węglem, który jest pojedynczo związany z innym atomem lub grupą atomów). Podczas gdy dostarczony kod do podstawiania pseudoatomów jest specyficzny dla zastąpienia trójwartościowego fosforu podstrukturą alaniny, użytkownicy mogą dostosować kod do pracy z różnymi pseudoatomami lub podstrukturami zastępczymi, potencjalnie używając wielu pseudoatomów podczas początkowego generowania struktury, a następnie zastępując każdy pseudoatom większą podstrukturą molekularną.
Metody generowania struktur podobne do tych stosowanych w ramach projektu MAYGEN (i inne metody, takie jak sieci neuronowe) są już wykorzystywane w odkrywaniu leków do generowania bibliotek związków do badań przesiewowych in silico ; W niedawnym przeglądzie4 omówiono te metody bardziej szczegółowo. Ponieważ metody te są przeznaczone przede wszystkim do tworzenia cząsteczek podobnych do leków, istnieją pewne ograniczenia ich zdolności do generowania cząsteczek, takie jak wykorzystanie właściwości biologicznych lub farmaceutycznych w celu ograniczenia tworzonych struktur (odwrotne QSPR/QSAR) lub tworzenie struktur z ustalonej liczby bloków budulcowych podstruktury. Ponieważ astrobiologia skupia się bardziej na mnogości związków organicznych, które mogą tworzyć się abiotycznie, a mniej na jakichkolwiek produktach końcowych lub ich właściwościach, wyczerpujące generowanie struktur w ramach projektu MAYGEN jest idealne do tworzenia bibliotek struktur w celu odpowiedzi na pytania astrobiologiczne. Opisane tutaj podejście do filtrowania podkonstrukcji (wykonywane po wygenerowaniu struktury za pomocą zewnętrznego programu) różni się od konkurencyjnego programu MOLGEN tym, że filtrowanie podstruktury MOLGEN następuje podczas generowania konstrukcji. Ponieważ MAYGEN jest oprogramowaniem typu open source, jest nie tylko bardziej dostępny niż MOLGEN ze względu na koszt licencji MOLGEN, ale osoby fizyczne mogą wdrażać nowe funkcje, takie jak filtrowanie podstruktur podczas generowania struktury.
Jak napisano, opisany tutaj protokół koncentruje się na generowaniu i utrzymywaniu bibliotek stosunkowo małych alfa-aminokwasów. Aby wygenerować różne biblioteki, użytkownicy mogą nadać MAYGENOWI różne wzory molekularne, zmienić filtrowanie podstruktury, zmieniając maksymalny dozwolony rozmiar pierścienia i wartościowość wiązania, lub edytować pliki goodlist i badlist, aby dodać lub usunąć wzorce podstruktury. Modyfikacje protokołu, które obejmują zmianę sposobu dodawania lub zastępowania atomów i podstruktur (podstawienie pseudoatomów i ograniczenie molekularne) są wykonalne, ale będą wymagały zwrócenia większej uwagi na ograniczenia walencyjne, aby uniknąć błędów RDKit dotyczących nieprawidłowych wartościowości w zmodyfikowanych strukturach.
Opisany powyżej protokół jest przeznaczony dla małych alfa-aminokwasów. Jednak ogólny format (kompleksowe generowanie struktury za pomocą pseudoatomów, a następnie filtrowanie podstruktury i modyfikacje molekularne) jest bardzo elastyczny dla związków innych niż małe aminokwasy. Nawet w astrobiologii podobną niedawną procedurę z użyciem MOLGEN zastosowano do zbadania konstytucyjnych izomerów kwasów nukleinowych26. Oprócz narzędzi opisanych powyżej, MAYGEN można połączyć z innymi narzędziami chemicznymi typu open source, aby tworzenie i analizowanie nowych struktur chemicznych było przystępne cenowo i dostępne dla szerokiego wachlarza dziedzin badawczych.

Generowanie struktur molekularnych służy jako praktyczne zastosowanie ogólnego problemu wyczerpującego generowania grafów; biorąc pod uwagę kilka węzłów (atomów) i ograniczenia ich łączności (np. wartościowości, wielokrotności wiązań, pożądane/niepożądane podstruktury), ile połączonych grafów (molekuł) jest możliwych? Generatory struktur znalazły szerokie zastosowanie w odkrywaniu leków i rozwoju farmaceutycznym, gdzie mogą tworzyć ogromne biblioteki nowych struktur do badań przesiewowych in silico1.
Pierwszy generator struktur, CONGEN, został opracowany dla pierwszego projektu sztucznej inteligencji w chemii organicznej, DENDRAL2 (skrót od DENDRitic ALgorithm). W literaturze opisano kilka następców oprogramowania DENDRAL; Jednak nie wszystkie z nich były utrzymywane lub wydajne. Obecnie MOLGEN3 jest najnowocześniejszym generatorem struktur molekularnych. Niestety dla większości potencjalnych użytkowników jest to oprogramowanie o zamkniętym kodzie źródłowym i wymaga opłaty licencyjnej. W związku z tym zaistniała potrzeba wydajnego generatora struktur typu open source, który można łatwo dostosować do określonych zastosowań. Jednym z wyzwań stojących przed wydajnym generatorem struktur jest zarządzanie eksplozją kombinatoryczną; Wraz ze wzrostem rozmiaru wzoru cząsteczkowego wykładniczo zwiększa się rozmiar chemicznej przestrzeni wyszukiwania. Niedawny przegląd dokładniej bada historię i wyzwania związane z generowaniem struktur molekularnych4.
Przed 2021 rokiem, Parallel Molecule Generator (PMG)5 był najszybszym generatorem struktur open-source, ale nadal był wolniejszy od MOLGEN o rzędy wielkości. MAYGEN6 jest około 47 razy szybszy niż PMG i około 3 razy wolniejszy niż MOLGEN, co czyni MAYGEN najszybszym i najbardziej wydajnym generatorem struktur open-source dostępnym na rynku. Bardziej szczegółowe porównania i testy porównawcze można znaleźć w artykule wprowadzającym MAYGEN6. Kluczową cechą programu jest test oparty na porządkowaniu leksykograficznym dla struktur kanonicznych, uporządkowana metoda generowania grafów oparta na algorytmie Schreier-Sims7. Oprogramowanie można łatwo zintegrować z innymi projektami i ulepszyć zgodnie z potrzebami użytkowników.
Podobnie jak MOLGEN i PMG, MAYGEN bierze zdefiniowany przez użytkownika wzór molekularny i generuje wszystkie możliwe struktury dla tego wzoru. Na przykład, jeśli użytkownik uruchomi MAYGEN ze wzorem C5H12, MAYGEN wygeneruje wszystkie możliwe struktury zawierające pięć atomów węgla i dwanaście atomów wodoru. W przeciwieństwie do swojego odpowiednika PMG o otwartym kodzie źródłowym, MAYGEN może również obsługiwać "rozmyte" wzory molekularne, które używają interwałów zamiast dyskretnych liczb do liczenia każdego pierwiastka. Na przykład, jeśli użytkownik uruchomi MAYGEN ze wzorem C5-7H12-15, MAYGEN wygeneruje wszystkie możliwe struktury, które zawierają od pięciu do siedmiu atomów węgla i dwanaście do piętnastu atomów wodoru, co pozwala na proste generowanie struktur o szerokim zakresie składu atomowego.
Astrobiologia jest jedną z takich dziedzin, która może skorzystać z generatorów struktur molekularnych. Popularnym tematem w astrobiologii jest ewolucja alfabetu aminokwasów wspólnego dla całego istniejącego życia na Ziemi. Jedną z cech charakterystycznych Ostatniego Uniwersalnego Wspólnego Przodka (LUCA) jest użycie dwudziestu genetycznie kodowanych aminokwasów do budowy białek8,9. Na podstawie metaanaliz pracy w wielu dziedzinach10,11,12, około 10 z tych aminokwasów (Gly, Ala, Val, Asp, Glu, Ser, Thr, Leu, Ile, Pro) łatwo tworzy się w warunkach abiotycznych i prawdopodobnie tworzy alfabet aminokwasów organizmów pre-Luca. Z biegiem czasu ten "wczesny" alfabet został rozszerzony w odpowiedzi na różne potrzeby strukturalne i funkcjonalne. Na przykład, niedawna recenzja z Moosmann13 twierdzi, że dodanie nowszych członków genetycznie kodowanych aminokwasów (mianowicie Met, Tyr i Trp) pozwoliło na przetrwanie w środowisku bogatym w tlen, zapobiegając wewnątrzkomórkowej proliferacji reaktywnych form tlenu.
Stale rosnący zestaw technik chemii analitycznej pozwala na wgląd w struktury aminokwasów, które mogą powstawać w warunkach abiotycznych. Niedawny przegląd14 autorstwa Simkusa i innych szczegółowo opisuje metody używane do wykrywania licznych związków organicznych w meteorytach, a także związków organicznych z symulacji in vitro wczesnych środowisk Ziemi15,16,17. Systematyczne generowanie struktur chemicznych pozwala naukowcom na eksplorację poza związki organiczne wykryte za pomocą oprzyrządowania, wypełniając przestrzeń strukturalną wokół strukturalnych "wysp" zidentyfikowanych przez chemię analityczną. W przypadku "wczesnych" aminokwasów, to systematyczne generowanie struktur pokazuje możliwe chemie białek dostępne dla wczesnego życia, bez ograniczania eksploracji do struktur, które zostały eksperymentalnie wykryte w warunkach syntezy abiotycznej. Dzięki zestawom narzędzi chemicznych i wydajnym generatorom struktur, takim jak MAYGEN, tworzenie i eksplorowanie nowatorskich bibliotek struktur chemicznych jest teraz łatwiejsze niż kiedykolwiek wcześniej i może prowadzić bardziej szczegółowe badania nad alternatywnymi substancjami chemicznymi życia.

<table><tbody><tr><td>conda v. 4.10.3</td><td></td><td></td><td>https://www.anaconda.com/products/individual</td></tr><tr><td>Java 17</td><td></td><td></td><td>https://java.com/en/download/help/download_options.html</td></tr><tr><td>MAYGEN v. 1.8</td><td></td><td></td><td>https://github.com/MehmetAzizYirik/MAYGEN/releases</td></tr><tr><td>PaDEL-Descriptor v. 2.21</td><td></td><td></td><td>http://www.yapcwsoft.com/dd/padeldescriptor/</td></tr><tr><td>python v. 3.7.11</td><td></td><td></td><td>zawarty w środowisku Anaconda</td></tr><tr><td>RDKit v. 2020.09.1.0</td><td></td><td></td><td>https://www.rdkit.org/docs/Install.html, lub zainstalowany przez conda: https://anaconda.org/rdkit/rdkit</td></tr><tr><td>*Te konkretne wersje zostały użyte w tym manuskrypcie; użytkownik może uzyskać nowsze wersje, jeśli są dostępne.</td><td></td><td></td><td></td></tr></tbody></table>

kuratorstwo obliczeniowych bibliotek chemicznych zademonstrowanych za pomocą alfa-aminokwasów

Curation of Computational Chemical Libraries Demonstrated with Alpha-Amino Acids

Wyczerpujące generowanie struktur molekularnych ma liczne zastosowania chemiczne i biochemiczne, takie jak projektowanie leków, budowa molekularnych baz danych, eksploracja alternatywnych biochemii i wiele innych. Matematycznie rzecz biorąc, są to generatory grafów z ograniczeniami chemicznymi. W tej dziedzinie najbardziej wydajnym obecnie generatorem (MOLGEN) jest produkt komercyjny, co ogranicza jego zastosowanie. Alternatywnie do tego, inny generator struktury molekularnej, MAYGEN, jest najnowszym narzędziem typu open source o wydajności porównywalnej z MOLGEN i zdolnością użytkowników do zwiększania jego wydajności poprzez dodawanie nowych funkcji. Jedną z dziedzin badań, które mogą skorzystać na tym rozwoju, jest astrobiologia; Generatory struktur umożliwiają naukowcom uzupełnianie danych eksperymentalnych o możliwości obliczeniowe dla alternatywnej biochemii. Protokół ten szczegółowo opisuje jeden przypadek użycia do generowania struktur w astrobiologii, a mianowicie generowania i selekcjonowania bibliotek alfa-aminokwasów. Korzystając z generatorów struktur typu open source i narzędzi cheminoinformatycznych, opisane tutaj praktyki mogą być wdrażane poza astrobiologię w celu taniego tworzenia i zarządzania bibliotekami struktur chemicznych dla dowolnego pytania badawczego.

<table border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always">
	<tbody >
		<tr >
			<td rowspan="2"></td>
			<td rowspan="2">Biblioteka</td>
			<td colspan="2" rowspan="2">Wzór</td>
			<td colspan="2" rowspan="2">Dodatkowe ograniczenia</td>
			<td rowspan="2">Aminokwasy kodowane jako „wczesne”</td>
			<td rowspan="2">Czas generowania (ms)</td>
			<td colspan="2">Struktury</td>
		</tr>
		<tr >
			<td>Początkowe</td>
			<td>Końcowe</td>
		</tr>
		<tr >
			<td>1</td>
			<td>Gly</td>
			<td colspan="2">C2H5NO2</td>
			<td colspan="2">zawierać podstrukturę Gly</td>
			<td>Gly</td>
			<td>192</td>
			<td>84</td>
			<td>1</td>
		</tr>
		<tr >
			<td>2</td>
			<td>VAIL</td>
			<td colspan="2">PC0-3H3-9</td>
			<td colspan="2"></td>
			<td>Val, Ala, Ile, Leu</td>
			<td>172</td>
			<td>70</td>
			<td>22</td>
		</tr>
		<tr >
			<td>3</td>
			<td>DEST</td>
			<td colspan="2">PC0-3O1-2H3-5</td>
			<td colspan="2"></td>
			<td>Asp, Glu, Ser, Thr</td>
			<td>481</td>
			<td>1928</td>
			<td>254</td>
		</tr>
		<tr >
			<td>4</td>
			<td>Pro</td>
			<td colspan="2">C2-5NO2H7-11</td>
			<td colspan="2">Zawierać podstrukturę N-meGly lub N-meAla</td>
			<td>Pro</td>
			<td>4035</td>
			<td>79777</td>
			<td>16</td>
		</tr>
		<tr >
			<td>5</td>
			<td>VAIL_S</td>
			<td colspan="2">PSC0-2H3-7</td>
			<td colspan="2"></td>
			<td></td>
			<td>122</td>
			<td>65</td>
			<td>31</td>
		</tr>
		<tr >
			<td>6</td>
			<td>DEST_S</td>
			<td colspan="2">PSC0-2O1-2H3</td>
			<td colspan="2"></td>
			<td></td>
			<td>349</td>
			<td>1075</td>
			<td>79</td>
		</tr>
		<tr >
			<td>7</td>
			<td>Pro_S</td>
			<td colspan="2">C2-4SNO2H7-9</td>
			<td colspan="2">Zawierać podstrukturę N-meGly lub N-meAla</td>
			<td></td>
			<td>3999</td>
			<td>75734</td>
			<td>10</td>
		</tr>
	</tbody>
</table>
Tabela 1: Biblioteki związków użyte w tym przykładzie. Biblioteki zbudowane na podstawie wzorów 1-4 (Gly, VAIL, DEST i Pro) opierają się na wcześniej opublikowanych wzorach rozmytych (fuzzy formulae) dla „wczesnych” kodowanych aminokwasów21, natomiast biblioteki zbudowane na podstawie wzorów 5-7 (VAIL_S, DEST_S i Pro_S) opierają się na wariantach wzorów 2-4, w których jeden z atomów węgla zostaje zastąpiony dwuwartościową siarką. Liczba struktur odzwierciedla liczbę cząsteczek wygenerowanych przez MAYGEN dla każdego wzoru („Wstępna”) oraz liczbę cząsteczek pozostałych po odfiltrowaniu tych z niepożądanymi podstrukturami („Końcowa”). Skróty: VAIL = walina, alanina, izoleucyna, leucyna; DEST = kwas asparaginowy, kwas glutaminowy, seryna, treonina; X_S = dwuwartościowa siarka zastępuje jeden z atomów węgla w bibliotece X; N-meX = N-metyloX.
Powyższe metody ogólne zastosowano do formuł opartych na „wczesnych” kodowanych aminokwasach, zgodnie z procedurą Meringera i wsp.21 Struktury z tzw. „czarnej listy” (badlist) pobrano z tego samego źródła i przekształcono w ciągi SMARTS, aby łatwo reprezentować wzorce podstrukturalne. W tym przykładzie nie wykorzystano dwóch podstruktur z czarnej listy: struktura 018 (CH3-CH-N) odpowiadała niemal izomerom proliny, które same w sobie nie były niestabilne; struktura 106 (R-C-C-OH, gdzie R=podstruktura alaniny przyłączona w węglu beta) odpowiadała kwasowi glutaminowemu, który jest kodowanym aminokwasem. Oprócz tych formuł chemicznych stworzono warianty, w których dwuwartościowa siarka zajmuje miejsce jednego atomu węgla i dwóch atomów wodoru. Ze względów wydajnościowych w kilku z tych formuł zastosowano trójwartościowy atom fosforu (np. „pseudoatom”) jako zamiennik węgla beta w podstrukturze alaniny. Tabela 1 przedstawia biblioteki wygenerowane w tym przykładzie, formuły użyte do ich utworzenia oraz liczbę zawartych w nich związków. Nazwy bibliotek opierają się na kodowanych aminokwasach, z których pochodzą: przy użyciu skrótu 3-literowego (Gly = glicyna, Pro = prolina) lub skrótu 1-literowego (VAIL = walina, alanina, izoleucyna, leucyna; DEST = kwas asparaginowy, kwas glutaminowy, seryna, treonina). Sufiks „_S” oznacza, że w formule oryginalnej biblioteki węgiel zastąpiono siarką (np. VAIL_S jest zbudowany na podstawie tej samej formuły przybliżonej co VAIL, ale z dwuwartościową siarką zastępującą jeden z atomów węgla).
Po wygenerowaniu struktur za pomocą MAYGEN, wynikające z tego biblioteki przefiltrowano w celu usunięcia związków zawierających co najmniej jedną podstrukturę znajdującą się na liście wykluczeń (badlist). Po tym filtrowaniu wszystkie atomy fosforu zastąpiono podstrukturą alaniny. Następnie utworzono wersje „zablokowane” (capped) wszystkich struktur, dodając grupę acetylową do N-końca oraz grupę N-metyloamidową do C-końca. Wykonano to, aby wyeliminować wpływ wolnych grup aminowych i karboksylowych w szkielecie alfa-aminokwasów na hydrofobowość. Program PaDEL-Descriptor został wykorzystany do obliczenia XLogP dla wszystkich struktur zablokowanych oraz obliczenia objętości van der Waalsa (VABC) dla wszystkich struktur niezablokowanych.
Rysunek 2 przedstawia przestrzeń chemiczną przefiltrowanych bibliotek, zdefiniowaną przez deskryptory VABC i XLogP. Tutaj zakres możliwych wartości logP rośnie wraz z objętością cząsteczki, nawet w bibliotekach pozbawionych wyraźnych łańcuchów bocznych hydrofilowych (np. VAIL, Pro). Aminokwasy kodowane z hydrokarbonowymi łańcuchami bocznymi były bardziej hydrofobowe niż większość innych aminokwasów o porównywalnej objętości z ich odpowiedniej biblioteki. Wydaje się, że tak samo jest w przypadku Met i Cys w porównaniu z innymi członkami biblioteki VAIL_S o podobnych objętościach. Aminokwasy kodowane z hydroksylowymi łańcuchami bocznymi (Ser i Thr) należały do najmniejszych członków biblioteki DEST, przy czym Asp był tylko nieco większy od Thr.
Rysunek 3 oraz Rysunek 4 przedstawiają wpływ zastąpienia węgla dwuwartościową siarką w łańcuchu bocznym alfa-aminokwasu na objętość i logP. Podstawienie siarką doprowadziło do niewielkiego wzrostu objętości cząsteczkowej we wszystkich bibliotekach (Rysunek 3). Wpływ podstawienia siarką na logP nie jest tak jednolity jak w przypadku objętości (Rysunek 4). Średni logP biblioteki VAIL_S jest nieco niższy niż w bibliotece VAIL, jednak efekt ten nie występuje w żadnej z pozostałych par bibliotek (DEST i DEST_S, Pro i Pro_S).
Rysunek 5 ilościowo określa wpływ na generowanie struktur zastosowania pseudoatomu zastępującego powszechną podstrukturę; w tym przypadku trójkoordynowany P zastąpił resztę alaniny podczas generowania struktury. Wykorzystanie pseudoatomu w procesie generowania struktur znacznie zmniejszyło liczbę wygenerowanych struktur o ok. 3 rzędy wielkości (Rysunek 5A) oraz całkowity czas potrzebny na wygenerowanie tych struktur o 1-2 rzędy wielkości (Rysunek 5B).
<img alt="Wykres rozrzutu właściwości chemicznych; xLogP w funkcji objętości Van der Waalsa; wykres porównawczy aminokwasów." class="xfigimg" src="/files/ftp_upload/63632/63632fig02.jpg" data-altupdatedat="1788282455000" data-alt="Rycina 2"> 
Rycina 2: Przestrzeń chemiczna wszystkich przefiltrowanych bibliotek aminokwasów. Czarne markery reprezentują aminokwasy z bibliotek bez siarki; żółte markery reprezentują aminokwasy z bibliotek wzbogaconych o siarkę. Koła: VAIL i VAIL_S; kwadraty: DEST i DEST_S; trójkąty: Pro i Pro_S; gwiazdki: aminokwasy kodowane. Należy zauważyć, że dwa kodowane aminokwasy zawierające siarkę (Met i Cys) nie są uznawane za aminokwasy „wczesne”, lecz znajdują się w bibliotece VAIL_S. Skróty: XLogP = współczynnik podziału; VAIL = walina, alanina, izoleucyna, leucyna; DEST = kwas asparaginowy, kwas glutaminowy, seryna, treonina; X_S = dwuwartościowa siarka zastępuje jeden z atomów węgla w bibliotece X. <a href="https://www.jove.com/files/ftp_upload/63632/63632fig02large.jpg" target="_blank">Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.</a>
<img alt="Wykres słupkowy porównujący objętości van der Waalsa (Å³) dla VAIL, DEST i Pro z słupkami błędów." class="xfigimg" src="/files/ftp_upload/63632/63632fig03.jpg" data-altupdatedat="1788282455000" data-alt="Rysunek 3"> 
Rysunek 3: Średnie objętości van der Waalsa (w Å3) bibliotek z siarką i bez siarki. Czarne słupki reprezentują średnie objętości bibliotek bez siarki (VAIL, DEST, Pro), natomiast żółte słupki reprezentują średnie objętości wersji tych bibliotek z podstawioną siarką (VAIL_S, DEST_S, Pro_S). Słupki błędów przedstawiają odchylenie standardowe. Skróty: VAIL = walina, alanina, izoleucyna, leucyna; DEST = kwas asparaginowy, kwas glutaminowy, seryna, treonina; X_S = dwuwartościowa siarka zastępuje jeden z atomów węgla w bibliotece X. <a href="https://www.jove.com/files/ftp_upload/63632/63632fig03large.jpg" target="_blank">Kliknij tutaj, aby wyświetlić większą wersję tego rysunku.</a>
<img alt="Wykres słupkowy porównujący wartości XLogP dla VAIL, DEST, Pro z słupkami błędów, przedstawiający hydrofobowość." class="xfigimg" src="/files/ftp_upload/63632/63632fig04.jpg" data-altupdatedat="1788282455000" data-alt="Rysunek 4"> 
Rysunek 4: Średnie wartości XLogP bibliotek z siarką i bez siarki. Czarne słupki reprezentują biblioteki bez siarki (VAIL, DEST, Pro), natomiast żółte słupki reprezentują wersje tych bibliotek z podstawnikiem siarkowym (VAIL_S, DEST_S, Pro_S). Słupki błędów przedstawiają odchylenie standardowe. Skróty: XLogP = współczynnik podziału; VAIL = walina, alanina, izoleucyna, leucyna; DEST = kwas asparaginowy, kwas glutaminowy, seryna, treonina; X_S = dwuwartościowa siarka zastępuje jeden z atomów węgla w bibliotece X. <a href="https://www.jove.com/files/ftp_upload/63632/63632fig04large.jpg" target="_blank">Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.</a>
<img alt="Wykresy słupkowe porównujące liczbę wygenerowanych struktur i czas wykonania dla VAIL i DEST, analiza w skali logarytmicznej." class="xfigimg" src="/files/ftp_upload/63632/63632fig05.jpg" data-altupdatedat="1788282455000" data-alt="Rysunek 5"> 
Rysunek 5: Wpływ trójwalentnego pseudoatomu na generowanie struktur w MAYGEN. Wszystkie testy przeprowadzono na komputerze PC z procesorem Intel i7-7700HQ o taktowaniu 2.8 GHz i 16 GB pamięci RAM, bez zapisywania struktur do pliku oraz z opcją -m w celu wykorzystania wielowątkowości. W testach z użyciem pseudoatomu zastosowano formuły rozmyte opisane w Tabeli 1. W testach bez pseudoatomu zastosowano te same formuły rozmyte, co w Tabeli 1, z następującymi zmianami: P zastąpiono N; liczbę atomów węgla zwiększono o 3; liczbę atomów wodoru zwiększono o 7; liczbę atomów tlenu zwiększono o 2. Czarne słupki przedstawiają biblioteki wygenerowane z pseudoatomem; szare słupki przedstawiają biblioteki wygenerowane bez pseudoatomu. (A) Liczba struktur wygenerowanych przy użyciu formuł rozmytych zastosowanych do budowy bibliotek VAIL i DEST z trójwalentnym fosforem zastępującym podstrukturę alaniny oraz bez niego. (B) Czas (w ms) potrzebny do zbudowania bibliotek VAIL i DEST z trójwalentnym fosforem zastępującym podstrukturę alaniny oraz bez niego. Skróty: VAIL = walina, alanina, izoleucyna, leucyna; DEST = kwas asparaginowy, kwas glutaminowy, seryna, treonina. <a href="https://www.jove.com/files/ftp_upload/63632/63632fig05large.jpg" target="_blank">Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.</a>
Plik uzupełniający 1: Notatnik przesiewowy podstruktur. <a href="https://www.jove.com/files/ftp_upload/63632/supple_1.zip" target="_blank">Aby pobrać ten plik, kliknij tutaj.</a>
Plik uzupełniający 2: Przykładowa lista błędów (badlist). <a href="https://www.jove.com/files/ftp_upload/63632/supple_2.zip" target="_blank">Kliknij tutaj, aby pobrać ten plik.</a>
Plik uzupełniający 3: Przykładowa lista prawidłowych próbek (goodlist). <a href="https://www.jove.com/files/ftp_upload/63632/supple_3.zip" target="_blank">Kliknij tutaj, aby pobrać ten plik.</a>
Plik uzupełniający 4: Notatnik zastępowania pseudoatomów. <a href="https://www.jove.com/files/ftp_upload/63632/supple_4.zip" target="_blank">Kliknij tutaj, aby pobrać ten plik.</a>
Plik uzupełniający 5: Notatnik dotyczący blokowania grup aminowych. <a href="https://www.jove.com/files/ftp_upload/63632/supple_5.zip" target="_blank">Kliknij tutaj, aby pobrać ten plik.</a>

Celem tego protokołu jest efektywne generowanie i zarządzanie bibliotekami struktur małych cząsteczek za pomocą oprogramowania open-source.

Kuratorstwo obliczeniowych bibliotek chemicznych zademonstrowanych za pomocą alfa-aminokwasów

Exhaustive generation of molecular structures has numerous chemical and biochemical applications such as drug design, molecular database construction, ...

Nasz protokół pokazuje, w jaki sposób oprogramowanie open source może pozwolić każdemu badaczowi na stworzenie i zarządzanie biblioteką struktur obliczeniowych. Atrakcyjność tego protokołu wynika z jego otwartości i elastyczności. Każdy może z niego korzystać i modyfikować go tak, aby pasował do jego konkretnego pytania badawczego.Wersje tego protokołu mogą być stosowane w aplikacjach związanych z odkrywaniem leków, szybko tworząc biblioteki specyficznych struktur do badań przesiewowych in silico. Chociaż protokół jest wyjaśniany krok po kroku, jeśli użytkownicy nie są zaznajomieni z Javą lub podstawowym kodowaniem, mogą najpierw zapoznać się z nimi przed wdrożeniem protokołu. Zacznij od utworzenia nowego katalogu dla projektu.Umieść wszystkie pliki i pliki wykonywalne w tym katalogu, aby mieć do nich łatwy dostęp. Pobierz najnowszą wersję Maygen jako plik jar i oprogramowanie do zarządzania pakietami Anaconda. W systemach Windows wyszukaj monit Anaconda i kliknij wynikowy skrót, aby uruchomić.Aby utworzyć środowisko RDKit w programie Anaconda i pobrać narzędzie RDKit do środowiska, wpisz polecenie wyświetlane na ekranie, naciśnij Enter, aby uruchomić i odpowiedzieć tak na wszelkie pytania, które pojawią się podczas instalacji. Następnie pobierz notesy Jupyter Notebook i pliki tekstowe wzorców podłoża z plików uzupełniających, od jednego do pięciu. W wierszu polecenia przejdź do katalogu zawierającego maygen.plik wykonywalny jar. Dla każdego wzoru chemicznego, który Cię interesuje, użyj polecenia pokazanego na ekranie, aby uruchomić Maygen. Jeśli formuła jest formułą rozmytą, a nie formułą dyskretną, zastąp flagę łącznika F flagą łącznika rozmytego i umieść wszystkie interwały elementów w nawiasach kwadratowych.W wierszu polecenia programu Anaconda przejdź do folderu zawierającego notesy Jupyter Notebook i aktywuj środowisko RDKit. Pobrane notesy wymagają narzędzia RDKit. Tak więc każde przyszłe użycie ich w tym protokole będzie wymagało ich otwarcia w środowisku RDKit.Następnie otwórz Jupyter Notebook w celu filtrowania podstruktury i zamknij nazwę pliku w cudzysłowie, jeśli zawiera spacje. W wyznaczonej komórce na początku notesu wprowadź pełną ścieżkę wejściowego pliku sdf. Pełna ścieżka do żądanego pliku wyjściowego sdf i ścieżka do pliku z listą uszkodzonych plików jako ciągi.Jeśli niektóre struktury podrzędne w filtrowanej bibliotece lub dobrej liście muszą zostać zachowane, utwórz plik txt z wzorcami SMARTS dla tych struktur podrzędnych i umieść ścieżkę do pliku dobrej listy w wyznaczonym wierszu na początku notesu. Z menu u góry wybierz pozycję jądro, uruchom ponownie i uruchom wszystko, aby ponownie uruchomić jądro notesu i uruchomić wszystkie komórki. Plik sdf o żądanej nazwie zostanie utworzony w określonym folderze wyjściowym.Powtórz te kroki dla każdego pliku struktury wygenerowanego przez Maygen. W celu zastąpienia pseudoatomu otwórz wiersz polecenia Anaconda, przejdź do folderu zawierającego notesy Jupyter Notebook i aktywuj środowisko RDKit. Następnie otwórz Jupyter Notebook w celu zastąpienia pseudoatomu.W wyznaczonej komórce na początku notesu wprowadź pełną ścieżkę wejściowego pliku sdf i pełną ścieżkę pliku wyjściowego sdf jako ciągi. Uruchom ponownie jądro notesu i uruchom wszystkie komórki, aby uzyskać plik sdf o żądanej nazwie w określonym folderze wyjściowym. Podobnie otwórz monit Anaconda dla limitowania końców aminokwasów N i C.Przejdź do folderu zawierającego notesy Jupyter Notebook i aktywuj środowisko RDKit. Otwórz Jupyter Notebook w celu ograniczenia aminokwasów. W wyznaczonej komórce na początku notesu wprowadź pełną ścieżkę wejściowego pliku sdf i pełną ścieżkę pliku wyjściowego sdf jako ciągi.Uruchom ponownie jądro notesu i uruchom wszystkie komórki, aby uzyskać plik sdf o żądanej nazwie w określonym folderze wyjściowym. Do generowania deskryptorów należy umieścić wszystkie pliki sdf, dla których mają być obliczane deskryptory, w jednym folderze. Następnie pobierz deskryptor PaDEL, rozpakuj go i rozpakuj do tego folderu.Otwórz wiersz polecenia, przejdź do folderu zawierającego plik jar deskryptora PaDEL i uruchom deskryptor PaDEL dla zebranych plików sdf. Przestrzeń chemiczna wszystkich filtrowanych bibliotek aminokwasów jest pokazana tutaj. Czarne znaczniki reprezentują aminokwasy z bibliotek bez siarki, a żółte znaczniki reprezentują aminokwasy z bibliotek wzbogaconych w siarkę.Tutaj biblioteki VAIL i VAIL_S są reprezentowane przez okręgi. Biblioteki DEST i DEST_S są reprezentowane przez kwadraty. Biblioteki Proline i Pro S są reprezentowane przez trójkąty, a gwiazdy reprezentują zakodowane aminokwasy.Zakres możliwych wartości logarytmu P zwiększa się wraz z objętością cząsteczkową nawet w bibliotekach, które wyraźnie nie mają hydrofilowych łańcuchów bocznych. Kodowane aminokwasy z węglowodorowymi łańcuchami bocznymi są bardziej hydrofobowe niż większość innych aminokwasów o porównywalnej objętości z odpowiedniej biblioteki. Dotyczy to również nalegania na metioninę w porównaniu z innymi członkami biblioteki VAILS o podobnych objętościach.Kodowane aminokwasy z hydroksylowymi łańcuchami bocznymi należały do najmniejszych członków biblioteki DEST, przy czym kwas asparaginowy był tylko nieznacznie większy niż trzy aniny. Przedstawiony obraz przedstawia średnie objętości bibliotek Van der Waala z siarką i bez siarki. Substytucja siarki doprowadziła do nieznacznego wzrostu objętości molekularnej we wszystkich bibliotekach.W tym miejscu pokazane są średnie wartości współczynników partycji bibliotek z siarką i bez siarki. Wpływ podstawienia siarki na log P nie jest tak jednorodny jak w przypadku objętości. Reprezentatywny obraz przedstawia wpływ trójwartościowego pseudoatomu na generowanie struktury Maygena.Użycie pseudoatomu w generowaniu struktur zmniejszyło liczbę generowanych struktur o około trzy rzędy wielkości w całkowitym czasie potrzebnym do wygenerowania tych struktur o jeden do dwóch rzędów wielkości. Zgodnie z tym protokołem, w przyszłości można zintegrować dodatkowe funkcjonalności w oparciu o potrzeby badaczy. Na przykład można zintegrować filtry podstruktury z Maygen, aby uniknąć etapu przetwarzania końcowego.Generowanie, selekcja i modyfikacja bibliotek. Ten ogólny proces może pomieścić inne struktury molekularne i modyfikacje z pewną wiedzą na temat kodowania, co pozwoli naukowcom badać biblioteki obliczeniowe wykraczające poza te dotyczące aminokwasów alfa. Protokół ten pomoże naukowcom usprawnić ich prace obliczeniowe w dziedzinie początków życia.Zestawy narzędzi open source znacznie pomogą w tych wysiłkach.

curation-computational-chemical-libraries-demonstrated-with-alpha

Research

JoVE Journal

Chemistry

Chemia

Synteza i charakterystyka trzeciorzędowej biblioteki amidów peptydów

Biblioteki wyświetlania bakterii w biopanowaniu to sprawdzona technika odkrywania odczynników powinowactwa peptydowego, solidna alternatywa dla przeciwciał. Półautomatyczna metoda sortowania usprawniła biopanoramowanie w celu zmniejszenia występowania wyników fałszywie dodatnich. Poniżej ilustrujemy proces myślowy i techniki stosowane w ocenie kandydatów i minimalizacji dalszej analizy.

Półautomatyczne biopanowanie bakteryjnych bibliotek wyświetlających w celu odkrywania odczynników powinowactwa peptydowego i analizy otrzymanych izolatów

Tworzenie chemicznie indukowanych systemów dimeryzacji białek o pożądanym powinowactwie i specyficzności dla dowolnego ligandu małych cząsteczek miałoby wiele zastosowań w biologicznym wykrywaniu i uruchamianiu. W tym miejscu opisujemy wydajną, uogólnioną metodę inżynierii de novo chemicznie indukowanych systemów dimeryzacji poprzez krokowy wybór kombinatorycznej biblioteki przeciwciał jednodomenowych wyświetlanej przez fagi.

Tworzenie wysoce specyficznych, chemicznie indukowanych systemów dimeryzacji białek poprzez stopniową selekcję fagów w kombinatorycznej bibliotece przeciwciał jednodomenowych

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

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.

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

Tutaj opisujemy protokół optymalizacji i parametryzacji reszt aminokwasowych modyfikowanych reaktywnymi gatunkami karbonylowymi, które można dostosować do systemów białkowych. Etapy protokołu obejmują projektowanie i optymalizację struktury, przypisywanie ładunków, konstrukcję parametrów i przygotowanie układów białkowych.

Synteza aminokwasów modyfikowanych reaktywnymi karbonylkami in silico w celu oceny efektów strukturalnych przy użyciu symulacji dynamiki molekularnej

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

Projektowanie białek polega na budowie sekwencji aminokwasów i włączaniu określonych motywów w celu stworzenia funkcjonalnych wariantów. Takie podejście ma kluczowe znaczenie dla opracowania peptydów przeciwdrobnoustrojowych (AMP) do zwalczania patogenów opornych na antybiotyki. W artykule przedstawiono procedurę budowy białek z wykorzystaniem różnych narzędzi bioinformatycznych.

Zastosowanie I TASSER, trRosetta, UCSF Chimera, serwera HADDOCK i lorii HEX do projektowania białek De Novo i In Silico

W tym artykule opisano AMOS, internetową bazę danych Analytical Methods and Open Spectra, aplikację cheminoinformatyczną zaprojektowaną w celu zapewnienia naukowcom łatwego dostępu do metod analitycznych i danych spektralnych.

Zastosowanie chemininformatyki do opracowania bazy danych metod analitycznych z możliwością przeszukiwania struktury

Poznaj metody kuratorowania obliczeniowych bibliotek chemicznych z wykorzystaniem oprogramowania open-source. Protokół umożliwia generowanie struktur alfa-aminokwasów i ich przesiewanie in silico.

Kuratorstwo obliczeniowych bibliotek chemicznych zademonstrowanych za pomocą alfa-aminokwasów

Przeglądaj więcej filmów