An algorithm to identify and count coplanar isomeric molecules formed

of complexities of intermediates as a comparative measure for the efficiency of synthesis.1 In comparisons of syntheses for the same target, however, ...
3 downloads 0 Views 296KB Size
J. Chem. ZnJ: Comput. Sci. 1987, 27,67-69

for the two orders are shown in Figure 5 (functional groups and the minor symmetry terms are not included). The results show a clear preference for the convergent order, as expected from other considerations.6 This procedure points up another interesting conclusion about synthesis. In a real convergent synthesis with added refunctionalization reactions, the total complexity of intermediates will be much less if these added reactions precede the final coupling of the two skeletal halves. The conclusion for synthesis planning is clear: the final coupling of intermediates should come near the end rather than the beginning of a convergent synthetic route. In more general terms, any step that exhibits a large increase in complexity is more efficiently positioned near the end rather than the beginning of a synthetic sequence. The calculation of molecular complexity is rendered much easier to carry out by hand when the variations in eq 1, 3, 5, and 7 are used and the values from Tables I and I1 are applied. The method presented here, used by hand or computer, yields the same complexity values as the Bertz method. Also, it is less liable to error than identifying and counting bond connectivities, when used by hand, and much more amenable to computerization in this way as well. The program CPXCAL~ is available to anyone with an interest in comparing molecular complexities in any molecular families of interest.

llnnnr 2

i”

P‘ 2, eZL&J n - 4

2

4

n-6

4

16

n-8

6

31

n-1P

8

48

10

66

n n

-

P

2

2

n - 4

86

12

Z

-

I:

251

.Intermediates .are l i n e a r skeletons of n atoms

(I)

-

n

-

67

-

ACKNOWLEDGMENT

78

2)

We are grateful to the National Science Foundation for a grant (CHE-9 102972) that has generously supported this work.

Figure 5. Complexity comparison of linear and convergent syntheses.

lessened, Le., QT lg 1,with 2 > Q > 1. The application of complexity calculations to assess synthesis efficiency has been illustrated by Bertz, who has used the sum of complexities of intermediates as a comparative measure for the efficiency of synthesis.’ In comparisons of syntheses for the same target, however, this will generally simply prefer the synthesis with the fewest steps, hence fewest intermediates to sum. However, in a comparison of two syntheses with the same starting fragments and the same number of steps, this measure should show a preference for one. This amounts to two different orders of assembly of the same units and is characteristic of the difference between a convergent and a linear synthesis? If we apply this method to the assembly of a linear target skeleton of 16 atoms from 8 2-carbon starting units, the results

REFERENCES AND NOTES (1) (a) Bertz, S. H. Chem. Commun. 1981, 818. (b) Bertz, S. H. J . Am. Chem. Sac. 1981, 103, 3599; 1982, 104, 5801. (c) Bertz, S. H. Bull. Math. B i d . 1983,45,849. (d) Bertz, S. H. In Chemical Applications of Topologv and Graph Theory; King, R. B., Ed.; Elsevier: New York, 1983; p 206. (2) Hendrickson, J. B.; Grier, D. L.; Toczko, A. G. J. Am. Chem. SOC. 1985, 107, 5228. (3) The three-atom skeleton extracted may be of any three linked, non-

hydrogen atoms, Le., C C - C , C-0-C, C-N-S, etc.

(4) In eq 4, h is the number of attached hydrogens plus the number of

unshared electron pairs, regarded as the conjugate bases of potentially attached hydrogens. (5) Hendrickson, J. B.; Toczko, A. G.J . Chem. l n j Comput. Sci. 1983,

23, 171. (6) Hendrickson, J. B. J . Am. Chem. Sac. 1977, 99, 5439. (7) CPXCAL is written for a DEC VAX computer with an input module for

direct graphic input from a Tektronix graphics terminal.

An Algorithm To Identify and Count Coplanar Isomeric Molecules Formed by the Linear Fusion of Cyclopentane Modules SEYMOUR B. ELK Elk Technical Associates, New Milford, New Jersey 07646 Received September 18, 1985 Because each of the various possible isomers formed by the coplanar linear fusion of cyclopentane modules may be represented by a binary sequence, the reverse technique of examining each binary sequence of a specified length underlies the formation of an algoithm to identify and count such isomeric molecules. T h e algorithm involves the specification of a set of three binary operations-which correspond to allowable physical transformations. This may be expressed in the form of a formal algebraic table of operations. Application of this algorithm, with Patterson’s drawing convention, produces a canonical representation for each such isomer.

