Name: curation of computational chemical libraries demonstrated with alpha-amino acids
Uploaded: 2022-04-13T13:30:00
Duration: 8 min 21 s
Description: Watch this Scientific Journal Video about Curation of Computational Chemical Libraries Demonstrated with Alpha-Amino Acids at JoVE.com

MAY признает финансирование со стороны Фонда Carl-Zeiss. Все рисунки были сгенерированы с помощью Microsoft Excel.

ПРИМЕЧАНИЕ: См. Рисунок 1 для краткого описания протокола и Таблицу материалов для получения подробной информации об используемом программном обеспечении.
<img alt="Figure 1" class="xfigimg" src="/files/ftp_upload/63632/63632fig01.jpg"> Рисунок 1: Сводная блок-схема протокола. <a href="https://www.jove.com/files/ftp_upload/63632/63632fig01large.jpg" target="_blank">Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.</a>
1. Загрузка программного обеспечения и файлов
ПРИМЕЧАНИЕ: Все программы бесплатны для индивидуального использования и могут быть запущены на персональном компьютере.
<ol><li>Создайте новый каталог для этого проекта. Разместите файлы и исполняемые файлы здесь для легкого доступа.</li>
	<li>Загрузите и установите необходимые пакеты программного обеспечения.<ol>
<li>Загрузите последнюю версию MAYGEN в виде файла .jar. ПРИМЕЧАНИЕ: MAYGEN находится в свободном доступе в виде файла .jar от https://github.com/MehmetAzizYirik/MAYGEN/releases</li>
		<li>Загрузите и установите программное обеспечение для управления пакетами Conda и инструментарий хеминформатики RDKit18. ПРИМЕЧАНИЕ: RDKit будет фильтровать молекулярные структуры, производимые MAYGEN, и лучше всего работает в среде Conda. Инструкции по загрузке платформы Conda можно найти на https://conda.io/projects/conda/en/latest/user-guide/install/index.html. Инструкции по установке и настройке среды RDKit можно найти в https://www.rdkit.org/docs/Install.html.<ol>
<li>Установите RDKit в основной среде Conda вместо отдельной среды RDKit с помощью приглашения Anaconda. В системах Windows найдите «Приглашение Anaconda» и нажмите на полученный ярлык для запуска. В системах MacOS и Linux взаимодействуйте с Conda через терминал без запуска каких-либо дополнительных программ. Затем введите следующую команду и нажмите клавишу ВВОД для запуска и ответьте «да» на все вопросы, возникающие во время установки: conda install -c rdkit rdkit. Хотя существует множество свободно доступных программ для вычисления дескрипторов, в этом примере используется PaDEL-Descriptor19, бесплатный и быстрый калькулятор для молекулярных дескрипторов и отпечатков пальцев.</li>
		</ol>
</li>
		<li>Загрузите и сохраните .jar файл в папке проекта. ПРИМЕЧАНИЕ: PaDEL-Дескриптор можно скачать бесплатно с http://www.yapcwsoft.com/dd/padeldescriptor/.</li>
	</ol>
</li>
	<li>Загрузите записные книжки Jupyter и текстовые файлы шаблонов подструктур из Дополнительных файлов 1-5. ПРИМЕЧАНИЕ: Записные книжки Jupyter также можно загрузить со следующей страницы GitHub: https://github.com/cmayerb1/AA-structure-manip.</li>
</ol>2. Генерация структуры с помощью MAYGEN
<ol><li>В командной строке перейдите в каталог, содержащий исполняемый файл MAYGEN .jar.</li>
	<li>Для каждой интересующей химической формулы запустите MAYGEN с помощью следующей команды: java -jar [MAYGEN .jar имя файла] -f [химическая формула] -v -o [папка для вывода MAYGEN] -m -sdf. ПРИМЕЧАНИЕ: Это сохранит файл .sdf в указанной папке, названной в честь используемой формулы.<ol>
<li>Если формула представляет собой нечеткую формулу, а не дискретную формулу, замените флаг -f флагом -f на флаг -fuzzy и заключите любые интервалы элементов в скобки (например, используйте C[5-7]H[12-15], чтобы гарантировать, что все генерируемые структуры имеют от 5 до 7 атомов углерода и от 12 до 15 атомов водорода).</li>
	</ol>
</li>
</ol>3. Фильтрующие соединения с нежелательными основаниями
<ol><li>Откройте приглашение Anaconda (см. шаг 1.2.2.1) и перейдите в папку, содержащую записные книжки Jupyter, загруженные из дополнительного файла 1.</li>
	<li>Откройте записную книжку Jupyter для фильтрации подструктуры с помощью следующей команды: jupyter notebook [имя файла записной книжки]
</li>
	<li>В назначенной ячейке в начале записной книжки введите полный путь к входному .sdf файлу (сгенерированному MAYGEN), полный путь к нужному выходному файлу .sdf и путь к файлу "badlist" в виде строк (в кавычках). Пример плохого списка см. в дополнительном файле 2 .<ol>
<li>Если некоторые подструктуры в отфильтрованной библиотеке (goodlist) должны быть сохранены, создайте .txt файл шаблонов SMARTS20 для этих подструктур (goodlist) и поместите путь к файлу goodlist в указанную строку в начале записной книжки. Пример хорошего списка см. в дополнительном файле 3 .</li>
	</ol>
</li>
	<li>Перезапустите ядро записной книжки и запустите все ячейки (из меню вверху выберите Kernel, Restart & Run All), чтобы получить файл .sdf с нужным именем в указанной выходной папке.</li>
	<li>Повторите предыдущие два шага для каждого файла структуры, созданного MAYGEN на шаге 2.</li>
</ol>4. (Необязательно) Дополнительные изменения структуры
ПРИМЕЧАНИЕ: Они выполняются в этом примере, но могут не понадобиться для курирования других библиотек.
<ol><li>Замена псевдоатома. ПРИМЕЧАНИЕ: Здесь псевдоатом является уникальным атомом, используемым для представления более крупной подструктуры, разделяемой всеми генерируемыми структурами, тем самым уменьшая время генерации MAYGEN. Пример замены псевдоатома см. в дополнительном файле 4 .<ol>
<li>Откройте приглашение Anaconda (см. шаг 1.2.2.1) и перейдите в папку, содержащую записные книжки Jupyter.</li>
		<li>Откройте записную книжку Jupyter для замены псевдоатома: jupyter notebook [имя файла записной книжки]
