Chemical graph generator

A chemical graph generator is a software package to generate computer representations of chemical structures adhering to certain boundary conditions. The development of such software packages is a research topic of cheminformatics. Chemical graph generators are used in areas such as virtual library generation in drug design, in molecular design with specified properties, called inverse QSAR/QSPR, as well as in organic synthesis design, retrosynthesis or in systems for computer-assisted structure elucidation (CASE). CASE systems again have regained interest for the structure elucidation of unknowns in computational metabolomics, a current area of computational biology.

History

Overlapping substructures of caffeine. Two substructures of a caffeine molecule are given, (A) and (B). The overlap of these substructures is highlighted in green in the caffeine structure (C).

Molecular structure generation is a branch of graph generation problems. Molecular structures are graphs with chemical constraints such as valences, bond multiplicity and fragments. These generators are the core of CASE systems. In a generator, the molecular formula is the basic input. If fragments are obtained from the experimental data, they can also be used as inputs to accelerate structure generation. The first structure generators were versions of graph generators modified for chemical purposes. One of the first structure generators was CONGEN,[1] originally developed for the DENDRAL project, the first artificial intelligence project in organic chemistry.[2] DENDRAL was developed as a part of the Mariner program launched by the NASA to search for life on Mars.[3][4] CONGEN dealt well with overlaps in substructures. The overlaps among substructures rather than atoms were used as the building blocks. For the case of stereoisomers, symmetry group calculations were performed for duplicate detection.

After DENDRAL, another mathematical method, MASS,[5] a tool for mathematical synthesis and analysis of molecular structures, was reported. As with CONGEN, the MASS algorithm worked as an adjacency matrix generator. Many mathematical generators are descendants of efficient branch-and-bound methods from Igor Faradjev[6] and Ronald C. Read's orderly generation method.[7] Although their reports are from the 1970s, these studies are still the fundamental references for structure generators. In the orderly generation method, specific order-check functions are performed on graph representatives, such as vectors. For example, MOLGEN[8] performs a descending order check while filling rows of adjacency matrices. This descending order check is based on an input valence distribution. The literature classifies generators into two major types: structure assembly and structure reduction. The algorithmic complexity and the run time are the criteria used for comparison.

Structure assembly

The generation process starts with a set of atoms from the molecular formula. In structure assembly, atoms are combinatorically connected to consider all possible extensions. If substructures are obtained from the experimental data, the generation starts with these substructures. These substructures provide known bonds in the molecule. One of the earliest attempts was made by Hidetsugu Abe in 1975 using a pattern recognition-based structure generator.[9] The algorithm had two steps: first, the prediction of the substructure from low-resolution spectral data; second, the assembly of these substructures based on a set of construction rules. Hidetsugu Abe and the other contributors published the first paper on CHEMICS,[10] which is a CASE tool comprising several structure generation methods. The program relies on a predefined non-overlapping fragment library. CHEMICS generates different types of component sets ranked from primary to tertiary based on component complexity. The primary set contains atoms, i.e., C, N, O and S, with their hybridization. The secondary and tertiary component sets are built layer-by-layer starting with these primary components. These component sets are represented as vectors and are used as building blocks in the process.

Substantial contributions were made by Craig Shelley and Morton Munk, who published a large number of CASE papers in this field. The first of these papers reported a structure generator, ASSEMBLE.[11] The algorithm is considered one of the earliest assembly methods in the field. As the name indicates, the algorithm assembles substructures with overlaps to construct structures. ASSEMBLE overcomes overlapping by including a “neighbouring atom tag”. The generator is purely mathematical and does not involve the interpretation of any spectral data. Spectral data are used for structure scoring and substructure information. Based on the molecular formula, the generator forms bonds between pairs of atoms, and all the extensions are checked against the given constraints. If the process is considered as a tree, the first node of the tree is an atom set with substructures if any are provided by the spectral data. By extending the molecule with a bond, an intermediate structure is built. Each intermediate structure can be represented by a node in the generation tree. ASSEMBLE was developed with a user-friendly interface to facilitate use. The second version of ASSEMBLE was released in 2000.[12] Another assembly method is GENOA.[13] Compared to ASSEMBLE and many other generators, GENOA is a constructive substructure search-based algorithm, and it assembles different substructures by also considering the overlaps.