Despite the presence of a certain amount of noncoplanarity in the cyclopentane molecule (which is caused by the relieving of the strain that would result if the hydrogen atoms were

allowed to remain eclipsed),’ the simplified geometrical skeleton model formed by the successive “straight line”* concatenation of regular pentagonal modules gives a fairly good

0095-2338/87/1627-0067$01.50/00 1987 American Chemical Society

ELK

68 J . Chem. In$ Comput. Sci., Vol. 27, No. 2, 1987 1

2

*w* 1

7

3

7

Y

5

Figure 4.

...

... =

Figure 1.

...oJ

Figure 5.

Figure 6.

Figure 2. Figure 7.

Figure 3.

description of an interesting class of fused pentagonal rings. The first two members of this set (Figure 1) are bicyclo[ 3.3 .O]octane tricyclo [6.3 .0.02,6] undecane The description of this class of molecules is e n h a n d through the use of the drawing convention prescribed by Patterson3 for the pentagonal module-inasmuch as such a convention assigns a unique orientation that is readily machine codable. In fact, with this convention, the idea of “straight” is fulfilled by an alternation of pentagons in the “up” vs. “down” position, whereas the intuitive alignment (represented by a straight line in the drawing along the top or bottom of a connected set of pentagons) produces a molecule that is, in actuality, circular (Figure 2). The representation of polyfused pentagonal rings using Patterson’s format yields a sequence of rings pointing up or down-which may be machine coded as “0” or “l”, respectively. The algorithm to manipulate each string of 0’s and 1’s is based on the following three principles: (A) The two end digits of each binary number are immaterial; Le., 0’s and 1’s may be interchanged at will. From a chemistry perspective, this is equivalent to saying that there is a plasticity to each end ring (Figure 3); Le., an end ring may be represented in either format without changing the incidence relations (vertices 1 and 5 have incidence = 3, while vertices 2-4 have incidence = 2). (B) Each binary number is equal to its 1’s complement. From a chemistry perspective, this is equivalent to a horizontal reflection (Figure 4). (C) Each binary number is equal to its reverse-order number. From a chemistry perspective, this is equivalent to a vertical reflection (Figure 5). From the vantage point of a chemist, all three operations must be examined inasmuch as this is the standard method for drawing such molecules. However, from the viewpoint of an information scientist, Ryon4has demonstrated that the final form of the algorithm may be simplified by noting that an arc

Figure 8.

through each sequence of three rings will point in the same direction as the center ring. In other words, for the five-ring sequence shown in Figure 2, instead of having to examine 25 = 32 numbers and then ignoring the two end rings (Le., the two end digits), by considering the direction of bend of the arc through each sequence of three rings there is the need to examine only 23 = 8 binary numbers-rather than 32-in order to find the set of all of the different pentacycloheptadecane isomers. Even with this set of eight numbers, however, there is still multiple redundancy; namely, there are only three distinct isomers (Figures 6-8). In other words, the final form of the algorithm will need to consider the set of r - 2 rings (without the end rings) and the use of operations B and C only. Furthermore, we may note that the application of first B and then C is equal to applying first C and then B; also, both are geometrically equivalent to rotating the configuration through 180’. This latter transformation may be represented by the combined operation BC. Also, applying either B or C twice gives the original configuration. This is expressed algebraically as BB = CC = I, where I is the identity operation ( M= XI = X , for every X). Since geometrically there are three operations (two reflections and one rotation), there exists a corresponding algebraic table of operations-containing four rows (and columns), one for each operation plus an additional one for the identity operation-which describes the manipulations that may be performed on each sequence of the r - 2 digit^:^

BC C BC

C BC

BC

I

B

C

B

I

ALGORITHM FOR ISOMERS FROM FUSIONOF CYCLOPENTANES

J . Chem. InJ Comput. Sei., Vol. 27, No. 2, 1987 69

Table I" comment tricyclo

0

tetracyclo

00 = 0 01 = 1

pentacyclo

