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 erkänner finansiering från Carl-Zeiss-Foundation. Alla siffror genererades med hjälp av Microsoft Excel.

OBS: Se figur 1 för en sammanfattning av protokollet och materialförteckningen för detaljer om den programvara som används.
<img alt="Figure 1" class="xfigimg" src="/files/ftp_upload/63632/63632fig01.jpg"> Bild 1: Sammanfattning av protokollet. <a href="https://www.jove.com/files/ftp_upload/63632/63632fig01large.jpg" target="_blank"></a>
1. Nedladdning av programvara och filer
OBS: Alla program är gratis för individuell användning och kan köras på en persondator.
<ol><li>Skapa en ny katalog för det här projektet. Placera filer och körbara filer här för enkel åtkomst.</li>
	<li>Ladda ner och installera nödvändiga programvarupaket.<ol>
<li>Ladda ner den senaste versionen av MAYGEN som en .jar fil. OBS: MAYGEN är fritt tillgänglig som en .jar fil från https://github.com/MehmetAzizYirik/MAYGEN/releases</li>
		<li>Ladda ner och installera pakethanteringsprogrammet Conda och cheminformatics toolkit RDKit18. OBS: RDKit filtrerar de molekylära strukturerna som produceras av MAYGEN och fungerar bäst i en Conda-miljö. Instruktioner för nedladdning av Conda-plattformen finns på https://conda.io/projects/conda/en/latest/user-guide/install/index.html. Installations- och miljöinstallationsinstruktioner för RDKit finns på https://www.rdkit.org/docs/Install.html.<ol>
<li>Installera RDKit i Conda-huvudmiljön istället för en separat RDKit-miljö via Anaconda-prompten. På Windows-system söker du efter "Anaconda prompt" och klickar på den resulterande genvägen för att köra. På MacOS- och Linux-system, interagera med Conda via terminalen utan att köra några ytterligare program. Skriv sedan följande kommando och tryck på Enter för att köra och svara ja på alla frågor som kommer upp under installationen: conda installera -c rdkit rdkit. Även om det finns många fritt tillgängliga deskriptorberäkningsprogram, använder detta exempel PaDEL-Descriptor19, en gratis och snabb kalkylator för molekylära deskriptorer och fingeravtryck.</li>
		</ol>
</li>
		<li>Ladda ned och spara den .jar filen i projektmappen. PADEL-Descriptor kan laddas ner gratis från http://www.yapcwsoft.com/dd/padeldescriptor/.</li>
	</ol>
</li>
	<li>Ladda ned Jupyter Notebooks och textfiler med understrukturmönster från Tilläggsfiler 1–5. Jupyter Notebooks kan också laddas ned från följande GitHub sida: https://github.com/cmayerb1/AA-structure-manip.</li>
</ol>2. Strukturgenerering med MAYGEN
<ol><li>I en kommandotolk navigerar du till katalogen som innehåller MAYGEN .jar körbar fil.</li>
	<li>För varje kemisk formel av intresse kör du MAYGEN med följande kommando: java -jar [MAYGEN .jar filnamn] -f [kemisk formel] -v -o [mapp för MAYGEN-utdata] -m -sdf. OBS: Detta sparar en .sdf-fil i den angivna mappen, uppkallad efter den använda formeln.<ol>
<li>Om formeln är en luddig formel istället för en diskret formel, ersätt flaggan -f med en -fuzzy flagga och bifoga eventuella elementintervall inom parentes (t.ex. använd C[5-7]H[12-15] för att säkerställa att alla genererade strukturer har mellan 5 och 7 kolatomer och mellan 12 och 15 väteatomer).</li>
	</ol>
</li>
</ol>3. Filtrera föreningar med oönskade understrukturer
<ol><li>Öppna en Anaconda-prompt (se steg 1.2.2.1) och navigera till mappen som innehåller Jupyter Notebooks som hämtats från tilläggsfil 1.</li>
	<li>Öppna Jupyter Notebook för filtrering av understrukturer med följande kommando: jupyter notebook [anteckningsbokens filnamn]
</li>
	<li>I den angivna cellen i början av anteckningsboken anger du den fullständiga filsökvägen för indatafilen .sdf (genererad av MAYGEN), fullständig filsökväg för önskad .sdf utdatafil och filsökväg för "badlist" -filen som strängar (inom citat). Se Tilläggsfil 2 för ett exempel på en dålig lista.<ol>
<li>Om vissa understrukturer i det filtrerade biblioteket (en goodlist) ska behållas skapar du en .txt-fil med SMARTS-mönster20 för dessa understrukturer (en goodlist) och placerar goodlist-filsökvägen på den angivna raden i början av anteckningsboken. Se Tilläggsfil 3 för ett exempel på en bra lista.</li>
	</ol>
</li>
	<li>Starta om notebook-kärnan och kör alla celler (från menyn högst upp, välj Kernel, Restart &Run All) för att få en .sdf fil med önskat namn i den angivna utdatamappen.</li>
	<li>Upprepa de två föregående stegen för varje strukturfil som genereras av MAYGEN i steg 2.</li>
