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 erkent financiering door de Carl-Zeiss-Foundation. Alle cijfers zijn gegenereerd met behulp van Microsoft Excel.

OPMERKING: Zie figuur 1 voor een samenvatting van het protocol en de tabel met materialen voor meer informatie over de gebruikte software.
<img alt="Figure 1" class="xfigimg" src="/files/ftp_upload/63632/63632fig01.jpg"> Figuur 1: Samenvatting stroomschema van het protocol. <a href="https://www.jove.com/files/ftp_upload/63632/63632fig01large.jpg" target="_blank">Klik hier om een grotere versie van deze figuur te bekijken.</a>
1. Software en bestandsdownloads
OPMERKING: Alle programma's zijn gratis voor individueel gebruik en kunnen op een pc worden uitgevoerd.
<ol><li>Maak een nieuwe map voor dit project. Plaats de bestanden en uitvoerbare bestanden hier voor eenvoudige toegang.</li>
	<li>Download en installeer de benodigde softwarepakketten.<ol>
<li>Download de nieuwste versie van MAYGEN als .jar-bestand. OPMERKING: MAYGEN is vrij beschikbaar als .jar bestand van https://github.com/MehmetAzizYirik/MAYGEN/releases</li>
		<li>Download en installeer de pakketbeheersoftware Conda en de cheminformatics toolkit RDKit18. OPMERKING: RDKit filtert de moleculaire structuren die door MAYGEN worden geproduceerd en werkt het beste in een Conda-omgeving. Instructies voor het downloaden van het Conda-platform zijn te vinden op https://conda.io/projects/conda/en/latest/user-guide/install/index.html. RdKit-installatie- en omgevingsinstellingen zijn te vinden op https://www.rdkit.org/docs/Install.html.<ol>
<li>Installeer RDKit in de belangrijkste Conda-omgeving in plaats van een afzonderlijke RDKit-omgeving via de Anaconda-prompt. Zoek op Windows-systemen naar "Anaconda-prompt" en klik op de resulterende snelkoppeling om uit te voeren. Op MacOS- en Linux-systemen kunt u via de terminal met Conda communiceren zonder extra programma's uit te voeren. Typ vervolgens de volgende opdracht en druk op Enter om uit te voeren en antwoord ja op eventuele vragen die tijdens de installatie opkomen: conda installeer -c rdkit rdkit. Hoewel er veel vrij beschikbare descriptorberekeningsprogramma's zijn, gebruikt dit voorbeeld PaDEL-Descriptor19, een gratis en snelle rekenmachine voor moleculaire descriptoren en vingerafdrukken.</li>
		</ol>
</li>
		<li>Download en sla het .jar bestand op in de projectmap. LET OP: PaDEL-Descriptor is gratis te downloaden vanaf http://www.yapcwsoft.com/dd/padeldescriptor/.</li>
	</ol>
</li>
	<li>Download de Jupyter notitieblokken en tekstbestanden van substructuurpatronen uit Supplemental Files 1-5. OPMERKING: Jupyter-notebooks kunnen ook worden gedownload van de volgende GitHub-pagina: https://github.com/cmayerb1/AA-structure-manip.</li>
</ol>2. Structuurgeneratie met MAYGEN
<ol><li>Navigeer in een opdrachtprompt naar de map met het uitvoerbare bestand MAYGEN .jar.</li>
	<li>Voer voor elke chemische formule van belang MAYGEN uit met de volgende opdracht: java -jar [MAYGEN .jar bestandsnaam] -f [chemische formule] -v -o [map voor MAYGEN output] -m -sdf. OPMERKING: Hiermee wordt een .sdf bestand opgeslagen in de aangewezen map, vernoemd naar de gebruikte formule.<ol>
<li>Als de formule een vage formule is in plaats van een discrete formule, vervang dan de -f-vlag door een -fuzzy flag en plaats eventuele elementintervallen tussen haakjes (gebruik bijvoorbeeld C[5-7]H[12-15] om ervoor te zorgen dat alle gegenereerde structuren tussen 5 en 7 koolstofatomen en tussen 12 en 15 waterstofatomen hebben).</li>
	</ol>
</li>
</ol>3. Filterverbindingen met ongewenste onderconstructies
<ol><li>Open een Anaconda-prompt (zie stap 1.2.2.1) en navigeer naar de map met de Jupyter-notitieblokken die zijn gedownload van Supplemental File 1.</li>
	<li>Open het Jupyter-notitieblok voor substructuurfiltering met de volgende opdracht: jupyter notebook [notebook bestandsnaam]
</li>
	<li>Voer in de aangewezen cel aan het begin van het notitieblok het volledige bestandspad van het invoer- .sdf bestand (gegenereerd door MAYGEN), het volledige bestandspad van het gewenste .sdf uitvoerbestand en het bestandspad van het "badlist"-bestand in als tekenreeksen (binnen aanhalingstekens). Zie Aanvullend bestand 2 voor een voorbeeld van een badlist.<ol>
<li>Als sommige substructuren in de gefilterde bibliotheek (een goodlist) moeten worden behouden, maakt u een .txt bestand met SMARTS-patronen20 voor die substructuren (een goodlist) en plaatst u het goodlist-bestandspad in de aangewezen regel aan het begin van het notitieblok. Zie Aanvullend bestand 3 voor een voorbeeld van een goodlist.</li>
	</ol>
</li>
	<li>Start de notebookkernel opnieuw op en voer alle cellen uit (selecteer in het menu bovenaan Kernel, Opnieuw opstarten en Alles uitvoeren) om een .sdf bestand met de gewenste naam in de opgegeven uitvoermap te krijgen.</li>
	<li>Herhaal de vorige twee stappen voor elk structuurbestand dat door MAYGEN is gegenereerd in stap 2.</li>
