Relationship between physicochemical properties and herbicidal activity of 1,2,5-oxadiazole N-oxide derivatives.

The relationship between the herbicidal activity of a number of novel 1,2,5-oxadiazole N-oxides and some physicochemical properties potentially related with this bioactivity, such as polarity, molecular volume, proton acceptor ability, lipophilicity, and reduction potential were studied. The semiempirical molecular orbital method AM1 was used to calculate theoretical descriptors such as dipolar moment, molecular volume, Mulliken's charge and the octanol/water partition coefficients (log P(o/w)). The values of the reduction potentials (E(r)) were obtained by cyclic voltammetry. In addition, the retention factors (log (k'w)) on a reversed-phase high-performance liquid chromatography(RP-HPLC) column in pure aqueous mobile phases were measured for several N-oxide derivatives. The log (k'w) values show good correlation with the calculated values of log P(o/w), showing that the chromatographic parameter can be used as lipophilicity descriptor for these compounds. The multiple regression analysis between the descriptors for the N-oxide derivatives and the herbicide activity indicate that the variance in the biological activity can be explained by changes in the lipophilicity and in the reduction potential.

Introduction

The bio-response of a compound can be the result of several types of interactions between the bioactive compound and the receptor, such as hydrophobic and electrostatic forces, hydrogen bonding and electron donor acceptor complex [1]. These interactions are closely related to physicochemical and structural properties of its component molecules. Therefore it is possible the prediction or explanation of the biological behaviour of molecules from their physicochemical properties. The basic assumption underlying this field of research called quantitative structure-activity relationships (QSAR) is that the structure of a molecule determines its performance [2]. This paradigm can be expressed by P = f (S) [3], where P is any physical, agrochemical, biomedical, toxicological or environmental activity of interest and S may represent either an empirical property of the total molecular structure, a relevant substructure fragment or a theoretical structural descriptor (or a set of descriptors) quantifying some aspects of molecular structure.

Drug lipophilicity is a physicochemical property very useful for bioactivity predictions [4,5]. Transport and distribution processes within biological systems are to a large extent controlled by lipophilicity of the system components. The highly hydrophobic interior of a bilayer membrane enables or facilitates the passage of lipophilic substances and prevents the free diffusion of polar molecules except water in and out of cells [6].

Lipophilicity is usually measured by determining the equilibrium concentration of the compound in two immiscible polar/non polar liquids and expressed as the logarithm of the partition coefficient. The octanol-water partition coefficient (log Po/w) is the parameter most widely used to measure hydrophobicity [7] because it has been shown that this partition system is a good model for many biological processes [6]. However, experimental log Po/w measurements are time-consuming and are limited to a certain range, e.g. –3< log Po/w <3 [8]. It is possible to obtain log Po/w values from fragmental contributions of the different atoms using computational calculations. [9]. Chromatographic methods have been also been used successfully to assess lipophilicity of organic compounds [10, 11]. Thus, the retention factor in RP-HPLC with pure aqueous mobile phases (log k’w) is commonly used as a lipophilic descriptor [8, 12]. Several advantages are attributed to log k’w. A priori, it reflects polar/no polar partitioning in a manner similar to shake-flask measurements and it is dependent on the solute structure and polar functionalities. However, it is difficult to measure directly, because of the prohibitively long retention times of organic solutes in pure water as mobile phase.

On the other hand, electrostatic properties play a crucial role in the receptor-bioactive compound recognition process. A great number of electrostatic descriptors have been described. Basically, semi-empirical and ab initio methods have to be distinguished in the electrostatic descriptors determinations. The most frequently used descriptor includes the highest occupied and lower unoccupied molecular orbital, frontier orbital electron densities, Mulliken´s population charge distribution and dipole moments [13,14].

In this sense, some electrostatic properties i.e. electron affinity, LUMO’s and HOMO’s energies can be related to experimental properties such as redox potentials. Furthermore, the latter is known to be related to bio-reduction or bio-oxidation processes of bioactive compounds [15]. The electrochemical response can be experimentally studied, i.e. by polarographic or by cyclic voltammetry techniques [16].

In this work we studied the relationship between the phytotoxic activity of a number of novel 1,2,5-oxadiazole N-oxide with a variegated set of substituents (Scheme I) and some physicochemical properties potentially related with such activity, such as lipophilicity descriptors, dipolar moment, molecular volume, hydrogen bond acceptor ability and reduction potential. Previous studies [17] have shown that these 1,2,5-oxadiazole N-oxide derivatives exhibit moderate to good herbicidal activity against Triticum aestivum, indicating that the N-oxide moiety could be participating in the phytophore, since in absence of it, a sharp loss of activity was observed (derivative 19, Scheme 1).

Scheme 1

Results and Discussion

The following parameters have used: the chromatographic log k’w and the calculated log Po/w parameters for lipophilicity; the calculated molecular volume, V, as a geometrical (steric) descriptor; charge density, qo, on the O atom of the N→O group to estimate the hydrogen-bond acceptor ability [18,19]; reduction potential, Er, and the LUMO’s energy as accounting for electron acceptor interactions or electronic affinity; and the calculated dipolar moment, μ as an electrostatic descriptor. All the measured and calculated parameters for those compounds are reported in Table 1. In the following sections the details the methodology used to obtain this descriptors is described.

Table 1

Herbicidal activity indexes for the N-oxide derivatives and some of their experimental and calculated descriptors.

Compounda	Rwb	Awb	Lflb	Erc	log k’wd	log Po/we	Vf	μf	qoN→O, f
1	-19	18	-4	-1.20	-	4.95	534	2.06	0.238
2	-48	16	-22	-1.19	3.34	3.76	742	4.52	0.304
3	-57	8	-30	-1.03	3.24	3.69	723	3.94	0.304
4	-15	-7	-4	-0.99	4.24	4.25	573	4.7	0.292
5	0	0	0	-1.22	-	4.03	943	4.07	0.298
6	-11	-5	-5	-1.62	0.22	1.88	375	2.17	0.302
7	97	-50	-50	-2.08	1.47	2.33	575	5.67	0.312
8	-75	-85	-50	-1.59	4.39	4.01	859	7.91	0.316
9	-86	-95	--	--	1.18	2.25	568	1.38	0.309
10	-69	-58	-58	-2.00	2.18	3.42	451	4.76	0.274
11	-4	-5	-30	-0.93	2.40	4.27	454	1.39	0.274
12	-23	-9	-16	--	-	2.66	475	2.02	0.280
13	-17	13	-6	-1.17	2.02	3.08	660	5.20	0.276
14	--	--	--	-1.10	2.31	2.60	809	8.8	0.278
15	-11	-2	-23	-1.54	1.46	2.20	685	3.59	0.281
16	-57	-52	-55	-2.06	1.27	4.89	675	2.59	0.278
17	-51	-66	-71	-1.54	2.74	5.75	706	4.93	0.269
18	-91	-73	-50	-1.43	-	6.74	751	3.65	0.265
19	10	20	35	--	2.21	1.80	635	3.36	c

a Compounds 1-9 oxadiazole derivatives, 10-19 benzooxadiazole derivatives. b Pre-emergent herbicidal activity indexes(change of the corresponding parameters respect to untreated reference, (-) denote decreasing parameters, (+) denote increasing parameters : radicular weight (rw), aerial weight (aw), and length of the leaf (lfl); c Reduction Potential (V, vs Calomel saturated electrode); d The logarithm of the chromatographic retention factor extrapolated to pure water; e Computationally calculated parameters: Po/w is the octanol-water partition coefficient, Vf is the molecular volume, μ the dipolar moment and qoN→O the atomic charge on the O atom of N-oxide moiety respectively;

Lipophilicity determinations using an octadecylsilane (ODS) RP-HPLC column.

Thus, log k'w is most often estimated by extrapolating a linear plot of log k' vs. volume percent of organic modifier (Φ) to 100% water [20], equation (1): log k'= log k'

Kaliszan [21] discussed the shortcomings of using log k’w as a lipophilic descriptor, that includes the actual dependence of log k’ with Φ (linear or quadratic) and, in some instances, with the organic modifier which results in a questionable physical meaning for log k’w. Neither the log Po/w nor any of the chromatographic hydrophobicity scales is unique, pharmacologically distinguished [22]. In light of this, we believe that the best chromatographic or partitioning system for a particular case under study is that which better mimics the molecular interactions between the drug and the biological system, but up to the present, finding that scale is an empirical process.

An alternative method to estimate the log k'w values uses the ET(30) solvatochromic scale, an empirical measure of solvent polarity (equation 2): log k'

Here log k'w is obtained by replacing the ET(30) value of pure water in equation 2 and the a and b coefficients are obtained by linear extrapolation of log k'w vs ET(30) values for the different mobile phases employed [22]. The logarithm of the solute retention factor, k’ = [tr – t0]/t0, (tr and to being the solute retention time and the dead time, respectively) in an ODS column obtained for each mobile phase composition was extrapolated to pure water (by using equation 2) to obtain the lipophilic descriptor log k’w . We used equation 2 since no linearity was observed with equation 1 for most of the compounds. The relative standard deviations for three replicates of k’ were, in each instance, fairly low (under 0.5%). The regression equation parameters (equation 2) are reported in Table 2.

Table 2

Retention model parameters from equation 2 for the studied N-oxide derivatives.

Compound	log k’w	-a	b	r2,a
1c	-	-	-	-
2	3.34	25.1	0.45	0.9992
3	3.24	25.2	0.45	0.9996
4	4.24	29.1	0.53	0.9932
5c	-	-	-	-
6	0.22	8.5	0.14	0.9722
7	1.47	13.5	0.24	0.9956
8	4.39	31.8	0.57	0.9914
9	1.18	11.9	0.21	0.9954
10	2.18	13.5	0.25	0.9920
11	2.40	19.1	0.34	0.9565
12b	-	-	-	-
13	2.02	16.7	0.29	0.9980
14	2.31	21.1	0.37	0.9978
15	1.46	15.6	0.27	0.9972
16	1.27	10.3	0.18	0.9851
17	2.74	18.47	0.34	0.9958
18c	-	-	-	-
19	2.21	19.2	0.34	0.9864

a Square regression coefficient; b Appropriate elution was not possible for the whole range

According to structural and polarity distribution analogies between a typical solvated HPLC reverse phase and the octanol-water system, it is expected for the theoretically calculated log Po/w parameter to be linearly related to the partition chromatographic parameter, log k'w (equation 3) [23]. log P where, m and n are the regression constants.

In Figure 1 the extrapolated log k´w for the ODS column are plotted against the calculated log Po/w data. A fair correlation (equation 4) is obtained: log k' r

Figure 1

Relationship between log k’w and log Po/w for ODS column.

Thus, considering the approximations of the calculated log Po/w data and assuming that this parameter as mostly reflecting hydrophobic interactions (as in an octanol-water partition system), the obtained linear regression means that the molecular interactions between these N-oxide derivatives and the ODS column are mainly of hydrophobic nature. Consequently, both descriptors can be used alternatively to model lipophilic interactions.

Cyclic Voltammetry and ESR

Table 1 lists the values of the first cathodic peaks, the reduction potentials, Er, obtained by cyclic voltammetry. These compounds displayed comparable voltammetric behaviour, showing one, two or three well-defined reduction waves both in DMSO and in DMF. The first wave corresponds to a quasi-reversible process that could be assigned to the following process:

To confirm the free radical production in the first wave we performed ESR spectroscopy. Thus, the N-oxide free radicals were prepared in situ by electrochemical reductions in DMSO, applying a potential corresponding to the first wave for the N-oxides as obtained from the cyclic voltammetric experiments. The interpretation of the ESR spectra by means of a simulation process led to the determination of the coupling constants for all the magnetic nuclei, confirmed by theoretical calculations. [26] In Figure 2 the ESR spectrum of compound 18 is shown as an example. In this case, two triplets could be assigned to the two nitrogen atoms of the oxadiazole system (aN1= 7.10 G; aN2= 2.50 G) and one triplet was accounted for two aromatic ring hydrogen atoms (aH= 0.70 G).

Figure 2

ESR spectra of compound 18 in DMSO.

Calculation of quantum-chemical descriptors

The calculation of the quantum-chemical descriptorsμ, and q0, for the different compounds were performed using a semiempirical molecular orbital method (AM1) [31], starting from standard bond lengths and bond angles. All geometries were fully optimised by minimising the energy with respect to geometrical variables without symmetry constraints, using a 0.01 Kcal/mol gradient and the Polak-Ribierie algorithm as convergence criteria. The log Po/w and V values were calculated by using the ChemPlus extension.

Complexes of pyridine N-oxides with some acids are used in agriculture to regulate plant growth rates. In order to analyse the mechanism of this biological activity in [18] the structural and electronic characteristic of the hydrogen bonding interaction between pyridine derivatives and phenols by experimental and semiempirical methods (AM1 and PM3) was studied. They concluded that the proton acceptor ability of pyridine N-oxide derivatives correlates with the charge (qo) on the oxygen atom of the N→O group. Changes in the charge of the oxygen atom are reflected in the equilibrium constant of complex formation. In the present work, qo was utilized as hydrogen bond acceptor ability index.

In the oxadiazole derivatives (compounds 1-9) the different functionalities are included in the same cycle than the N→O. Consequently, the changes observed in the Mulliken´s charge (qo) on the oxygen atom of the N→O group for different functionalities, are more significant than in the benzo-oxadiazole derivatives (compounds 10-18, Table 1).

Contributing factors to the biological response

The choice of the descriptors utilized in this work is based in previous studies that have shown that electronic effects and lipophilic-hydrophilic balance of the compounds may explain some differences in biological activity [24]. As multilinearity diagnostic the correlation matrix for the used physicochemical properties descriptors was performed (Table 3). The values obtained indicate a small cross-correlation between most of the descriptors, except to logk´w and V. The analysis of the correlation matrix allowed one to safely use these factors in the multi-parameter regression [24]. However, calculated log Po/w and qo show certain correlation when they are analyzed by family of compounds ( r2 = 0.924 for compounds 1-5 and r2 = 0.775 for compounds 10-18) such was described by [19].

Table 3

Correlation matrix for the used physicochemical descriptors

	Er	logk´w	logPo/w	V	μ	qo
Er	1	0.25	0.03	0.02	0.05	0.01
Logk´w	–	1	0.92 (4 outliers)	0.78	0.37	0.02
logPo/w	–	–	1	0.13	0.01	0.23
V	–	–	–	1	0.29	0.13
μ	–	–	–	–	1	0.17
qo	–	–	–	–	–	1

Multi-parameter regressions between the different herbicide bioactivity indexes and the calculated and experimental descriptors were performed. The adjust regression coefficient (AdjR-squared) has been utilized as evaluation criteria of the regression sub-models. The regression model quality has been validated in function of the analysis of residuals. The best results are summarized in Table 4. No good regressions were obtained with the calculated quantum-chemical parameters V, qo, and μ.

Table 4

Multi-parameter regressions between the biological indexes and physico-chemical properties of N-oxide derivatives

N°	Multi-parameter regression	R-squared	Adj-R-squared	Prob> F	n b
I	aw = (167 ± 26) + (95±12) Er – (18±3) log k’w	0.8866	0.8542	0.005	10
II	aw = (132 ± 24) + (66±10)Er – (15±3) log P o/w	0.7490	0.7074	0.0005	13
III	lfl = (78 ± 20) + (56 ± 9) Er – (9±3) log k’w	0.8270	0.7780	0.002	10
IV	lfl = (74 ± 17) + (48±9) Er – (8±2) log P o/w	0.8445	0.8000	0.0015	10

a number of compounds

However, acceptable regression was obtained with the experimental descriptors Er and log k’w, and the biological index. (Table 4, I). The use of the calculated lipophilic descriptor logPo/w values for correlation (Table 4, II) yield comparable results that as in I. That demonstrates than log k’w and log Po/w descriptors can be used alternatively for these compounds to model lipophilic interactions. There are similar correlations between physicochemical descriptors and both (Table 4, I and II) and biological indexes (Table 4, III and IV). These results are not surprising since and are related parameters as is shown in the correlation matrix given in Table 5.

Table 5

The correlation matrix between herbicidal activity indexes of 1,2,5-oxadiazole N-oxide derivatives.

	rw	aw	lfl
rw	1	0.19	0.13
aw	–	1	0.67
lfl	–	–	1

On the other hand, very poor correlations were obtained with the index. Other physicochemical descriptors could be related to the oxadiazole’s effects over radicular weight, presumably because the mechanism of action at this level is different from that who affects the aerial parts of the plant.

Also, a relevant relationship was observed between the variance in the biological activity produced by changes in the different functionalities of the N-oxide derivatives and both the reduction potential and lipophilicity (Table 4). We propose that the interaction between the drug and the receptor is electronic in nature, where the N-oxide moiety participates in the redox process, but also it depends on the hydrophobic balance of the compounds.

Conclusions

The results obtained show that log k’w and log Po/w descriptors can be used alternatively to model lipophilic interactions for the N-oxide series studied. An irrelevant relationship was observed between the variance in the biological activity produced by changes in the different functionalities in the N-oxide derivatives and the proton acceptor ability, measured by the Mülliken´s charge on the oxygen atom of the N-oxide (qo), probably due to the small variation observed in this descriptor for the studied compounds. The obtained physicochemical descriptors Er and log k’w allow a qualitative, reasonable explanation for the bioactivity data of these N-oxide derivatives with a variegated set of substituents.

Experimental

General

The 1,2,5-oxadiazole N-oxide derivatives studied, summarized in Scheme 1, were prepared as described previously [17]. Acetonitrile (ACN) was HPLC grade from Sintorgan and used as received except for the chromatographic measurements for which was filtered through a Nylon 66, 0.45 μm membrane. For the spectroscopic studies, the UV cut-off point in 10 mm cell against air was used as purity criteria [25]. Ultra-pure water was obtained using LABCONCO model 90901-01 equipment and then filtered through a cellulose nitrate 0.45 μm filter (MSF). A Varian Model 5000 liquid chromatograph with manual loop injector connected to a model 2550 UV-visible detector and to a Varian 4270/4290 integrator was used for liquid chromatography. A Micropack MCH-10 column, 150 x 0.46 mm i.d. (ODS) thermostated at 25°C was used. The dead time (to) was measured using an aqueous solution of NaNO2 at 210 nm. The UV-visible spectroscopic measurements were made in a Shimadzu UV 2401 PC spectrophotometer thermostated cells. For cyclic voltammetry studies, dimethylsulfoxide (DMSO) and dimethylformamide (DMF) (spectroscopy grade) were obtained from Aldrich. Tetrabutylammonium perchlorate (TBAP) used as supporting electrolyte was obtained from Fluka. Cyclic voltammetry was carried out using a Weenking POS 88 instrument with a Kipp Zonen BD93 recorder, in DMSO or DMF (ca 1.0 x 10-3 mol dm-3), under a nitrogen atmosphere, with TBAP (ca. 0.1 mol dm-3), using three-electrode cells. A mercury-dropping electrode was used as the working electrode, a platinum wire as the auxiliary electrode, and saturated calomel as the reference electrode [26]. To analyze the products obtained in the cyclic voltammetry first waves (one electro reduction product), electron spin resonance (ESR) spectroscopy was employed [27,28]. The N-oxide radicals were generated by in situ electrolytic reduction at 25°C and the ESR spectra were recorded in the X band (9.85 GHz) using a Bruker ECS 106 spectrometer with a rectangular cavity and 50 KHz field modulation. The hyperfine splitting constants were estimated to be accurate within 0.05 G [29,30]. The quantum chemical descriptors (dipolar moment, μ, and the Mulliken´s charge (qo) on the oxygen atom of the N→O group of the N-oxide derivatives) were calculated by semiempirical molecular orbital methods AM1 [31] using the HyperChem software, version 5.0. The log Po/w , and V values were calculated by using the ChemPlus module of HyperChem 5.0. The biological activity, evaluated in pre-emergence herbicide assays, was expressed as the decrease in the following parameters: length of first leaf (lfl), aerial weight (aw), radicular weight (rw), compared with an untreated control. The herbicide evaluation was described in detail in reference [17].