</li>
		<li>В назначенной ячейке в начале записной книжки введите полный путь к входному .sdf файлу и полный путь к файлу нужного .sdf выходного файла в виде строк (в кавычках).</li>
		<li>Перезапустите ядро записной книжки и запустите все ячейки, чтобы получить файл .sdf с нужным именем в указанной выходной папке.</li>
	</ol>
</li>
	<li>Аминокислоты N- и C-термины укупорочные ПРИМЕЧАНИЕ: Эта процедура специфична для альфа-аминокислот, добавляя молекулярные колпачки к N- и C-концам альфа-аминокислотных магистралей. Пример укупорки аминокислот см. в дополнительном файле 5 .<ol>
<li>Откройте приглашение Anaconda (см. шаг 1.2.2.1) и перейдите в папку, содержащую записные книжки Jupyter.</li>
		<li>Откройте блокнот Jupyter для укупорки аминокислот: jupyter notebook [имя файла записной книжки]
</li>
		<li>В назначенной ячейке в начале записной книжки введите полный путь к входному .sdf файлу и полный путь к файлу нужного .sdf выходного файла в виде строк (в кавычках).</li>
		<li>Перезапустите ядро записной книжки и запустите все ячейки, чтобы получить файл .sdf с нужным именем в указанной выходной папке.</li>
	</ol>
</li>
</ol>5. Генерация дескрипторов
<ol><li>Перед созданием дескриптора поместите все файлы .sdf, для которых дескрипторы должны быть рассчитаны, в одну папку. ПРИМЕЧАНИЕ: Если это еще не сделано, дайте этим файлам описательные имена для легкой фильтрации после генерации дескриптора.</li>
	<li>Откройте командную строку и перейдите в папку, содержащую файл .jar PaDEL-Descriptor.</li>
	<li>Запустите PaDEL-Descriptor для собранных файлов .sdf с помощью следующей команды: java -jar PaDEL-Descriptor.jar -dir [каталог .sdf файлов] -файл [путь к файлу .csv для результатов] -2d -retainorder -usefilenameasmolname ПРИМЕЧАНИЕ: Файл результатов будет иметь имя молекулы в первом столбце и каждый дескриптор в последующих столбцах.</li>
	<li>Экспортируйте эти данные в любое программное обеспечение для работы с электронными таблицами для дальнейшего анализа.</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/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Enumeration+of+166+billion+organic+small+molecules+in+the+chemical+universe+database+GDB-17.">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., Webber, B. L., Nilsson, N. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Dendral+and+Meta-Dendral:+their+applications+dimension.">Dendral and Meta-Dendral: their applications dimension.</a> Readings in Artificial Intelligence. , 313-322 (1981).</li><li>Gugisch, R., Basak, S. C., Restrepo, G., Villaveces, J. L., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=MOLGEN+5.0,+A+Molecular+Structure+Generator.">MOLGEN 5.0, A Molecular Structure Generator.</a> Advances in Mathematical Chemistry and Applications. , 113-138 (2015).</li><li>Yirik, M. A., Steinbeck, C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Chemical+graph+generators.">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/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=PMG:+multi-core+metabolite+identification.">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/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=MAYGEN:+an+open-source+chemical+structure+generator+for+constitutional+isomers+based+on+the+orderly+generation+principle.">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., Leech, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Computational+methods+in+the+study+of+permutation+groups.">Computational methods in the study of permutation groups.</a> Computational Problems in Abstract Algebra. , 169-183 (1970).</li><li>Mat, W. -. K., Xue, H., Wong, J. T. -. F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+genomics+of+LUCA.">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/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Ancestral+reconstruction+of+a+pre-LUCA+aminoacyl-tRNA+synthetase+ancestor+supports+the+late+addition+of+Trp+to+the+genetic+code.">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/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+Thermodynamic+basis+for+prebiotic+amino+acid+synthesis+and+the+nature+of+the+first+genetic+code.">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/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Consensus+temporal+order+of+amino+acids+and+evolution+of+the+triplet+code.">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/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+origin+of+the+biologically+coded+amino+acids.">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/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Redox+biochemistry+of+the+genetic+code.">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/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Methodologies+for+analyzing+soluble+organic+compounds+in+extraterrestrial+samples:+amino+acids,+amines,+monocarboxylic+acids,+aldehydes,+and+ketones.">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&#237;a-Ruiz, J. M., Saladino, R., Di Mauro, E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+role+of+borosilicate+glass+in+Miller-Urey+experiment.">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/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Primordial+synthesis+of+amines+and+amino+acids+in+a+1958+Miller+H2S-rich+spark+discharge+experiment.">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/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=New+insights+into+prebiotic+chemistry+from+Stanley+Miller's+spark+discharge+experiments.">New insights into prebiotic chemistry from Stanley Miller's spark discharge experiments.</a> Chemical Society Reviews. 42 (5), 2186-2196 (2013).</li><li>Yap, C. W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=PaDEL-descriptor:+An+open+source+software+to+calculate+molecular+descriptors+and+fingerprints.">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/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Beyond+terrestrial+biology:+charting+the+chemical+universe+of+&#945;-amino+acid+structures.">Beyond terrestrial biology: charting the chemical universe of &#945;-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/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Speciation+of+organosulfur+compounds+in+carbonaceous+chondrites.">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/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+hydrophobic+effect+and+the+organization+of+living+matter.">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/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Amino+acid+difference+formula+to+help+explain+protein+evolution.">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/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Covalent+radii+revisited.">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/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=One+among+millions:+the+chemical+space+of+nucleic+acid-like+molecules.">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>

У авторов нет конфликта интересов для раскрытия.