</ol>4. (Optioneel) Aanvullende structuurwijzigingen
OPMERKING: Deze worden in dit voorbeeld uitgevoerd, maar zijn mogelijk niet nodig voor het beheren van andere bibliotheken.
<ol><li>Pseudoatom vervanging. OPMERKING: Hier is een pseudoatom een uniek atoom dat wordt gebruikt om een grotere substructuur weer te geven die wordt gedeeld door alle gegenereerde structuren, waardoor de generatietijd van MAYGEN wordt verkort. Zie Aanvullend bestand 4 voor een voorbeeld van pseudoatom-vervanging.<ol>
<li>Open een Anaconda-prompt (zie stap 1.2.2.1) en navigeer naar de map met de Jupyter-notitieblokken.</li>
		<li>Open het Jupyter-notitieblok voor pseudoatom-vervanging: jupyter notebook [notebook bestandsnaam]
</li>
		<li>Voer in de aangewezen cel aan het begin van het notitieblok het volledige bestandspad van het invoer- .sdf bestand en het volledige bestandspad van het gewenste .sdf uitvoerbestand in als tekenreeksen (tussen aanhalingstekens).</li>
		<li>Start de notebookkernel opnieuw op en voer alle cellen uit om een .sdf bestand met de gewenste naam in de opgegeven uitvoermap te krijgen.</li>
	</ol>
</li>
	<li>Aminozuur N- en C-termini capping OPMERKING: Deze procedure is specifiek voor alfa-aminozuren, waarbij moleculaire doppen worden toegevoegd aan de N- en C-termini van alfa-aminozuurruggegraten. Zie Supplemental File 5 voor een voorbeeld van aminozuur capping.<ol>
<li>Open een Anaconda-prompt (zie stap 1.2.2.1) en navigeer naar de map met de Jupyter-notitieblokken.</li>
		<li>Open het Jupyter-notitieblok voor aminozuurafdekking: jupyter notebook [notebook bestandsnaam]
</li>
		<li>Voer in de aangewezen cel aan het begin van het notitieblok het volledige bestandspad van het invoer- .sdf bestand en het volledige bestandspad van het gewenste .sdf uitvoerbestand in als tekenreeksen (tussen aanhalingstekens).</li>
		<li>Start de notebookkernel opnieuw op en voer alle cellen uit om een .sdf bestand met de gewenste naam in de opgegeven uitvoermap te krijgen.</li>
	</ol>
</li>
</ol>5. Descriptor generatie
<ol><li>Plaats voorafgaand aan het genereren van descriptor alle .sdf bestanden waarvoor descriptoren moeten worden berekend in één map. OPMERKING: Als dit nog niet is gebeurd, geef deze bestanden dan beschrijvende namen voor eenvoudige filtering na het genereren van descriptor.</li>
	<li>Open een opdrachtprompt en navigeer naar de map met het paDEL-descriptorbestand .jar.</li>
	<li>Voer PaDEL-Descriptor uit voor de verzamelde .sdf bestanden met de volgende opdracht: java -jar PaDEL-Descriptor.jar -dir [directory of the .sdf files] -file [bestandspad van een .csv bestand voor resultaten] -2d -retainorder -usefilenameasmolname OPMERKING: Het resultatenbestand heeft de molecuulnaam in de eerste kolom en elke descriptor in de volgende kolommen.</li>
	<li>Exporteer deze gegevens naar spreadsheetsoftware voor verdere analyse.</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>

De auteurs hebben geen belangenconflicten te onthullen.