</ol>4. (Valfritt) Ytterligare strukturändringar
Dessa utförs i det här exemplet men kanske inte behövs för att kurera andra bibliotek.
<ol><li>Pseudoatom ersättning. OBS: Här är en pseudoatom en unik atom som används för att representera en större understruktur som delas av alla genererade strukturer, vilket minskar MAYGEN: s generationstid. Se Kompletterande fil 4 för ett exempel på pseudoatombyte.<ol>
<li>Öppna en Anaconda-prompt (se steg 1.2.2.1) och navigera till mappen som innehåller Jupyter Notebooks.</li>
		<li>Öppna Jupyter Notebook för pseudoatombyte: jupyter notebook [anteckningsbokens filnamn]
</li>
		<li>I den angivna cellen i början av anteckningsboken anger du den fullständiga filsökvägen för indatafilen .sdf och den fullständiga filsökvägen för önskad .sdf utdatafilen som strängar (inom citat).</li>
		<li>Starta om notebook-kärnan och kör alla celler för att hämta en .sdf-fil med önskat namn i den angivna utdatamappen.</li>
	</ol>
</li>
	<li>Aminosyra N- och C-termini capping OBS: Denna procedur är specifik för alfa-aminosyror och lägger till molekylära kepsar till N- och C-termini av alfa-aminosyra ryggrad. Se Kompletterande fil 5 för ett exempel på aminosyrakapsling.<ol>
<li>Öppna en Anaconda-prompt (se steg 1.2.2.1) och navigera till mappen som innehåller Jupyter Notebooks.</li>
		<li>Öppna Jupyter Notebook för aminosyrakapsling: jupyter notebook [anteckningsbokens filnamn]
</li>
		<li>I den angivna cellen i början av anteckningsboken anger du den fullständiga filsökvägen för indatafilen .sdf och den fullständiga filsökvägen för önskad .sdf utdatafilen som strängar (inom citat).</li>
		<li>Starta om notebook-kärnan och kör alla celler för att hämta en .sdf-fil med önskat namn i den angivna utdatamappen.</li>
	</ol>
</li>
</ol>5. Generering av deskriptorer
<ol><li>Innan du genererar deskriptorer placerar du alla .sdf filer som deskriptorerna ska beräknas för i en enda mapp. OM det inte redan är gjort, ge dessa filer beskrivande namn för enkel filtrering efter deskriptorgenerering.</li>
	<li>Öppna en kommandotolk och navigera till mappen som innehåller Filen PaDEL-Descriptor .jar.</li>
	<li>Kör PaDEL-Descriptor för de insamlade .sdf filerna med följande kommando: java -jar PaDEL-Descriptor.jar -dir [katalog över .sdf filer] -fil [filsökväg för en .csv fil för resultat] -2d -retainorder -usefilenameasmolname Resultatfilen kommer att ha molekylnamnet i den första kolumnen och varje deskriptor i de efterföljande kolumnerna.</li>
	<li>Exportera dessa data till valfritt kalkylprogram för vidare analys.</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>

Författarna har inga intressekonflikter att avslöja.