Одной из особенностей «ранних» аминокислот является недостаток серы. Мета-анализы, упомянутые ранее, обычно считают, что серосодержащие кодированные аминокислоты (Cys и Met) были относительно поздними дополнениями к генетическому коду, выводы, подкрепленные отсутствием серосодержащих аминокислот в метеоритах и экспериментах с искровыми трубками. Однако сероорганические соединения легко обнаруживаются в кометах и метеоритах22, и повторный анализ экспериментов с искровыми трубками с использованием газа H2S обнаружил аминокислоты и другие органические соединения, содержащие серу16. При рассмотрении альтернативного аминокислотного алфавита стоит изучить тот, который обогащен серой.
В приведенном выше протоколе генерация структуры и фильтрация подструктур считаются критическими этапами; в зависимости от состава готовой библиотеки структуры, исследователю может потребоваться только выполнить эти два шага. Инструкции и программное обеспечение для дополнительных действий (замена псевдоатомов и добавление подструктур (в данном случае укупорка аминокислот)) включены для более релевантного расчета дескриптора (укупорка гарантирует, что на расчеты XLogP влияет боковая цепь, а не аминные или карбоксильные группы) и более быстрой генерации структуры с помощью псевдоатома, который более подробно обсуждается ниже. Кроме того, здесь выполняется расчет дескриптора как простой способ визуализации разнообразия генерируемых структур и сравнения эффектов обогащения серы в готовых библиотеках.
В то время как PaDEL-Дескриптор может вычислять тысячи молекулярных свойств, молекулярный объем (как рассчитанный объем Ван-дер-Ваальса) и коэффициент разделения (как XLogP) использовались здесь по двум различным причинам. Во-первых, эти два дескриптора измеряют молекулярные свойства (размер и гидрофобность соответственно), которые знакомы большинству химиков и биологов. Во-вторых, в случае аминокислот эти два свойства являются значительными. В течение десятилетий было известно, что размер аминокислот и гидрофобность влияют на термодинамику сворачивания белка23. Эти два свойства помогают объяснить частоты замещения аминокислот, которые были неотъемлемой частью понимания эволюции белка24.
Приведенный выше пример показывает, что в двух изученных дескрипторах (молекулярный объем и гидрофобность) замена двухвалентной серы на углерод и два водорода не дает существенных изменений. Незначительное, незначительное увеличение среднего молекулярного объема от замещения серы (рисунок 3) может быть связано с большим ковалентным радиусом серы (~ 103 пм) по сравнению с sp3 (~ 75 pm) илиsp2 (~ 73 pm) углерода25. Аналогичным образом, замещение серы оказывает минимальное влияние на среднее значение XLogP (рисунок 4). Наибольший эффект был между библиотеками VAIL и VAIL_S, вероятно, из-за того, что сочетание библиотеки VAIL было особенно гидрофобным (боковые цепи являются только углеводородами), а сульфгидрильные группы были гораздо более кислыми, чем метильные группы, которые они заменили бы. Минимальный эффект замещения серы очевиден на рисунке 2, где библиотеки с замещением серы занимают то же химическое пространство, что и аналогичные библиотеки без замещения серы.
Уменьшение количества структур (рисунок 5A) и времени, необходимого для генерации этих структур (рисунок 5B) при использовании псевдоатома, неудивительно. Использование псевдоатома уменьшает количество тяжелых атомов, которые необходимо включить в химический граф, уменьшая количество узлов графа и приводя к экспоненциальному сокращению времени генерации и числа структур. Здесь выбор трехвалентного фосфора в качестве псевдоатома проистекает из основной биохимии (отсутствие посттрансляционного добавления фосфатных групп, ни одна генетически закодированная аминокислота не содержит фосфора) и валентность атома, который бы его заменил (трехвалентный фосфор может быть легко заменен четырехвалентным углеродом, который отдельно связан с другим атомом или группой атомов). Хотя предоставленный код для замещения псевдоатомов специфичен для замены трехвалентного фосфора аланиновой субструктурой, пользователи могут настроить код для работы с различными псевдоатомами или замещающими подструктурами, потенциально используя несколько псевдоатомов во время первоначальной генерации структуры с последующей заменой каждого псевдоатома более крупной молекулярной подструктурой.
Методы генерации структуры, аналогичные тем, которые используются MAYGEN (и другими методами, такими как нейронные сети), уже используются в открытии лекарств для создания библиотек соединений для скрининга in silico ; в недавнем обзоре4 эти методы рассматриваются более подробно. Поскольку эти методы предназначены в первую очередь для создания лекарственно-подобных молекул, существуют некоторые ограничения на их способность генерировать молекулы, такие как использование биологических или фармацевтических свойств для ограничения создаваемых структур (обратный QSPR / QSAR) или создание структур из заданного числа строительных блоков подструктуры. Поскольку астробиология больше сосредоточена на множестве органических соединений, которые могут образовываться абиотически, и меньше на любых конечных продуктах или их свойствах, исчерпывающая генерация структур MAYGEN идеально подходит для создания структурных библиотек для решения астробиологических вопросов. Описанный здесь подход к фильтрации подструктур (выполняемый после генерации структуры с помощью внешней программы) отличается от программы конкурента MOLGEN тем, что фильтрация подструктуры MOLGEN происходит во время генерации структуры. Поскольку MAYGEN является открытым исходным кодом, он не только более доступен, чем MOLGEN, из-за стоимости лицензирования MOLGEN, но и отдельные лица могут внедрять новые функции, такие как фильтрация подструктур во время генерации структуры.
Как написано, протокол, описанный здесь, ориентирован на генерацию и курирование библиотек относительно небольших альфа-аминокислот. Для создания различных библиотек пользователи могут давать различные молекулярные формулы MAYGEN, изменять фильтрацию подструктуры, изменяя максимально допустимый размер кольца и валентность связи, или редактировать файлы goodlist и badlist для добавления или удаления шаблонов подструктуры. Модификации протокола, которые включают изменение способа добавления или замены атомов и подструктур (замена псевдоатомов и молекулярное ограничение), возможны, но потребуют большего внимания к валентным ограничениям, чтобы избежать ошибок RDKit о неправильных валентностях в модифицированных структурах.
Протокол, описанный выше, предназначен для малых альфа-аминокислот. Однако общий формат (комплексная генерация структуры с использованием псевдоатомов с последующей фильтрацией субструктур и молекулярными модификациями) является очень гибким для соединений, выходящих за рамки небольших аминокислот. Даже в астробиологии аналогичная недавняя процедура с использованием MOLGEN использовалась для исследования конституциональных изомеров нуклеиновых кислот26. В дополнение к инструментам, описанным выше, MAYGEN может быть сопряжен с другими инструментами хеминформатики с открытым исходным кодом, чтобы сделать создание и анализ новых химических структур доступными и доступными для широкого спектра областей исследований.