Een kenmerk van de "vroege" aminozuren is een gebrek aan zwavel. De eerder genoemde meta-analyses beschouwen de zwavelhoudende gecodeerde aminozuren (Cys en Met) over het algemeen als relatief late toevoegingen aan de genetische code, conclusies ondersteund door het ontbreken van zwavelhoudende aminozuren in meteorieten en vonkbuisexperimenten. Organozwavelverbindingen worden echter gemakkelijk gedetecteerd in kometen en meteorieten22, en heranalyse van vonkbuisexperimenten met behulp van H2S-gas vond aminozuren en andere organische verbindingen die zwavelbevatten 16. Bij het overwegen van een alternatief aminozuuralfabet, is een verrijkt met zwavel het verkennen waard.
In het bovenstaande protocol worden structuurgeneratie en substructuurfiltering beschouwd als kritieke stappen; afhankelijk van de samenstelling van de voltooide structuurbibliotheek, hoeft een onderzoeker mogelijk alleen die twee stappen uit te voeren. Instructies en software voor aanvullende acties (pseudoatom vervanging en toevoeging van substructuren (in dit geval aminozuur capping)) zijn opgenomen voor meer relevante descriptorberekening (capping zorgt ervoor dat XLogP-berekeningen worden beïnvloed door de sidechain en niet door de backbone-amine- of carboxylgroepen) en snellere structuurgeneratie via het gebruik van een pseudoatom, die hieronder in meer detail wordt besproken. Bovendien wordt hier descriptorberekening uitgevoerd als een eenvoudige manier om de diversiteit van de gegenereerde structuren te visualiseren en de effecten van zwavelverrijking in de voltooide bibliotheken te vergelijken.
Terwijl PaDEL-Descriptor duizenden moleculaire eigenschappen kan berekenen, werden moleculair volume (zoals berekend van der Waals volume) en verdelingscoëfficiënt (als XLogP) hier om twee verschillende redenen gebruikt. Ten eerste meten deze twee descriptoren moleculaire eigenschappen (respectievelijk grootte en hydrofobiciteit) die bekend zijn bij de meeste chemici en biologen. Ten tweede, in het geval van aminozuren, zijn deze twee eigenschappen significant. Decennialang was bekend dat aminozuurgrootte en hydrofobiciteit de thermodynamica van eiwitvouwing beïnvloeden23. Deze twee eigenschappen helpen bij het verklaren van aminozuursubstitutiefrequenties die een integraal onderdeel zijn geweest van het begrijpen van eiwitevolutie24.
Het bovenstaande voorbeeld laat zien dat in de twee bestudeerde descriptoren (moleculair volume en hydrofobiciteit) het vervangen van een tweewaardige zwavel door een koolstof en twee waterstofatomen geen significante veranderingen oplevert. De lichte, niet-significante toename van het gemiddelde moleculaire volume van zwavelsubstitutie (figuur 3) kan worden toegeschreven aan de grotere covalente straal van zwavel (~ 103 pm) in vergelijking met sp3 (~ 75 pm) of sp2 (~ 73 pm) koolstof25. Evenzo heeft zwavelsubstitutie een minimaal effect op de gemiddelde XLogP (figuur 4). Het grootste effect was tussen de VAIL- en VAIL_S bibliotheken, waarschijnlijk als gevolg van een combinatie van de VAIL-bibliotheek die vooral hydrofoob is (de sidechains zijn alleen koolwaterstoffen) en sulfhydrylgroepen die veel zuurder zijn dan de methylgroepen die ze zouden vervangen. Het minimale effect van zwavelsubstitutie is duidelijk in figuur 2, waar bibliotheken met zwavelsubstitutie dezelfde chemische ruimte innemen als analoge bibliotheken zonder zwavelsubstitutie.
De afname van het aantal structuren (figuur 5A) en de tijd die nodig is om die structuren te genereren (figuur 5B) bij het gebruik van een pseudoatom is niet verrassend. Het gebruik van een pseudoatom vermindert het aantal zware atomen dat in een chemische grafiek moet worden opgenomen, waardoor het aantal grafiekknooppunten wordt verminderd en exponentiële afnames in de generatietijd en het aantal structuren worden gegenereerd. Hier komt de keuze voor driewaardig fosfor als pseudoatom voort uit de basisbiochemie (afwezige posttranslationele toevoeging van fosfaatgroepen, geen genetisch gecodeerde aminozuren bevatten fosfor) en de valentie van het atoom dat het zou vervangen (een driewaardig fosfor kan gemakkelijk worden vervangen door een tetravalente koolstof die afzonderlijk is gebonden aan een ander atoom of een groep atomen). Hoewel de verstrekte code voor pseudoatomsubstitutie specifiek is voor het vervangen van een driewaardig fosfor door een alaninesubstructuur, kunnen gebruikers de code aanpassen om te werken met verschillende pseudoatomen of vervangende substructuren, mogelijk met behulp van meerdere pseudoatomen tijdens de eerste structuurgeneratie, gevolgd door het vervangen van elke pseudoatom door een grotere moleculaire substructuur.
Methoden voor het genereren van structuren die vergelijkbaar zijn met die van MAYGEN (en andere methoden zoals neurale netwerken) worden al gebruikt bij het ontdekken van geneesmiddelen om samengestelde bibliotheken te genereren voor in silico-screening ; een recente review4 bespreekt deze methoden in meer detail. Aangezien deze methoden voornamelijk bedoeld zijn voor het creëren van medicijnachtige moleculen, zijn er enkele beperkingen op hun vermogen om moleculen te genereren, zoals het gebruik van biologische of farmaceutische eigenschappen om de gecreëerde structuren te beperken (inverse QSPR / QSAR) of het creëren van structuren uit een vooraf ingesteld aantal substructuurbouwstenen. Omdat astrobiologie meer gericht is op de veelheid aan organische verbindingen die zich abiotisch kunnen vormen en minder op eindproducten of hun eigenschappen, is de uitputtende structuurgeneratie van MAYGEN ideaal voor het maken van structuurbibliotheken om astrobiologische vragen aan te pakken. De hier beschreven benadering van substructuurfiltering (uitgevoerd na het genereren van structuren via een extern programma) verschilt van het concurrerende programma MOLGEN doordat molgen's substructuurfiltering plaatsvindt tijdens het genereren van structuren. Omdat MAYGEN open-source is, is het niet alleen toegankelijker dan MOLGEN vanwege de licentiekosten van MOLGEN, maar individuen kunnen ook nieuwe functies implementeren, zoals substructuurfiltering tijdens het genereren van structuren.
Zoals geschreven, is het hier beschreven protocol gericht op het genereren en beheren van bibliotheken van relatief kleine alfa-aminozuren. Om verschillende bibliotheken te genereren, kunnen gebruikers verschillende molecuulformules aan MAYGEN geven, de substructuurfiltering wijzigen door de maximaal toegestane ringgrootte en bindingsvalentie te wijzigen, of de goodlist- en badlist-bestanden bewerken om substructuurpatronen toe te voegen of te verwijderen. Protocolwijzigingen waarbij de manier waarop atomen en substructuren worden toegevoegd of vervangen (pseudoatom-substitutie en moleculaire capping) worden gewijzigd, zijn haalbaar, maar vereisen meer aandacht voor valentiebeperkingen om RDKit-fouten over onjuiste valenties in gewijzigde structuren te voorkomen.
Het hierboven beschreven protocol is ontworpen voor kleine alfa-aminozuren. Het algemene formaat (uitgebreide structuurgeneratie met behulp van pseudoatomen, gevolgd door substructuurfiltering en moleculaire modificaties) is echter zeer flexibel voor verbindingen die verder gaan dan kleine aminozuren. Zelfs in de astrobiologie werd een vergelijkbare recente procedure met MOLGEN gebruikt om constitutionele isomeren van nucleïnezuren te onderzoeken26. Naast de hierboven beschreven tools, kan MAYGEN worden gecombineerd met andere open-source cheminformatics-tools om het maken en analyseren van nieuwe chemische structuren betaalbaar en toegankelijk te maken voor een breed scala aan onderzoeksgebieden.