En egenskap hos de "tidiga" aminosyrorna är brist på svavel. De metaanalyser som nämnts tidigare anser i allmänhet att de svavelhaltiga kodade aminosyrorna (Cys och Met) har varit relativt sena tillägg till den genetiska koden, slutsatser som stöds av bristen på svavelhaltiga aminosyror i meteoriter och gniströrsexperiment. Organosulfurföreningar detekteras emellertid lätt i kometer och meteoriter22, och omanalys av gniströrsexperiment med användning avH2S-gashittade aminosyror och andra organiska föreningar innehållande svavel16. När man överväger ett alternativt aminosyraalfabet är ett berikat med svavel värt att utforska.
I ovanstående protokoll betraktas strukturgenerering och understrukturfiltrering som kritiska steg; beroende på sammansättningen av det färdiga strukturbiblioteket kan en forskare bara behöva utföra dessa två steg. Instruktioner och programvara för ytterligare åtgärder (pseudoatombyte och tillägg av understrukturer (i detta fall aminosyrakapsling)) ingår för mer relevant deskriptorberäkning (capping säkerställer att XLogP-beräkningar påverkas av sidokedjan och inte ryggradsamin- eller karboxylgrupperna) och snabbare strukturgenerering via användning av en pseudoatom, som diskuteras mer detaljerat nedan. Dessutom görs deskriptorberäkning här som ett enkelt sätt att visualisera mångfalden av de genererade strukturerna och jämföra effekterna av svavelanrikning i de färdiga biblioteken.
Medan PaDEL-Descriptor kan beräkna tusentals molekylära egenskaper, användes molekylvolym (som beräknad van der Waals-volym) och fördelningskoefficient (som XLogP) här av två olika skäl. För det första mäter dessa två deskriptorer molekylära egenskaper (storlek respektive hydrofobicitet) som är bekanta för de flesta kemister och biologer. För det andra, när det gäller aminosyror, är dessa två egenskaper signifikanta. I årtionden var aminosyrastorlek och hydrofobicitet kända för att påverka termodynamiken för proteinveckning23. Dessa två egenskaper hjälper till att förklara aminosyrasubstitutionsfrekvenser som har varit integrerade för att förstå proteinutveckling24.
Ovanstående exempel visar att i de två studerade deskriptorerna (molekylvolym och hydrofobicitet) ger ersättning av ett tvåvärt svavel för ett kol och två väten inte signifikanta förändringar. Den lilla, icke-signifikanta ökningen av medelmolekylvolymen från svavelsubstitution (figur 3) kan hänföras till svavels större kovalenta radie (~ 103 pm) jämfört med antingen sp3 (~ 75 pm) eller sp2 (~ 73 pm) kol25. På samma sätt har svavelsubstitution minimal effekt på medelvärdet XLogP (figur 4). Den största effekten var mellan VAIL- och VAIL_S-biblioteken, troligen på grund av att en kombination av VAIL-biblioteket var särskilt hydrofobt (sidokedjorna är bara kolväten) och sulfhydrylgrupper är mycket surare än de metylgrupper de skulle ersätta. Den minimala effekten av svavelsubstitution framgår av figur 2, där bibliotek med svavelsubstitution upptar samma kemiska utrymme som analoga bibliotek utan svavelsubstitution.
Minskningen av antalet strukturer (figur 5A) och den tid som krävs för att generera dessa strukturer (figur 5B) vid användning av en pseudoatom är inte förvånande. Att använda en pseudoatom minskar antalet tunga atomer som måste införlivas i en kemisk graf, vilket minskar antalet grafnoder och ger exponentiella minskningar i generationstid och antal strukturer. Här härrör valet av trivalent fosfor som pseudoatom från grundläggande biokemi (frånvarande posttranslationell tillsats av fosfatgrupper, inga genetiskt kodade aminosyror innehåller fosfor) och valensen hos atomen som skulle ersätta den (en trivalent fosfor kan lätt ersättas med ett tetravalent kol som är enskilt bundet till en annan atom eller grupp av atomer). Medan den tillhandahållna koden för pseudoatomsubstitution är specifik för att ersätta en trivalent fosfor med en alaninunderstruktur, kan användare anpassa koden för att arbeta med olika pseudoatomer eller ersättningsunderstrukturer, potentiellt med hjälp av flera pseudoatomer under den första strukturgenereringen följt av att ersätta varje pseudoatom med en större molekylär understruktur.
Strukturgenereringsmetoder som liknar de som används av MAYGEN (och andra metoder som neurala nätverk) används redan vid läkemedelsupptäckt för att generera sammansatta bibliotek för in silico-screening ; en ny granskning4 diskuterar dessa metoder mer detaljerat. Eftersom dessa metoder främst är avsedda för att skapa läkemedelsliknande molekyler finns det vissa begränsningar för deras förmåga att generera molekyler, såsom att använda biologiska eller farmaceutiska egenskaper för att begränsa de skapade strukturerna (invers QSPR / QSAR) eller skapa strukturer från ett förinställt antal byggstenar för understruktur. Eftersom astrobiologi fokuserar mer på de många organiska föreningar som kan bildas abiotiskt och mindre på några slutprodukter eller deras egenskaper, är MAYGEN: s uttömmande strukturgenerering idealisk för att skapa strukturbibliotek för att ta itu med astrobiologiska frågor. Det tillvägagångssätt för understrukturfiltrering som beskrivs här (utfört efter strukturgenerering via ett externt program) skiljer sig från konkurrentprogrammet MOLGEN genom att MOLGEN:s understrukturfiltrering sker under strukturgenerering. Eftersom MAYGEN är öppen källkod är det inte bara mer tillgängligt än MOLGEN på grund av MOLGEN: s licenskostnad, utan individer kan implementera nya funktioner som understrukturfiltrering under strukturgenerering.
Som skrivet är protokollet som beskrivs här inriktat på att generera och kurera bibliotek med relativt små alfa-aminosyror. För att generera olika bibliotek kan användare ge olika molekylformler till MAYGEN, ändra understrukturfiltreringen genom att ändra den maximala tillåtna ringstorleken och bindningsvalensen, eller redigera goodlist- och badlist-filerna för att lägga till eller ta bort understrukturmönster. Protokolländringar som innebär att man ändrar hur atomer och understrukturer läggs till eller ersätts (pseudoatomsubstitution och molekylär capping) är genomförbara men kommer att kräva mer uppmärksamhet åt valensbegränsningar för att undvika RDKit-fel om felaktiga valenser i modifierade strukturer.
Protokollet som beskrivs ovan är utformat för små alfa-aminosyror. Det allmänna formatet (omfattande strukturgenerering med pseudoatomer, följt av understrukturfiltrering och molekylära modifieringar) är emellertid mycket flexibelt för föreningar bortom små aminosyror. Även inom astrobiologi användes ett liknande nyligen genomfört förfarande med MOLGEN för att undersöka konstitutionella isomerer av nukleinsyror26. Förutom de verktyg som beskrivs ovan kan MAYGEN paras ihop med andra kemininformatikverktyg med öppen källkod för att göra det möjligt att skapa och analysera nya kemiska strukturer överkomliga och tillgängliga för ett brett spektrum av forskningsområden.