Генерация молекулярной структуры служит практическим применением общей проблемы генерации исчерпывающих графов; Учитывая несколько узлов (атомов) и ограничения на их связность (например, валентности, кратности связей, желаемые/нежелательные подструктуры), сколько связанных графов (молекул) возможно? Структурные генераторы нашли широкое применение в открытии лекарств и фармацевтических разработках, где они могут создавать обширные библиотеки новых структур для скрининга in silico 1.
Первый генератор структуры, CONGEN, был разработан для первого проекта искусственного интеллекта в органической химии, DENDRAL2 (сокращение от DENDRitic ALgorithm). В литературе сообщалось о нескольких программных преемниках DENDRAL; однако не все из них были сохранены или эффективны. В настоящее время MOLGEN3 является современным генератором молекулярной структуры. К сожалению для большинства потенциальных пользователей, он является закрытым исходным кодом и требует лицензионного сбора. Таким образом, возникла потребность в эффективном генераторе структуры с открытым исходным кодом, который может легко адаптироваться к конкретным приложениям. Одной из проблем для эффективного генератора структуры является управление комбинаторным взрывом; по мере увеличения размера молекулярной формулы размер пространства химического поиска увеличивается экспоненциально. Недавний обзор дополнительно исследует историю и проблемы молекулярной структуры поколения4.
До 2021 года генератор параллельных молекул (PMG)5 был самым быстрым генератором структуры с открытым исходным кодом, но он все еще был медленнее, чем MOLGEN на порядки. MAYGEN6 примерно в 47 раз быстрее, чем PMG и примерно в 3 раза медленнее, чем MOLGEN, что делает MAYGEN самым быстрым и эффективным генератором структуры с открытым исходным кодом. Более подробные сравнения и сравнительные тесты можно найти в статье, представляющей MAYGEN6. Ключевой особенностью программы является ее лексикографический тест на основе упорядочения канонических структур, метод упорядоченной генерации графов, основанный на алгоритме Шрайера-Симса7 . Программное обеспечение может быть легко интегрировано в другие проекты и улучшено для нужд пользователей.
Как и MOLGEN и PMG, MAYGEN принимает определяемую пользователем молекулярную формулу и генерирует все структуры, возможные для этой формулы. Например, если пользователь запустит MAYGEN с формулой C5H12, MAYGEN сгенерирует все возможные структуры, содержащие пять атомов углерода и двенадцать атомов водорода. В отличие от своего аналога с открытым исходным кодом PMG, MAYGEN также может вместить «нечеткие» молекулярные формулы, которые используют интервалы вместо дискретных чисел для подсчета каждого элемента. Например, если пользователь запустит MAYGEN с формулой C5-7H12-15, MAYGEN сгенерирует все возможные структуры, которые содержат от пяти до семи атомов углерода и двенадцать и пятнадцать атомов водорода, что позволит просто генерировать структуры с широким диапазоном атомных составов.
Астробиология является одной из таких областей, которая может извлечь выгоду из генераторов молекулярной структуры. Популярной темой в астробиологии является эволюция аминокислотного алфавита, разделяемого всей существующей жизнью на Земле. Одной из определяющих особенностей последнего универсального общего предка (LUCA) является использование двадцати генетически закодированных аминокислот для построения белка 8,9. Основываясь на мета-анализе работы в нескольких областях 10,11,12, примерно 10 из этих аминокислот (Gly, Ala, Val, Asp, Glu, Ser, Thr, Leu, Ile, Pro) легко образуются в абиотических условиях и, вероятно, составляют аминокислотный алфавит организмов до LUCA. Со временем этот «ранний» алфавит был расширен в ответ на различные структурные и функциональные потребности. Например, в недавнем обзоре Moosmann13 утверждается, что добавление более поздних членов генетически закодированных аминокислот (а именно Met, Tyr и Trp) позволило выжить в богатых кислородом средах, предотвращая внутриклеточную пролиферацию активных форм кислорода.
Постоянно растущий набор методов аналитической химии позволяет понять аминокислотные структуры, которые могут образовываться в абиотических условиях. В недавнем обзоре14 Симкуса и других подробно описываются методы, используемые для обнаружения многочисленных органических соединений в метеоритах, а также органических соединений из моделирования in vitro ранних сред Земли 15,16,17. Систематическая генерация химических структур позволяет исследователям исследовать за пределами органических соединений, обнаруженных с помощью приборов, заполняя структурное пространство вокруг структурных «островов», идентифицированных аналитической химией. В случае «ранних» аминокислот эта систематическая генерация структуры показывает возможную химию белка, доступную для ранней жизни, не ограничивая исследование структурами, которые были экспериментально обнаружены в условиях абиотического синтеза. С инструментами хеминформатики с открытым исходным кодом и эффективными генераторами структур, такими как MAYGEN, создание и изучение новых библиотек химических структур теперь проще, чем когда-либо прежде, и может направлять более подробные исследования альтернативных химических веществ жизни.

<table border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always">
	<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>included in Anaconda environment</td>
		</tr>
		<tr>
			<td>RDKit v. 2020.09.1.0</td>
			<td></td>
			<td></td>
			<td>https://www.rdkit.org/docs/Install.html, or installed via conda: https://anaconda.org/rdkit/rdkit</td>
		</tr>
		<tr>
			<td>*These specific versions were used for this manuscript; user can obtain more recent versions if available.</td>
			<td></td>
			<td></td>
			<td></td>
		</tr>
	</tbody>
</table>

curation of computational chemical libraries demonstrated with alpha-amino acids