Moleculaire structuurgeneratie dient als een praktische toepassing van het algemene probleem van uitputtende grafiekgeneratie; gezien verschillende knooppunten (atomen) en beperkingen op hun connectiviteit (bijv. Valenties, bindingsmultipliciteiten, gewenste / ongewenste substructuren), hoeveel verbonden grafieken (moleculen) zijn mogelijk? Structuurgeneratoren hebben uitgebreide toepassing gezien in medicijnontdekking en farmaceutische ontwikkeling, waar ze enorme bibliotheken met nieuwe structuren kunnen creëren voor in silico-screening 1.
De eerste structuurgenerator, CONGEN, werd ontwikkeld voor het eerste kunstmatige intelligentieproject in de organische chemie, DENDRAL2 (afkorting van DENDRitic ALgorithm). Verschillende software-opvolgers van DENDRAL werden in de literatuur vermeld; ze waren echter niet allemaal onderhouden of efficiënt. Momenteel is MOLGEN3 de state-of-the-art moleculaire structuurgenerator. Helaas voor de meeste potentiële gebruikers is het closed-source en vereist het een licentievergoeding. Er is dus behoefte aan een efficiënte open-source structuurgenerator die zich gemakkelijk kan aanpassen aan specifieke toepassingen. Een uitdaging voor een efficiënte structuurgenerator is het beheren van combinatorische explosies; naarmate de grootte van een molecuulformule toeneemt, neemt de grootte van de chemische zoekruimte exponentieel toe. Een recent overzicht onderzoekt verder de geschiedenis en uitdagingen van moleculaire structuurgeneratie4.
Vóór 2021 was de Parallel Molecule Generator (PMG)5 de snelste open-source structuurgenerator, maar deze was nog steeds langzamer dan MOLGEN in ordes van grootte. MAYGEN6 is ongeveer 47 keer sneller dan PMG en ongeveer 3 keer langzamer dan MOLGEN, waardoor MAYGEN de snelste en meest efficiënte open-source structuurgenerator is die beschikbaar is. Meer gedetailleerde vergelijkingen en benchmarkingtests zijn te vinden in de paper die MAYGEN6 introduceert. Een belangrijk kenmerk van het programma is de lexicografische ordeningstest voor canonieke structuren, een geordende methode voor het genereren van grafieken op basis van het Schreier-Sims7-algoritme . De software kan eenvoudig worden geïntegreerd in andere projecten en worden verbeterd voor de behoeften van de gebruikers.
Net als MOLGEN en PMG neemt MAYGEN een door de gebruiker gedefinieerde molecuulformule en genereert alle mogelijke structuren voor die formule. Als een gebruiker bijvoorbeeld MAYGEN uitvoert met de formule C5H12, genereert MAYGEN alle mogelijke structuren die vijf koolstofatomen en twaalf waterstofatomen bevatten. In tegenstelling tot zijn open-source tegenhanger PMG, kan MAYGEN ook "fuzzy" moleculaire formules bevatten die intervallen gebruiken in plaats van discrete getallen voor de telling van elk element. Als een gebruiker bijvoorbeeld MAYGEN uitvoert met de formule C5-7H12-15, genereert MAYGEN alle mogelijke structuren die tussen de vijf en zeven koolstofatomen en twaalf en vijftien waterstofatomen bevatten, waardoor eenvoudige generatie van structuren met een breed scala aan atomaire samenstellingen mogelijk is.
Astrobiologie is zo'n veld dat kan profiteren van moleculaire structuurgeneratoren. Een populair onderwerp in de astrobiologie is de evolutie van het aminozuuralfabet dat wordt gedeeld door al het bestaande leven op aarde. Een van de bepalende kenmerken van de Last Universal Common Ancestor (LUCA) is het gebruik van twintig genetisch gecodeerde aminozuren voor eiwitconstructie 8,9. Op basis van meta-analyses van werk op meerdere gebiedenvormen 10,11,12, ongeveer 10 van deze aminozuren (Gly, Ala, Val, Asp, Glu, Ser, Thr, Leu, Ile, Pro) gemakkelijk onder abiotische omstandigheden en vormen waarschijnlijk het aminozuuralfabet van pre-LUCA-organismen. In de loop van de tijd werd dit "vroege" alfabet uitgebreid als reactie op verschillende structurele en functionele behoeften. Een recente review van Moosmann13 beweert bijvoorbeeld dat de toevoeging van meer recente leden van de genetisch gecodeerde aminozuren (namelijk Met, Tyr en Trp) het mogelijk maakte om te overleven in zuurstofrijke omgevingen door de intracellulaire proliferatie van reactieve zuurstofsoorten te voorkomen.
Een steeds groter wordende reeks analytische chemietechnieken biedt inzicht in de aminozuurstructuren die zich onder abiotische omstandigheden kunnen vormen. Een recente review14 door Simkus en anderen beschrijft de methoden die worden gebruikt om talrijke organische verbindingen in meteorieten te detecteren, evenals organische verbindingen uit in vitro simulaties van vroege aardse omgevingen 15,16,17. Systematische generatie van chemische structuren stelt onderzoekers in staat om verder te gaan dan de organische verbindingen die via instrumentatie worden gedetecteerd, waarbij de structurele ruimte rond structurele "eilanden" wordt bevolkt die door analytische chemie zijn geïdentificeerd. In het geval van de "vroege" aminozuren toont deze systematische structuurgeneratie mogelijke eiwitchemie die beschikbaar is voor het vroege leven zonder de exploratie te beperken tot structuren die experimenteel zijn gedetecteerd onder abiotische syntheseomstandigheden. Met open-source cheminformatics toolkits en efficiënte structuurgeneratoren zoals MAYGEN is het maken en verkennen van nieuwe chemische structuurbibliotheken nu eenvoudiger dan ooit tevoren en kan het meer gedetailleerde onderzoeken naar alternatieve chemische stoffen van het leven begeleiden.

