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

מאי מכירה במימון של קרן קרל-צייס. כל הנתונים נוצרו באמצעות 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 ופועלת בצורה הטובה ביותר בסביבת קונדה. הוראות להורדת פלטפורמת קונדה ניתן למצוא בכתובת 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 prompt" ולחץ על קיצור הדרך שנוצר כדי לפעול. במערכות MacOS ו-Linux, יש לקיים אינטראקציה עם Conda דרך המסוף מבלי להפעיל תוכניות נוספות. לאחר מכן, הקלד את הפקודה הבאה והקש Enter כדי לפעול, וענה בחיוב על כל שאלה שעולה במהלך ההתקנה: conda install -c rdkit rdkit. בעוד שישנן תוכניות רבות לחישוב מתארים הזמינות באופן חופשי, דוגמה זו משתמשת ב- PaDEL-Descriptor19, מחשבון חופשי ומהיר לתיאורים מולקולריים וטביעות אצבע.</li>
		</ol>
</li>
		<li>הורד ושמור את קובץ .jar בתיקיית הפרוייקט. הערה: ניתן להוריד את PaDEL-Descriptor בחינם 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 בדגל -מטושטש , וסגור את כל מרווחי היסוד בסוגריים (לדוגמה, השתמש ב-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 [שם קובץ מחברת]
</li>
	<li>בתא המיועד בתחילת המחברת, הזן את נתיב הקובץ המלא של קובץ .sdf הקלט (שנוצר על-ידי MAYGEN), את נתיב הקובץ המלא של קובץ הפלט הרצוי .sdf, ואת נתיב הקובץ של קובץ "badlist" כמחרוזות (בתוך מרכאות). ראה קובץ משלים 2 לדוגמה של רשימה גרועה.<ol>
<li>אם יש לשמור כמה מבני משנה בספריה המסוננת (רשימה טובה), צור קובץ .txt של תבניות SMARTS20 עבור אותם מבני משנה (רשימה טובה) והצב את נתיב הקובץ של הרשימה הטובה בשורה המיועדת בתחילת המחברת. ראה קובץ משלים 3 לדוגמה של רשימה טובה.</li>
	</ol>
</li>
	<li>הפעל מחדש את ליבת המחברת והפעל את כל התאים (מהתפריט בחלק העליון, בחר ליבה, הפעל מחדש והפעל הכל) כדי לקבל קובץ .sdf עם השם הרצוי בתיקיית הפלט שצוינה.</li>
	<li>חזור על שני השלבים הקודמים עבור כל קובץ מבנה שנוצר על-ידי MAYGEN בשלב 2.</li>
</ol>4. (אופציונלי) שינויים נוספים במבנה
הערה: אלה מבוצעים בדוגמה זו, אך ייתכן שאין בהם צורך באוצרות ספריות אחרות.
<ol><li>החלפת פסאודואטום. הערה: כאן, פסאודו-אטום הוא אטום ייחודי המשמש לייצוג תת-מבנה גדול יותר המשותף לכל המבנים הנוצרים, ובכך מקצר את זמן הדור של MAYGEN. ראה קובץ משלים 4 לדוגמה של החלפת פסאודו-אטום.<ol>
<li>פתח הנחיית Anaconda (ראה שלב 1.2.2.1) ונווט אל התיקיה המכילה את מחברות Jupyter.</li>
		<li>פתח את מחברת Jupyter להחלפה פסאודו-אטומית: מחברת jupyter [שם קובץ מחברת]
</li>
		<li>בתא המיועד בתחילת המחברת, הזן את נתיב הקובץ המלא של קובץ .sdf הקלט ואת נתיב הקובץ המלא של קובץ הפלט הרצוי .sdf כמחרוזות (בתוך מרכאות).</li>
		<li>הפעל מחדש את ליבת המחברת והפעל את כל התאים כדי לקבל קובץ .sdf עם השם הרצוי בתיקיית הפלט שצוינה.</li>
	</ol>
</li>
	<li>חומצת אמינו N- ו-C-termini מכסה הערה: הליך זה הוא ספציפי לחומצות אלפא-אמינו, ומוסיף מכסים מולקולריים ל-N ו-C-termini של עמודי השדרה של חומצות אלפא-אמינו. ראה קובץ משלים 5 לדוגמה של הגבלת חומצות אמינו.<ol>
<li>פתח הנחיית Anaconda (ראה שלב 1.2.2.1) ונווט אל התיקיה המכילה את מחברות Jupyter.</li>
		<li>פתח את מחברת Jupyter למכסה חומצות אמינו: מחברת jupyter [שם קובץ מחברת]
</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 [directory of the .sdf files] -file [נתיב קובץ של קובץ .csv לתוצאות] -2d -keeporder -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 and Met) כתוספות מאוחרות יחסית לקוד הגנטי, מסקנות הנתמכות על ידי היעדר חומצות אמינו המכילות גופרית במטאוריטים ובניסויים בצינורות ניצוץ. עם זאת, תרכובות אורגנוסולפור מזוהות בקלות בשביטים ובמטאוריטים22, ואנליזה מחודשת של ניסויים בצינורות ניצוץ באמצעות גז H2S מצאה חומצות אמינו ותרכובות אורגניות אחרות המכילות גופרית16. כאשר שוקלים אלפבית של חומצות אמינו חלופיות, כדאי לחקור אחת מועשרת בגופרית.
בפרוטוקול הנ"ל, יצירת מבנים וסינון תת-מבנה נחשבים לשלבים קריטיים; בהתאם להרכב של ספריית המבנה המוגמרת, ייתכן שחוקר יצטרך לבצע רק את שני השלבים הללו. הוראות ותוכנה לפעולות נוספות (החלפת פסאודו-אטום והוספת תת-מבנים (במקרה זה, כיפה של חומצות אמינו)) נכללות לחישוב מתאר רלוונטי יותר (ההצמדה מבטיחה שחישובי XLogP יושפעו מהצד ולא מקבוצות האמין או הקרבוקסיל של עמוד השדרה) ויצירת מבנה מהיר יותר באמצעות שימוש בפסאודואטום, הנדון ביתר פירוט בהמשך. בנוסף, חישוב תיאור נעשה כאן כדרך קלה להמחיש את מגוון המבנים שנוצרו ולהשוות את ההשפעות של העשרת גופרית בספריות המוגמרות.
בעוד ש-PaDEL-Descriptor יכול לחשב אלפי תכונות מולקולריות, נפח מולקולרי (כפי שחושב נפח ואן דר ואלס) ומקדם חלוקה (כ-XLogP) שימשו כאן משתי סיבות שונות. ראשית, שני מתארים אלה מודדים תכונות מולקולריות (גודל והידרופוביות, בהתאמה) המוכרות לרוב הכימאים והביולוגים. שנית, במקרה של חומצות אמינו, שתי תכונות אלה הן משמעותיות. במשך עשרות שנים, גודל חומצות אמינו והידרופוביות היו ידועים כמשפיעים על התרמודינמיקה של קיפול חלבונים23. שתי תכונות אלה מסייעות להסביר את תדירות ההחלפה של חומצות אמינו שהיו חלק בלתי נפרד מהבנת האבולוציה של חלבונים24.
הדוגמה לעיל מראה כי בשני התיאורים שנחקרו (נפח מולקולרי והידרופוביות), החלפת גופרית דיוולנטית בפחמן ושני מימן אינה מניבה שינויים משמעותיים. את העלייה הקלה והלא משמעותית בנפח המולקולרי הממוצע כתוצאה מהחלפת גופרית (איור 3) ניתן לייחס לרדיוס הקוולנטי הגדול יותר של הגופרית (כ-22:00) בהשוואה לפחמן sp3 (~75 pm) או sp2 (~73 pm)פחמן 25. באופן דומה, לתחליפי גופרית יש השפעה מינימלית על ה-XLogP הממוצע (איור 4). ההשפעה הגדולה ביותר הייתה בין ספריות ה-VAIL וה-VAIL_S, ככל הנראה בשל שילוב של ספריית ה-VAIL כהידרופובית במיוחד (ה-sidechains הם רק פחמימנים) וקבוצות סולפהידריל חומציות הרבה יותר מקבוצות המתיל שהן היו מחליפות. ההשפעה המינימלית של החלפת גופרית ניכרת באיור 2, שבו ספריות עם תחליפי גופרית תופסות את אותו מרחב כימי כמו ספריות מקבילות ללא החלפת גופרית.
הירידה במספר המבנים (איור 5A) ובזמן הדרוש ליצירת המבנים האלה (איור 5B) בעת שימוש בפסאודואטום אינה מפתיעה. שימוש בפסאודואטום מקטין את מספר האטומים הכבדים שיש לשלב בגרף כימי, מקטין את מספר צמתי הגרפים ומניב ירידות אקספוננציאליות בזמן הייצור ובמספר המבנים. כאן, הבחירה בזרחן טריוולנטי כפסאודו-אטום נובעת מביוכימיה בסיסית (בהיעדר תוספת פוסט-טרנסלציונית של קבוצות פוספטים, אין חומצות אמינו מקודדות גנטית המכילות זרחן) וערכיות האטום שתחליף אותו (זרחן טריוולנטי יכול בקלות להיות מוחלף בפחמן טטרוולנטי שנקשר בנפרד לאטום אחר או לקבוצת אטומים אחרים). בעוד שהקוד המסופק להחלפת פסאודו-אטומים הוא ספציפי להחלפת זרחן טריוולנטי בתת-מבנה אלנין, משתמשים יכולים להתאים אישית את הקוד כך שיעבוד עם פסאודו-אטומים שונים או תת-מבנים חלופיים שונים, וייתכן שהם משתמשים במספר פסאודו-אטומים במהלך יצירת המבנה הראשוני ולאחר מכן החלפת כל פסאודו-אטום בתת-מבנה מולקולרי גדול יותר.
שיטות ליצירת מבנה הדומות לאלה שבהן משתמשת MAYGEN (ושיטות אחרות כגון רשתות עצביות) כבר משמשות לגילוי תרופות כדי ליצור ספריות מורכבות לסינון סיליקו ; סקירה אחרונה4 דנה בשיטות אלה ביתר פירוט. מכיוון ששיטות אלה מיועדות בעיקר ליצירת מולקולות דמויות תרופה, ישנן מגבלות מסוימות על יכולתן לייצר מולקולות, כגון שימוש בתכונות ביולוגיות או פרמצבטיות כדי להגביל את המבנים שנוצרו (QSPR/QSAR הפוך) או יצירת מבנים ממספר מוגדר מראש של אבני בניין תת-מבניות. ככל שהאסטרוביולוגיה מתמקדת יותר בריבוי התרכובות האורגניות שיכולות להיווצר באופן אביוטי ופחות במוצרים סופיים כלשהם או בתכונותיהם, יצירת המבנה הממצה של MAYGEN אידיאלית ליצירת ספריות מבנים כדי לענות על שאלות אסטרוביולוגיות. הגישה לסינון תת-מבנה המתוארת כאן (המבוצעת לאחר יצירת מבנה באמצעות תוכנית חיצונית) שונה מהתוכנית המתחרה MOLGEN בכך שסינון תת-המבנה של MOLGEN מתרחש במהלך יצירת המבנה. מכיוון ש- MAYGEN הוא קוד פתוח, לא רק שהוא נגיש יותר מ- MOLGEN בשל עלות הרישוי של MOLGEN, אלא שאנשים פרטיים יכולים ליישם תכונות חדשות כגון סינון תת-מבנה במהלך יצירת מבנה.
כפי שנכתב, הפרוטוקול המתואר כאן מתמקד ביצירה ואוצרות של ספריות של חומצות אלפא-אמינו קטנות יחסית. כדי ליצור ספריות שונות, משתמשים יכולים לתת נוסחאות מולקולריות שונות ל- MAYGEN, לשנות את סינון תת-המבנה על-ידי שינוי גודל הטבעת המרבי המותר ואת ערכיות הקשר, או לערוך את קבצי הרשימה הטובה והרשימה הרעה כדי להוסיף או להסיר תבניות תת-מבנה. שינויים בפרוטוקולים הכוללים שינוי אופן ההוספה או ההחלפה של אטומים ותתי-מבנים (החלפה פסאודו-אטומית וכיפה מולקולרית) הם אפשריים, אך ידרשו תשומת לב רבה יותר למגבלות ערכיות כדי למנוע שגיאות RDKit לגבי ערכיות שגויה במבנים שהשתנו.
הפרוטוקול המפורט לעיל מיועד לחומצות אלפא-אמינו קטנות. עם זאת, הפורמט הכללי (יצירת מבנה מקיף באמצעות פסאודו-אטומים, ואחריו סינון תת-מבנה ושינויים מולקולריים) גמיש מאוד עבור תרכובות מעבר לחומצות אמינו קטנות. אפילו באסטרוביולוגיה, הליך דומה שבוצע לאחרונה באמצעות MOLGEN שימש לחקר איזומרים חוקתיים של חומצות גרעין26. בנוסף לכלים שתוארו לעיל, ניתן לשלב את MAYGEN עם כלים אחרים של כימינופורמטיקה בקוד פתוח כדי להפוך את היצירה והניתוח של מבנים כימיים חדשניים למשתלמים ונגישים למגוון רחב של תחומי מחקר.