Исчерпывающая генерация молекулярных структур имеет многочисленные химические и биохимические применения, такие как разработка лекарств, построение молекулярных баз данных, исследование альтернативных биохимий и многое другое. Математически говоря, это графовые генераторы с химическими ограничениями. В полевых условиях наиболее эффективным генератором в настоящее время (MOLGEN) является коммерческий продукт, ограничивающий его использование. Альтернативой этому является другой генератор молекулярной структуры, MAYGEN, является недавним инструментом с открытым исходным кодом с эффективностью, сопоставимой с MOLGEN, и способностью пользователей повышать его производительность за счет добавления новых функций. Одной из областей исследований, которые могут извлечь выгоду из этого развития, является астробиология; Структурные генераторы позволяют исследователям дополнять экспериментальные данные вычислительными возможностями для альтернативной биохимии. Этот протокол подробно описывает один из вариантов использования генерации структуры в астробиологии, а именно генерацию и курирование библиотек альфа-аминокислот. Используя генераторы структур с открытым исходным кодом и инструменты хеминформатики, описанные здесь практики могут быть реализованы за пределами астробиологии для недорогого создания и курирования библиотек химических структур для любого исследовательского вопроса.

<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">Библиотека</td>
			<td colspan="2" rowspan="2">Формула</td>
			<td colspan="2" rowspan="2">Дополнительные ограничения</td>
			<td rowspan="2">«Ранние» кодированные аминокислоты</td>
			<td rowspan="2">Время генерации (мс)</td>
			<td colspan="2">Структуры</td>
		</tr>
<tr>
<td>Начальный</td>
			<td>Последний</td>
		</tr>
<tr>
<td>1</td>
			<td>Гли</td>
			<td colspan="2">C2Ч5NO2
</td>
			<td colspan="2">включить подструктуру Gly</td>
			<td>Гли</td>
			<td>192</td>
			<td>84</td>
			<td>1</td>
		</tr>
<tr>
<td>2</td>
			<td>СКЛОНЯТЬ</td>
			<td colspan="2">ПК0-3Ч3-9
</td>
			<td colspan="2"></td>
			<td>Валь, Ала, Иль, Лей</td>
			<td>172</td>
			<td>70</td>
			<td>22</td>
		</tr>
<tr>
<td>3</td>
			<td>ДЕСТ</td>
			<td colspan="2">ПК0-3O1-2Ч3-5
</td>
			<td colspan="2"></td>
			<td>Жерех, Клей, Сер, Тр</td>
			<td>481</td>
			<td>1928</td>
			<td>254</td>
		</tr>
<tr>
<td>4</td>
			<td>Профессионал</td>
			<td colspan="2">С2-5NO2Ч7-11
</td>
			<td colspan="2">Включить N-meGly или N-meAla подструктуру</td>
			<td>Профессионал</td>
			<td>4035</td>
			<td>79777</td>
			<td>16</td>
		</tr>
<tr>
<td>5</td>
			<td>VAIL_S</td>
			<td colspan="2">ЦОН0-2Ч3-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">ЦОН0-2О1-2Ч3
</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-4СНО2Ч7-9
</td>
			<td colspan="2">Включить N-meGly или N-meAla подструктуру</td>
			<td></td>
			<td>3999</td>
			<td>75734</td>
			<td>10</td>
		</tr>
