Quantitative structure-activity relationships in insecticidal oxathiolane

Jan 1, 1987 - Roy T. McBurney and John C. Walton. Journal of the American Chemical Society 2013 135 (19), 7349-7354. Abstract | Full Text HTML | PDF ...
0 downloads 0 Views 780KB Size
J. Agric. Food Chem. 1987, 35, 115-121

115

Quantitative Structure-Activity Relationships in Insecticidal Oxathiolane and Dithiolane Oxime Carbamates and Related Compounds A. Peter Kurtz* and John A. Durden, Jr.

The insecticidal and acaricidal activities and the mammalian toxicity of a series of 28 oxathia and dithia heterocyclic oxime methylcarbamates have been related to certain of their physicochemical properties. Multiple regression studies have shown that fly head acetylcholinesterase inhibition and housefly, buckthorn aphid, two-spotted spider mite, southern armyworm, Mexican bean beetle, and rat toxicities are enhanced by increased steric requirement of ring substituents (Taft E&. Increasing lipophilicity (log P) uniformly diminishes activity against the arthropods and the rat while the presence of a @-sulfur as compared to a @-oxygenis associated with increased activity against all test organisms. These correlation studies suggest that further synthesis in this series is unlikely to yield effective insecticides in which the arthropod/mammal selectivity is improved over that of compounds already in hand.

The syntheses and the insecticidal properties of carbamate esters of 4-oximino-1,3-oxathiolanes (Ia), 3-oximino-1,4-oxathianes (Ib; Kurtz et al., 1987), 4-oximino-1,3dithiolanes (IC),and 3-oximino-1,4-dithianes (Id; D'Silva et al., 1985) have been reported separately (Chart I). These materials are broad-spectrum insecticides having excellent activity against the buckthorn aphid, Mexican bean beetle, southern armyworm, and the housefly. The dithia compounds and selected members of the oxathia series in addition exhibit good acaricidal potency. The more insecticidally active monomethylcarbamates in both series of compounds, however, possess significant levels of toxicity to warm-blooded animals as measured by acute toxicity to the rat. As part of the examination of these compounds we sought to define which structural and physicochemical features of the monomethylcarbamates in the series are significant in controlling toxicity to both insect and mammalian species. With the role of such physicochemicalfeatures defined, the activities of synthesis candidates could be predicted. Should features leading to biological selectivity be identified, structures could be devised to improve the insect/mammalian toxicity ratio. This report is concerned with development of quantitative structureactivity relationships (QSAR)required to reach these goals. Multiple-regression analyses following the approach developed by Hansch (Tute, 1971) are presented. An examination of the data given for I in companion reports (D'Silva et al., 1985; Kurtz et al., 1987) reveals that the N-monomethylcarbamate compounds (I: R1 = CH3, R2 = H) are the most potent toxicants in the series. In order to focus on the role of various structural features in the oxime moiety, the compound data bank was limited to N-monomethylcarbamates in the present study. In Table I are presented chemical structures, dose/response activity data, and physicochemical parameters for the 15 oxathia (compounds 1-15) and 13 dithia (compounds S 164328) heterocyclic oxime methylcarbamates examined in this work. (The prefix S on a compound number designates the material as a member of the dithia series.) Compounds in these series can exist as syn (oximino oxygen cis to the thiohydroximate sulfur) or anti isomers. All data reviewed in the present report are for syn isomers. Methods used to establish configuration are given in the companion reports (DSilva et al., 1985; Kurtz et al., 1987). Research and Development Department, Union Carbide Agricultural Products Company, Research Triangle Park, North Carolina 27709. 002 1-85611871 1435-0 115$0 1.5010

Chart I

f:

/R'

R4Rw-0-c-N\

R2

-x -n la

0

1

Ib

0

2

IC

s

1

Id

S

2

ESTIMATION/DETERMINATIONOF PHYSICOCHEMICAL PARAMETERS

1. Octanol/Water log P. log P, a measure of the overall lipophilicity/hydrophilicity in molecules, was determined experimentally according to the procedure of Fujita et al. (1964). The data are given in Table I where experimental determinations are noted by asterisks. The other log P values given in the table were estimated as detailed in the footnotes of Table I. The lipophilic character of an entire molecule is important to its transport through and between cells and membranes. Of equal importance is the relative lipophilicity/hydrophilicity of specific groups within a molecule-i.e., side chains. For instance, a lipophilic side chain may promote binding in a lipophilic "pocket" in the active site of a target enzyme. However, the role of localized side-chain lipophilicity could not be independently studied in this work since ET parameters at the a and y ring positions (defined below) are highly correlated with the corresponding Taft CE, steric parameters for the group of compounds under study here (correlation coefficients 0.79 and 0.80, respectively). Thus, the E E 8 parameters may, in this case, express not only the role of substituent bulk in the biological activity of these compounds but also the role of regional lipophilicity. 2. Taft Steric Parameters (CE,",CEsy). To quantitate the role of the size of substituents in the biological potency of these compounds, the Taft steric parameter, E, (Craig, 1971),of substituents in the position a (adjacent to) the oximino carbon were summed to yield XESawhile E, for the substituents at the position adjacent to the distal heteroatom (@to the oximino group) were summed to give 0 1987 American Chemical Society

116 J. Agric. Food Chem., Vol. 35, No. 1, 1987

Kurtz a n d Durden

Table I. Physicochemical Parameters and Biological Data for 15 Oxathia and 13 Dithia Heterocyclic Oxime Methylcarbarnates U

II “OCNHCH3

xP

no. X I O 2 0 3 0 4 0 5 0 6 0

7 0 8

0

9 10 11 12 13

0 0 0

0 0 14 0 15 0

S16 S S17 S 518

S

S19 S S20 S S21 S S22 S S23 S S24 S S25 S S26 S S27 S S28 S

ring size 5 5 5 5 5 5 5 5 5 5 5 5 5 5 6 6 6 6 6 6 6 6 6 5 5 5 5 5

a subst

n-C3H7, H i-C3H7, H CH30CH2CH2, H CHnOCHI, H

CH,, CH3

y subst

1%

CEsaa ZESTb P 2.48 1.24 2.48 2.48 1.17 0.88 1.24 0.77 1.10 1.05 1.24 1.17 1.24 0.00 2.48 2.48 1.24 0.00 2.48 1.24 0.00 0.88 0.77 2.48 1.24 0.00 2.48 1.24

2.48 2.48 1.24 1.17 2.48 2.48 1.24 2.48 2.48 2.48 1.17 1.24 0.88 2.48 2.48 2.48 2.48 2.48 1.24 1.24 1.24 2.48 2.48 2.48 2.48 2.48 0.00 0.00

0.06* 0.46* 0.56 1.06 0.96 1.52* 0.96* 1.48* 0.51* 0.01 1.46 1.46 1.96 0.82* O.OO* 0.49* 0.89 1.51* 0.99 1.39 2.01 1.89 1.89 0.51* 0.91* 1.31* 1.33* 1.73

aphid

mite

SAW

MBB

HF

25.0 6.0 33.0 45.0 25.0 500.0 15.0 40.0

150.0 340.0 200.0 200.0 180.0 500.0 85.0 500.0 500.0 200.0 10.0 85.0 500.0 40.0 500.0 6.0 6.0 4.0 12.0 30.0 7.0 500.0 250.0 55.0 15.0 2.0 13.0 7.0

150.0 180.0 14.0 14.0 500.0 70.0 500.0 500.0 110.0 25.0 500.0 500.0 25.0 17.0 500.0 60.0 500.0 170.0 inact 10.0 140.0 70.0 140.0 60.0 500.0 160.0 11.0 14.0 100.0 38.0 40.0 52.0 11.0 60.0 140.0 80.0 100.0 35.0 110.0 65.0 50.0 1000.0 500.0 150.0 500.0 150.0 28.0 100.0 10.0 22.0 23.0 16.0 35.0 250.0 21.0 75.0

3.0 5.0 3.0 18.0 3.0 4.0 11.0 9.0 5.0 50.0 60.0 35.0 50.0 16.0 2.0 0.8 4.0 12.0 10.0 13.0 45.0 60.0 27.0 0.8 1.5 4.0 6.0 6.0

inact 44.0 10.0 13.0 65.0 4.0 500.0 5.0 4.0 0.8 8.0 3.0 1.5 75.0 40.0 6.0 5.0 6.0 2.0 3.0

zapmg/kg 2.50 0.80 1.00 0.40 0.50 0.05 0.80

49.2 11.7 12.3 59.5 23.6 83.0 4.5 16.2 8.0 15.0

1.50 1.00 0.07 0.10 4.00 0.15 0.05 0.20 3.50 3.00 2.00 0.25 0.20

1.6 37.0 15.4 2.2 0.27 5.10

13.0 1.0 0.7 2.0 1.0

MW 176 190 190 204 204 218 204 218 234 220 218 218 232 204 190 206 220 234 220 234 248 248 248 192 206 220 220 234

C of Taft steric parameters (Craig, 1971) for Substituents at the a-position. The E, of methoxyethyl was taken as -0.14, twice that of of Taft steric parameters for ethyl (-0.07); the E, of methoxymethyl, -0.19, was taken as half that of ethoxymethyl (-0.37; Talvik, 1971). substituents at the y-position. ‘log of octanol/water partition coefficient determined spectrophotometrically according to the method of Fujita et al. (1964) or calculated from log P of nearest analogue and K values derived from experiment or the literature (Craig, 1971). Values marked with an asterisk were determined experimentally; others were calculated from data for the nearest analogue as follows (compound, calculation): 3,log P (1)+ 0.50; 4,log P (3)+ 0.50; 5, log P (2)+ 0.50; 10,log P (9)- 0.50; 11, log P (2)+ 1.00; 12,log P (2) + 1.00; 13,log P (11)+ 0.50; 517,log P (S16)+ 0.40 (A between 1 and 2);S19,log P (516)+ 0.50;520, log P (516)+ 0.50 + 0.40; 521,log P (Sl8)+ 0.50; 522,log P (517)+ 1.00; 523,log P (517)+ 1.00; 528,log P (527)+ 0.40. 0.50 and 1.00 are the literature ?r values for, respectively, CH3 and CH3CHzincrements. dBean plants sprayed to runoff infested with buckthorn aphid, two-spotted spider mite, southern armyworm larvae, or Mexican bean beetle larvae. Housefly test: cotton fibers impregnated with solution of appropriate concentration (Payne et al., 1966). ‘Determined with Warburg apparatus (Moorefield and Tefft, 1958). ’Intubation in corn oil; acute oral LD, determinations were carried out by C. S. Carpenter and colleagues at the Carnegie-Mellon Institute of Research, Carnegie Hygiene Fellowship, Pittsburgh, PA.

CE,?. Data and estimation details are given in the footnotes to Table I. 3. Indicator Variable for X = S . A major difference between the dithia and oxathia compounds in this series is the presence or absence of an oxidizable sulfur atom. An indicator variable, Ia-s was created to account for this difference in correlation studies. For oxathia compounds, was set equal to 0; for dithia compounds, 1.0. 4. Ring Size. Table I contains compounds having either five- or six-membered rings. Preliminary regression work indicated no significance to an indicator variable designed to examine the role of ring size on activity. BIOLOGICAL DATA

The methods of synthesis and biological evaluation of compounds reviewed in this study are described in companion papers (D’Silva et al., 1985; Kurtz et al., 1987) and references cited therein. Toxicity data (LC,o, ppm) from foliage spray tests against buckthorn aphid adults, two-spotted spider mite (all stages), southern armyworm larvae, and Mexican bean

beetle larvae and from a bait test involving the adult housefly (Payne et al., 1966) are presented in Table I. Mammalian toxicity data from acute oral tests in the rat (LD,,, mg/kg, stomach intubation suspended in corn oil) and the molar Imfor fly head acetylcholinesterase inhibition (Moorefield and Tefft, 1958) are also given in Table I. All data are from a single series of dose/response experiments. In general, the 95% confidence interval for individual LCm’sin these tests is LCm/2 to LCw X 2 (e.g. the CI for an LCm of 50 is 25-100.) In foliage spray tests the arthropods are exposed to the toxicant by both contact and feeding routes as the animals forage on treated leaves. For the regression trials, biological activity is expressed on a log scale in which larger numbers indicate greater activity as follows (MW = molecular weight): aphid, mite, SAW, MBB, HF, rat I507

log (1/C) = -log (LC,/MW) log (1/C) = -log Im

MULTIPLE REGRESSION ANALYSIS

The data presented in Table I were analyzed by standard multiple linear regression and other statistics routines

J. Agric. FoodChem., Voi. 35, No. 1, 1987

SAR in Oxathiolane and Dithioiane Oxime Carbamates

Table 11. Range of Data for the Mixed Data Set of 28 Comnounds variable mean S min max range

XES" XES7 log P [(log P - log P)']"

MW 18-S aphidb miteb

SAWb MBBb HFb IwC ratd

1.28 1.89 1.11 0.30 214.9 0.50 1.15 0.68 0.26 0.68 1.44 0.30 1.48

0.90 0.0 0.80 0.0 0.56 0.0 0.33 0.0009 18.78 176.0 0.51 0.0 -0.42 0.77 0.84 -0.42 -0.61 0.53 0.48 -0.39 0.56 0.51 0.63 -0.60 0.70 0.42

2.48 2.48 2.01 1.188 248.0 1.0 2.47 2.04 1.27 1.44 2.41 1.30 2.89

2.48 2.48 2.01 1.187 72.0 1.0 2.89 2.89 1.88 1.83 1.85 1.90 1.47

The log P used in the regression analysis was 1.09 for the mixed and oxathia seta; log P = 1.31 was used for the dithia set. "Variable: -log (LC,/MW). 'Variable: -log Im Range of data given for Zm is for the set of 20 mixed compounds. See Table IV for a listing of the exclusions. dRange of data given for the rat is for the set of 21 mixed compounds. See Table IV for a listing of the exclusions.

on an Il3M-370 M V S batch system. The means and ranges covered by the physicochemical parameters that serve as independent variables in these regression studies are presented in Table 11. A percent covariance matrix (r2) reviewing one-on-one intercorrelation between all variables is given in Table 111. For the -mixed" data set of 28 compounds, the range in log P (Table 11)is limited (range 2.01); in this situation, log P is highly correlated with its square (P = 0.96), making specious the examination of (log P)2 terms in regression trials. To break this internal correlation, the mean log P (log P = 1.09) for the oxathia data set was subtracted from each log P prior to squaring for regression studies -of the oxathia and mixed data sets; the mean log P (log P = 1.31) for the dithia data set was subtracted prior to squaring for studies on the dithia wries. Raw regressions in the form [log (l/C) = a + b(1og P) + c(1og P - logP)2 + d(CE,*) + e(CE,Y) + f(l,+s)] were algebraically expanded to the form [log (l/C) = a + b(1og P ) + c(1og P ) 2 + d(C,*) + e(CE,Y) + f(I,+s)]before extracting the coefficient data recorded for each equation in Table IV. In Table IV are presented data for equations correlating activity with physicochemical effects for the four insects, the mite, the rat, and fly head AChE for the oxathia set (maximum of 15 compounds), the dithia set (maximum of 13 compounds), and the mixed set (maximum of 28 com-

117

pounds). Where n in Table IV is less than 28 for the mixed data set, less than 15 for the oxathia data set, or less than 13 for the dithia data set, specific compounds are missing from the regression trial either because data were missing or compounds proved to be serious outliers in preliminary regression runs. Details on outlier exclusions are clearly noted in Table IV. In no case were more than two compounds removed from regression data sets because of lack of fit (see Table IV); where required, such exclusions were made to assure that equations used to predict activity for new structures had coefficients as reliable as possible. The sense and size of equation coefficients were not significantly affected in any case by the noted removals of outliers. RESULTS AND DISCUSSION

Intervariable Correlations. The correlations presented in Table I11 provide insights into which parameters exert discernible effects on biological activity when examined one-on-one and also provide an opportunity to examine the structural and physicochemical parameters for covariance, a property to be avoided where possible in the independent variables in multiple-regression work. The following generalizations can be made from the data in Table 111: 1. There is little one-on-one correlation between Im (a measure of intrinsic activity) and the activity of the compounds studied against the four insects, the mite, or the rat. This suggests that factors other than activity at the target enzyme are important, e.g., efficiency of transport in vivo or susceptibility to detoxifying metabolism. 2. Rat toxicity is highly correlated with activity against the aphid and the mite. 3. Mite and aphid activities are highly correlated as are bean beetle and armyworm toxicities. 4. Housefly activity is poorly correlated with activity against any of the other species. 5. XES"may be important to both Im and rat toxicity. 6. 150correlates somewhat with log P. It is important to note, however, that the more detailed multivariate analyses presented in Table IV clearly show the XESaand CE,Y play an important role in activity against both the enzyme and in the activities in most of the whole-animal tests. The role of CE,* and CE,Y in enzyme activity would not be apparent in the correlation matrix since that matrix only evaluates one-on-onecorrelations. A special regression run showed that log P is covariant with ( X E S a CE,Y) for the 150data set. For this data set, log P or X E S a and CE,Y express in slightly different ways that ring substituents promote enzyme inhibition. More detailed analysis

+

Table 111. Percent Covariance (r2X 100)

aphid mite SAW MBB HF I50

rat

aphid

mite

+100.0 +55.1 +7.9 +17.5 -0.3 +14.6b +65.4b -6.2 -8.8 +7.3 -13.0 +31.8

+55.1 +100.0 +15.5 +29.8 +2.7 +15.7b +62.4b -9.0 -1.9 +1.2 -15.4 +45.0

SAW +7.9 +15.5 +100.0 +56.6 +12.8 O.Ob

+28.4b -5.9 +7.0 -19.6 -12.2 +0.5

MBB

HF

+17.5 +29.8 +56.6 +100.0 +12.0 +2.7b +38.gb -9.0 +1.1 -7.4 -10.0 +8.9

-0.3 +2.7 +12.8 +12.0 +100.0 -19.5b +0.8b +13.2 +12.5 -44.2 -0.8 +4.0

I50

+14.6b +15.7b +O.Ob

+2.7b -19.5b +100.0 +21.2b -38.1 -23.3 +43.2d +25.3 -11.6

rat

XE,*

+65.4b +62.4b +28.4b +38.gb +0.8b +21.2b +100.0 -31.3 -2.5 +0.18 -16.3 +43.7

-6.2 -9.0 -5.9 -9.0 +13.2 -38.1 -31.3 +100.0 -8.0 -20.6 +3.0 -0.9

EE*7c -8.8 -1.9 +7.0 +1.1 +12.5 -23.3 -2.56 -8.0 +100.0 -13.8 +1.3 -0.2

log P

(loa .- P log P)2

ID-S

+7.3 +1.2 -19.6 -7.4 -44.2 +43.2d +0.18 -20.6 -13.8 +100.0 -1.1 +8.9

-13.0 -15.4 -12.2 -10.0 -0.8 +25.3 -16.3 +3.0 +1.3 -1.1 +100.0 -1.3

+31.8 +45.0 +0.5 +8.9 +4.0 -11.6 +43.7 -0.9 -0.2 +8.9 -1.3 +100.0

"Except where otherwise noted, the number of compounds involved in each correlation is as described in the notes to Table II; signs noted are for r , the correlation coefficient. bNumber of compounds limited to those where all biological variables were available: n = 20 for Iw/insect intercorrelation; n = 16 for rat/Z, intercorrelation; rat/insect, n = 21. (-) sign indicates biological potency increases with increasingly negative XES,i.e. greater steric bulk. dPossibly a function of overall molecular size. r2 for log P vs. MW was 0.841.

118

J. Agric. Food Chem., Vol. 35, No. 1, 1987

Kurtz and Durden

I lk? 00

22 m

cu 0

Y

."6

g 0:

0

A

6

w

m

Y

9 B 0 Y

+ I

n 3 Ll

ET 3 m m

SAR in

Oxathiolane and Dithiolane Oxime Carbamates

detailed below indicates that higher levels of lipophilicity actually detract from activity in the whole-animal tests. 7. Housefly activity is inversely correlated with log P. The negative role of log P in activity in essentially all the whole-animal tests developed below is most significant in the housefly. It appears substantial in the correlation matrix only for housefly because only with housefly does log P alone explain a high percentage of variation in the biological variable. 8. Much of the variation in biological activity for the aphid, mite, and the rat in the mixed data sets is accounted for by the presence or absence of a @sulfur atom, i.e. an oxidizable sulfur. 9. Notwithstanding the log [P/(CE,“+ CE,T)] functional relationship noted above, there are no serious oneon-one intercorrelations of the physicochemical independent variables. Regression Equations (Table IV). Role of Ring Substitution in Acetylcholinesterase Inhibition. One of the more fundamental aspects of the toxicology of oxime methylcarbamates is their activity as inhibitors of acetylcholinesterase, the target enzyme. Molar fly head acetylcholinesterase 150data (molar concentrations giving 50% inhibition of the enzyme) are a measure of relative activity in this regard (Moorefield and Tefft, 1958). Data were available for 8 oxathia and 12 dithia compounds. Despite the relatively low r2 values that were achieved, eq 2 and 3 for the oxathia and dithia sets demonstrate the importance to enzyme inhibition of steric bulk or branching in both the CY- and y-positions of these compounds. Since groups with greater steric requirements have increasingly negative E, values, the coefficients of these E, terms suggest that potency should increase with increasing bulk. The similarity of coefficients for eq 2 and 3 prompted examination of a combined data set; eq 1was the result. The equation could not be improved by forcing in log P terms or the 18-sparameter. Compounds 517, 519, and 523 are major outliers in eq 1. A rationale for their lack of fit is not readily apparent. It is worth noting that their removal from the data set gave an equation [-log 150= 1.96 - 0.4CEa“ (*0.08) - O.48CEsY (f0.15); n = 17, s = 0.28, r2 = 0.781 essentially identical with eq 1in which nearly 80% of the variation in data for 17 compounds is explained by the two steric terms. The critical role of the E, terms to Im found in all of eq 1-3 extends into the equations for the insects, mite, and rat as described below. Because of the high correlation between log P and (CE; + ZEST)noted earlier, the present equations do not permit a distinction between lipophilicity and steric bulk as operative factors in the enzyme inhibition effectiveness of these materials. The insensitivity of the Im’s to I p s is not surprising since sulfoxidation, which occurs in vivo, should not occur during the in vitro test. The advantageous role of branching and the size of groups a t the CY- or y-positions in these compounds is reminiscent of the branched quaternary group in the natural substrate acetylcholine (Payne et al., 1966). Indeed, Cohen et al. (1985) have assembled good evidence that the “anionic site” in acetylcholinesterase might better be considered a “branched alkyl” site: acetylcholine mimics having branching but no quaternary center bind only 1.5 kcal/mol more weakly than quaternary analogues. Role of the Size of Substituents, Lipophilicity, and the Presence of a n Oxidizable Sulfur in Insecticidal, Acaricidal, and Mammalian Toxicity. Regression Equations (Table I V ) . Increased steric requirement (branching) at

J. Agric. Food Chem., Vol. 35, No. 1, 1987

119

the a- and/or y-positions promotes activity in all three series in every case except the housefly equations where variation in potency is too limited to support both E, and log P terms: the log P terms in the housefly equations probably “include” a favorable steric component that is masked by an unfavorable role of the lipophilicity of substituents as discussed below. Two new factors, however, show up in the QSAR for the insects, the mite, and the rat, all in vivo tests, which were not present in the 150equations: 1. As log P increases above about 0.25, potency decreases. 2. Potency increases with replacement of the oxygen in the oxathia series with a readily oxidizable sulfur atom. The coefficients for the log P and (log P ) 2terms given in Table IV yield three basic patterns: 1. In eight equations (6, 11, 13, 15, 16, 18-20) activity decreases linearly over the range of log P studied (0.0-2.01) according to a function having a slope averaging -0.97 (range of slopes -0.79 to -1.23). A steeper linear function (slope -2.04) was noted for the dithia mite eq 9 on the basis of the remarkable inactivity of compounds 522 and 523; a much more gradual slope (-0.56) was observed in the oxathia housefly equation (17). 2. The Mexican bean beetle oxathia equation (14) and the rat dithia equation (21) show no‘dependence on log P. In the case of Mexican bean beetle, there is insufficient variation in biological data to support a strong regression analysis. The dithia rat equation shows no overt role of log P because any possible real effect of log P is hidden in the CE, terms (log P is a linear function of the X E S terms as noted a t the bottom of Table IV). 3. For all other cases (eq 4,5, 7, 10, and 12), log P and (log P)2terms described parabolas over the range of data in log P present in these data sets. The shapes of these parabolas are fairly uniform, with log P optima occurring between -0.11 and +0.76. Data sets having a point a t which a dramatic loss in activity occurred at high log P values generally resulted in a parabolic log P function. Two curves are plotted in Figure 1: one is a line for the log P function of eq 19 having a slope of -0.79 representing the family of equations having linear log P functions; the second is a parabola for the log P function of eq 4 having optima at 0.33 and representing the family of equations having parabolic log P functions. It is important to note that, within the range of data studied, the parabolic function follows the linear function closely, certainly within the precision of the data and the standard errors of the regression equations. The general conclusion can be drawn that, for the biological activity parameters examined in this study, activity decreases monotonously with increase in log P above 0.5 with a slope approximating -1.0. Further definition within this generalization requires more data, particularly data for compounds having log P substantially below zero. Finally, the la.s terms in each mixed-set equation clearly express where and to what extent activity is promoted by the presence of an oxidizable sulfur atom. Aphid activity, mite activity, and acute toxicity to the rat are all increased about 1 log unit when the oxygen of the oxathia compounds is replaced with sulfur. Activity against the southern armyworm, Mexican bean beetle, and housefly is much less affected by the oxygen/sulfur change. Predictions from the Regression Equations. In Table V are given predicted vs. observed data for eight compounds illustrating the predictive utility of the QSAR regression equations. Two examples (section a, Table V) are given illustrating predicted vs. observed data for compounds which were part of the data set underlying the

J. Agric. Food Chem., Vol. 35, No. 1 , 1987

120

Kurtz and Durden

Table V. Regression Predictions (Mixed Equations) predicted/observed

no. XES” x E B y 18.9

LCm, PPm bean two-spotted southern Mexican LD,, m g / k l o g P MW aphid mite armyworm bean beetle housefly acute oral, rat

8

0.77

a. Illustrative (Predicted/Observed) from Regression Set 2.48 0.0 1.48 218 57/40 437/500 227/500 137/60

30/9

25/16.2

14

0.0

2.48

8/16

1.6/1.6

T-1

1.24

b. Test Cases (Predicted/Observed) - New Compounds 1.17 1.0 1.89’ 248 15/60 35/45 580/1000 88/35

23/90

30/ ...

T-2

1.24 2.48 0.0

T-3

0.0

T-4

P-1

0.0

0.82

204 5/4

38/40

15/11

14/14

CH3”0c0NHCH3

bS CHZCH~

0.452 204 22/12

128/112

38/8

27/6

4/3

8.9/5.6

2.48 0.0

0.763

218 5/8

38/18

16/17

14/3

8/18

1.7/0.6

1.17

2.48 0.0

0.954

218 32/20

207/100

88/1000

73/1000

11/11

20.2/ 19.0

0.45

2.48

0.0

425

mNoCoNHCH3

114

50

4

P-2 -0.72

2.48

1.0

507

359

5

CH~HsNOCONHCH~

0

u

CH3

NOCONHCH3

c. Other Predictions (Compounds Not Made) 2.125 258 423 4080 3700

v

3.316

276 10,000

66,000

313,000

NOCONHCH3

S V S

“Key: 1, 520 + 0.50; 2, 15

+ 2.00.

+ 0.45; 3, 14 - 1 + 15; 4, T-2 + 0.50; 5 , 1 + 2.51 (cyclohexyl; Tute, 1971) - 0.45 (ring-ring

regression equations, compounds 8 and 14. The validity of the models for compounds not present in the regression data sets is tested with compounds T-1-T-4 (section b, Table V), compounds made subsequent to the development of the model. The data for T-1-T-4 indicate how well the contributions of branching, log P, and Ia.s calculated for the four test compounds correctly predict biological activity. Finally, section c of the table illustrates how activity levels predicted by QSAR equations can discourage the synthesis of disadvantageous analogues, p-1 and P-2 (P as a compound label prefix designates a compound proposed for synthesis but not made). C0N CLU S ION S

The potency of insecticidal oxathia and dithia heterocyclic oxime methylcarbamates is promoted by an increase in steric requirement at the a and y sites on the molecule. This advantageous role of increased alkyl bulk is, however, modified by the negative contribution to activity in the insects, the mite, and the rat of increased lipophilicity. Changing the oxygen atom in the oxathia compounds to sulfur uniformly enhances activity. It is interesting to speculate that the improvement in activity with a and y branching is related to an improvement in how well the inhibitor mimics the natural target enzyme substrate acetylcholine. It is further interesting to attribute the negative role of increased lipophilicity to

attachment); 6, 526

the role of increased lipophilicity in promoting metabolic detoxification. There is precedent (Gaudette and Brodie, 1959) for the premise that metabolism in vivo is promoted by increased lipophilicity. Finally, the advantageous role of sulfur vs. oxygen may lie in the ease with which the former is oxidized in vivo; not only are sulfoxide-bearing oxime carbamates known to be excellent acetylcholinesterase inhibitors (Metcalf et al., 1965; Payne et al., 1966; Weiden, 1968), but sulfoxidation substantially lowers log P, a factor shown by the present work to be advantageous. The roles of substituent size a t the a- and y-positions and log P are essentially identical for the oxathia and dithia series; the individual QSAR for the oxathia series may be superimposed on the individual QSAR for the dithia series by the simple expedient of expressing the presence or absence of the oxidizable sulfur using the Ip-s term. Selective improvement in insect vs. mammalian toxicity is not likely to be achieved by further synthesis in this series. The compounds are toxic to the rat. Minimization of that toxicity by decreasing the size of groups a t the aor y-positions or increasing log P is accompanied by parallel reduction in insecticidal activity. The single defined factor making possible reduced mammalian toxicity without equal loss of insect potency is the reduction in toxicity associated with removal of the sulfur. Unfortu-

J. Agric. FoodChem., Vol. 35, No. 1, 1987

SAR In Oxathiolane and Dithiolane Oxime Carbamates

121

expense of maximum insecticidal activity by some combination of limitation of branching LY and y to the oxime, avoidance of the oxidizable sulfur, or use of log P levels higher than those optimum for maximum insecticide performance. ACKNOWLEDGMENT

Range in log p

Thanks to T. D. D’Silva and colleagues for data for the compounds in the dithia series. Insect test data were developed by A. A. Sousa of our Biological Evaluations Department, and Is0determinations were carried out by M. H. J. Weiden of our Biological Fundamentals Group. We are indebted to C. S. Weil and colleagues of the Carnegie-Mellon Hygiene Fellowship at Bushey Run, PA, for carrying out the acute oral rat toxicology studies. 3

4

5

log P

Figure 1. Biological potency vs. log P.

nately, reduction in bean aphid and two-spotted mite potency is realized concomitantly. The optimum compounds in the series have probably been svnthesized. Thev are analormes in which somewhat atteniated mammalia; toxicity cas been achieved a t the

LITERATURE CITED Cohen, S. G.; Chrishti, S. B.; Elkind, J. L.; Reese, H.; Cohen, J. B. J. Med. Chem. 1985,28, 1309. Craig, P. N. J. Med. Chem. 1971, 14, 680. D’Silva, T. D. J.; Durden, J. A., Jr.; Sousa, A. A.; Weiden, M. H. J. J.Agric. Food Chern. 1985, 33, 110. 1964,86,5175. Fujita, T.; Isawa, J.; Hansch, C. J. Am. Chem. SOC. Gaudette, L. E.; Brodie, B. B. Biochem. Pharmacol. 1959,2,89. Kurtz, A. P.; Durden, J. A., Jr.; Sousa, A. A.; Weiden, M. H. J. J. Agric. Food Chem. 1987, preceding paper in this issue. Metcalf, R. L.; Fukuto, T. R.; Frederickson, M.; Peak, L. J. Agric. Food Chem. 1965, 13, 473. Moorefield, H. H.; Tefft, E. R. Contrib. Boyce Thompson Inst. 1958, 19, 295. Payne, L. K., Jr.; Stansbury, H. A., Jr.; Weiden, M. H. J. J. Agric. Food Chem. 1966, 14, 356. Talvik, I.; Palm, V. Org. React. (USSR) 1971, 8, 445. Tute, M. S. Advances in Drug Research; Harper, N. J., Simmonds, A. B., Eds.; Academic: New York, 1971; Vol. 6, Chapter 1,pp 1-78. Weiden, M. H. J. J. Sci. Food Agric. Suppl. 1968, 19. Received for review March 14,1986. Revised September 11,1986. Accepted October 22, 1986. Presented in part a t the 172nd National Meeting of the American Chemical Society, San Francisco, CA, Aug 29-Sept 3, 1976; PEST 14.