יצירת מבנה מולקולרי משמשת כיישום מעשי של הבעיה הכללית של יצירת גרפים ממצה; בהינתן מספר צמתים (אטומים) ואילוצים על הקישוריות שלהם (למשל, ערכיות, ריבוי קשרים, תת-מבנים רצויים/לא רצויים), כמה גרפים מחוברים (מולקולות) אפשריים? מחוללי מבנים ראו יישום נרחב בגילוי תרופות ובפיתוח תרופות, שם הם יכולים ליצור ספריות עצומות של מבנים חדשים לסינון סיליקו 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, Thr, Leu, Ile, Pro) נוצרות בקלות בתנאים אביוטיים וככל הנראה מהוות את האלפבית של חומצות האמינו של אורגניזמים טרום-LUCA. עם הזמן, האלפבית "המוקדם" הזה הורחב בתגובה לצרכים מבניים ותפקודיים שונים. לדוגמה, סקירה עדכנית של Moosmann13 טוענת כי תוספת של חברים חדשים יותר של חומצות האמינו המקודדות גנטית (כלומר Met, Tyr ו- Trp) אפשרה הישרדות בסביבות עשירות בחמצן על ידי מניעת התפשטות תוך-תאית של מיני חמצן תגובתי.
חבילה הולכת וגדלה של טכניקות כימיה אנליטית מאפשרת תובנה לגבי מבני חומצות האמינו שיכולים להיווצר בתנאים אביוטיים. סקירה אחרונה14 של סימקוס ואחרים מפרטת את השיטות ששימשו לאיתור תרכובות אורגניות רבות במטאוריטים, כמו גם תרכובות אורגניות מהדמיות במבחנה של סביבות כדור הארץ הקדומות 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">זמן דור (ms)</td>
			<td colspan="2">מבנים</td>
		</tr>