The efficiency and exhaustivity of generators are also related to the data structures. Unlike previous methods, AEGIS[14] was a list-processing generator. Compared to adjacency matrices, list data requires less memory. As no spectral data was interpreted in this system, the user needed to provide substructures as inputs. Structure generators can also vary based on the type of data used, such as HMBC, HSQC and other NMR data. LUCY is an open-source structure elucidation method based on the HMBC data of unknown molecules,[15] and involves an exhaustive 2-step structure generation process where first all combinations of interpretations of HMBC signals are implemented in a connectivity matrix, which is then completed by a deterministic generator filling in missing bond information. This platform could generate structures with any arbitrary size of molecules; however, molecular formulas with more than 30 heavy atoms are too time consuming for practical applications. This limitation highlighted the need for a new CASE system. SENECA was developed to eliminate the shortcomings of LUCY.[16] To overcome the limitations of the exhaustive approach, SENECA was developed as a stochastic method to find optimal solutions. The systems comprise two stochastic methods: simulated annealing and genetic algorithms. First, a random structure is generated; then, its energy is calculated to evaluate the structure and its spectral properties. By transforming this structure into another structure, the process continues until the optimum energy is reached. In the generation, this transformation relies on equations based on Jean-Loup Faulon's rules.[17] LSD (Logic for Structure Determination)[18] is an important contribution from French scientists. The tool uses spectral data information such as HMBC and COSY data to generate all possible structures. LSD is an open source structure generator released under the General Public License (GPL). A well-known commercial CASE system, StrucEluc,[19] also features a NMR based generator. This tool is from ACD Labs and, notably, one of the developers of MASS, Mikhail Elyashberg. COCON[20] is another NMR based structure generator, relying on theoretical data sets for structure generation. Except J-HMBC and J-COSY, all NMR types can be used as inputs.

In 1994, Chinese scientists reported an integer partitioning-based structure generator.[21] The decomposition of the molecular formula into fragments, components and segments was performed as an application of integer partitioning. These fragments were then used as building blocks in the structure generator. This structure generator was part of a CASE system, ESESOC.[22]

Breadth-first search generation. Molecular structure generation is explained step by step. Starting from a set of atoms, bonds are added between atom pairs until reaching saturated structures.

A series of stochastic generators was reported by Jean-Loup Faulon. The software, MOLSIG,[23] was integrated into this stochastic generator for canonical labelling and duplicate checks.[24] As for many other generators, the tree approach is the skeleton of Jean-Loup Faulon's structure generators. However, considering all possible extensions leads to a combinatorial explosion. Orderly generation is performed to cope with this exhaustivity. Many assembly algorithms, such as OMG,[25] MOLGEN and Jean-Loup Faulon's structure generator,[26] are orderly generation methods. Jean-Loup Faulon's structure generator relies on equivalence classes over atoms. Atoms with the same interaction type and element are grouped in the same equivalence class. Rather than extending all atoms in a molecule, one atom from each class is connected with other atoms. Similar to the former generator, Julio Peironcely's structure generator, OMG, takes atoms and substructures as inputs and extends the structures using a breadth-first search method. This tree extension terminates when all the branches reach saturated structures.

OMG generates structures based on the canonical augmentation method from Brendan McKay's NAUTY package. The algorithm calculates canonical labelling and then extends structures by adding one bond. To keep the extension canonical, canonical bonds are added.[27] Although NAUTY is an efficient tool for graph canonical labelling, OMG is approximately 2000 times slower than MOLGEN.[28] The problem is the storage of all the intermediate structures. OMG has since been parallelized, and the developers released PMG (Parallel Molecule Generator).[29] MOLGEN outperforms PMG using only 1 core; however, PMG outperforms MOLGEN by increasing the number of cores to 10.

A constructive search algorithm is a branch-and-bound method, such as Igor Faradjev's algorithm, and an additional solution to memory problems. Branch-and-bound methods are matrix generation algorithms. In contrast to previous methods, these methods build all the connectivity matrices without building intermediate structures. In these algorithms, canonicity criteria and isomorphism checks are based on automorphism groups from mathematical group theory. MASS, SMOG[30] and Ivan Bangov's algorithm[31] are good examples in the literature. MASS is a method of mathematical synthesis. First, it builds all incidence matrices for a given molecular formula. The atom valences are then used as the input for matrix generation. The matrices are generated by considering all the possible interactions among atoms with respect to the constraints and valences. The benefit of constructive search algorithms is their low memory usage. SMOG is a successor of MASS.

