Physicochemical Evidence on Sublethal Neonicotinoid Imidacloprid

Apr 2, 2017 - Zhejiang Provincial Key Laboratory of Biometrology and Inspection & Quarantine, College of Life Sciences, China Jiliang University, Hang...
0 downloads 0 Views 1MB Size
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