Molekylär strukturgenerering fungerar som en praktisk tillämpning av det allmänna problemet med uttömmande grafgenerering; givet flera noder (atomer) och begränsningar för deras anslutning (t.ex. valenser, bindningsmultiplikationer, önskade /oönskade understrukturer), hur många anslutna grafer (molekyler) är möjliga? Strukturgeneratorer har sett omfattande tillämpning inom läkemedelsupptäckt och läkemedelsutveckling, där de kan skapa stora bibliotek med nya strukturer för in silico screening1.
Den första strukturgeneratorn, CONGEN, utvecklades för det första artificiella intelligensprojektet inom organisk kemi, DENDRAL2 (förkortning för DENDRitic ALgorithm). Flera programvaruföljare av DENDRAL rapporterades i litteraturen; men inte alla var underhållna eller effektiva. För närvarande är MOLGEN3 den senaste molekylära strukturgeneratorn. Tyvärr för de flesta potentiella användare är det sluten källa och kräver en licensavgift. Således har det funnits behov av en effektiv strukturgenerator med öppen källkod som enkelt kan anpassas till specifika applikationer. En utmaning för en effektiv strukturgenerator är att hantera kombinatorisk explosion; när storleken på en molekylformel ökar ökar storleken på det kemiska sökutrymmet exponentiellt. En nyligen genomförd granskning utforskar ytterligare historien och utmaningarna med molekylär struktur generation4.
Före 2021 var Parallel Molecule Generator (PMG)5 den snabbaste strukturgeneratorn med öppen källkod, men den var fortfarande långsammare än MOLGEN i storleksordningar. MAYGEN6 är ungefär 47 gånger snabbare än PMG och cirka 3 gånger långsammare än MOLGEN, vilket gör MAYGEN till den snabbaste och mest effektiva open source-strukturgeneratorn som finns tillgänglig. Mer detaljerade jämförelser och benchmarkingtester finns i dokumentet som introducerar MAYGEN6. Ett viktigt inslag i programmet är dess lexikografiska beställningsbaserade test för kanoniska strukturer, en ordnad grafgenereringsmetod baserad på Schreier-Sims7-algoritmen . Programvaran kan enkelt integreras i andra projekt och förbättras för användarnas behov.
Liksom MOLGEN och PMG tar MAYGEN en användardefinierad molekylformel och genererar alla möjliga strukturer för den formeln. Till exempel, om en användare kör MAYGEN med formelnC5H12, kommer MAYGEN att generera alla möjliga strukturer som innehåller fem kolatomer och tolv väteatomer. Till skillnad från sin motsvarighet med öppen källkod PMG kan MAYGEN också rymma "fuzzy" molekylformler som använder intervall istället för diskreta tal för antalet av varje element. Till exempel, om en användare kör MAYGEN med formeln C5-7H12-15, kommer MAYGEN att generera alla möjliga strukturer som innehåller mellan fem och sju kolatomer och tolv och femton väteatomer, vilket möjliggör enkel generering av strukturer med ett brett spektrum av atomkompositioner.
Astrobiologi är ett sådant område som kan dra nytta av molekylära strukturgeneratorer. Ett populärt ämne inom astrobiologi är utvecklingen av aminosyraalfabetet som delas av allt existerande liv på jorden. En av de definierande egenskaperna hos Den sista universella gemensamma förfadern (LUCA) är dess användning av tjugo genetiskt kodade aminosyror för proteinkonstruktion 8,9. Baserat på metaanalyser av arbete inom flera fält 10,11,12, bildas cirka 10 av dessa aminosyror (Gly, Ala, Val, Asp, Glu, Ser, Thr, Leu, Ile, Pro) lätt under abiotiska förhållanden och utgjorde sannolikt aminosyraalfabetet för pre-LUCA-organismer. Med tiden utvidgades detta "tidiga" alfabet som svar på olika strukturella och funktionella behov. Till exempel hävdar en ny granskning från Moosmann13 att tillsatsen av nyare medlemmar av de genetiskt kodade aminosyrorna (nämligen Met, Tyr och Trp) möjliggjorde överlevnad i syrerika miljöer genom att förhindra intracellulär spridning av reaktiva syrearter.
En ständigt växande svit av analytiska kemitekniker ger insikt i de aminosyrastrukturer som kan bildas under abiotiska förhållanden. En nyligen genomförd granskning14 av Simkus och andra beskriver de metoder som används för att detektera många organiska föreningar i meteoriter, liksom organiska föreningar från in vitro-simuleringar av tidiga jordmiljöer 15,16,17. Systematisk generering av kemiska strukturer gör det möjligt för forskare att utforska bortom de organiska föreningar som detekteras via instrumentering och befolka det strukturella utrymmet runt strukturella "öar" som identifierats av analytisk kemi. När det gäller de "tidiga" aminosyrorna visar denna systematiska strukturgenerering möjliga proteinkemier tillgängliga för tidigt liv utan att begränsa utforskningen till strukturer som experimentellt har detekterats under abiotiska syntesbetingelser. Med verktygssatser för kemininformatik med öppen källkod och effektiva strukturgeneratorer som MAYGEN är det nu enklare än någonsin att skapa och utforska nya kemiska strukturbibliotek och kan vägleda mer detaljerade undersökningar av livets alternativa kemier.

<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

Uttömmande generering av molekylära strukturer har många kemiska och biokemiska tillämpningar som läkemedelsdesign, molekylär databaskonstruktion, utforskning av alternativa biokemister och många fler. Matematiskt sett är dessa grafgeneratorer med kemiska begränsningar. På fältet är den mest effektiva generatorn för närvarande (MOLGEN) en kommersiell produkt som begränsar dess användning. Alternativt till det, en annan molekylär strukturgenerator, MAYGEN, är ett nytt verktyg med öppen källkod med effektivitet jämförbar med MOLGEN och kapaciteten för användare att öka dess prestanda genom att lägga till nya funktioner. Ett av de forskningsområden som kan dra nytta av denna utveckling är astrobiologi; strukturgeneratorer gör det möjligt för forskare att komplettera experimentella data med beräkningsmöjligheter för alternativ biokemi. Detta protokoll beskriver ett användningsfall för strukturgenerering inom astrobiologi, nämligen generering och kurering av alfa-aminosyrabibliotek. Med hjälp av strukturgeneratorer med öppen källkod och kemininformatikverktyg kan de metoder som beskrivs här implementeras utöver astrobiologi för lågkostnadsskapande och kurering av kemiska strukturbibliotek för alla forskningsfrågor.