<tr>
<td>הראשונית</td>
			<td>הסופי</td>
		</tr>
<tr>
<td>1</td>
			<td>גלי</td>
			<td colspan="2">ג2ח5לא2
</td>
			<td colspan="2">כולל תת-מבנה גלי</td>
			<td>גלי</td>
			<td>192</td>
			<td>84</td>
			<td>1</td>
		</tr>
<tr>
<td>2</td>
			<td>וייל</td>
			<td colspan="2">מחשב0-3H3-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-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>מקצוען</td>
			<td colspan="2">ג2-5לא2ח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">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">ג2-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 עבור כל נוסחה ("ראשונית") ואת מספר המולקולות שנותרו לאחר סינון אלה עם תת-מבנים לא רצויים ("סופי"). קיצורים: VAIL = ואלין, אלנין, איזולאוצין, לאוצין; DEST = חומצה אספרטית, חומצה גלוטמית, סרין, תראונין; X_S = גופרית דיוולנטית מחליפה את אחד הפחמנים בספריה X; N-meX = N-מתילX.
השיטות הכלליות לעיל יושמו על נוסחאות המבוססות על חומצות האמינו המקודדות "המוקדמות", בעקבות ההליך של Meringer et al.21 מבנים Badlist נלקחו מאותו מקור והומרו למחרוזות SMARTS כדי לייצג בקלות תבניות תת-מבניות. שני תתי-מבנים של רשימות רעות לא שימשו בדוגמה זו: מבנה 018 (CH 3-CH-N) התאים כמעט לאיזומרים של פרולין שלא היו בעצמם לא יציבים; מבנה 106 (R-C-C-OH, שבו R=תת-מבנה אלנין מתחבר בבטא-פחמן) התאים לחומצה גלוטמית, חומצת אמינו מקודדת. בנוסף לנוסחאות הכימיות הללו, נוצרו גרסאות עם גופרית דיוואלנטית התופסת את מקומו של אטום פחמן ושני אטומי מימן. מטעמי ביצועים, כמה מהנוסחאות הללו משתמשות באטום זרחן טריוולנטי (למשל, "פסאודו-אטום") כתחליף לבטא-פחמן של תת-מבנה אלנין. טבלה 1 מפרטת את הספריות שנוצרו בדוגמה זו, את הנוסחאות ששימשו ליצירתן ואת מספר התרכובות הכלולות בהן. שמות הספריות מבוססים על חומצות האמינו המקודדות שמהן הם נגזרים: או באמצעות קיצור בן 3 אותיות (Gly = glycine, Pro = proline) או קיצור של אותיות בודדות (VAIL = Valine, Alanine, איזולאוצין, לאוצין; DEST = חומצה אספרטית, חומצה גלוטמית, סרין, תראונין). הסיומת "_S" מצביעה על כך שגופרית הוחלפה בפחמן בנוסחה המקורית של הספרייה (למשל, VAIL_S בנוי עם אותה נוסחה מטושטשת כמו VAIL, אך עם גופרית דיאלוונטית שמחליפה את אחד הפחמנים).
לאחר יצירת מבנים עם MAYGEN, הספריות שהתקבלו סוננו מתרכובות המכילות לפחות תת-מבנה אחד הכלול ברשימה הרעה. בעקבות סינון זה, כל אטומי הזרחן הוחלפו בתת-מבנה אלנין. לאחר מכן, נוצרו גרסאות "מכוסות" של כל המבנים, עם קבוצת אצטיל שנוספו ל-N-terminus וקבוצת N-מתיל אמיד שנוספה ל-C-terminus. זה נעשה כדי להסיר את ההשפעה על ההידרופוביות של קבוצות האמין החופשי והחומצה הקרבוקסילית בעמוד השדרה של חומצת האלפא-אמינו. PaDEL-Descriptor שימש לחישוב XLogP עבור כל המבנים המכוסים וחישב את נפח ואן דר ואלס (VABC) עבור כל המבנים הלא מכוסים.
איור 2 מציג את המרחב הכימי של הספריות המסוננות, כפי שהוגדר על ידי מתארי VABC ו-XLogP. כאן, טווח ערכי logP אפשריים גדל עם נפח מולקולרי, אפילו בתוך ספריות שאין בהן צדדים הידרופיליים במפורש (למשל, VAIL, Pro). חומצות אמינו מקודדות עם ציד פחמימנים היו הידרופוביות יותר מרוב חומצות האמינו האחרות בנפח דומה מהספרייה שלהן. נראה שזה המקרה גם עבור Met ו- Cys בהשוואה לחברים אחרים בספריית VAIL_S עם כרכים דומים. חומצות אמינו מקודדות עם שרשראות צד הידרוקסיל (Ser ו-Thr) היו בין החברים הקטנים ביותר בספריית DEST, כאשר Asp היה גדול רק במעט מ-Thr.
איור 3 ואיור 4 מראים את ההשפעות על נפח ולוגפ כאשר גופרית דיואלנט מחליפה פחמן בשרשרת צד של חומצת אלפא-אמינו. תחליפי גופרית הובילו לעלייה קלה בנפח המולקולרי בכל הספריות (איור 3). ההשפעה של החלפת גופרית על logP אינה הומוגנית כמו נפח (איור 4). הלוג הממוצע של ספריית VAIL_S נמוך מעט מזה של ספריית VAIL, אך אפקט זה אינו נראה באף אחד מזוגות הספריות האחרים (DEST ו- DEST_S, Pro ו- Pro_S).
איור 5 מכמת את ההשפעות על יצירת מבנה של פסאודו-אטום העומד בפני תת-מבנה משותף; כאן, P טריוולנטי החליף את המואטי של אלנין במהלך יצירת מבנה. שימוש בפסאודואטום ביצירת מבנים הפחית מאוד את מספר המבנים הנוצרים על ידי סדרי גודל של כ-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_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 GHz, 16 GB של זיכרון RAM, ללא מבני שמירה לקובץ, ואפשרות -m לשימוש בריבוי הליכי משנה. בדיקות באמצעות פסאודו-אטום השתמשו בנוסחאות המטושטשות כמתואר בטבלה 1. עבור בדיקות ללא פסאודואטום, הנוסחאות המטושטשות שבהן נעשה שימוש היו זהות לאלה המתוארות בטבלה 1 עם השינויים הבאים: P הוחלף ב-N; ספירת הפחמן הוגדלה ב-3; ספירת המימן הוגדלה ב-7; ספירת החמצן עלתה ב-2. פסים שחורים מציגים ספריות שנוצרו עם פסאודו-אטום; פסים אפורים מציגים ספריות שנוצרו ללא פסאודו-אטום. (A) מספר המבנים שנוצרו באמצעות הנוסחאות המטושטשות ששימשו לבניית ספריות VAIL ו-DEST עם ובלי החלפת זרחן טריוולנטית עבור תת-מבנה אלנין. (B) זמן (ב-ms) הדרוש לבניית ספריות 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, ...