Unlike previous methods, MOLGEN is the only maintained efficient generic structure generator, developed as a closed-source platform by a group of mathematicians as an application of computational group theory. MOLGEN is an orderly generation method. Many different versions of MOLGEN have been developed, and they provide various functions. Based on the users' needs, different types of inputs can be used. For example, MOLGEN-MS[32] allows users to input mass spectrometry data of an unknown molecule. Compared to many other generators, MOLGEN approaches the problem from different angles. The key feature of MOLGEN is generating structures without building all the intermediate structures and without generating duplicates.

In the field, the recent studies are from Kimito Funatsu's research group. As a type of assembly method, building blocks, such as ring systems and atom fragments, are used in the structure generation.[33] Every intermediate structure is extended by adding building blocks in all possible ways. To reduce the number of duplicates, Brendan McKay's canonical path augmentation method is used. To overcome the combinatorial explosion in the generation, applicability domain and ring systems are detected based on inverse QSPR/QSAR analysis.[34] The applicability domain, or target area, is described based on given biological as well as pharmaceutical activity information from QSPR/QSAR.[35] In that study, monotonically changed descriptors (MCD) are used to describe applicability domains. For every extension in intermediate structures, the MCDs are updated. The usage of MCDs reduces the search space in the generation process. In the QSPR/QSAR based structure generation, there is the lack of synthesizability of the generated structures. Usage of retrosynthesis paths in the generation makes the generation process more efficient. For example, a well-known tool called RetroPath[36] is used for molecular structure enumeration and virtual screening based on the given reaction rules.[37] Its core algorithm is a breadth-first method, generating structures by applying reaction rules to each source compound. Structure generation and enumeration are performed based on Brendan McKay's canonical augmentation method. RetroPath 2.0 provides a variety of workflows such as isomer transformation, enumeration, QSAR and metabolomics.

Besides these mathematical structure generation methods, the implementations of neural networks, such as generative autoencoder models,[38][39] are the novel directions of the field.

Structure reduction

Unlike these assembly methods, reduction methods make all the bonds between atom pairs, generating a hypergraph. Then, the size of the graph is reduced with respect to the constraints. First, the existence of substructures in the hypergraph is checked. Unlike assembly methods, the generation tree starts with the hypergraph, and the structures decrease in size at each step. Bonds are deleted based on the substructures. If a substructure is no longer in the hypergraph, the substructure is removed from the constraints. Overlaps in the substructures were also considered due to the hypergraphs. The earliest reduction-based structure generator is COCOA,[40] an exhaustive and recursive bond-removal method. Generated fragments are described as atom-centred fragments to optimize storage, comparable to circular fingerprints[41] and atom signatures.[42] Rather than storing structures, only the list of first neighbours of each atom is stored. The main disadvantage of reduction methods is the massive size of the hypergraphs. Indeed, for molecules with unknown structures, the size of the hyper structure becomes extremely large, resulting in a proportional increase in the run time.

The structure generator GEN[43] by Simona Bohanec combines two tasks: structure assembly and structure reduction. Like COCOA, the initial state of the problem is a hyper structure. Both assembly and reduction methods have advantages and disadvantages, and the GEN tool avoids these disadvantages in the generation step. In other words, structure reduction is efficient when structural constraints are provided, and structure assembly is faster without constraints. First, the useless connections are eliminated, and then the substructures are assembled to build structures. Thus, GEN copes with the constraints in a more efficient way by combining these methods. GEN removes the connections creating the forbidden structures, and then the connection matrices are filled based on substructure information. The method does not accept overlaps among substructures. Once the structure is built in the matrix representation, the saturated molecule is stored in the output list. The COCOA method was further improved and a new generator was built, HOUDINI.[44] It relies on two data structures: a square matrix of compounds representing all bonds in a hyper structure is constructed, and second, substructure representation is used to list atom-centred fragments. In the structure generation, HOUDINI maps all the atom-centred fragments onto the hyper structure.

Mathematical basis