<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">Bibliotek</td>
			<td colspan="2" rowspan="2">Formel</td>
			<td colspan="2" rowspan="2">Ytterligare begränsningar</td>
			<td rowspan="2">"Tidiga" kodade aminosyror</td>
			<td rowspan="2">Generationstid (ms)</td>
			<td colspan="2">Strukturer</td>
		</tr>
<tr>
<td>Initial</td>
			<td>Sist</td>
		</tr>
<tr>
<td>1</td>
			<td>Gly</td>
			<td colspan="2">C2H5NEJ2
</td>
			<td colspan="2">inkludera Gly-understruktur</td>
			<td>Gly</td>
			<td>192</td>
			<td>84</td>
			<td>1</td>
		</tr>
<tr>
<td>2</td>
			<td>VAIL</td>
			<td colspan="2">PC0-3H3-9
</td>
			<td colspan="2"></td>
			<td>Val, Ala, Ile, Leu</td>
			<td>172</td>
			<td>70</td>
			<td>22</td>
		</tr>
<tr>
<td>3</td>
			<td>DEST</td>
			<td colspan="2">PC0-3O1-2H3-5
</td>
			<td colspan="2"></td>
			<td>Asp, Glu, Ser, Thr</td>
			<td>481</td>
			<td>1928</td>
			<td>254</td>
		</tr>
<tr>
<td>4</td>
			<td>PRO</td>
			<td colspan="2">C2-5NEJ2H7-11
</td>
			<td colspan="2">Inkludera N-meGly- eller N-meAla-understruktur</td>
			<td>PRO</td>
			<td>4035</td>
			<td>79777</td>
			<td>16</td>
		</tr>
<tr>
<td>5</td>
			<td>VAIL_S</td>
			<td colspan="2">Psc0-2H3-7
</td>
			<td colspan="2"></td>
			<td></td>
			<td>122</td>
			<td>65</td>
			<td>31</td>
		</tr>
<tr>
<td>6</td>
			<td>DEST_S</td>
			<td colspan="2">Psc0-2O1-2H3
</td>
			<td colspan="2"></td>
			<td></td>
			<td>349</td>
			<td>1075</td>
			<td>79</td>
		</tr>
<tr>
<td>7</td>
			<td>Pro_S</td>
			<td colspan="2">C2-4SNO2H7-9
</td>
			<td colspan="2">Inkludera N-meGly- eller N-meAla-understruktur</td>
			<td></td>
			<td>3999</td>
			<td>75734</td>
			<td>10</td>
		</tr>