</tbody></table>Таблица 1: Составные библиотеки, используемые в этом примере. Библиотеки, построенные из формул 1-4 (Gly, VAIL, DEST и Pro), основаны на ранее опубликованных нечетких формулах «ранних» кодированных аминокислот21, в то время как библиотеки, построенные из формул 5-7 (VAIL_S, DEST_S и Pro_S), основаны на вариантах формул 2-4, которые представляют собой двухвалентную серу, заменяющую один из атомов углерода. Количество структур отражает количество молекул, генерируемых MAYGEN для каждой формулы («Initial»), и количество молекул, оставшихся после фильтрации молекул с нежелательными подструктурами («Final»). Сокращения: VAIL = валин, аланин, изолейцин, лейцин; DEST = аспарагиновая кислота, глутаминовая кислота, серин, треонин; X_S = Двухвалентная сера заменяет один из атомов углерода в библиотеке X; N-meX = N-метилX.
Общие методы, приведенные выше, были применены к формулам, основанным на «ранних» кодированных аминокислотах, следуя процедуре Meringer et al.21 Структуры Badlist были взяты из этого же источника и преобразованы в строки SMARTS для легкого представления субструктурных паттернов. В этом примере не использовались две подструктуры badlist: структура 018 (CH 3-CH-N) соответствовала почти изомерам пролина, которые сами по себе не были нестабильными; структура 106 (R-C-C-OH, где R=аланиновая субструктура, присоединяющаяся к бета-углероду) соответствовала глутаминовой кислоте, кодированной аминокислоте. В дополнение к этим химическим формулам были созданы варианты с двухвалентной серой, занимающей место атома углерода и двух атомов водорода. По соображениям производительности некоторые из этих формул используют трехвалентный атом фосфора (например, «псевдоатом») в качестве заменителя бета-углерода подструктуры аланина. В таблице 1 перечислены библиотеки, созданные в этом примере, формулы, используемые для их создания, и количество соединений, содержащихся внутри. Названия библиотек основаны на кодированных аминокислотах, из которых они получены: либо с использованием 3-буквенной аббревиатуры (Gly = глицин, Pro = пролин), либо однобуквенной аббревиатуры (VAIL = Valine, Alanine, Isoleucine, Leucine; DEST = аспарагиновая кислота, глутаминовая кислота, серин, треонин). Суффикс «_S» указывает на то, что сера была заменена на углерод в формуле оригинальной библиотеки (например, VAIL_S построена с той же нечеткой формулой, что и VAIL, но с двухвалентной серой, заменяющей один из атомов углерода).
После генерации структуры с помощью MAYGEN полученные библиотеки фильтровали соединения, содержащие, по меньшей мере, одну подструктуру, содержащуюся в плохом списке. После этой фильтрации любые атомы фосфора были заменены подструктурой аланина. Затем были созданы «закрытые» версии всех структур, с ацетильной группой, добавленной к N-концу, и N-метиламидной группой, добавленной к С-концу. Это было сделано для устранения влияния на гидрофобность групп свободных аминов и карбоновых кислот в альфа-аминокислотной основе. PaDEL-Descriptor использовался для расчета XLogP для всех закрытых структур и расчета объема Ван-дер-Ваальса (VABC) для всех безрамочных структур.
На рисунке 2 показано химическое пространство отфильтрованных библиотек, как определено дескрипторами VABC и XLogP. Здесь диапазон возможных значений logP увеличивается с молекулярным объемом, даже в библиотеках, в которых отсутствуют явно гидрофильные боковые цепи (например, VAIL, Pro). Закодированные аминокислоты с углеводородными боковыми цепями были более гидрофобными, чем большинство других аминокислот сопоставимого объема из их соответствующей библиотеки. Это также, по-видимому, относится к Met and Cys по сравнению с другими членами библиотеки VAIL_S с аналогичными томами. Кодированные аминокислоты с гидроксильными боковыми цепями (Ser и Thr) были одними из самых маленьких членов библиотеки DEST, причем Asp лишь немного больше, чем Thr.
На рисунках 3 и 4 показано влияние на объем и logP, когда двухвалентная сера заменяет углерод в боковой цепи альфа-аминокислоты. Замещение серы привело к небольшому увеличению молекулярного объема во всех библиотеках (рисунок 3). Влияние замещения серы на logP не так однородно, как для объема (рисунок 4). Средний logP библиотеки VAIL_S немного ниже, чем у библиотеки VAIL, но этот эффект не наблюдается ни в одной из других пар библиотек (DEST и DEST_S, Pro и Pro_S).
На рисунке 5 количественно определено влияние на структурную генерацию псевдоатома, стоящего за общей подструктурой; здесь трехвалентный Р заменяет фрагмент аланина во время генерации структуры. Использование псевдоатома в генерации структур значительно уменьшило количество структур, генерируемых ~3 порядками величины (рисунок 5A), и общее время, необходимое для генерации этих структур, на 1-2 порядка (рисунок 5B).
<img alt="Figure 2" class="xfigimg" src="/files/ftp_upload/63632/63632fig02.jpg"> Рисунок 2: Химическое пространство всех фильтрованных библиотек аминокислот. Черные маркеры представляют аминокислоты из библиотек без серы; желтые маркеры представляют аминокислоты из библиотек, обогащенных серой. Круги: VAIL и VAIL_S; квадраты: DEST и DEST_S; треугольники: Pro и Pro_S; звезды: кодированные аминокислоты. Обратите внимание, что две серосодержащие кодированные аминокислоты (Met и Cys) не считаются «ранними» аминокислотами, но присутствуют в библиотеке VAIL_S. Сокращения: XLogP = коэффициент разбиения; VAIL = валин, аланин, изолейцин, лейцин; DEST = аспарагиновая кислота, глутаминовая кислота, серин, треонин; X_S = Двухвалентная сера заменяет один из атомов углерода в библиотеке X. <a href="https://www.jove.com/files/ftp_upload/63632/63632fig02large.jpg" target="_blank">Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.</a>
<img alt="Figure 3" class="xfigimg" src="/files/ftp_upload/63632/63632fig03.jpg"> Рисунок 3: Среднее значение томов ван дер Ваальса (в Å3) библиотек с серой и без нее. Черные полосы представляют средние объемы библиотек без серы (VAIL, DEST, Pro), в то время как желтые полосы представляют средние объемы серозамещенных версий этих библиотек (VAIL_S, DEST_S Pro_S). На панелях ошибок отображается стандартное отклонение. Сокращения: VAIL = валин, аланин, изолейцин, лейцин; DEST = аспарагиновая кислота, глутаминовая кислота, серин, треонин; X_S = Двухвалентная сера заменяет один из атомов углерода в библиотеке X. <a href="https://www.jove.com/files/ftp_upload/63632/63632fig03large.jpg" target="_blank">Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.</a>
<img alt="Figure 4" class="xfigimg" src="/files/ftp_upload/63632/63632fig04.jpg"> Рисунок 4: Средние значения XLogP библиотек с серой и без нее. Черные полосы представляют библиотеки без серы (VAIL, DEST, Pro), в то время как желтые полосы представляют серозамещенные версии этих библиотек (VAIL_S, DEST_S, Pro_S). На панелях ошибок отображается стандартное отклонение. Сокращения: XLogP = коэффициент разбиения; VAIL = валин, аланин, изолейцин, лейцин; DEST = аспарагиновая кислота, глутаминовая кислота, серин, треонин; X_S = Двухвалентная сера заменяет один из атомов углерода в библиотеке X. <a href="https://www.jove.com/files/ftp_upload/63632/63632fig04large.jpg" target="_blank">Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.</a>
<img alt="Figure 5" class="xfigimg" src="/files/ftp_upload/63632/63632fig05.jpg"> Рисунок 5: Влияние трехвалентного псевдоатома на генерацию структуры MAYGEN. Все тесты проводились на ПК с процессором Intel i7-7700HQ с тактовой частотой 2,8 ГГц, 16 ГБ оперативной памяти, отсутствием сохранения структур в файл и опцией -m для использования многопоточности. Тесты с использованием псевдоатома использовали нечеткие формулы, как описано в таблице 1. Для испытаний без псевдоатома использовались нечеткие формулы, аналогичные описанным в таблице 1 со следующими изменениями: P был заменен на N; количество углерода было увеличено на 3; количество водорода было увеличено на 7; количество кислорода было увеличено на 2. Черные полосы показывают библиотеки, сгенерированные псевдоатомом; серые полосы показывают библиотеки, созданные без псевдоатома. (A) Число структур, созданных с использованием нечетких формул, используемых для построения библиотек VAIL и DEST с трехвалентным фосфором, заменяющим подструктуру аланина, и без нее. (B) Время (в мс), необходимое для создания библиотек VAIL и DEST с трехвалентным фосфором, заменяющим подструктуру аланина, и без нее. Сокращения: VAIL = валин, аланин, изолейцин, лейцин; DEST = аспарагиновая кислота, глутаминовая кислота, серин, треонин. <a href="https://www.jove.com/files/ftp_upload/63632/63632fig05large.jpg" target="_blank">Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.</a>
Дополнительный файл 1: Блокнот для скрининга подструктуры. <a href="https://www.jove.com/files/ftp_upload/63632/supple_1.zip" target="_blank">Пожалуйста, нажмите здесь, чтобы загрузить этот файл.</a>
Дополнительный файл 2: Образец плохого списка. <a href="https://www.jove.com/files/ftp_upload/63632/supple_2.zip" target="_blank">Пожалуйста, нажмите здесь, чтобы загрузить этот файл.</a>
Дополнительный файл 3: Пример списка товаров. <a href="https://www.jove.com/files/ftp_upload/63632/supple_3.zip" target="_blank">Пожалуйста, нажмите здесь, чтобы загрузить этот файл.</a>
Дополнительный файл 4: Блокнот для замены псевдоатомов. <a href="https://www.jove.com/files/ftp_upload/63632/supple_4.zip" target="_blank">Пожалуйста, нажмите здесь, чтобы загрузить этот файл.</a>
Дополнительный файл 5: Блокнот для укупорки аминокислот. <a href="https://www.jove.com/files/ftp_upload/63632/supple_5.zip" target="_blank">Пожалуйста, нажмите здесь, чтобы загрузить этот файл.</a>