Chemical graphs

Graph representation of the serotonin molecule. (A) Molecular structure of serotonin. (B) Graph representation of the molecule.

In a graph representing a chemical structure, the vertices and edges represent atoms and bonds, respectively. The bond order corresponds to the edge multiplicity, and as a result, chemical graphs are vertex and edge-labelled graphs. A vertex and edge-labelled graph ${\displaystyle G=(V,E)}$ is described as a chemical graph where ${\displaystyle V}$ is the set of vertices, i.e., atoms, and ${\displaystyle E}$ is the set of edges, which represents the bonds.

In graph theory, the degree of a vertex is its number of connections. In a chemical graph, the maximum degree of an atom is its valence, and the maximum number of bonds a chemical element can make. For example, carbon's valence is 4. In a chemical graph, an atom is saturated if it reaches its valence. A graph is connected if there is at least one path between each pair of vertices. Although chemical mixtures[45] are one of the main interests of many chemists, due to the computational explosion, many structure generators output only connected chemical graphs. Thus, the connectivity check is one of the mandatory intermediate steps in structure generation because the aim is to generate fully saturated molecules. A molecule is saturated if all its atoms are saturated.

Symmetry groups for molecular graphs

Molecular Symmetry. (A) The initial labelling of 2,4-Dimethyl-2-pentene. (B) and (C) are symmetries of the same molecule with different labels.

For a set of elements, a permutation is a rearrangement of these elements.[46] An example is given below:

Table 1: Permutation of set of integers.

${\displaystyle x}$	1	2	3	4	5	6	7	8	9	10	11
${\displaystyle f(x)}$	4	2	11	6	1	5	8	9	7	10	3

The second line of Table 1 shows a permutation of the first line. The multiplication of permutations, ${\displaystyle a}$ and ${\displaystyle b}$, is defined as a function composition, as shown below.

${\displaystyle (a*b)(x)=a(b(x))}$

The combination of two permutations is also a permutation. A group, ${\displaystyle G}$, is a set of elements together with an associative binary operation ${\displaystyle *}$ defined on ${\displaystyle G}$ such that the following are true:

• There is an element ${\displaystyle I}$ in ${\displaystyle G}$ satisfying ${\displaystyle g*I=g}$, for all elements ${\displaystyle g}$ of ${\displaystyle G}$.

• For each element of G, there is an element ${\displaystyle g^{-1}}$ such that ${\displaystyle g*g^{-1}}$ is equal to the identity element.

The order of a group is the number of elements in the group. Let us assume ${\displaystyle X}$ is a set of integers. Under the function composition operation, ${\displaystyle Sym(X)}$ is a symmetry group, the set of all permutations over X. If the size of ${\displaystyle X}$ is ${\displaystyle n}$, then the order of ${\displaystyle Sym(X)}$ is ${\displaystyle n!}$. Set systems consist of a finite set ${\displaystyle X}$ and its subsets, called blocks of the set. The set of permutations preserving the set system is used to build the automorphisms of the graph. An automorphism permutes the vertices of a graph; in other words, it maps a graph onto itself. This action is edge-vertex preserving. If ${\displaystyle (u,v)}$ is an edge of the graph, ${\displaystyle G=(V,E)}$, and ${\displaystyle a}$ is a permutation of ${\displaystyle V}$, then

${\displaystyle a({u,v})=(a(u),a(v))}$

A permutation ${\displaystyle a}$ of ${\displaystyle V}$ is an automorphism of the graph ${\displaystyle G=(E,V)}$ if

${\displaystyle a(u,v)}$ is an element of ${\displaystyle E}$, if ${\displaystyle {(u,v)}}$ is an element of ${\displaystyle E}$.

The automorphism group of a graph ${\displaystyle G}$, denoted ${\displaystyle Aut(G)}$, is the set of all automorphisms on ${\displaystyle V}$. In molecular graphs, canonical labelling and molecular symmetry detection are implementations of automorphism groups. Although there are well known canonical labelling methods in the field, such as InChI[47] and ALATIS,[48] NAUTY is a commonly used software package for automorphism group calculations and canonical labelling.

List of available structure generators

The available software packages and their links are listed below.

Table 2: List of Available Structure Generators.