<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

Uitputtende generatie van moleculaire structuren heeft tal van chemische en biochemische toepassingen zoals medicijnontwerp, moleculaire databaseconstructie, verkenning van alternatieve biochemici en nog veel meer. Wiskundig gezien zijn dit grafiekgeneratoren met chemische beperkingen. In het veld is de meest efficiënte generator momenteel (MOLGEN) een commercieel product, waardoor het gebruik ervan wordt beperkt. Een andere moleculaire structuurgenerator, MAYGEN, is een recente open-source tool met een efficiëntie die vergelijkbaar is met MOLGEN en de capaciteit voor gebruikers om de prestaties te verbeteren door nieuwe functies toe te voegen. Een van de onderzoeksgebieden die baat kunnen hebben bij deze ontwikkeling is astrobiologie; structuurgeneratoren stellen onderzoekers in staat om experimentele gegevens aan te vullen met computationele mogelijkheden voor alternatieve biochemie. Dit protocol beschrijft één use case voor het genereren van structuren in de astrobiologie, namelijk het genereren en cureren van alfa-aminozuurbibliotheken. Met behulp van open-source structuurgeneratoren en cheminformatics-tools kunnen de hier beschreven praktijken buiten de astrobiologie worden geïmplementeerd voor het goedkoop maken en cureren van chemische structuurbibliotheken voor elke onderzoeksvraag.

<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">Bibliotheek</td>
			<td colspan="2" rowspan="2">Formule</td>
			<td colspan="2" rowspan="2">Aanvullende beperkingen</td>
			<td rowspan="2">"Vroege" gecodeerde aminozuren</td>
			<td rowspan="2">Generatietijd (ms)</td>
			<td colspan="2">Structuren</td>
		</tr>
<tr>
<td>Initiaal</td>
			<td>Finaal</td>
		</tr>
<tr>
<td>1</td>
			<td>Gly</td>
			<td colspan="2">C2H5NR2
</td>
			<td colspan="2">inclusief Gly substructuur</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-5NEE2H7-11
</td>
			<td colspan="2">Inclusief N-meGly of N-meAla substructuur</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">Inclusief N-meGly of N-meAla substructuur</td>
			<td></td>
			<td>3999</td>
			<td>75734</td>
			<td>10</td>
		</tr>