000 = 0 010 = 2 001 = 1 011 Bc-001 = 1 100 through 11 1 % 3 through 0, respectively

exactly three isomers (0,1, and 2) for five rings-see

0000 = 0 0 1 0 0 2 0010 = 2 0001 = 1 0101 = 5 0010 = 2 0110 = 6 0011 = 3 O l l l ~ o o 0 1= 1 1000 through 11 11 % 7 through 0, respectively

exactly six isomers (0,1, 2, 3, 5 , and 6) for six rings-see

hexacyclo

1 5 0

exactly one isomer for three rings-see

1 0 % 01 = 1 11 JLOO = 0

exactly two isomers (0 and 1) for four rings-see

heptacyclo canonical 0, 1, 2, 3, 4, 5, 6, 9, 10, and 14 C transformed 8 (=2) and 12 (=6) BC transformed 7 (=3), 11 ( = 5 ) , 13 (=9), and I5 (=1) exactly ten isomers for seven rings "Note that the second half of each set of numbers is always equivalent to the first half in the reverse order-by omitted. as it cannot describe a new isomer.

v

3

Figure 9 Figure 10

Figure 11

operation B-and

Figure 12

thus may be

.

5

tricyclo [6.3 .0.02*6] undecane

8

10

Figure 9.

tetracycl0[6.6.0.O~~~,O~J~]tetradecane tetracycl0[6.6.0.0~~~.O'~*'~]tetradecane

Figure 11. xoooOx = hexacyclo[9.9.0.02*9.03~'.013~20.01s~19]ic~~ane

Figure 10.

xOOOl x = hexacyclo[9.9.0.02~9.03~7.0'3~20.0'4~L8]icosane xOOlOx = hexacyclo[9.9.0.@~9.03~7.0'2~'9.0~4~18]icosane

In order to describe each of the different ways for forming isomers having r rings, consider the set of binary numbers from 0 to 2r2 - 1. Let n designate the given binary number without the end digits, B(n) be the binary number formed by taking the 1's complement of n, C(n) be the binary number formed by writing n in the reverse order, and BC(n) be the binary number formed by taking the 1's complement of the digits of n in the reverse order. The smallest of the numbers [n,B(n), C(n),and BC(n)J is the canonical representation of each of the different isomers. In practice, however, if we are cataloging the entire set of r rings, rather than just examining a single arbitrary configuration, the process may be greatly simplified. The formal algorithm is as follows: Apply the three operations B(n), C(n), and BC(n) to each of the binary numbers formed starting with a string of r - 2 0's and progressing through r - 2 1's (see Table I),5 as follows: (1) If the first digit of n is a "0", prooeed to step 2; otherwise, apply B(n) to get a smaller equivalent binary representation (i.e., set this new number equal to n) and then proceed to step 2. (2) If the final digit is a "0", apply C(n). The canonical representation is min [n,C(n)]. Otherwise, the final digit must

xOOl1x = hexacyclo[9.9.0.02~9.03~7.012J9.0'3J7]icosane xOlOl x = hexacyclo[9.9.0.02~9.0'~8.0~2~19.014J8]ic~ane x0 11Ox = hexacyclo[9.9.0.02~9.04~8.013~20.0149~8] icosane

Figure 12.

be "1". Apply BC(n). The canonical representation is min [n,BC(n)l. Note that applying any of the other operations will give a larger number than the starting one or else an isomer already included in the set. REFERENCES AND NOTES Greenberg, A.; Liebman, J. F."Strained Organic Molecules"; Academic Press: New York, 1978; p 65. Elk, S.B.'Refinement of Systematic Nomenclature for Polybenzenes and Its Extension to Systems of General Arenes". Moth. Chem. 1982, 13, 261. Patterson, A. M. 'Proposed Intemational Rules for Numbering Organic Ring Systems". J. Am. Chem. Soc. 1925, 47, 547. Ryon, J. W., Professor of Computer Science, New Jersey Institute of Tochnology, personal communication. Figures 9-12 give the ictorial representations and the corresponding IUPAC nomenclature for each of the different tricyclo through hexacyclo compounds. International Union of Pure and Applied Chemistry 'Nomenclature of Organic Chemistry: Section A"; Pergamon Press: Oxford, England, 1979.

f