Watch this Scientific Journal Video about Curation of Computational Chemical Libraries Demonstrated with Alpha-Amino Acids at JoVE.com

Curation of Computational Chemical Libraries Demonstrated with Alpha-Amino Acids

Целью этого протокола является эффективное создание и курирование библиотек структуры малых молекул с использованием программного обеспечения с открытым исходным кодом.

Курирование вычислительных химических библиотек, демонстрируемых с помощью альфа-аминокислот

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

Наш протокол демонстрирует, как программное обеспечение с открытым исходным кодом может позволить любому исследователю создать и курировать библиотеку вычислительной структуры. Эта апелляция к протоколам исходит из его открытости и гибкости. Любой может использовать его и модифицировать в соответствии со своим конкретным исследовательским вопросом.Версии этого протокола могут быть применены к приложениям для обнаружения лекарств, быстро создавая специальные структурные библиотеки для скрининга in silico. Хотя протокол объясняется шаг за шагом, если пользователи не знакомы с Java или базовым кодированием, они могут сначала взглянуть на них перед реализацией протокола. Начните с создания нового каталога для проекта.Поместите все файлы и исполняемые файлы в этот каталог для легкого доступа. Загрузите последнюю версию Maygen в виде jar-файла и программного обеспечения для управления пакетами Anaconda. В системах Windows найдите подсказку Anaconda и нажмите на полученный ярлык для запуска.Чтобы создать среду RDKit в Anaconda и загрузить RDKit в среду, введите команду, показанную на экране, нажмите клавишу ВВОД, чтобы запустить и ответить «да» на любые вопросы, возникающие во время установки. Затем загрузите записные книжки Jupyter и текстовые файлы шаблонов подложки из дополнительных файлов, от одного до пяти. В командной строке перейдите в каталог, содержащий maygen.исполняемый файл jar. Для каждой интересующей химической формулы используйте команду, показанную на экране, чтобы запустить Maygen. Если формула является нечеткой формулой, а не дискретной формулой, замените флаг дефиса F флагом нечеткого дефиса и заключите любые интервалы элементов в скобки.В командной строке Anaconda перейдите в папку, содержащую записные книжки Jupyter, и активируйте среду RDKit. Для загруженных ноутбуков требуется RDKit. Таким образом, любое будущее использование их в этом протоколе потребует их открытия в среде RDKit.Затем откройте записную книжку Jupyter для фильтрации подструктуры и закройте имя файла в кавычках, если оно содержит пробелы. В назначенной ячейке в начале записной книжки введите полный путь к входному sdf-файлу. Полный путь к нужному выходному файлу SDF и путь к файлу плохого списка в виде строк.Если необходимо сохранить некоторые подструктуры в отфильтрованной библиотеке или хороший список, создайте txt-файл шаблонов SMARTS для этих подструктур и поместите путь к файлу хорошего списка в назначенную строку в начале записной книжки. В меню вверху выберите ядро, перезапустите и запустите все, чтобы перезапустить ядро ноутбука и запустить все ячейки. SDF-файл с нужным именем будет создан в указанной выходной папке.Повторите эти шаги для каждого файла структуры, созданного Maygen. Для замены псевдоатома откройте приглашение Anaconda, перейдите в папку, содержащую записные книжки Jupyter, и активируйте среду RDKit. Затем откройте записную книжку Jupyter для замены псевдоатома.В назначенной ячейке в начале записной книжки введите полный путь к входному sdf-файлу и полный путь к нужному выходному файлу sdf в виде строк. Перезапустите ядро записной книжки и запустите все ячейки, чтобы получить SDF-файл с нужным именем в указанной выходной папке. Аналогичным образом, откройте подсказку Анаконды для аминокислот N и C терминов.Перейдите в папку, содержащую записные книжки Jupyter, и активируйте среду RDKit. Откройте блокнот Jupyter для укупорки аминокислот. В назначенной ячейке в начале записной книжки введите полный путь к входному sdf-файлу и полный путь к нужному выходному файлу sdf в виде строк.Перезапустите ядро записной книжки и запустите все ячейки, чтобы получить SDF-файл с нужным именем в указанной выходной папке. Для генерации дескриптора поместите все SDF-файлы, для которых дескрипторы должны быть вычислены, в одну папку. Затем загрузите дескриптор PaDEL, распакуйте его и извлеките в эту папку.Откройте командную строку, перейдите в папку, содержащую JAR-файл дескриптора PaDEL, и запустите дескриптор PaDEL для собранных SDF-файлов. Здесь показано химическое пространство всех фильтрованных библиотек аминокислот. Черные маркеры представляют аминокислоты из библиотек без серы, а желтые маркеры представляют аминокислоты из библиотек, обогащенных серой.Здесь библиотеки VAIL и VAIL_S представлены кругами. Библиотеки DEST и DEST_S представлены квадратами. Библиотеки Proline и Pro S представлены треугольниками, а звезды представляют закодированные аминокислоты.Диапазон возможных значений log P увеличивается с молекулярным объемом даже в библиотеках, в которых явно отсутствуют гидрофильные боковые цепи. Кодированные аминокислоты с углеводородными боковыми цепями более гидрофобны, чем большинство других аминокислот сопоставимого объема из их соответствующей библиотеки. Это также относится к метионину, настаивающему на сравнении с другими членами библиотеки VAILS с аналогичными томами.Кодированные аминокислоты с гидроксильными боковыми цепями были одними из самых маленьких членов библиотеки DEST с аспарагиновой кислотой, лишь немного превышающей три анина. На представленном изображении показаны средние тома библиотек Ван-дер-Ваала с серой и без серы. Замещение серы привело к небольшому увеличению молекулярного объема во всех библиотеках.Здесь показаны средние значения коэффициента разбиения библиотек с серой и без нее. Влияние замещения серы на log P не так однородно, как для объема. Репрезентативное изображение показывает влияние трехвалентного псевдоатома на генерацию структуры Майгена.Использование псевдоатома в генерации структур уменьшило количество структур, генерируемых примерно на три порядка за общее время, необходимое для создания этих структур, на один-два порядка. Следуя этому протоколу, дополнительные функциональные возможности могут быть интегрированы в будущем на основе потребностей исследователей. Например, можно интегрировать фильтры подструктуры в Maygen, чтобы избежать этапа постобработки.Генерация, курирование и модификация библиотек. Этот общий процесс может вместить другие молекулярные структуры и модификации с некоторыми знаниями кодирования, что позволит исследователям исследовать вычислительные библиотеки, выходящие за рамки библиотек альфа-аминокислот. Этот протокол поможет исследователям улучшить свою вычислительную работу в области происхождения жизни.Инструментарий с открытым исходным кодом будет в значительной степени способствовать этим усилиям.

