Subscriber access provided by HACETTEPE UNIVERSITESI KUTUPHANESI
Article
Physicochemical evidences on the sublethal neonicotinoid imidacloprid interacting with an odorant-binding protein from the tea geometrid moth, Ectropis obliqua Hongliang LI, Lei Zhao, Xiaobin Fu, Xinmi Song, Fan Wu, Mingzhu Tang, Hongchun Cui, and Jizhong Yu J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.7b00597 • Publication Date (Web): 02 Apr 2017 Downloaded from http://pubs.acs.org on April 3, 2017
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
Journal of Agricultural and Food Chemistry is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 35
Journal of Agricultural and Food Chemistry
1
Physicochemical evidences on the sublethal neonicotinoid
2
imidacloprid interacting with an odorant-binding protein from
3
the tea geometrid moth, Ectropis obliqua
4 Hongliang Li, *, †, Lei Zhao †, Xiaobin Fu †, Xinmi Song †, Fan Wu †, Mingzhu Tang †,
5
,‡
Hongchun Cui ‡, Jizhong Yu *
6 7
†
8
Sciences, China Jiliang University, Hangzhou 310018, China;
9
‡
Zhejiang Provincial Key Laboratory of Biometrology and Inspection & Quarantine, College of Life
Tea Research Institute, Hangzhou Academy of Agricultural Sciences, Hangzhou 310024, China.
10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
1
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
27
ABSTRACT: Nowadays the excessive usage of neonicotinoid insecticides always
28
results in the residues in Chinese tea fields. It is not clear whether the insecticide
29
residue at the sublethal level influence on the physiological processes of tea pests.
30
Here, we provided the interacting evidences between neonicotinoid imidacloprid and
31
a general odorant-binding protein, EoblGOBP2, from the tea geometrid moth,
32
Ectropis obliqua. The interacting process was demonstrated through multiple
33
fluorescence spectra, UV absorption spectra, circular dichroism (CD) spectra and
34
molecular docking, etc. The binding mode was determined to be static (from 300 to
35
310 K) and dynamic quenching (from 290 to 300 K). The binding distance was
36
calculated to be 6.9 nm based on the FRET theory. According to the thermodynamic
37
analysis, the process was mainly driven by enthalpy (∆H97 % purity) and E-2-hexenal (>98 % purity) were both
113
purchased from TCI co. Japan, and dissolved in spectroscopic pure grade methanol
114
(Tedia, USA) to prepare 1.0 ×10
115
dark. Milli-Q water (18.2 MΩ, Millipore, US) was used, and all the other solvents and
116
chemicals used in this study were of analytical reagent grade.
−3
mol/L stock solutions, and stored at 4 °C in the
117 118
Insects and Extract of Total Moth Antennae RNA
119
The caterpillars of tea geometrid Ectropis obliqua Prout were collected from tea
120
plantation in HangZhou Daqinggu, and then kept with fresh tea leaves until eclosion
121
in the laboratory at the College of Life Sciences, China Jiliang University at 16 h
122
light : 8 h dark. The antennae of female or male moths were excised, and immediately
123
frozen in liquid nitrogen, and stored at –70 °C until use.
124
According to the kits manual of TRIzol® Reagent (Invitrogen, USA), the total
125
RNA was extracted from the antennae of 100 male and female moths, respectively.
126
The UV A260/280 values for all total antennal RNA were ranging from 1.90 to 2.0.
127
Using 2 µg antennal total RNA, the first-strand cDNA was synthesized by PrimeScript
128
Reverse Transcriptase System (Takara, Japan) according to the kits manual.
129 130 131
Preparation of Recombinant EoblGOBP2 Protein The recombinant pET32-EobGOBP2 plasmid, constructed and obtained in the 24
132
previous study,
was transformed into BL21(DE3) E. coli competent cell, then the
133
positive clone was inoculated into BL medium with ampicillin. The recombinant
134
EoblGOBP2 protein was induced in the presence of 1.0 ×10
135
thiogalactoside (IPTG) with at 37 °C with 200 rpm for 4 hours. After ultrasonication, 6
ACS Paragon Plus Environment
−3
mol/L isopropyl-β-d-
Page 6 of 35
Page 7 of 35
Journal of Agricultural and Food Chemistry
136
from broken cells, the recombinant EoblGOBP2 protein was purified by the Ni2+-NTA
137
column with gradient concentration imidazole washing. When the primary purified
138
protein was dialyzed by PBS buffer (pH 7.4) until imidazole free, the final purified
139
protein including N-terminal tag was digested by thrombin, and then obtained by
140
using general gel filtration method. After detected by SDS-PAGE, the protein was
141
diluted into 1.0×10
142
until further analysis.
−6
mol/L stock solutions in PBS (pH 7.4), and stored at −20 °C
143 144
Absolute Quantification of EoblGOBP2 in Moth Antennae
145
For the assessment of EoblGOBP2 amounts in E. oblique, the absolute
146
quantification of EoblGOBP2 in total mRNA of E. oblique was calculated based on
147
the standard curve of pMD-18T/EoblGOBP2 plasmid. The primers of qRT-PCR were
148
designed as sense primer is: 5’-CAGAGGCTCGTGAGGATG-3’, and antisense
149
primer is: 5’- CCTGACCATTAGGGAAGC-3’, and synthesized by Corporation
150
Sangon Bietech (Shanghai, CN).
151
Real time PCR was performed in Bio-Rad iQ5 (Bio-Rad, USA). The 20 µL
152
reaction system consisted of 10 µL of 2 ×SYBR Premix EX Taq, 0.4 µL each primer
153
(10 µM), 2 µL cDNA sample (total 12 tissues), 7.2 µL ddH2O. The PCR program was
154
at 95℃ for 30 sec, 40 cycles of 95℃ for 10 sec, and 59℃ for 30 sec. The negative
155
controls were set, and the interaction was replaced with ddH2O instead of cDNA
156
samples. Three experimental technical replications were performed for each sample.
157 158
Fluorescence Spectra Measurements
159
The fluorescence spectral data of EoblGOBP2 with imidacloprid were recorded
160
using a RF-5301 PC spectrofluorimeter (Shimadzu, Japan) with a 1.0 cm quartz cell at
161
different temperatures (290, 300, and 310 K). An electronic thermostat water bath 7
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 8 of 35
162
(9012, PolyScience, USA) was used to precisely control the temperature of the
163
reaction. For the fluorescence quenching spectra, the stock solution of EoblGOBP2
164
(1.0 × 10 −6 mol/L) was titrated with the working solution of imidacloprid (1.0 × 10 −3
165
mol/L). Three individual experiments included:
166
(1). Fluorescence quenching spectra. The excitation and emission slit widths were set
167
at 5.0 nm. The optimum excitation wavelength was 282 nm, determined by the
168
scanning of fluorescence excitation spectra of protein alone (Fig. S1), and then the
169
fluorescence quenching / emission spectra were recorded ranging from 290 to 550 nm.
170
(2). Synchronous fluorescence spectra. The spectra were recorded when ∆λ (the
171
difference between excitation and emission wavelength) was set as 15 nm and 60 nm,
172
respectively.
173
(3). Competitive fluorescence spectra. E-2-hexenal (1.0 × 10
174
into the unique EoblGOBP2 stock solution (1.0 × 10 −6 mol/L) and the complex of the
175
equal concentration of EoblGOBP2 with imidacloprid (1.0 × 10
176
respectively. The spectra scanning parameters are the same to the fluorescence
177
quenching spectra, and then the both dissociation constants of E-2-hexenal with both
178
systems were calculated and compared.
−3
mol/L) was titrated
−6
mol/L),
179 180
UV Measurements
181
The UV absorption spectra of imidacloprid and EoblGOBP2 were measured
182
using a Shimadzu UV-1800 UV spectrophotometer (Shimadzu, Japan) in the
183
wavelength range of 190~400 nm with a 1.0 cm quartz cell at room temperature.
184
Based on the results of fluorescence quenching, the UV absorption spectra of the
185
recombinant EoblGOBP2 protein alone, imidacloprid alone, and the corresponding
186
ratio between imidacloprid and EoblGOBP2 were recorded and obtained. All the UV
8
ACS Paragon Plus Environment
Page 9 of 35
Journal of Agricultural and Food Chemistry
187
measurements were carried out in PBS buffer (pH 7.4) at room temperature.
188 189
Circular Dichroism (CD)
190
The CD spectra were measured using a CD spectrometer (Jasco-815, Japan) with
191
a 1.0 cm path length quartz cuvette (200 µL of total volume). In the presence of
192
imidacloprid, the CD spectra of recombinant EoblGOBP2 (1.0 × 10
193
recorded in the range of 190~270 nm at room temperature. The molar ratios of
194
imidacloprid to EoblGOBP2 were varied as 0:1, 3:1, and 5:1. All the observed CD
195
spectra were baseline-subtracted for PBS buffer (pH 7.4), and the results were taken
196
as CD ellipticity. The contents of the secondary conformation forms of EoblGOBP2,
197
e.g., α-helix, β-sheet, β-turn, and random coil, were analyzed from the CD
198
spectroscopic data using the Raussens’s method. 29
−6
mol⋅L−1) were
199 200 201
Molecular Docking The binding docking analysis of EoblGOBP2 with imidacloprid was performed 30
202
by the Molegro Virtual Docker (MVD 4.1, free trial) software.
203
dimensional (3D) conformer structure of imidacloprid was downloaded from the
204
chemical compound database of PubChem.
205
EoblGOBP2 was obtained from the Swiss-Model Workspace
206
crystal template in Protein Data Bank (PDB). The potential binding sites of
207
EoblGOBP2 were determined using the developed grid-based cavity algorithm. The
208
binding pose of EoblGOBP2-imidacloprid complex was obtained according to the
209
searching algorithm of MolDock Optimizer and the energetic evaluation of MolDock
210
Score. By default, the torsions in the imidacloprid ligand are set flexible during the
211
docking simulation. The optimized binding pose and hydrogen bonds were then
212
analyzed by Ligplot+ software 33 and displayed by Pymol software.34
31
The three-
The predicted 3D crystal structure of
9
ACS Paragon Plus Environment
32
based on the 3D
Journal of Agricultural and Food Chemistry
213
RESULTS AND DISCUSSION
214
Protein Purification and Fluorescence Quenching Assay
215
After the recombinant EoblGOBP2 protein induced and purified, as seen in the
216
SDS-PAGE gel in Fig. 1(A), from the expressional bacteria including the pET32-
217
EoblGOBP2 plasmid, lane 2 and 3 were the uninduced and induced full lysis proteins
218
of bacteria, respectively. As the predicted molecular weight of the recombinant
219
EoblGOBP2 protein are about 17 kDa, and lane 4 and 5 showed the purified target
220
recombinant EoblGOBP2 protein with and without N-terminal tags, respectively. It
221
showed that the recombinant protein was successfully obtained.
222
Since the EoblGOBP2 protein contains a tryptophan residue (Trp37) which
223
consists of a fluorescence group of indole ring, the EoblGOBP2 protein can emit the
224
fluorescence when excited by appropriative ultraviolet light. 35 After diluted into 1.0 ×
225
10
226
EoblGOBP2 in the absence and presence of an increasing amount of imidacloprid are
227
shown in Fig. 1(B). The molecular structure of imidacloprid is shown in the inset of
228
Fig. 1(B). Upon excitation at 282 nm, EoblGOBP2 exhibited a strong fluorescence
229
emission signal and the maximum emission spectrum was observed at 340 nm. The
230
fluorescence intensity of EoblGOBP2 gradually decreased upon increasing the
231
imidacloprid concentration, without changing the maximum emission wavelength and
232
the overall shape of the peaks. This latter result is a clear indication for interactions
233
between imidacloprid and EoblGOBP2. Therefore, the fluorescence quenching effect
234
can most likely be attributed to the generation of non-fluorescent complexes.
235
−6
mol·L −1 stock solutions in PBS (pH 7.4), the fluorescence emission spectra of
[Fig.1 insert here]
236 237 238
Mechanism of Binding Interaction The fluorescence quenching of biological macromolecules by small-molecule 10
ACS Paragon Plus Environment
Page 10 of 35
Page 11 of 35
Journal of Agricultural and Food Chemistry
239
fluorescent probes can be explained through dynamic quenching and static quenching
240
phenomena. The fluorescence quenching mechanism can be commonly analyzed by
241
using the following Stern-Volmer36 and Lineweaver-Burk equation:37
242
243 244
F0 = 1 + K qτ o [Q ] = 1 + K sv [Q ] F
(1)
1 1 KD = + F0 − F F0 F0 [ Q ]
(2)
245
Here F0 and F are the fluorescence intensity in the absence and presence of a
246
quenching species at a given concentration [Q]. Kq is the quenching rate constant of
247
the biomolecule, τ0 is the average lifetime of the molecule with a value of 10 −8 s and
248
without a quenching species, Ksv is the Stern-Volmer dynamic quenching constant and
249
KD is the dissociation constant.38, 39 Both the Stern-Volmer and Lineweaver-Burk plots
250
can be found in Fig. 2(A) and (B), respectively; Ksv, Kq and KD are listed in Table 1.
251
When the quenching process follows a dynamic process, the value of Ksv, i.e. the
252
slope of the plot, increases with increasing temperature, on the contrary, being a static
253
process if Ksv decrease with increasing temperature. 40 As shown Fig. 2-A and Table 1,
254
interestingly, the binding of imidacloprid with EoblGOBP2 seems to be a complex
255
process. For example, the value of Ksv at 300 K was higher than those at both 290 K
256
and 300 K, while that at 290 K was closer to that at 300 K than 310 K (Fig.2-A). It
257
indicates that the interaction between imidacloprid and EoblGOBP2 is a probably
258
dynamic process from 290 to 300 K, while clearly static ranging from 300 to 310 K
259
(Table 1). The latter is most notably, due to the fact that the static quenching is
260
generated by a new complex, and the stability of the complex decreases with
261
increasing temperature. 41
262
[Fig.2 insert here] 11
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 12 of 35
263 264
In addition, UV spectra of EoblGOBP2 in the absence and presence of
265
imidacloprid were recorded at room temperature to further elucidate the quenching
266
mechanism. As showed in Fig. 3, when imidacloprid was added to EoblGOBP2 ([c] =
267
1.0 × 10
268
curve of EoblGOBP2 with imidacloprid (the both [c] = 1.0 × 10
269
maximum peak of UV spectra of EoblGOBP2 evidently increased. This also indicates
270
that the binding interaction between imidacloprid and EoblGOBP2 is static quenching.
271
Because the static quenching usually results from the ground state of the fluorescent
272
molecule, while the overall shape of the spectrum does not change as it is the case in
273
dynamic quenching. 42
274
−6
mol·L
−1
) at equal concentration, compared with the theoretical merging −6
mol·L
−1
), the
[Fig. 3 insert here]
275 276
All the above data obtained demonstrate that fluorescence quenching between
277
EoblGOBP2 and imidacloprid follows a stable static quenching mode at or above
278
room temperature. Nevertheless, as other OBP type, ASP2 specifically expressed in
279
the antennae of honeybee, shows a binding mode that is static upon binding with
280
imidacloprid in the range of 290 to 300 K.
281
natural floral volatile β-ionone proves to follow a dynamic mode in the range of 300
282
to 310 K. 44 Therefore, the detailed binding mode between OBP and the ligand might
283
result from the difference of molecular polarity.
43
However, the binding of ASP2 to the
284
To sum up, as a water-soluble systemic insecticide, imidacloprid can be easily
285
absorbed and transported via a binding protein like EoblGOBP2, which is widely
286
distributed in the hemolymph of the insect. This binding in the range of 300 to 310 K
287
is similar to the static binding interaction between serum albumin and numerous drug
288
molecules or other bioactive components,
45
requiring highly stable temperatures
12
ACS Paragon Plus Environment
Page 13 of 35
Journal of Agricultural and Food Chemistry
289
when transported in vivo. This result is also accordance with the characteristics that
290
the pesticide effect of imidacloprid was greater at the higher temperatures in the range
291
of 30 to 39 ℃ (303 to 312 K). 46
292 293
Determination of Interaction Force
294
It has already been shown that a variety of interactions exist between proteins
295
and organic molecules, e.g. electrostatic forces, hydrophobic interactions, hydrogen
296
bonds, and van der Waals interactions.
297
EoblGOBP2 in this study, the interaction forces can be calculated using the
298
thermodynamic equations as follows:
47
For the interaction of imidacloprid with
299
∆ G = − RT ln K = ∆ H − T ∆ S
300
∆H =
301
∆S = ( ∆H − ∆G ) T
(2)
RT1T2 ln( K 0 , 2 K 0 ,1 ) T2 − T1
(3) (4)
302
where ∆H, ∆S, and ∆G are enthalpy change, entropy change, and the free energy
303
change, respectively. The values of ∆G, ∆H, and ∆S for imidacloprid and EoblGOBP2
304
are calculated. As listed in Table 1, the interactions between EoblGOBP2 and
305
imidacloprid could be divided into two cases, one is the temperature range of 290 to
306
300 K and another is from 300 to 310 K. According to the characteristics of the
307
theoretical parameters, ∆H0 indicate a typical hydrophobic
309
interaction; ∆H0 indicate hydrophobic and electrostatic interactions.
310
suggested that these two cases were spontaneous (due to ∆G