הפרוטוקול שלנו מדגים כיצד תוכנת קוד פתוח יכולה לאפשר לכל חוקר ליצור ולאצור ספריית מבנה חישובית. ערעור הפרוטוקולים הזה נובע מהפתיחות והגמישות שלו. כל אחד יכול להשתמש בו ולשנות אותו כך שיתאים לשאלת המחקר הספציפית שלו.ניתן ליישם גרסאות של פרוטוקול זה על יישומי גילוי תרופות, תוך יצירת ספריות מבנה ספציפיות במהירות לסינון סיליקו. למרות שהפרוטוקול מוסבר צעד אחר צעד, אם משתמשים אינם מכירים את Java או קידוד בסיסי הם יכולים תחילה להסתכל על אלה לפני יישום הפרוטוקול. התחל על ידי יצירת ספרייה חדשה עבור הפרויקט.מקם את כל הקבצים והתקני ההפעלה בספריה זו לגישה נוחה. הורד את הגרסה האחרונה של Maygen כקובץ צנצנת ואת תוכנת ניהול החבילות Anaconda. במערכות Windows חפש את הנחיית Anaconda ולחץ על קיצור הדרך שנוצר כדי לפעול.כדי ליצור סביבת RDKit באנקונדה ולהוריד את ה- RDKit לסביבה, הקלד את הפקודה המוצגת על המסך, לחץ על Enter כדי לפעול וענה כן על כל שאלה שעולה במהלך ההתקנה. לאחר מכן הורד את מחברות Jupyter וקבצי הטקסט של תבניות המצע מהקבצים המשלימים, אחד עד חמישה. בשורת הפקודה, נווט אל הספריה המכילה את המייגן.קובץ הפעלה של צנצנת. עבור כל נוסחה כימית מעניינת השתמש בפקודה המוצגת על המסך כדי להפעיל את Maygen. אם הנוסחה היא נוסחה מטושטשת במקום נוסחה בדידה, החלף את דגל המקף F בדגל מקף מטושטש והקף את כל מרווחי הרכיבים בסוגריים.בשורת אנקונדה נווט אל התיקיה המכילה את מחברות 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, פתח אותו וחלץ אותו לתיקיה זו.פתח שורת פקודה, נווט אל התיקיה המכילה את קובץ צנצנת המתאר PaDEL והפעל את מתאר PaDEL עבור קבצי sdf שנאספו. המרחב הכימי של כל ספריות חומצות האמינו המסוננות מוצג כאן. סמנים שחורים מייצגים חומצות אמינו מהספריות ללא גופרית וסמנים צהובים מייצגים חומצות אמינו מספריות מועשרות בגופרית.כאן, VAIL וספריות VAIL_S מיוצגות על ידי מעגלים. ספריות DEST ו- DEST_S מיוצגות על ידי ריבועים. ספריות Proline ו-Pro S מיוצגות על ידי משולשים וכוכבים מייצגים חומצות אמינו מקודדות.טווח ערכי הלוג P האפשריים גדל עם הנפח המולקולרי אפילו בתוך הספריות שחסרות במפורש שרשראות צד הידרופיליות. חומצות אמינו מקודדות עם שרשראות צד פחמימנים הן הידרופוביות יותר מרוב חומצות האמינו האחרות בנפח דומה מהספרייה המתאימה להן. זה גם המקרה של מתיונין מתעקש בהשוואה לחברים אחרים בספריית VAILS עם כרכים דומים.חומצות אמינו מקודדות עם שרשראות צד הידרוקסיל היו בין החברים הקטנים ביותר בספריית DEST עם חומצה אספרטית רק מעט גדולה יותר משלוש אנין. התמונה המיוצגת מציגה את הנפחים הממוצעים של ואן דר ואל של ספריות עם גופרית וללא גופרית. תחליפי גופרית הובילו לעלייה קלה בנפח המולקולרי בכל הספריות.ערכי מקדם החלוקה הממוצעים של ספריות עם ובלי גופרית מוצגים כאן. ההשפעה של החלפת גופרית על לוג 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

זיהוי איכותי של חומצות קרבוקסיליות, חומצות boronic, ואמיני שימוש 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

ההפרדה סולפט ידי סלקטיבי התגבשות עם ליגנד Bis-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

קביעת ההרכב הכימי של ממשקי קורוזיה מונעים / מתכת עם 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; ה</em&gt;) - Dienamides על ידי Olefination של אלקנים אלקטרון לקוי עם Allyl אצטט

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

פיצי הפקה של תרכובות אורגניות נדיפות בשילוב עם לחץ אטמוספרי כימי יוניזציה Quadropole מסה ספקטרומטריית

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 עם חומצה - Nitrobenzoic p

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. ...