Research

JoVE Journal

Chemistry

Free Radicals in Chemical Biology: from Chemical Behavior to Biomarker Development

Свободные радикалы в химической биологии: от химического поведения в биомаркеров развития

The involvement of free radicals in life sciences has constantly increased with time and has been connected to several physiological and pathological ...

Qualitative Identification of Carboxylic Acids, Boronic Acids, and Amines Using Cruciform Fluorophores

Качественная идентификация карбоновые кислоты, борной кислоты и аминов с использованием крестообразных Флуорофоры

Molecular cruciforms are X-shaped systems in which two conjugation axes intersect at a central core. If one axis of these molecules is substituted with ...

Sulfate Separation by Selective Crystallization with a Bis-iminoguanidinium Ligand

Сульфат Разделение селективной кристаллизацией с Бис-iminoguanidinium-лигандом

A simple and effective method for selective sulfate separation from aqueous solutions by crystallization with a bis-guanidinium ligand, ...

Insights into the Interactions of Amino Acids and Peptides with Inorganic Materials Using Single-Molecule Force Spectroscopy

Исследования в взаимодействий аминокислот и пептидов с использованием неорганических материалов с использованием одномолекулярной силовой микроскопии

The interactions between proteins or peptides and inorganic materials lead to several interesting processes. For example, combining proteins with minerals ...

Determining the Chemical Composition of Corrosion Inhibitor/Metal Interfaces with XPS: Minimizing Post Immersion Oxidation

Определение химического состава ингибитора коррозии / Metal Интерфейсы с XPS: Минимизация Сообщение Погружение Окисление

An approach for acquiring more reliable X-ray photoelectron spectroscopy data from corrosion inhibitor/metal interfaces is described. More specifically, ...

Facile Preparation of (2Z,4E)-Dienamides by the Olefination of Electron-deficient Alkenes with Allyl Acetate

Facile Preparation of 2Z,4E-Dienamides by the Olefination of Electron-deficient Alkenes with Allyl Acetate

Подготовка фаски (2<em&gt; Z</em&gt; 4<em&gt; E</em&gt;) - Диенамиды путем олефинирования электронно-дефицитных алкенов с аллилацетатом

Direct cross-coupling between two alkenes via vinylic C-H bond activation represents an efficient strategy for the synthesis of butadienes with high ...

Fizzy Extraction of Volatile Organic Compounds Combined with Atmospheric Pressure Chemical Ionization Quadrupole Mass Spectrometry

Физзи-экстракция летучих органических соединений в сочетании с химическим ионизацией при атмосферном давлении. Квадрупольная масс-спектрометрия

Chemical analysis of volatile and semivolatile compounds dissolved in liquid samples can be challenging. The dissolved components need to be brought to ...

Continuous Flow Chemistry: Reaction of Diphenyldiazomethane with p-Nitrobenzoic Acid

Continuous Flow Chemistry: Reaction of Diphenyldiazomethane with p-Nitrobenzoic Acid

Непрерывный поток химия: Реакция Diphenyldiazomethane с p- Nitrobenzoic кислота

Continuous flow technology has been identified as instrumental for its environmental and economic advantages leveraging superior mixing, heat transfer and ...

Chemical Precipitation Method for the Synthesis of Nb2O5 Modified Bulk Nickel Catalysts with High Specific Surface Area

Chemical Precipitation Method for the Synthesis of Nb2O5 Modified Bulk Nickel Catalysts with High Specific Surface Area

Метод химического осаждения для синтеза Nb2O5 изменение массовых никель катализаторов с высокой удельной площади поверхности

We demonstrate a method for the synthesis of NixNb1-xO catalysts with sponge-like and fold-like nanostructures. By varying the Nb:Ni ratio, a series of ...

A Two-Step Pyrolysis-Gas Chromatography Method with Mass Spectrometric Detection for Identification of Tattoo Ink Ingredients and Counterfeit Products

Двухступенчатый пиролиза-газ хроматографии метод с масс спектрометрическое обнаружение для идентификации татуировки чернила ингредиенты и контрафактной продукции

Tattoo inks are complex mixtures of ingredients. Each of them possesses different chemical properties which have to be addressed upon chemical analysis. ...