</tbody></table>Tabell 1: Sammansatta bibliotek som används i det här exemplet. Bibliotek byggda av formlerna 1-4 (Gly, VAIL, DEST och Pro) är baserade på tidigare publicerade fuzzy formler av de "tidiga" kodade aminosyrorna21, medan bibliotek byggda av formlerna 5-7 (VAIL_S, DEST_S och Pro_S) är baserade på varianter av formler 2-4 som föreställer sig att ett tvåvärt svavel ersätter en av kolatomerna. Strukturantal återspeglar antalet molekyler som genereras av MAYGEN för varje formel ("Initial") och antalet molekyler som återstår efter filtrering av de med oönskade understrukturer ("Final"). Förkortningar: VAIL = valin, alanin, isoleucin, leucin; DEST = asparaginsyra, glutaminsyra, serin, treonin; X_S = Tvåvärt svavel ersätter ett av kolen i bibliotek X; N-meX = N-metylX.
De allmänna metoderna ovan tillämpades på formler baserade på de "tidiga" kodade aminosyrorna, enligt proceduren för Meringer et al.21 Badliststrukturer togs från samma källa och konverterades till SMARTS-strängar för att enkelt representera substrukturella mönster. Två badlist-understrukturer användes inte i det här exemplet: struktur 018 (CH3-CH-N) matchade nära isomerer av prolin som inte själva var instabila; struktur 106 (R-C-C-OH, där R = alaninunderstruktur som fäster vid beta-kolet) matchade glutaminsyra, en kodad aminosyra. Förutom dessa kemiska formler skapades varianter med tvåvärt svavel som tog platsen för en kolatom och två väteatomer. Av prestandaskäl använder flera av dessa formler en trivalent fosforatom (t.ex. en "pseudoatom") som ersättning för beta-kolet i en alaninunderstruktur. Tabell 1 visar de bibliotek som genereras i det här exemplet, formlerna som används för att generera dem och antalet föreningar som finns i. Biblioteksnamn är baserade på de kodade aminosyrorna från vilka de härrör: antingen med 3-bokstavsförkortningen (Gly = glycin, Pro = prolin) eller förkortning med en bokstav (VAIL = Valine, Alanin, Isoleucin, Leucin; DEST = Asparaginsyra, glutaminsyra, serin, treonin). Suffixet "_S" indikerar att ett svavel ersattes med ett kol i det ursprungliga bibliotekets formel (t.ex. VAIL_S är byggt med samma fuzzy formel som VAIL, men med ett tvåvärt svavel som ersätter ett av kolen).
Efter strukturgenerering med MAYGEN filtrerades de resulterande biblioteken av föreningar som innehöll minst en understruktur som finns i badlisten. Efter denna filtrering ersattes alla fosforatomer med en alaninunderstruktur. Därefter skapades "begränsade" versioner av alla strukturer, med en acetylgrupp tillsatt till N-terminalen och en N-metylamidgrupp tillsatt till C-terminalen. Detta gjordes för att avlägsna effekten på hydrofobiciteten hos de fria amin- och karboxylsyragrupperna i alfa-aminosyrans ryggrad. PaDEL-Descriptor användes för att beräkna XLogP för alla begränsade strukturer och beräknade van der Waals volym (VABC) för alla obegränsade strukturer.
Figur 2 visar det kemiska utrymmet för de filtrerade biblioteken, enligt definitionen av VABC- och XLogP-deskriptorer. Här ökar utbudet av möjliga logP-värden med molekylär volym, även inom bibliotek som uttryckligen saknar hydrofila sidokedjor (t.ex. VAIL, Pro). Kodade aminosyror med kolväte-sidokedjor var mer hydrofoba än de flesta andra aminosyror med jämförbar volym från deras respektive bibliotek. Detta verkar också vara fallet för Met och Cys jämfört med andra medlemmar i VAIL_S bibliotek med liknande volymer. Kodade aminosyror med hydroxyl sidokedjor (Ser och Thr) var bland de minsta medlemmarna i DEST-biblioteket, med Asp bara något större än Thr.
Figur 3 och figur 4 visar påverkan på volym och logP när ett tvåvärt svavel ersätter ett kol i en alfa-aminosyra sidokedja. Svavelsubstitution ledde till en liten ökning av molekylvolymen i alla bibliotek (Figur 3). Effekten av svavelsubstitution på logP är inte lika homogen som för volym (figur 4). Den genomsnittliga logP för VAIL_S-biblioteket är något lägre än för VAIL-biblioteket, men den här effekten ses inte i något av de andra biblioteksparen (DEST och DEST_S, Pro och Pro_S).
Figur 5 kvantifierar effekterna på strukturgenerering av en pseudoatom som står för en gemensam understruktur; här ersatte en trivalent P en alanindel under strukturgenerering. Att använda en pseudoatom i strukturgenerering minskade kraftigt antalet strukturer som genererades av ~ 3 storleksordningar (figur 5A) och den totala tiden som behövs för att generera dessa strukturer med 1-2 storleksordningar (figur 5B).
<img alt="Figure 2" class="xfigimg" src="/files/ftp_upload/63632/63632fig02.jpg"> Figur 2: Kemiskt utrymme för alla filtrerade aminosyrabibliotek. Svarta markörer representerar aminosyror från bibliotek utan svavel; gula markörer representerar aminosyror från svavelberikade bibliotek. Cirklar: VAIL och VAIL_S; rutor: DEST och DEST_S; trianglar: Pro och Pro_S; stjärnor: kodade aminosyror. Observera att de två svavelhaltiga kodade aminosyrorna (Met och Cys) inte anses vara "tidiga" aminosyror utan finns i VAIL_S bibliotek. Förkortningar: XLogP = fördelningskoefficient; VAIL = valin, alanin, isoleucin, leucin; DEST = asparaginsyra, glutaminsyra, serin, treonin; X_S = Tvåvärt svavel ersätter ett av kolen i bibliotek X. <a href="https://www.jove.com/files/ftp_upload/63632/63632fig02large.jpg" target="_blank">Klicka här för att se en större version av denna figur.</a>
<img alt="Figure 3" class="xfigimg" src="/files/ftp_upload/63632/63632fig03.jpg"> Figur 3: Medelvärden för van der Waals volymer (i Å3) av bibliotek med och utan svavel. Svarta staplar representerar medelvolymerna av bibliotek utan svavel (VAIL, DEST, Pro), medan gula staplar representerar medelvolymer av de svavelsubstituerade versionerna av dessa bibliotek (VAIL_S, DEST_S, Pro_S). Felstaplar visar standardavvikelse. Förkortningar: VAIL = valin, alanin, isoleucin, leucin; DEST = asparaginsyra, glutaminsyra, serin, treonin; X_S = Tvåvärt svavel ersätter ett av kolen i bibliotek X. <a href="https://www.jove.com/files/ftp_upload/63632/63632fig03large.jpg" target="_blank">Klicka här för att se en större version av denna figur.</a>
<img alt="Figure 4" class="xfigimg" src="/files/ftp_upload/63632/63632fig04.jpg"> Bild 4: Genomsnittliga XLogP-värden för bibliotek med och utan svavel. Svarta staplar representerar bibliotek utan svavel (VAIL, DEST, Pro), medan gula staplar representerar svavelsubstituerade versioner av dessa bibliotek (VAIL_S, DEST_S, Pro_S). Felstaplar visar standardavvikelse. Förkortningar: XLogP = fördelningskoefficient; VAIL = valin, alanin, isoleucin, leucin; DEST = asparaginsyra, glutaminsyra, serin, treonin; X_S = Tvåvärt svavel ersätter ett av kolen i bibliotek X. <a href="https://www.jove.com/files/ftp_upload/63632/63632fig04large.jpg" target="_blank">Klicka här för att se en större version av denna figur.</a>
<img alt="Figure 5" class="xfigimg" src="/files/ftp_upload/63632/63632fig05.jpg"> Figur 5: Effekter av en trivalent pseudoatom på MAYGEN-strukturgenerering. Alla tester gjordes på en dator med en Intel i7-7700HQ-processor vid 2, 8 GHz, 16 GB RAM, inga spara strukturer i en fil och -m-alternativet att använda multithreading. Tester med en pseudoatom använde de luddiga formlerna som beskrivs i tabell 1. För tester utan pseudoatom var de fuzzy formler som användes desamma som beskrivs i tabell 1 med följande ändringar: P ersattes med N; kolantalet ökades med 3; väteantalet ökades med 7; syreantalet ökades med 2. Svarta staplar visar bibliotek som genereras med en pseudoatom; grå staplar visar bibliotek som genereras utan pseudoatom. (A) Antal strukturer som genereras med hjälp av de fuzzy formler som används för att bygga VAIL- och DEST-biblioteken med och utan en trivalent fosfor som ersätter en alaninunderstruktur. (B) Tid (i ms) som behövs för att bygga VAIL- och DEST-biblioteken med och utan en trivalent fosfor som ersätter en alaninunderstruktur. Förkortningar: VAIL = valin, alanin, isoleucin, leucin; DEST = asparaginsyra, glutaminsyra, serin, treonin. <a href="https://www.jove.com/files/ftp_upload/63632/63632fig05large.jpg" target="_blank">Klicka här för att se en större version av denna siffra.</a>
Tilläggsfil 1: Anteckningsbok för screening av understruktur. <a href="https://www.jove.com/files/ftp_upload/63632/supple_1.zip" target="_blank">Klicka här för att ladda ner den här filen.</a>
Kompletterande fil 2: Exempel på badlist. <a href="https://www.jove.com/files/ftp_upload/63632/supple_2.zip" target="_blank">Klicka här för att ladda ner den här filen.</a>
Kompletterande fil 3: Exempel på bra lista. <a href="https://www.jove.com/files/ftp_upload/63632/supple_3.zip" target="_blank">Klicka här för att ladda ner den här filen.</a>
Kompletterande fil 4: Pseudoatom ersättning anteckningsbok. <a href="https://www.jove.com/files/ftp_upload/63632/supple_4.zip" target="_blank">Klicka här för att ladda ner den här filen.</a>
Kompletterande fil 5: Aminosyra capping anteckningsbok. <a href="https://www.jove.com/files/ftp_upload/63632/supple_5.zip" target="_blank">Klicka här för att ladda ner den här filen.</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