</tbody></table>Tabel 1: Samengestelde bibliotheken die in dit voorbeeld worden gebruikt. Bibliotheken gebouwd op formule 1-4 (Gly, VAIL, DEST en Pro) zijn gebaseerd op eerder gepubliceerde vage formules van de "vroege" gecodeerde aminozuren21, terwijl bibliotheken gebouwd uit formule 5-7 (VAIL_S, DEST_S en Pro_S) zijn gebaseerd op varianten van formule 2-4 die zich voorstellen dat een tweewaardige zwavel een van de koolstofatomen vervangt. Structuurtellingen weerspiegelen het aantal moleculen dat door MAYGEN wordt gegenereerd voor elke formule ("Initieel") en het aantal moleculen dat overblijft na het filteren van moleculen met ongewenste substructuren ("Final"). Afkortingen: VAIL = valine, alanine, isoleucine, leucine; DEST = asparaginezuur, glutaminezuur, serine, threonine; X_S = Tweewaardige zwavel vervangt een van de koolstofatomen in bibliotheek X; N-meX = N-methylX.
De bovenstaande algemene methoden werden toegepast op formules op basis van de "vroege" gecodeerde aminozuren, volgens de procedure van Meringer et al.21 Badlist-structuren werden uit dezelfde bron genomen en omgezet in SMARTS-strings om gemakkelijk substructurele patronen weer te geven. Twee badlist-substructuren werden in dit voorbeeld niet gebruikt: structuur 018 (CH 3-CH-N) kwam overeen met near-isomeren van proline die zelf niet onstabiel waren; structuur 106 (R-C-C-OH, waarbij R=alanine substructuur hecht aan de bèta-koolstof) kwam overeen met glutaminezuur, een gecodeerd aminozuur. Naast deze chemische formules werden varianten gemaakt met tweewaardige zwavel die de plaats innam van een koolstofatoom en twee waterstofatomen. Om prestatieredenen gebruiken verschillende van deze formules een driewaardig fosforatoom (bijvoorbeeld een "pseudoatom") als vervanging voor de bèta-koolstof van een alanine-substructuur. Tabel 1 bevat de bibliotheken die in dit voorbeeld zijn gegenereerd, de formules die zijn gebruikt om ze te genereren en het aantal verbindingen dat erin is opgenomen. Bibliotheeknamen zijn gebaseerd op de gecodeerde aminozuren waaruit ze zijn afgeleid: ofwel met behulp van de 3-letterige afkorting (Gly = glycine, Pro = proline) of eenletterige afkorting (VAIL = Valine, Alanine, Isoleucine, Leucine; DEST = Asparaginezuur, Glutaminezuur, Serine, Threonine). Het achtervoegsel "_S" geeft aan dat een zwavel werd vervangen door een koolstof in de formule van de oorspronkelijke bibliotheek (bijvoorbeeld VAIL_S is gebouwd met dezelfde wazige formule als VAIL, maar met een tweewaardige zwavel die een van de koolstofatomen vervangt).
Na het genereren van structuren met MAYGEN werden de resulterende bibliotheken gefilterd op verbindingen die ten minste één substructuur in de badlist bevatten. Na deze filtering werden alle fosforatomen vervangen door een alanine-substructuur. Vervolgens werden "afgetopte" versies van alle structuren gemaakt, waarbij een acetylgroep werd toegevoegd aan de N-terminus en een N-methylamidegroep werd toegevoegd aan de C-terminus. Dit werd gedaan om het effect op de hydrofobiciteit van de vrije amine- en carbonzuurgroepen in de alfa-aminozuurruggegraat te verwijderen. PaDEL-Descriptor werd gebruikt om XLogP te berekenen voor alle afgedekte structuren en berekende van der Waals volume (VABC) voor alle niet-afgedekte structuren.
Figuur 2 toont de chemische ruimte van de gefilterde bibliotheken, zoals gedefinieerd door VABC- en XLogP-descriptoren. Hier neemt het bereik van mogelijke logP-waarden toe met het moleculaire volume, zelfs binnen bibliotheken die expliciet hydrofiele sidechains missen (bijv. VAIL, Pro). Gecodeerde aminozuren met koolwaterstof sidechains waren meer hydrofoob dan de meeste andere aminozuren met een vergelijkbaar volume uit hun respectieve bibliotheek. Dit lijkt ook het geval te zijn voor Met en Cys in vergelijking met andere leden van de VAIL_S bibliotheek met vergelijkbare volumes. Gecodeerde aminozuren met hydroxyl zijketens (Ser en Thr) behoorden tot de kleinste leden van de DEST-bibliotheek, met Asp slechts iets groter dan Thr.
Figuur 3 en figuur 4 tonen de effecten op volume en logP wanneer een tweewaardige zwavel een koolstof in een alfa-aminozuur zijketen vervangt. Zwavelsubstitutie leidde tot een lichte toename van het moleculair volume in alle bibliotheken (figuur 3). Het effect van zwavelsubstitutie op logP is niet zo homogeen als voor volume (figuur 4). De gemiddelde logP van de VAIL_S bibliotheek is iets lager dan die van de VAIL-bibliotheek, maar dit effect wordt niet gezien in een van de andere bibliotheekparen (DEST en DEST_S, Pro en Pro_S).
Figuur 5 kwantificeert de effecten op structuurgeneratie van een pseudoatom die invalt voor een gemeenschappelijke substructuur; hier verving een driewaardige P een alaninegroep tijdens het genereren van de structuur. Het gebruik van een pseudoatom bij het genereren van structuren verminderde het aantal structuren gegenereerd door ~ 3 ordes van grootte (figuur 5A) en de totale tijd die nodig was om die structuren te genereren met 1-2 ordes van grootte (figuur 5B).
<img alt="Figure 2" class="xfigimg" src="/files/ftp_upload/63632/63632fig02.jpg"> Figuur 2: Chemische ruimte van alle gefilterde aminozuurbibliotheken. Zwarte markers vertegenwoordigen aminozuren uit bibliotheken zonder zwavel; gele markers vertegenwoordigen aminozuren uit met zwavel verrijkte bibliotheken. Cirkels: VAIL en VAIL_S; pleinen: DEST en DEST_S; driehoeken: Pro en Pro_S; sterren: gecodeerde aminozuren. Merk op dat de twee zwavelhoudende gecodeerde aminozuren (Met en Cys) niet als "vroege" aminozuren worden beschouwd, maar aanwezig zijn in de VAIL_S bibliotheek. Afkortingen: XLogP = verdelingscoëfficiënt; VAIL = valine, alanine, isoleucine, leucine; DEST = asparaginezuur, glutaminezuur, serine, threonine; X_S = Tweewaardige zwavel vervangt een van de koolstofatomen in bibliotheek X. <a href="https://www.jove.com/files/ftp_upload/63632/63632fig02large.jpg" target="_blank">Klik hier om een grotere versie van deze figuur te bekijken.</a>
<img alt="Figure 3" class="xfigimg" src="/files/ftp_upload/63632/63632fig03.jpg"> Figuur 3: Gemiddelde van der Waals-volumes (in Å3) van bibliotheken met en zonder zwavel. Zwarte balken vertegenwoordigen de gemiddelde volumes van bibliotheken zonder zwavel (VAIL, DEST, Pro), terwijl gele balken gemiddelde volumes van de door zwavel vervangende versies van die bibliotheken vertegenwoordigen (VAIL_S, DEST_S, Pro_S). Foutbalken geven de standaarddeviatie weer. Afkortingen: VAIL = valine, alanine, isoleucine, leucine; DEST = asparaginezuur, glutaminezuur, serine, threonine; X_S = Tweewaardige zwavel vervangt een van de koolstofatomen in bibliotheek X. <a href="https://www.jove.com/files/ftp_upload/63632/63632fig03large.jpg" target="_blank">Klik hier om een grotere versie van deze figuur te bekijken.</a>
<img alt="Figure 4" class="xfigimg" src="/files/ftp_upload/63632/63632fig04.jpg"> Figuur 4: Gemiddelde XLogP-waarden van bibliotheken met en zonder zwavel. Zwarte balken vertegenwoordigen bibliotheken zonder zwavel (VAIL, DEST, Pro), terwijl gele balken zwavelvervangende versies van die bibliotheken vertegenwoordigen (VAIL_S, DEST_S, Pro_S). Foutbalken geven de standaarddeviatie weer. Afkortingen: XLogP = verdelingscoëfficiënt; VAIL = valine, alanine, isoleucine, leucine; DEST = asparaginezuur, glutaminezuur, serine, threonine; X_S = Tweewaardige zwavel vervangt een van de koolstofatomen in bibliotheek X. <a href="https://www.jove.com/files/ftp_upload/63632/63632fig04large.jpg" target="_blank">Klik hier om een grotere versie van deze figuur te bekijken.</a>
<img alt="Figure 5" class="xfigimg" src="/files/ftp_upload/63632/63632fig05.jpg"> Figuur 5: Effecten van een driewaardige pseudoatom op maygen-structuurgeneratie. Alle tests werden uitgevoerd op een pc met een Intel i7-7700HQ-processor op 2,8 GHz, 16 GB RAM, geen opslagstructuren in een bestand en de -m-optie om multithreading te gebruiken. Bij proeven met een pseudoatom werden de fuzzy formules gebruikt zoals beschreven in tabel 1. Voor tests zonder pseudoatom waren de gebruikte fuzzy formules dezelfde als beschreven in tabel 1 met de volgende veranderingen: P werd vervangen door N; het aantal koolstofcellen werd verhoogd met 3; het aantal waterstof werd met 7 verhoogd; het zuurstofgehalte werd met 2 verhoogd. Zwarte balken tonen bibliotheken die zijn gegenereerd met een pseudoatom; grijze balken tonen bibliotheken die zijn gegenereerd zonder een pseudoatom. (A) Aantal structuren gegenereerd met behulp van de fuzzy formules die worden gebruikt om de VAIL- en DEST-bibliotheken te bouwen met en zonder een driewaardige fosfor die een alanine-substructuur vervangt. (B) Tijd (in ms) die nodig is om de VAIL- en DEST-bibliotheken te bouwen met en zonder een driewaardige fosfor die een alanine-substructuur vervangt. Afkortingen: VAIL = valine, alanine, isoleucine, leucine; DEST = asparaginezuur, glutaminezuur, serine, threonine. <a href="https://www.jove.com/files/ftp_upload/63632/63632fig05large.jpg" target="_blank">Klik hier om een grotere versie van deze figuur te bekijken.</a>
Aanvullend bestand 1: Substructuur screening notebook. <a href="https://www.jove.com/files/ftp_upload/63632/supple_1.zip" target="_blank">Klik hier om dit bestand te downloaden.</a>
Aanvullend bestand 2: Voorbeeld badlist. <a href="https://www.jove.com/files/ftp_upload/63632/supple_2.zip" target="_blank">Klik hier om dit bestand te downloaden.</a>
Aanvullend bestand 3: Voorbeeld goodlist. <a href="https://www.jove.com/files/ftp_upload/63632/supple_3.zip" target="_blank">Klik hier om dit bestand te downloaden.</a>
Aanvullend bestand 4: Pseudoatom vervangende notebook. <a href="https://www.jove.com/files/ftp_upload/63632/supple_4.zip" target="_blank">Klik hier om dit bestand te downloaden.</a>
Aanvullend bestand 5: Aminozuur capping notebook. <a href="https://www.jove.com/files/ftp_upload/63632/supple_5.zip" target="_blank">Klik hier om dit bestand te downloaden.</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