Name	Link
ASSEMBLE	https://www.upstream.ch/main.html?src=%2Findex.html
COCON	https://cocon.nmr.de
DENDRAL CONGEN+GENOA	http://www.softwarepreservation.org/projects/AI/DENDRAL/DENDRAL-CONGEN_GENOA.zip/view
LSD	http://eos.univ-reims.fr/LSD/index_ENG.html
MOLGEN	https://www.molgen.de/
MOLSIG	http://molsig.sourceforge.net
OMG	https://sourceforge.net/p/openmg/
PMG	https://sourceforge.net/projects/pmgcoordination/
SENECA	https://github.com/steinbeck/seneca
SMOG	http://ccl.net/cca/software/MS-DOS/SMOG/index.shtml

References

This article was adapted from the following source under a CC BY 4.0 license (2021) (reviewer reports): "Chemical graph generators". PLOS Computational Biology. 17 (1): e1008504. 5 January 2021. doi:10.1371/JOURNAL.PCBI.1008504. ISSN 1553-734X. PMC 7785115. PMID 33400699. Wikidata Q104747658.

1. Bruccoleri RE; Karplus M (1 January 1987). "Prediction of the folding of short polypeptide segments by uniform conformational sampling". Biopolymers. 26 (1): 137–168. doi:10.1002/BIP.360260114. ISSN 0006-3525. PMID 3801593. Wikidata Q69715633.
2. Sutherland G. DENDRAL-a computer program for generating and filtering chemical structures. STANFORD UNIV CALIF DEPT OF COMPUTER SCIENCE; 1967 Feb 15.
3. Robert K. Lindsay; Bruce G. Buchanan; Edward A. Feigenbaum; Joshua Lederberg (June 1993). "DENDRAL: A case study of the first expert system for scientific hypothesis formation". Artificial Intelligence. 61 (2): 209–261. doi:10.1016/0004-3702(93)90068-M. ISSN 0004-3702. Wikidata Q29387651.
4. Karina A. Gulyaeva; Irina L. Artemieva (2020). "The Ontological Approach in Organic Chemistry Intelligent System Development". Advances in Intelligent Systems and Computing: 69–78. doi:10.1007/978-981-32-9343-4_7. ISSN 2194-5357. Wikidata Q105092432.
5. V.V. Serov; M.E. Elyashberg; L.A. Gribov (April 1976). "Mathematical synthesis and analysis of molecular structures". Journal of Molecular Structure. 31 (2): 381–397. doi:10.1016/0022-2860(76)80018-X. ISSN 0022-2860. Wikidata Q99232065.
6. Faradzev IA. Constructive enumeration of combinatorial objects, Colloq. Internat. CNRS, 1978, vol. 260, pp. 131–135.
7. Charles J. Colbourn; Ronald C. Read (1979). "Orderly algorithms for generating restricted classes of graphs". Journal of Graph Theory. 3 (2): 187–195. doi:10.1002/JGT.3190030210. ISSN 0364-9024. Wikidata Q99232279.
8. Grüner T, Laue R, Meringer M, Bayreuth U. Algorithms for group actions: Homomorphism principle and orderly generation applied to graphs. DIMACS Series in Discrete Mathematics and Theoretical Computer Science. 1997;28:113-22.
9. Hidetsugu. Abe; Peter C. Jurs (September 1975). "Automated chemical structure analysis of organic molecules with a molecular structure generator and pattern recognition techniques". Analytical Chemistry. 47 (11): 1829–1835. doi:10.1021/AC60361A007. ISSN 0003-2700. Wikidata Q99232471.
10. Shin-ichi Sasaki; Hidetsugu Abe; Yuji Hirota; Yoshiaki Ishida; Yoshihiro Kudo; Shukichi Ochiai; Keiji Saito; Tohru Yamasaki (1 November 1978). "CHEMICS-F: A Computer Program System for Structure Elucidation of Organic Compounds". Journal of Chemical Information and Computer Sciences. 18 (4): 211–222. doi:10.1021/CI60016A007. ISSN 1520-5142. Wikidata Q99233202.
11. Craig A. Shelley; Morton E. Munk (November 1981). "Case, a computer model of the structure elucidation process". Analytica Chimica Acta. 133 (4): 507–516. doi:10.1016/S0003-2670(01)95416-9. ISSN 0003-2670. Wikidata Q99233261.
12. Martin Badertscher; Andrew Korytko; Klaus-Peter Schulz; et al. (May 2000). "Assemble 2.0: a structure generator". Chemometrics and Intelligent Laboratory Systems. 51 (1): 73–79. doi:10.1016/S0169-7439(00)00056-3. ISSN 0169-7439. Wikidata Q99233839.
13. Raymond E. Carhart; Dennis H. Smith; Neil A. B. Gray; James G. Nourse; Carl Djerassi (April 1981). "Applications of artificial intelligence for chemical inference. 37. GENOA: a computer program for structure elucidation utilizing overlapping and alternative substructures". The Journal of Organic Chemistry. 46 (8): 1708–1718. doi:10.1021/JO00321A037. ISSN 0022-3263. Wikidata Q99233344.
14. H.J. Luinge; J.H. Van Der Maas (June 1990). "AEGIS, an algorithm for the exhaustive generation of irredundant structures". Chemometrics and Intelligent Laboratory Systems. 8 (2): 157–165. doi:10.1016/0169-7439(90)80131-O. ISSN 0169-7439. Wikidata Q99233812.
15. Christoph Steinbeck (20 September 1996). "LUCY—A Program for Structure Elucidation from NMR Correlation Experiments". Angewandte Chemie International Edition. 35 (17): 1984–1986. doi:10.1002/ANIE.199619841. ISSN 1433-7851. Wikidata Q50368945.
16. Christoph Steinbeck (November 2001). "SENECA:  A Platform-Independent, Distributed, and Parallel System for Computer-Assisted Structure Elucidation in Organic Chemistry". Journal of Chemical Information and Computer Sciences. 41 (6): 1500–1507. doi:10.1021/CI000407N. ISSN 1520-5142. PMID 11749575. Wikidata Q28837910.
17. Jean-Loup Faulon (January 1996). "Stochastic Generator of Chemical Structure. 2. Using Simulated Annealing To Search the Space of Constitutional Isomers". Journal of Chemical Information and Modeling. 36 (4): 731–740. doi:10.1021/CI950179A. ISSN 1549-9596. Wikidata Q28837961.
18. Jean-Marc Nuzillard; Massiot Georges (January 1991). "Logic for structure determination". Tetrahedron. 47 (22): 3655–3664. doi:10.1016/S0040-4020(01)80878-4. ISSN 0040-4020. Wikidata Q57818172.
19. K. A. Blinov; M. E. Elyashberg; S. G. Molodtsov; A. J. Williams; E. R. Martirosian (1 April 2001). "An expert system for automated structure elucidation utilizing 1H-1H, 13C-1H and 15N-1H 2D NMR correlations". Fresenius' Journal of Analytical Chemistry. 369 (7–8): 709–714. doi:10.1007/S002160100757. ISSN 0937-0633. PMID 11371077. Wikidata Q43616194.
20. Jochen Junker (28 July 2011). "Theoretical NMR correlations based Structure Discussion". Journal of Cheminformatics. 3 (1): 27. doi:10.1186/1758-2946-3-27. ISSN 1758-2946. PMC 3162559. PMID 21797997. Wikidata Q38264559.
21. Chang-Yu Hu; Lu Xu (November 1994). "Principles for structure generation of organic isomers from molecular formula". Analytica Chimica Acta. 298 (1): 75–85. doi:10.1016/0003-2670(94)90044-2. ISSN 0003-2670. Wikidata Q99233968.
22. Junfeng Hao; Lu Xu; Changyu Hu (October 2000). "Expert system for elucidation of structures of organic compounds (ESESOC)". Science in China. Series B: Chemistry. 43 (5): 503–515. doi:10.1007/BF02969496. ISSN 1006-9291. Wikidata Q105032775.
23. Jean-Loup Faulon (1 September 1994). "Stochastic Generator of Chemical Structure. 1. Application to the Structure Elucidation of Large Molecules". Journal of Chemical Information and Modeling. 34 (5): 1204–1218. doi:10.1021/CI00021A031. ISSN 1549-9596. Wikidata Q99233862.
24. Jean-Loup Faulon; Carla J Churchwell; Donald P Visco (1 May 2003). "The signature molecular descriptor. 2. Enumerating molecules from their extended valence sequences". Journal of Chemical Information and Computer Sciences. 43 (3): 721–734. doi:10.1021/CI020346O. ISSN 1520-5142. PMID 12767130. Wikidata Q52016182.
25. Julio E Peironcely; Miguel Rojas-Chertó; Davide Fichera; Theo Reijmers; Leon Coulier; Jean-Loup Faulon; Thomas Hankemeier (17 September 2012). "OMG: Open Molecule Generator". Journal of Cheminformatics. 4 (1): 21. doi:10.1186/1758-2946-4-21. ISSN 1758-2946. PMC 3558358. PMID 22985496. Wikidata Q27499209.
26. Jean Loup Faulon (1 July 1992). "On using graph-equivalent classes for the structure elucidation of large molecules". Journal of Chemical Information and Modeling. 32 (4): 338–348. doi:10.1021/CI00008A013. ISSN 1549-9596. Wikidata Q99233853.
27. Brendan D. McKay; Adolfo Piperno (January 2014). "Practical graph isomorphism, II". Journal of Symbolic Computation. 60: 94–112. doi:10.1016/J.JSC.2013.09.003. ISSN 0747-7171. Wikidata Q99301767.
28. M.A. Yirik, ‘The Benchmark for Structure Generators’, Blogger, 2020.
29. Mohammad Mahdi Jaghoori; Sung-Shik T.Q. Jongmans; Frank de Boer; Julio Peironcely; Jean-Loup Faulon; Theo Reijmers; Thomas Hankemeier (December 2013). "PMG: Multi-core Metabolite Identification". Electronic Notes in Theoretical Computer Science. 299: 53–60. doi:10.1016/J.ENTCS.2013.11.005. ISSN 1571-0661. Wikidata Q105032974.
30. M. S. Molchanova; V. V. Shcherbukhin; N. S. Zefirov (January 1996). "Computer Generation of Molecular Structures by the SMOG Program". Journal of Chemical Information and Modeling. 36 (4): 888–899. doi:10.1021/CI950393Z. ISSN 1549-9596. Wikidata Q99233768.
31. I. P. Bangov; K. D. Kanev (February 1988). "Computer-assisted structure generation from a gross formula: II. Multiple bond unsaturated and cyclic compounds. Employment of fragments". Journal of Mathematical Chemistry. 2 (1): 31–48. doi:10.1007/BF01166467. ISSN 0259-9791. Wikidata Q105033085.
32. Kerber A, Laue R, Meringer M, Varmuza K. MOLGEN-MS: Evaluation of low resolution electron impact mass spectra with MS classification and exhaustive structure generation. Adv Mass Spectrom. 2001;15(939-940):22.
33. Tomoyuki Miyao; Hiromasa Kaneko; Kimito Funatsu (14 June 2016). "Ring system-based chemical graph generation for de novo molecular design". Journal of Computer - Aided Molecular Design. 30 (5): 425–446. doi:10.1007/S10822-016-9916-1. ISSN 0920-654X. PMID 27299746. Wikidata Q50627884.
34. Tomoyuki Miyao; Hiromasa Kaneko; Kimito Funatsu (26 November 2014). "Ring-System-Based Exhaustive Structure Generation for Inverse-QSPR/QSAR". Molecular Informatics. 33 (11–12): 764–778. doi:10.1002/MINF.201400072. ISSN 1868-1743. PMID 27485423. Wikidata Q39092888.
35. Tomoyuki Miyao; Masamoto Arakawa; Kimito Funatsu (1 January 2010). "Exhaustive Structure Generation for Inverse-QSPR/QSAR". Molecular Informatics. 29 (1–2): 111–125. doi:10.1002/MINF.200900038. ISSN 1868-1743. PMID 27463853. Wikidata Q51758769.
36. Baudoin Delépine; Thomas Duigou; Pablo Carbonell; Jean-Loup Faulon (9 December 2017). "RetroPath2.0: A retrosynthesis workflow for metabolic engineers". Metabolic Engineering. 45: 158–170. doi:10.1016/J.YMBEN.2017.12.002. ISSN 1096-7176. PMID 29233745. Wikidata Q47256449.
37. Mathilde Koch; Thomas Duigou; Pablo Carbonell; Jean-Loup Faulon (19 December 2017). "Molecular structures enumeration and virtual screening in the chemical space with RetroPath2.0". Journal of Cheminformatics. 9 (1): 64. doi:10.1186/S13321-017-0252-9. ISSN 1758-2946. PMC 5736515. PMID 29260340. Wikidata Q47199780.
38. Artur Kadurin; Sergey I. Nikolenko; Kuzma Khrabrov; Alex Aliper; Alex Zhavoronkov (13 July 2017). "druGAN: An Advanced Generative Adversarial Autoencoder Model for de Novo Generation of New Molecules with Desired Molecular Properties in Silico". Molecular Pharmaceutics. 14 (9): 3098–3104. doi:10.1021/ACS.MOLPHARMACEUT.7B00346. ISSN 1543-8384. PMID 28703000. Wikidata Q38681438.
39. Thomas Blaschke; Marcus Olivecrona; Ola Engkvist; Jürgen Bajorath; Hongming Chen (13 December 2017). "Application of Generative Autoencoder in de Novo Molecular Design". Molecular Informatics. 37 (1–2): 1700123. arXiv:1711.07839. doi:10.1002/MINF.201700123. ISSN 1868-1743. PMC 5836887. PMID 29235269. Wikidata Q48127458.
40. Christie BD; Munk ME (1 May 1988). "Structure generation by reduction: a new strategy for computer-assisted structure elucidation". Journal of Chemical Information and Computer Sciences. 28 (2): 87–93. doi:10.1021/CI00058A009. ISSN 1520-5142. PMID 3392122. Wikidata Q38594392.
41. Robert C Glem; Andreas Bender; Catrin H Arnby; Lars Carlsson; Scott Boyer; James Smith (1 March 2006). "Circular fingerprints: flexible molecular descriptors with applications from physical chemistry to ADME". IDrugs: the Investigational Drugs Journal. 9 (3): 199–204. ISSN 1369-7056. PMID 16523386. Wikidata Q51947334.
42. Jean-Loup Faulon; Michael J Collins; Robert D Carr (1 March 2004). "The signature molecular descriptor. 4. Canonizing molecules using extended valence sequences". Journal of Chemical Information and Modeling. 44 (2): 427–436. doi:10.1021/CI0341823. ISSN 1549-9596. PMID 15032522. Wikidata Q45023689.
43. Simona Bohanec (1 May 1995). "Structure Generation by the Combination of Structure Reduction and Structure Assembly". Journal of Chemical Information and Modeling. 35 (3): 494–503. doi:10.1021/CI00025A017. ISSN 1549-9596. Wikidata Q99233866.
44. Korytko A; Schulz KP; Madison MS; Munk ME (1 September 2003). "HOUDINI: a new approach to computer-based structure generation". Journal of Chemical Information and Modeling. 43 (5): 1434–1446. doi:10.1021/CI034057R. ISSN 1549-9596. PMID 14502476. Wikidata Q52009004.
45. G. Massiot; J. M. Nuzillard (July 1992). "Computer-assisted elucidation of structures of natural products". Phytochemical Analysis. 3 (4): 153–159. doi:10.1002/PCA.2800030403. ISSN 0958-0344. Wikidata Q57818162.
46. Donald L. Kreher; Douglas R. Stinson (March 1999). "Combinatorial algorithms: generation, enumeration, and search". ACM SIGACT News. 30 (1): 33–35. doi:10.1145/309739.309744. ISSN 0163-5700. Wikidata Q105033277.
47. Stephen R. Heller; Alan McNaught; Igor Pletnev; Stephen Stein; Dmitrii Tchekhovskoi (2015). "InChI, the IUPAC International Chemical Identifier". Journal of Cheminformatics. 7 (1): 23. doi:10.1186/S13321-015-0068-4. ISSN 1758-2946. PMC 4486400. PMID 26136848. Wikidata Q21146620.
48. Hesam Dashti; William M Westler; John L Markley; Hamid R Eghbalnia (23 May 2017). "Unique identifiers for small molecules enable rigorous labeling of their atoms". Scientific Data. 4: 170073. doi:10.1038/SDATA.2017.73. ISSN 2052-4463. PMC 5441290. PMID 28534867. Wikidata Q33718167.