3774
Y. Hu et al. / Bioorg. Med. Chem. Lett. 24 (2014) 3772–3776
Table 3
similar to Molecule A. The difference between Molecule A and B
was that Molecule B was a cis–trans-isomer due to the disordered
state of O4 and C13, while Molecule A was a trans-isomer.
Gaussian 03W software was used to decide the most stable con-
PLS statistics of CoMFA and CoMSIA 3D-QSAR models
PLS statistics
CoMFA
CoMSIA
q2a
0.565
0.983
0.090
180.464
5
0.823
0.945
0.151
103.152
3
0.164
0.330
0.024
0.024
0.458
former of the compounds.29 The trans-isomer
a and cis-isomer b
r2b
sc
were shown in Figure 2. The compounds 7f and 8a were chosen
as the models for calculation. The HF/3-21G was used for prelimin-
ary optimization, and B3LYP/6-31G⁄ was applied for further
optimization. Single point energies of two compounds were calcu-
lated with DFT method at B3LYP/6-311++G⁄⁄ level. The solvent
effect of DMSO was also taken in account. The calculated results
Fd
ONCe
Stericf
0.498
0.502
Electrostaticg
Donorh
Acceptorh
Hydrophobici
showed that the relative energies of 7fa and 7fb were 0.00 and
a
b
c
1.44 kJ/mol, while those of 8a and 8ab were 0.00 and 1.07 kJ/
a
Cross-validated correlation coefficient from leave-one-out.
Noncross-validated r2.
Standard error of estimate.
F-test value.
Optimum number of principal components.
Steric field contribution.
mol. It should be speculated that trans-isomer was predominant.
This was consistent with the result of the X-ray diffraction study.
The ratios of trans-isomer were estimated at 55–60% according to
the intensities of the two pairs of peaks in the NMR spectrum.
All the target compounds were screened for antifungal activity
against Botrytis cinerea.30 The fungicide procymidone was used as
control drug. Procymidone and iprodione are the general fungi-
cides against B. cinerea. The reports gave the corresponding EC50
d
e
f
g
h
Electrostatic field contribution.
Donor and acceptor, of hydrogen bond fields contribution,
respectively.
i
Hydrophobic field contribution.
values of 0.32 lg/mL and 0.80 l
g/mL, respectively.31,32 These data
demonstrated the potential application value of the high active
title compounds. The results in Table 4 showed that most target
compounds exhibited remarkable fungicidal activity. Among them,
the compounds 7f–h and 8g were highly effective against B. cinerea
and their EC50 values were 0.241
and 0.294 g/mL, respectively, close to or lower than 0.240
of the control drug.
l
g/mL, 0.283
l
g/mL, 0.167
l
l
g/mL
g/mL
l
To analyze the structure–activity relationship, the CoMFA and
CoMSIA 3D-QSAR models with a total of twenty-five target com-
pounds were developed.33,34 The statistical parameters were given
in Table 3. The CoMFA model (q2 = 0.565, r2 = 0.983) was based on
the steric and electrostatic fields, and the CoMSIA model
(q2 = 0.823, r2 = 0.945) was based on the steric, electrostatic, hydro-
phobic, hydrogen bond donor and hydrogen bond acceptor fields.
Partial least-square (PLS) analysis was performed to establish a
linear relationship between the molecular fields and the activity
of molecules. Figure 3 showed the alignment of all target com-
pounds used in the training set. Contour maps for the CoMFA
and CoMSIA models were displayed in Figure 4. Experimental
and predicted pEC50 values for the training set and test set were
reported in Table 4, and the correlation plots of CoMFA and CoMSIA
models were shown in Figures. 5 and 6, respectively.
Steric CoMFA map (Fig. 4A) showed green contour around
4-position of phenyl ring indicating bulky groups were favored at
this position. It was confirmed that the compounds 7f–h, 8f–h
and 9b–c exhibited higher antifungal activity. The yellow contours
around 3-position of phenyl ring indicated that compounds with
bulky groups at this position were less potent. This explained that
the compounds 7e, 8e and 9a had not good fungicidal activity.
Electrostatic CoMFA contour map (Fig. 4B) was shown in red
around 4-position of phenyl ring indicating that negative charge
might play a favorable role on activity, such as the compounds
7f–h, 8f–h, 9b–c. Combined with the steric CoMFA map, the bulky
Figure 3. Alignment of all target compounds in the training set.
groups with negative charges were favored at 4-position of phenyl
ring. That also explained that the compounds 7c and 8c had worse
fungicidal activity than 7f and 8f, respectively, due to the smaller
atom radius of fluorine. The blue contours around almost whole
phenyl ring showed that electronegative atom was disfavored at
other positions except for 4-position. The compounds 7b, 7d–e,
8b and 8d–e were good examples.
The colors of the steric and electrostatic contour maps in the
CoMSIA model had the same meanings as those of the CoMFA
model. Similar to steric CoMFA map, steric CoMSIA contour map
(Fig. 4C) indicated bulky group was favored at 4-position of phenyl
ring. In addition, the green contour around 5-position of furandi-
one nucleus suggested that bulky group should be introduced at
this position. So the compounds 7 showed better bioactivity than
the compounds 8, roughly. In agreement with CoMFA model,
electrostatic CoMSIA map (Fig. 4D) indicated that the group with
negative charge was also favored at 4-position of phenyl ring.
Hydrophobic CoMSIA contour map (Fig. 4E) showed yellow con-
tour around 4-position of phenyl ring indicating that hydrophobic
group was favored at the position, such as the compounds 7f–h,
8f–h and 9b–c. The grey contours around 2- and 3-position of phe-
nyl ring indicated that hydrophobic group played a disfavored role
at these positions. The compounds 7b, 7d–e, 8b and 8d–e were
good examples.
R1
R1
O
O
O
H
N
H
N
R3
N
H
N
H
R2
R2
O
O
O
R3
β
α
R3 = H, CH3
Figure 2. The cis–trans-isomers of target compounds.