Het doel van dit protocol is om efficiënt kleine molecuulstructuurbibliotheken te genereren en te beheren met behulp van open-source software.

Curatie van computationele chemische bibliotheken gedemonstreerd met alfa-aminozuren

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

Ons protocol laat zien hoe open source software elke onderzoeker in staat kan stellen om een computationele structuurbibliotheek te maken en te beheren. Dit protocol-appèl komt voort uit de openheid en flexibiliteit. Iedereen kan het gebruiken en aanpassen aan hun specifieke onderzoeksvraag.Versies van dit protocol kunnen worden toegepast op toepassingen voor het ontdekken van geneesmiddelen, waardoor snel specifieke structuurbibliotheken worden gemaakt voor in silico-screening. Hoewel het protocol stap voor stap wordt uitgelegd, kunnen gebruikers, als ze niet bekend zijn met Java of basiscodering, eerst die bekijken voordat ze het protocol implementeren. Begin met het maken van een nieuwe map voor het project.Plaats alle bestanden en uitvoerbare bestanden in deze map voor eenvoudige toegang. Download de nieuwste versie van Maygen als jar-bestand en de pakketbeheersoftware Anaconda. Zoek op Windows-systemen naar de Anaconda-prompt en klik op de resulterende snelkoppeling om uit te voeren.Om een RDKit-omgeving in Anaconda te maken en de RDKit naar de omgeving te downloaden, typt u de opdracht op het scherm, drukt u op enter om uit te voeren en ja te antwoorden op eventuele vragen die tijdens de installatie opkomen. Download vervolgens de Jupyter Notebooks en tekstbestanden van de substraatpatronen uit de aanvullende bestanden, één tot vijf. Navigeer in de opdrachtprompt naar de map met de maygen.jar uitvoerbaar bestand. Gebruik voor elke chemische formule van belang de opdracht die op het scherm wordt weergegeven om Maygen uit te voeren. Als de formule een vage formule is in plaats van een discrete formule, vervangt u het koppelteken F door een vage teken met koppeltekens en plaatst u eventuele elementintervallen tussen haakjes.Navigeer in een Anaconda-prompt naar de map met de Jupyter Notebooks en activeer de RDKit-omgeving. Voor de gedownloade notebooks is RDKit vereist. Dus elk toekomstig gebruik ervan in dit protocol vereist dat ze worden geopend in de RDKit-omgeving.Open vervolgens het Jupyter-notitieblok voor substructuurfiltering en sluit de bestandsnaam tussen aanhalingstekens als deze spaties bevat. Voer in de aangewezen cel aan het begin van het notitieblok het volledige bestandspad van het sdf-invoerbestand in. Het volledige bestandspad van het gewenste sdf-uitvoerbestand en het bestandspad van het slechte lijstbestand als tekenreeksen.Als sommige substructuren in de gefilterde bibliotheek of een goede lijst moeten worden behouden, maakt u een txt-bestand met SMARTS-patronen voor die substructuren en plaatst u het pad naar het goede lijstbestand in de aangewezen regel aan het begin van het notitieblok. Selecteer in het menu bovenaan kernel, start opnieuw op en voer alles uit om de notebookkernel opnieuw op te starten en alle cellen uit te voeren. Een sdf-bestand met de gewenste naam wordt gemaakt in de opgegeven uitvoermap.Herhaal deze stappen voor elk structuurbestand dat door Maygen wordt gegenereerd. Voor pseudoatom-vervanging opent u een Anaconda-prompt, navigeert u naar de map met de Jupyter Notebooks en activeert u de RDKit-omgeving. Open vervolgens de Jupyter Notebook voor pseudoatom vervanging.Voer in de aangewezen cel aan het begin van het notitieblok het volledige bestandspad van het sdf-invoerbestand en het volledige bestandspad van het gewenste sdf-uitvoerbestand in als tekenreeksen. Start de notebookkernel opnieuw op en voer alle cellen uit om een sdf-bestand met de gewenste naam in de opgegeven uitvoermap te krijgen. Open op dezelfde manier een Anaconda-prompt voor aminozuur N- en C-termini-capping.Navigeer naar de map met de Jupyter Notebooks en activeer de RDKit-omgeving. Open de Jupyter Notebook voor aminozuur capping. Voer in de aangewezen cel aan het begin van het notitieblok het volledige bestandspad van het sdf-invoerbestand en het volledige bestandspad van het gewenste sdf-uitvoerbestand in als tekenreeksen.Start de notebookkernel opnieuw op en voer alle cellen uit om een sdf-bestand met de gewenste naam in de opgegeven uitvoermap te krijgen. Plaats voor het genereren van descriptoren alle sdf-bestanden waarvoor descriptoren moeten worden berekend in één map. Download vervolgens de PaDEL-descriptor, pak deze uit en pak deze uit naar die map.Open een opdrachtprompt, navigeer naar de map met het PaDEL-descriptor jar-bestand en voer de PaDEL-descriptor uit voor de verzamelde sdf-bestanden. De chemische ruimte van alle gefilterde aminozuurbibliotheken wordt hier weergegeven. Zwarte markers vertegenwoordigen aminozuren uit de bibliotheken zonder zwavel en gele markers vertegenwoordigen aminozuren uit met zwavel verrijkte bibliotheken.Hier worden VAIL- en VAIL_S bibliotheken weergegeven door cirkels. DEST- en DEST_S bibliotheken worden vertegenwoordigd door vierkanten. Proline- en Pro S-bibliotheken worden vertegenwoordigd door driehoeken en sterren vertegenwoordigen gecodeerde aminozuren.Het bereik van mogelijke log P-waarden neemt toe met het moleculaire volume, zelfs binnen de bibliotheken die expliciet hydrofiele zijketens missen. Gecodeerde aminozuren met koolwaterstofzijketens zijn meer hydrofoob dan de meeste andere aminozuren met een vergelijkbaar volume uit hun respectieve bibliotheek. Dit is ook het geval voor methionine in vergelijking met andere leden van de VAILS-bibliotheek met vergelijkbare volumes.Gecodeerde aminozuren met hydroxylzijketens behoorden tot de kleinste leden van de DEST-bibliotheek met asparaginezuur dat slechts iets groter was dan drie anine. De weergegeven afbeelding toont de gemiddelde Van der Waal-volumes van bibliotheken met zwavel en zonder zwavel. Zwavelsubstitutie leidde tot een lichte toename van het moleculair volume in alle bibliotheken.De gemiddelde verdelingscoëfficiëntwaarden van bibliotheken met en zonder zwavel worden hier weergegeven. Het effect van zwavelsubstitutie op log P is niet zo homogeen als voor volume. De representatieve afbeelding toont de effecten van een trivalent pseudoatom op maygenstructuurgeneratie.Het gebruik van een pseudoatom bij het genereren van structuren verminderde het aantal structuren gegenereerd met ongeveer drie ordes van grootte in de totale tijd die nodig was om die structuren te genereren met één tot twee ordes van grootte. Volgens dit protocol kunnen in de toekomst extra functionaliteiten worden geïntegreerd op basis van de behoeften van onderzoekers. Men zou bijvoorbeeld substructuurfilters in Maygen kunnen integreren om de nabewerkingsstap te vermijden.Bibliotheekgeneratie, curatie en modificatie. Dit algemene proces kan andere moleculaire structuren en modificaties accommoderen met enige codeerkennis, waardoor onderzoekers computationele bibliotheken kunnen verkennen die verder gaan dan die van alfa-aminozuren. Dit protocol zal onderzoekers helpen hun computationele werk in de oorsprong van het levensveld te verbeteren.Open source toolkits zullen enorm helpen bij deze inspanningen.

Research

JoVE Journal

Chemistry

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

Free Radicals in Chemical Biology: van chemische gedrag tot 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

Kwalitatieve Identificatie van carbonzuren, boronzuren en aminen gebruiken Cruciform Fluoroforen

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

Sulfaat Scheiden door selectieve kristallisatie met een 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

Inzicht in de interacties van aminozuren en peptiden met anorganische materialen met behulp van enkel-molecuul Force Spectroscopie

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

Het bepalen van de chemische samenstelling van Corrosion Inhibitor / Metal Interfaces met XPS: minimaliseren Bericht Immersion Oxidatie

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;) - Dienamiden door de olefinering van elektron-deficiënte alkenen met allylacetaat

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 Extractie van vluchtige organische verbindingen gecombineerd met atmosferische druk Chemische ionisatie 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

Continue stroom chemie: Reactie van Diphenyldiazomethane met 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

Chemische neerslag methode voor de synthese van Nb2O5 bewerkt Bulk nikkel katalysatoren met hoge specifieke oppervlakte

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

Een twee-stap pyrolyse-Gas Chromatografie methode met massa spectrometrische detectie voor de identificatie van tatoeage inkt ingrediënten en namaakproducten

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