Syftet med detta protokoll är att effektivt generera och kurera småmolekylära strukturbibliotek med hjälp av programvara med öppen källkod.

Kurering av beräkningskemiska bibliotek demonstrerade med alfa-aminosyror

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

Vårt protokoll visar hur programvara med öppen källkod kan göra det möjligt för alla forskare att skapa och kurera ett beräkningsstrukturbibliotek. Detta protokolls överklagande kommer från dess öppenhet och flexibilitet. Vem som helst kan använda den och ändra den så att den passar deras specifika forskningsfråga.Versioner av detta protokoll kan tillämpas på läkemedelsupptäcktsapplikationer, vilket snabbt skapar specifika strukturbibliotek för in silico-screening. Även om protokollet förklaras steg för steg, om användare inte är bekanta med Java eller grundläggande kodning kan de först titta på dem innan de implementerar protokollet. Börja med att skapa en ny katalog för projektet.Placera alla filer och körbara filer i den här katalogen för enkel åtkomst. Ladda ner den senaste versionen av Maygen som en jar-fil och pakethanteringsprogrammet Anaconda. På Windows-system söker du efter Anaconda-prompten och klickar på den resulterande genvägen för att köra.För att skapa en RDKit-miljö i Anaconda och för att ladda ner RDKit till miljön, skriv kommandot som visas på skärmen, tryck på enter för att köra och svara ja på alla frågor som kommer upp under installationen. Ladda sedan ned Jupyter Notebooks och textfilerna för substratmönstren från tilläggsfilerna, en till fem. I kommandotolken navigerar du till katalogen som innehåller maygen.jar körbar fil. För varje kemisk formel av intresse använd kommandot som visas på skärmen för att köra Maygen. Om formeln är en luddig formel i stället för en diskret formel ersätter du bindestreckets F-flagga med en fuzzy-flagga med bindestreck och omger eventuella elementintervall inom hakparenteser.I en Anaconda-prompt navigerar du till mappen som innehåller Jupyter Notebooks och aktiverar RDKit-miljön. De nedladdade anteckningsböckerna kräver RDKit. Så all framtida användning av dem i detta protokoll kräver att de öppnas i RDKit-miljön.Öppna sedan Jupyter Notebook för understrukturfiltrering och stäng filnamnet inom citattecken om det innehåller blanksteg. I den angivna cellen i början av anteckningsboken anger du den fullständiga filvägen för den inmatade sdf-filen. Den fullständiga filsökvägen för önskad sdf-utdatafil och filsökvägen för den dåliga listfilen som strängar.Om vissa understrukturer i det filtrerade biblioteket eller en bra lista behöver behållas skapar du en txt-fil med SMARTS-mönster för dessa understrukturer och placerar sökvägen till good list-filen på den angivna raden i början av anteckningsboken. Från menyn högst upp väljer du kärna, startar om och kör alla för att starta om notebook-kärnan och köra alla celler. En sdf-fil med önskat namn skapas i den angivna utdatamappen.Upprepa dessa steg för varje strukturfil som genereras av Maygen. För pseudoatombyte öppnar du en Anaconda-prompt, navigerar till mappen som innehåller Jupyter Notebooks och aktiverar RDKit-miljön. Öppna sedan Jupyter Notebook för pseudoatombyte.I den angivna cellen i början av anteckningsboken anger du den fullständiga filsökvägen för den inmatade sdf-filen och den fullständiga filsökvägen för önskad sdf-utdatafil som strängar. Starta om notebook-kerneln och kör alla celler för att hämta en sdf-fil med önskat namn i den angivna utdatamappen. På samma sätt öppnar du en Anaconda-prompt för aminosyra N och C-termini-kapning.Gå till mappen som innehåller Jupyter Notebooks och aktivera RDKit-miljön. Öppna Jupyter Notebook för aminosyrakapsling. I den angivna cellen i början av anteckningsboken anger du den fullständiga filsökvägen för den inmatade sdf-filen och den fullständiga filsökvägen för önskad sdf-utdatafil som strängar.Starta om notebook-kerneln och kör alla celler för att hämta en sdf-fil med önskat namn i den angivna utdatamappen. För deskriptorgenereringen placerar du alla sdf-filer för vilka deskriptorer ska beräknas i en enda mapp. Ladda sedan ner PaDEL-beskrivningen, packa upp den och extrahera den till den mappen.Öppna en kommandotolk, navigera till mappen som innehåller PaDEL-deskriptor jar-filen och kör PaDEL-deskriptorn för de insamlade sdf-filerna. Det kemiska utrymmet för alla filtrerade aminosyrabibliotek visas här. Svarta markörer representerar aminosyror från biblioteken utan svavel och gula markörer representerar aminosyror från svavelberikade bibliotek.Här representeras VAIL- och VAIL_S-bibliotek av cirklar. DEST- och DEST_S-bibliotek representeras av rutor. Proline- och Pro S-bibliotek representeras av trianglar och stjärnor representerar kodade aminosyror.Intervallet av möjliga log P-värden ökar med molekylvolymen även inom biblioteken som uttryckligen saknar hydrofila sidokedjor. Kodade aminosyror med kolväte sidokedjor är mer hydrofoba än de flesta andra aminosyror med en jämförbar volym från deras respektive bibliotek. Detta är också fallet för metionin som insisterar jämfört med andra medlemmar i VAILS-biblioteket med liknande volymer.Kodade aminosyror med hydroxyl sidokedjor var bland de minsta medlemmarna i DEST-biblioteket med asparaginsyra endast något större än tre Anine. Den representerade bilden visar de genomsnittliga Van der Waal-volymerna av bibliotek med svavel och utan svavel. Svavelsubstitution ledde till en liten ökning av molekylvolymen i alla bibliotek.De genomsnittliga fördelningskoefficientvärdena för bibliotek med och utan svavel visas här. Effekten av svavelsubstitution på log P är inte lika homogen som för volym. Den representativa bilden visar effekterna av en trivalent pseudoatom på Maygen-strukturgenerering.Att använda en pseudoatom i strukturgenerering minskade antalet strukturer som genererades med cirka tre storleksordningar under den totala tid som behövs för att generera dessa strukturer med en till två storleksordningar. Efter detta protokoll kan ytterligare funktioner integreras i framtiden baserat på forskarnas behov. Till exempel kan man integrera understrukturfilter i Maygen för att undvika efterbehandlingssteget.Biblioteksgenerering, kurering och modifiering. Denna allmänna process kan rymma andra molekylära strukturer och modifieringar med viss kodningskunskap, vilket gör det möjligt för forskare att utforska beräkningsbibliotek utöver alfaaminosyror. Detta protokoll kommer att hjälpa forskare att förbättra sitt beräkningsarbete inom livets ursprung.Verktygssatser med öppen källkod kommer i hög grad att hjälpa till med dessa ansträngningar.

