Virtual libraries of tetrapyrrole macrocycles. Combinatorics, isomers, product distributions, and data mining

A software program (PorphyrinViLiGe) has been developed to enumerate the type and relative amounts of substituted tetrapyrrole macrocycles in a virtual library formed by one of four different classes of reactions. The classes include (1) 4-fold reaction of n disubstituted heterocycles (e.g., pyrroles or diiminoisoindolines) to form β-substituted porphyrins, β-substituted tetraazaporphyrins, or α- or β-substituted phthalocyanines; (2) combination of m aminoketones and n diones to form m × n pyrroles, which upon 4-fold reaction give β-substituted porphyrins; (3) derivatization of an 8-point tetrapyrrole scaffold with n reagents, and (4) 4-fold reaction of n aldehydes and pyrrole to form meso-substituted porphyrins. The program accommodates variable ratios of reactants, reversible or irreversible reaction (reaction classes 1 and 2), and degenerate modes of formation. Pólya's theorem (for enumeration of cyclic entities) has also been implemented and provides validation for reaction classes 3 and 4. The output includes the number and identity of distinct reaction-accessible substituent combinations, the number and identity of isomers thereof, and the theoretical mass spectrum. Provisions for data mining enable assessment of the number of products having a chosen pattern of substituents. Examples include derivatization of an octa-substituted phthalocyanine with eight reagents to afford a library of 2,099,728 members (yet only 6435 distinct substituent combinations) and reversible reaction of six distinct disubstituted pyrroles to afford 2649 members (yet only 126 distinct substituent combinations). In general, libraries of substituted tetrapyrrole macrocycles occupy a synthetically accessible region of chemical space that is rich in isomers (>99% or 95% for the two examples, respectively).