Research

JoVE Journal

Chemistry

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

Fria radikaler i kemisk biologi: från Chemical beteende till biomarkörer utveckling

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

Kvalitativ Identifiering av karboxylsyror, boronsyror och aminer Använda Cruciform Fluoroforer

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

Sulfat Separation genom selektiv kristallisa med en Bis-iminoguanidinium ligand

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

Insikter Interaktioner av aminosyror och peptider med oorganiska material med hjälp av enda molekyl Force spektroskopi

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

Fastställande av kemiska sammansättningen av Korrosionsinhibitor / metall gränssnitt med XPS: Minimera Post Immersion Oxidation

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

Facile Preparation of (2<em&gt; Z</em&gt;, 4<em&gt; E</em&gt;) - Dienamider genom olefinering av elektron-bristande alkener med allylacetat

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

Fizzy Extraction av flyktiga organiska föreningar kombinerat med atmosfärstryck Kemisk ionisering 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

Kontinuerligt flöde kemi: Reaktion av Diphenyldiazomethane med p- Nitrobenzoic Acid

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

Kemisk fällning metod för syntesen av Nb2O5 modifierad Bulk Nickel katalysatorer med höga specifika yta

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

En två-stegs pyrolys-gaskromatografi metod med masspektrometrisk detektion för identifiering av Tattoo Ink ingredienser och förfalskade produkter

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