Original Research ARTICLE
Binding Specificity of Two PBPs in the Yellow Peach Moth Conogethes punctiferalis (Guenée)
- 1State Key Laboratory for Biology of Plant Disease and Insect Pest, Institute of Plant Protection, Chinese Academy of Agricultural Science, Beijing, China
- 2Department of Plant Protection, Henan Institute of Science and Technology, Xinxiang, China
- 3Bangladesh Sugarcrop Research Institute, Pabna, Bangladesh
Pheromone binding proteins (PBPs) play an important role in olfaction of insects by transporting sex pheromones across the sensillum lymph to odorant receptors. To obtain a better understanding of the molecular basis between PBPs and semiochemicals, we have cloned, expressed, and purified two PBPs (CpunPBP2 and CpunPBP5) from the antennae of Conogethes punctiferalis. Fluorescence competitive binding assays were used to investigate binding affinities of CpunPBP2 and CpunPBP5 to sex pheromone and volatiles. Results indicate both CpunPBP2 and CpunPBP5 bind sex pheromones E10-16:Ald, Z10-16:Ald and hexadecanal with higher affinities. In addition, CpunPBP2 and CpunPBP5 also could bind some odorants, such as 1-tetradecanol, trans-caryopyllene, farnesene, and β-farnesene. Homology modeling to predict 3D structure and molecular docking to predict key binding sites were used, to better understand interactions of CpunPBP2 and CpunPBP5 with sex pheromones E10-16:Ald and Z10-16:Ald. According to the results, Phe9, Phe33, Ser53, and Phe115 were key binding sites predicted for CpunPBP2, as were Ser9, Phe12, Val115, and Arg120 for CpunPBP5. Binding affinities of four mutants of CpunPBP2 and four mutants of CpunPBP5 with the two sex pheromones were investigated by fluorescence competitive binding assays. Results indicate that single nucleotides mutation may affect interactions between PBPs and sex pheromones. Expression levels of CpunPBP2 and CpunPBP5 in different tissues were evaluated using qPCR. Results show that CpunPBP2 and CpunPBP5 were largely amplified in the antennae, with low expression levels in other tissues. CpunPBP2 was expressed mainly in male antennae, whereas CpunPBP5 was expressed mainly in female antennae. These results provide new insights into understanding the recognition between PBPs and ligands.
Insects depend on a well-developed olfactory system to distinguish odorants and sex pheromones. Odorant binding proteins (OBPs), chemosensory proteins (CSPs), odorant receptors (ORs) and odorant degrading enzymes (ODEs) are involved in the selectivity and sensitivity of olfaction (Leal, 2005, 2013; Fan et al., 2011; Ahmed et al., 2014, 2017). OBPs are small, water-soluble proteins identified from the chemosensory organs, that are carriers between the external environment and chemoreceptors (Ishida et al., 2002; Leal, 2013). As a multi-genes family, OBPs usually are divided into PBPs, general odorant binding proteins (GOBPs) and antennal binding proteins (ABPs) in lepidopteran insects, based on their binding affinity with sex pheromone and odorant molecules (Vogt et al., 1991; Krieger et al., 1996). Actually, GOBPs and ABPs in many insect species also play roles in pheromone detection, because some of them were found to be expressed in long trichoid sensilla, which are known as pheromone-sensitive sensilla, and most of the main contributors to the ligand binding pocket are conserved (Feng and Prestwich, 1997; Maibeche-Coisne et al., 1998; Zhou et al., 2009; He et al., 2010; Liu et al., 2012). Surprisingly, some GOBP have higher binding affinities with sex pheromone than PBP (Zhou et al., 2009; Liu et al., 2012). PBPs are thought to bind and transport hydrophobic sex pheromone molecules across the aqueous sensillum-lymph to specific pheromone receptors on the dendritic membrane of olfactory neurons (Vogt and Riddiford, 1981; Leal et al., 2005; Forstner et al., 2006; Pelosi et al., 2006). In the earlier studies, PBPs are considered mostly male-specific, while other OBPs are expressed in both males and females (Pelosi et al., 2006). As the first step of pheromone recognition, when PBPs bind to different components of sex pheromones, they can lead to species specificity (Willett and Harrison, 1999).
So far, the 3D structure of PBPs in Bombyx mori (Sandler et al., 2000; Horst et al., 2001), Antheraea polyphemus (Mohanty et al., 2004), Leucophaea maderae (Lartigue et al., 2003), Amyelois transitella (Xu et al., 2010; di Luccio et al., 2013), Apis mellifera (Lartigue et al., 2004) have been elucidated both alone and in combination with various ligands. The structure of B. mori PBP (BmorPBP) with bombykol was the first to be studied by X-ray diffraction spectroscopy and nuclear magnetic resonance (NMR) techniques (Sandler et al., 2000; Horst et al., 2001). The binding pocket of BmorPBP was formed by four antiparallel helices (α1, α4, α5, and α6; Sandler et al., 2000), and the conformational transition in solution displayed pH-dependence (Horst et al., 2001). Stability of protein and ligands are maintained by amino acid residues. Some of these residues are critical for binding ligands (Sandler et al., 2000; Mohanty et al., 2004; Thode et al., 2008; Jiang et al., 2009; Yin et al., 2015; Tian and Zhang, 2016; Zhu et al., 2016; Ahmed et al., 2017; Zhang et al., 2017). Of the residues in BmorPBP, Met5, Phe12, Phe36, Trp37, Ile52, Ser56, Phe76, Val94, Glu98, Ala115, and Phe118 are more conserved and involved in binding to bombykol, which suggests they are interacting with ligands (Sandler et al., 2000; Klusák et al., 2003). Thr57, Ser52 and Thr48 in Drosophila melanogaster LUSH are involved in the binding of short-chain n-alcohols. Thr57 mutants had a significant decrease in ability to bind alcohol compounds compared with wild type, which indicates Thr57 is the key site of LUSH binding to small alcohol molecules (Kruse et al., 2003; Thode et al., 2008).
Insect pheromones play an important role in intra-species communication, sexual attraction, mating aggregation and oviposition host-marking. In many moth species, sex pheromones are usually blends of chemical compounds. Airborne pheromones of moths often consist of two or three chemical components, each of which is perceived by specific olfactory receptor neurons (Abraham et al., 2005).
The yellow peach moth, Conogethes punctiferalis (Guenée; Lepidoptera: Crambidae), is an important agricultural pest of peach, apple, chestnut, maize, and sorghum (Luo and Honda, 2015; Ge et al., 2016). The main sex pheromone compounds of yellow peach moth are (E)-10-hexadecenal (E10-16:Ald), along with the two minor components (Z)-10-hexadecenal (Z10-16:Ald) and hexadecenal (16:Ald; Konno et al., 1982; Liu et al., 1994; Kyungsaeng and Park, 2005). Field trials indicate that Z10-16:Ald and 16:Ald alone do not attract males. A blend of these compounds (two or three) was more attractive (Liu et al., 1994). A better understanding of the molecular mechanisms of sex pheromone perception would improve the use of pheromones to control this pest. In this study, two PBP genes, CpunPBP2 and CpunPBP5, which were identified as pheromone binding proteins, are cloned in the antennae of C. punctiferalis and successfully expressed in Escherichia coli. In order to better understand the function of these PBPs, fluorescence displacement binding assays of CpunPBP2 and CpunPBP5 and their mutants are carried out with sex pheromone components.
Materials and Methods
C. punctiferalis larvae were collected from the sunflower Helianthus annuus at Langfang Experimental Station of Chinese Academy of Agricultural Sciences, Hebei Province, China, and reared on fresh maize in an environmentally controlled room. Rearing conditions were 27 ± 1°C, 70–80% relative humidity (RH) and 16:8 light: dark (L:D). Adults were provided with 10% honey solution. After eclosion, the antennae from males and females (80 pairs of each sex) were immediately cut and processed for RNA extraction.
RNA Extraction and Reverse Transcription
Total RNA was isolated from the antennae using Trizol Reagent (Invitrogen, Carlsbad, CA, USA) following manufacturer's recommendations. The integrity of total RNA was assessed with 1.2% agarose gel electrophoresis and the concentration was determined on a NanoDrop 2000 spectrophotometer (Thermo, USA). One μg RNA was added for reverse transcription to cDNA according to product kit instructions (TransGen, Beijing, China).
Cloning and Sequencing
CpunPBP2 (GenBank accession number: GEDO010000019.1; Jia et al., 2016) and CpunPBP5 (GenBank accession number KP985227) of C. punctiferalis were obtained from the antennal cDNA library. The primers were designed to clone the coding region of CpunPBP2 and CpunPBP5 (Table S1; Underlined bases show restriction enzyme sites for forward and reverse primers, respectively). PCR products were separated by electrophoresis on 1% agarose gels in 1 × TAE buffer. Then the specific fragments were cut and purified by DNA gel extraction kit (Axygen, Hangzhou, China) following the manufacturer's protocol. The purified products were cloned into pGEM-T easy vector (TransGen, Beijing, China) and then transformed to TransT1 E. coli competent cells (TransGen, Beijing, China). Positive clones were selected by PCR using M13 primers and then sequenced.
Sequences obtained for alignment and phylogenetic tree construction were downloaded from NCBI database (https://www.ncbi.nlm.nih.gov/), the putative signal peptides were predicted with SignalP 4.1 server (http://www.cbs.dtu.dk/services/SignalP/). Sequence alignments were produced with DNAMAN software. The phylogenetic tree was constructed using the neighbor-joining method with the MEGA 5.2 program (bootstrapping with 1,000 replications; Tamura et al., 2011). Evolutionary distances were computed using the Poisson correction method.
Recombinant Protein Expression and Purification
Prokaryotic expression system (Gu et al., 2012) was used to express CpunPBP2 and CpunPBP5. First, the pGEM plasmid containing the positive clones were digested by Bam HI and Hind III enzymes (NEB, Beijing, China). The expected band was purified and cloned into the bacterial expression vector pET 30a(+) digested with the same enzymes. The pET 30a(+)-CpunPBP2 and pET 30a(+)-CpunPBP5 were transformed into the TransT1 competent cells and grown on LB solid medium with 10 μL kanamycin (10 mg/mL). Positive colonies were selected by PCR using T7 primers and transformed into BL21 (DE3) competent cells (TransGen, Beijing, China). The verified single colony was cultured overnight in 5 mL LB broth including 50 μg/mL kanamycin. LB broth (0.5 L) was inoculated with 5 mL overnight culture at 37°C for 3 h until the absorbance at OD600 reached to 0.6. Then the protein was induced with isopropyl-β-d-thiogalactoside (IPTG) in a final concentration of 1 mM at 37°C for 6 h (Prestwich, 1993). The induced bacterial cells were centrifuged at 4°C for 10 min (10,000 rpm) and resuspended in the PBS buffer (NaCl 137 mmol/L, KCl 2.7 mmol/L, Na2HPO4 10 mmol/L, KH2PO4 2 mmol/L, pH 7.4), agitated by ultrasonic waves (an interval of 5 s, 10 min) and centrifuged again (15,000 rpm, 20 min, 4°C). The supernatant and pellet were analyzed by sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE), which showed that CpunPBP2 and CpunPBP5 were expressed mainly in the precipitate. Precipitate was resolved in 8M carbamide and purified by 6 × His-Tagged Purification Kit (CWbio, Beijing, China). Refolded proteins were dialyzed within PBS buffer overnight at 4°C and then concentrated using Amicom 10 kDa cutoff concentrators (Millipore Billerica, MA, USA). The purity and size were checked by SDS-PAGE. The concentration was determined by the Bradford method using bovine serum albumin (BSA) as standard protein.
Fluorescence Displacement Binding Assay
Fluorescence binding assay was used to measure the affinity of the CpunPBP2 and CpunPBP5 to 3 sex pheromone and 21 volatile compounds (Konno et al., 1982; Kyungsaeng and Park, 2005). The fluorescence intensity was recorded on a FluoroMax-4 spectrophotometer (Horiba Scientific, USA) at room temperature using a 1 cm light path fluorimeter quartz cuvette. The fluorescent probe N-phenyl-1-naphthylamine (1-NPN) and all the tested chemicals were dissolved in HPLC purity methanol. The final concentration was prepared 1 mM. To measure the affinity of florescent ligand 1-NPN to each protein, a 2 μM solution of the protein in 50 mM Tris-HCl, pH 7.4, was titrated with aliquots of 1 mM ligand in methanol to final concentrations of 1–8 μM. The fluorescence of 1-NPN was excited at 337 nm and emission spectra were recorded between 300 and 450 nm. The affinity of other ligands was measured in competitive binding assays, using 1-NPN as the fluorescent reporter at 2 μM concentration and different concentrations of each ligands. The GraphPad Prism 5 (GraphPad Software, Inc.) was used to estimate the K1−NPN (KD of complex protein /1-NPN) values by nonlinear regression for a unique site of binding. It was assumed that the proteins were 100% active, with a stoichiometry of 1:1 (protein:ligand) at saturation. For other competitor ligands, the dissociation constants were calculated from the corresponding IC50 (concentrations of ligands halving the initial fluorescence value of 1-NPN) values using Microsoft Office Excel 2010, with the formula: KD = [IC50]/(1+[1-NPN]/K1−NPN). In the equation, [1-NPN] is the free concentration of 1-NPN, and K1−NPN is the dissociation constant of the complex protein /1-NPN.
Sequences of CpunPBP2 and CpunPBP5 were submitted to the SWISS-MODEL server (http://swissmodel.expasy.org/) for structural modeling with all known proteins to obtain template sequences. Then target and template sequences were aligned with ClustalW program. Finally, three dimensional models of CpunPBP2 and CpunPBP5 were generating using I-TASSER Protein Structure and Function Prediction web server (http://zhanglab.ccmb.med.umich.edu/I-TASSER/; Zhang, 2008; Yang et al., 2015). The 3D structure of E10-16:Ald and Z10-16:Ald were obtained from ChemOffice (http://www.cambridgesoft.com/Ensemble_for_Chemistry/ChemOffice/ChemOfficeProfessional/) and was further refined by the CHARMm force field (http://www.charmm.org/). The model was rendered in PyMol (http://www.pymol.org/). The energy minimization was used to refine the ligand poses. Based on the established homology model, the docking program CDOCKER was used to dock the sex pheromone compounds (E10-16:Ald and Z10-16:Ald) with CpunPBP2 and CpunPBP5 models, respectively. The binding energy included van der Waals energy (Evdw), electrostatic interaction energy (Eeie) and total interaction energy (Etotal). The energy required for interactions among sex pheromone and CpunPBP2 and CpunPBP5 were calculated to select key residues.
Preparation of Site-Directed Mutants
Four mutants of CpunPBP2 and four mutants of CpunPBP5 were developed using the QuikChange Lightning Site-Directed Mutagenesis Kit (Stratagene, USA). The mutational primers were designed manually. Mutation sites are underlined in Table S2. The CpunPBP2/pGEM-T Easy construct was used as a template. The PCR conditions were 95°C for 5 min, followed by 30 cycles of 95°C for 30 s, 58°C for 30 s and 68°C for 1 min, and final extension at 72°C for 10 min. The correct insertion of mutation was subcloned into pGEM-T Easy vector (TransGen, Beijing, China). The expression system and fluorescence binding assay were conducted as mentioned for wild type proteins.
Relative Expression Pattern of CpunPBP2 and CpunPBP5
Antennae, proboscises, maxillary palps, thoraxes, legs, abdomens, heads (without antennae, proboscises, and maxillary palps), and wings (50 pairs of each sex) were collected for total RNA extraction using Trizol reagent (Invitrogen, Carlsbad, CA, USA). The first strand cDNA template was synthetized with One-Step gDNA removal and cDNA Synthesis kit (TransGen, Beijing, China) including oligo dt-primer according to product manual recommendations. The primers of CpunPBP2, CpunPBP5 and reference gene (β-actin, accession number JX119014) for real-quantitative PCR (qPCR) were designed using Primer premier 5.0 program (Premier Biosoft International, Palo Alto, CA, USA; Table S1). qPCR were conducted on ABI 7500 fast real-time PCR system (Applied Biosysterm, USA). Each amplification reaction was performed with 20 μL volume using SYBR Premix Ex Taq II (Tli RNaseH Plus) master mix (Takara-Bio, Shiga, Japan) under the following conditions: 95°C for 30 s, followed by 40 cycles of 95°C for 3 s and 60°C for 30 s. To check reproducibility, each test sample was done in triplicate technical replicates and three biological replicates. Relative quantification was analyzed using the comparative 2−ΔΔCT method (Livak and Schmittgen, 2001). The relative expression levels in different tissues were calculated with the transcript level of the female antennae used as the calibrator.
Sequence Analysis of CpunPBP2 and CpunPBP5
Coding regions of CpunPBP2 and CpunPBP5 were obtained from the antennal cDNA library. Sequence analysis shows that the full-length cDNA encoding CpunPBP2 consists of 513 nucleotides that encode 170 amino acid residues. SignalP predicts that the signal peptide cleavage sites have 25 amino acids. On the other hand, CpunPBP5 contains 507 nucleotides for a polypeptide of 168 amino acids with 25 amino acids as signal peptide. The alignment of amino acid sequences shows that CpunPBP2 and CpunPBP5 have six conserved cysteines, which are typical of classic OBPs (Figure S1; Pelosi et al., 2006). Additionally, a few amino acids also are conserved in the aligned sequences. Compared with the other 81 Lepidopteran PBPs, the phylogenetic tree based on the amino acid sequences shows that CpunPBP2 and CpunPBP5 share closer ancestry with PBPs in Crambidae, Lepidoptera (Figure S2).
Recombinant Protein Expression and Fluorescence Displacement Binding Assay
Recombinant CpunPBP2 and CpunPBP5 (wild type) proteins expressed in E. coli occurred in inclusion bodies and were high yield. The precipitate was resuspended and purified by affinity chromatography (Figure 1) to produce ~1 mg/mL protein, which was used in the fluorescence displacement binding assay.
Figure 1. SDS-PAGE analyses of expressed recombinant of CpunPBP2 (A) and CpunPBP5 (B). M, marker protein; 1, purified fusion protein; 2, inclusion body of induced cells; 3, supernatant from ultrasound treated cells; 4, IPTG induced E. coli pET30a (+)/CpunPBPs transformed BL 21(DE3) cells; 5, Non-induced pET30a (+)/CpunPBPs transformed BL 21(DE3) cells.
Fluorescence of CpunPBP2/1-NPN and CpunPBP5/1-NPN complexes were excited at 337 nm, and the fluorescence peak was 390–410 nm. The dissociation constants (KD) of CpunPBP2/1-NPN and CpunPBP5/1-NPN complexes are 0.76 ± 0.10 μM and 0.58 ± 0.04 μM as measured by Scatchard plots (Figure 2). The IC50 values and the calculated KD of 21 volatiles and 3 sex pheromone analogs to CpunPBP2 and CpunPBP5 are shown in Table 1.
Table 1. IC50 values (μM) and calculated dissociation constants (KD) (μM) of CpunPBP2 and CpunPBP5 to different ligands at pH = 7.4.
Fluorescence intensity of CpunPBP2 and CpunPBP5 gradually declined with the increased concentrations of volatiles and sex pheromone (Figure 3). The results show that CpunPBP2 and CpunPBP5 have the highest binding ability to sex pheromones E10-16:Ald, Z10-16:Ald compared with hexadecanal and other volatiles. Compared with CpunPBP5, CpunPBP2 has a higher binding affinity to the sex pheromone E10-16:Ald and hexadecanal. CpunPBP2 also has a similar binding affinity between E10-16:Ald and Z10-16:Ald. This result indicates there is a definite apparent interaction between the sex pheromones and the two PBPs. Among the volatiles, the binding results indicate that 1-tetradecanol had the highest binding affinity with CpunPBP2 and CpunPBP5, followed by trans-caryopyllene, farnesene, β-farnesene. Interestingly, results also indicate that CpunPBP2 and CpunPBP5 could discriminate the chiral structure of chemical molecules. The two PBPs could bind to α-ionone better than its isomer β-ionone, while is counter to the isomer of pinene. Hexenal, cis-3-hexen-1-ol, α-pinene, β-pinene, trans-2-nonenal and linalool had the minimum binding affinities to CpunPBP2. The vanillin, heptanal, cis-3-hexen-1-ol, linalool had the minimum binding abilities to CpunPBP5.
Figure 3. Competitive binding curves of CpunPBP2 and CpunPBP5 to different ligands. (A1) Competitive binding curves of CpunPBP2 to aldehydes. (A2) Competitive binding curves of CpunPBP2 to alcohols. (A3) Competitive binding curves of CpunPBP2 to olefins. (A4) Competitive binding curves of CpunPBP2 to other compounds. (B1) Competitive binding curves of CpunPBP5 to aldehydes. (B2) Competitive binding curves of CpunPBP5 to alcohols. (B3) Competitive binding curves of CpunPBP5 to olefins. (B4) Competitive binding curves of CpunPBP5 to other compounds.
To predict the 3D structure of CpunPBP2 and CpunPBP5, sequences from other similar proteins were compared. The search suggests BmorPBP (PDB id: 1ls8) and AtraPBP1 (PDB id: 2 kph) were used to construct the 3D structure of CpunPBP2 and CpunPBP5 with high similarity (54.0 and 45.8%), respectively (Figure 4). The predicted 3D structure of CpunPBP2 and CpunPBP5 consists of seven α-helices, and the antiparallel helices converge to form the hydrophobic binding pocket (Figure 4). To further study the binding site of sex pheromone to CpunPBP2 and CpunPBP5, E10-16:Ald and Z10-16:Ald were docked with the predicted CpunPBP2 and CpunPBP5 models (Figure 5). The interaction energies between key residues and the ligands are predicted and listed in Tables 2, 3. Based on the interaction energy of docking models, several residues including Phe9, Phe33, Ser53, and Phe115 in CpunPBP2 and Ser9, Phe12, Val115, and Arg120 in CpunPBP5 seem to play crucial roles in the binding to E10-16:Ald and Z10-16:Ald.
Figure 4. 3D structure model of the CpunPBP2 and CpunPBP5. (A) Predicted 3D model of CpunPBP2 was built on structure of PBP1 from Bombyx mori. Seven α-helixes, N-terminal (Nt) and C-terminal (Ct) are marked. (B) Sequence alignment of CpunPBP2 and BmorPBP. In the alignment of the two proteins, identical residues are highlighted in blue. (C) Predicted 3D model of CpunPBP5 was built on structure of PBP1 from Amyelois transitella. Seven α-helixes, N-terminal (Nt) and C-terminal (Ct) are marked. (D) Sequence alignment of CpunPBP5 and AtraPBP1. In the alignment of the two proteins, identical residues are highlighted in blue.
Figure 5. Docking of CpunPBP2 and CpunPBP5 with pheromone compounds. (A) Molecular docking of predicted CpunPBP2 model with E10-16:Ald. (B) Molecular docking of predicted CpunPBP2 model with Z10-16:Ald. (C) Molecular docking of predicted CpunPBP5 model with E10-16:Ald. (D) Molecular docking of predicted CpunPBP5 model with Z10-16:Ald.
Table 2. Interaction energies (kcal/mol) between the key residues of CpunPBP2 and pheromone compounds.
Table 3. Interaction energies (kcal/mol) between the key residues of CpunPBP5 and pheromone compounds.
Fluorescence Displacement Binding Assay of Mutants
The recombinant mutant proteins were expressed and purified as described for wild type and analyzed by SDS-PAGE (Figure S3). The emission wave lengths of mutants with 1-NPN were 400–410 nm. The binding curve (Figure 6) of CpunPBP2 and CpunPBP5 mutants with 1-NPN complexes were made. The binding affinities of mutant between proteins and sex pheromones are listed in Table 4. The results showed that, compared with CpunPBP2, the mutant Cpun2-m4 likely lost the binding ability to the two sex pheromones (Figure 6). The binding abilities of the three remaining mutants show no significant differences with wild CpunPBP2. Compared with CpunPBP5, the binding ability of all CpunPBP5 mutants to sex pheromones are reduced by varying degrees (Figure 6). The binding affinity of mutant of CpunPBP5-m3 to E10-16:Ald decreased the most, and the binding capacity of CpunPBP5-m4 to Z10-16:Ald also decreased considerably.
Figure 6. Binding of 1-NPN and ligands to CpunPBP2 and CpunPBP5 mutants. (A) The binding curve and relative Scatchard plots of CpunPBP2 and mutants. (B) Competitive binding curves of CpunPBP2 and mutants to E10-16:Ald. (C) Competitive binding curves of CpunPBP2 and mutants to Z10-16:Ald. (D) The binding curve and relative Scatchard plots of CpunPBP5 and mutants. (E) Competitive binding curves of CpunPBP5 and mutants to E10-16:Ald. (F) Competitive binding curves of CpunPBP5 and mutants to Z10-16:Ald.
Table 4. IC50 values (μM) and calculated dissociation constants (KD) (μM) of CpunPBP2 and CpunPBP5 with their mutants to two pheromones.
Tissues-Specific Expression Pattern of CpunPBP2 and CpunPBP5
The expression levels of CpunPBP2 and CpunPBP5 in different tissues (male and female antennae, proboscises, maxillary palps, thoraxes, legs, abdomens, heads, and wings) were evaluated using qPCR. The target product was largely amplified in antennae, with low expression level in other tissues (Figure 7). CpunPBP5 is mainly expressed in the female antennae, which contrasts with CpunPBP2 and its male-specific expression. In general, expression levels of CpunPBP2 and CpunPBP5 in proboscises, maxillary palps, thoraxes, legs, abdomens, heads, and wings were very low or null.
Figure 7. Expression pattern analysis of CpunPBP2 and CpunPBP5 in different tissues. (A) Expression levels of CpunPBP2 in different tissues. (B) Expression levels of CpunPBP5 in different tissues. The different letters (a, b, c) above each bar represented significant differences (p < 0.05).
Odorant binding proteins are essential for insect olfactory perception because they are transporters between the external environment and insect chemoreceptors (Sun Y. L. et al., 2013). Fluorescence binding affinity has emerged as an important method to demonstrate binding capacity with ligands and help elucidate mechanisms of OBPs (Campanacci et al., 2001; Fan et al., 2011). Jia et al. (2015) cloned a PBP from C. punctiferalis and named as CpunPBP1 (GenBank accession number: KP027286), which is similar to CpunPBP2 we obtained. But in 2016 (Jia et al., 2016), they got the same sequence by transcriptome analysis and named as CpunPBP2 (GenBank accession number: GEDO010000019.1). In order to eliminate the confusion, we use the second name in our study. In this study, CpunPBP2 and CpunPBP5 had strong binding abilities with two sex pheromone compounds, indicating that the two PBPs may play important roles in transporting sex pheromones within the sensillar lymph. Furthermore, CpunPBP2 and CpunPBP5 also bind volatiles: 1-tetrodecanol, trans-caryopyllene, farnesene, and β-farnesene, which suggest CpunPBP2 and CpunPBP5, may share similar amino acid binding sites with GOBPs associated with the volatiles (Mao et al., 2016). Interestingly, CpunPBP2 and CpunPBP5 discriminate the chiral structure of chemical molecules, similar to AlinOBP5 results in Adelphocoris lineolatus (Wang et al., 2013). We speculate that the chiral structure of ligands may affect the binding constants and need to be further investigated.
Protein structure plays crucial roles in recognition and binding of ligand molecules. Studies of the interactions between proteins and ligands are necessary to better understand the binding mechanism. Structures of OBP and PBP in other lepidopteran insects, such as B. mori (Sandler et al., 2000; Horst et al., 2001), A. polyphemus (Mohanty et al., 2004) and A. transitella (Xu et al., 2010; di Luccio et al., 2013), were used to provide insights into our PBPs. In this study, the key residues were evaluated based on the energy values. After site-directed mutagenesis, four mutants of CpunPBP2 and four mutants of CpunPBP5 protein were purified and used to analyze the binding mechanism. Compared with CpunPBP2, the binding ability of CpunPBP2 mutants were not significantly reduced, expect for CpunPBP2-m4. We speculate that the amino acid substitution of the three mutants of CpunPBP2 had a slight effect of relaxing the compact structure of the binding site, similar to the loss of high specificity with Plutella xylostella mutants (Zhu et al., 2016). Because Phe115 in CpunPBP2 had a stronger hydrophobic interaction than other amino acids (Table 2) and the binding affinity between CpunPBP2-m4 and sex pheromone compounds sharply decreased, we speculate that Phe115 in CpunPBP2 are involved in sex pheromone recognition. The binding abilities of CpunPBP5 mutants with sex pheromones varied, which suggests that the small protein modifications may have affected the hydrogen bond between protein and sex pheromones. These results may be due to the change of hydrocarbon interactions or the stabilization of the hydrophobic binding pocket. This suggests that the conformation of PBP was influenced by the transformation of the single amino acid (Zhang et al., 2017). Further research using NMR or x-ray to analyze the protein structure may be necessary to better understand these changes.
The expression levels measured by qPCR showed that CpunPBP2 and CpunPBP5 were mainly expressed in antennae, with low expression in the other tissues. These results suggest that CpunPBP2 and CpunPBP5 play a crucial role in odorant chemoreception (including sex pheromone). CpunPBP2 gene was more abundantly expressed in male antennae than in female antennae, which is similar to results found in other insects, including Spodoptera exigua, P. xylostella, Agrotis ipsilon, Helicoverpa armigera, and Maruca vitrata (Xiu and Dong, 2007; Zhang et al., 2011; Gu et al., 2013; Sun M. J. et al., 2013; Mao et al., 2016). High expression of CpunPBP2 in male antennae may indicate that CpunPBP2 is involved in male-female recognition. Expression level of the CpunPBP5 gene in male antennae was lower than that of female antennae, which is similar to results found with M. vitrata and Sesamia inferens (Jin et al., 2014; Mao et al., 2016). Thus, these results suggest CpunPBP2 may be involved in the detection of conspecific sex pheromone and autodetection of sex pheromone compounds (Yang et al., 2009; Holdcraft et al., 2016; Mao et al., 2016).
In conclusion, our study provides key information about CpunPBP2 and CpunPBP5 in C. punctiferalis, which may be useful for developing effective pest management strategies for this pest.
TZ and ZW: Conceived and designed the experimental plan; XG: Preformed the experiments; XG, TZ, and SB: Analyzed the sequence and data; ZW and KH: Provided all the materials and lab facilities necessary for this work; TA and ZW: Revised the manuscript. All authors read and approved the final manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This work was supported by Special Fund for Agro-scientific Research in the Public Interest (201303026) and China Agriculture Research System (CARS-02). The authors gratefully acknowledge Dr. Richard Hellmich for comments and suggestions on this manuscript.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2018.00308/full#supplementary-material
Figure S1. Multiple sequence alignment of CpunPBP2 and CpunPBP5 with other Lepidopteran insect PBPs.
Figure S2. Phylogenetic tree of CpunPBPs amino acid sequences with other 81 PBPs from different insect species.
Figure S3. SDS-PAGE analyses of purified CpunPBP2 and CpunPBP5 with their mutants. M: marker protein. (A)1:CpunPBP2 original protein. (A)2: CpunPBP2-Phe9 mutant (CpunPBP2-m1). (A)3: CpunPBP2-Phe33 mutant (CpunPBP2-m2). (A)4: CpunPBP2-Ser53 mutant (CpunPBP2-m3). (A)5: CpunPBP2-Phe115 mutant (CpunPBP2-m4). (B)1:CpunPBP5 original protein. (B)2: CpunPBP5-Ser9 mutant (CpunPBP5-m1). (B)3: CpunPBP5-Phe12 mutant (CpunPBP5-m2). (B)4: CpunPBP5-Val115 mutant (CpunPBP5-m3). (B)5: CpunPBP5-Arg120 mutant (CpunPBP5-m4).
Table S1. Primers for expression and qPCR.
Table S2. Primers for mutants of CpunPBP2 and CpunPBP5.
Abraham, D., Löfstedt, C., and Picimbon, J. F. (2005). Molecular characterization and evolution of pheromone binding protein genes in Agrotis moths. Insect Biochem. Mol. Biol.35, 1100–1111. doi: 10.1016/j.ibmb.2005.05.002
Ahmed, T., Zhang, T. T., Wang, Z. Y., He, K. L., and Bai, S. X. (2014). Three amino acid residues bind corn odorants to McinOBP1 in the polyembryonic endoparasitoid of Macrocentrus cingulum brischke. PLoS ONE 12:e93501. doi: 10.1371/journal.pone.0093501
Ahmed, T., Zhang, T., Wang, Z., He, K., and Bai, S. (2017). C-terminus methionene specifically involved in binding corn odorants to odorant binding protein4 in Macrocentrus cingulum. Front. Physiol. 8:62. doi: 10.3389/fphys.2017.00062
Campanacci, V., Krieger, J., Bette, S., Sturgis, J. N., Lartigue, A., Cambillau, C., et al. (2001). Revisiting the specificity of Mamestra brassicae and Antheraea polyphemus pheromone-binding proteins with a fluorescence binding assay. J. Biol. Chem. 276, 20078–20084. doi: 10.1074/jbc.M100713200
di Luccio, E., Ishida, Y., Leal, W. S., and Wilson, D. K. (2013). Crystallographic observation of pH-induced conformational changes in the Amyelois transitella pheromone-binding protein AtraPBP1. PLoS ONE 8:e53840. doi: 10.1371/journal.pone.0053840
Fan, J., Francis, F., Liu, Y., Chen, J. L., and Cheng, D. F. (2011). An overview of odorant-binding protein functions in insect peripheral olfactory reception. Genet. Mol. Res. 10, 3056–3069. doi: 10.4238/2011.December.8.2
Ge, X., Zhang, T., Wang, Z., He, K., and Bai, S. (2016). Identification of putative chemosensory receptor genes from yellow peach moth Conogethes punctiferalis (Guenée) antennae transcriptome. Sci. Rep. 6:32636. doi: 10.1038/srep32636
Gu, S. H., Wang, S. Y., Zhang, X. Y., Ji, P., Liu, J. T., Wang, G. R., et al. (2012). Functional characterizations of chemosensory proteins of the alfalfa plant bug Adelphocoris lineolatus indicate their involvement in host recognition. PLoS ONE 7:e42871. doi: 10.1371/journal.pone.0042871
Gu, S. H., Zhou, J. J., Wang, G. R., Zhang, Y. J., and Guo, Y. Y. (2013). Sex pheromone recognition and immunolocalization of three pheromone binding proteins in the black cutworm moth Agrotis ipsilon. Insect Biochem. Mol. Biol. 43, 237–251. doi: 10.1016/j.ibmb.2012.12.009
He, X., Tzotzos, G., Woodcock, C., Pickett, J. A., Hooper, T., Field, L. M., et al. (2010). Binding of the general odorant binding protein of Bombyx mori BmorGOBP2 to the moth sex pheromone components. J. Chem. Ecol. 36, 1293–1305. doi: 10.1007/s10886-010-9870-7
Horst, R., Damberger, F., Luginbuhl, P., Guntert, P., Peng, G., Nikonova, L., et al. (2001). NMR structure reveals intramolecular regulation mechanism for pheromone binding and release. Proc. Natl. Acad. Sci.U.S.A. 98, 14374–14379. doi: 10.1073/pnas.251532998
Jia, X. J., Hao, S. D., Du, Y. L., Zhang, M. Z., Qin, X. C., Wang, J. Z., et al. (2015). cDNA cloning, expression profiling and binding affinity assay of the pheromone binding protein Cpun-PBP1 in the yellow peach moth, Conogethes punctiferalis (Lepidoptera: Crambidae). Acta Entomol. Sin. 58, 1167–1176. doi: 10.16380/j.kcxb.2015.11.003
Jia, X. J., Wang, H. X., Yan, Z. G., Zhang, M. Z., Wei, C. H., Qin, X. C., et al. (2016). Antennal transcriptome and differential expression of olfactory genes in the yellow peach moth, Conogethes punctiferalis (Lepidoptera: Crambidae). Sci. Rep. 6:29067. doi: 10.1038/srep29067
Jiang, Q. Y., Wang, W. X., Zhang, Z. D., and Zhang, L. (2009). Binding specificity of locust odorant binding protein and its key binding site for initial recognition of alcohols. Insect Biochem. Mol. Biol. 39, 440–447. doi: 10.1016/j.ibmb.2009.04.004
Jin, J. Y., Li, Z. Q., Zhang, Y. N., Liu, N. Y., and Dong, S. L. (2014). Different roles suggested by sex-biased expression and pheromone binding affinity among three pheromone binding proteins in the pink rice borer, Sesamia inferens (Walker) (Lepidoptera: Noctuidae). J. Insect Physiol. 66, 71–79. doi: 10.1016/j.jinsphys.2014.05.013
Klusák, V., Havlas, Z., Rulísek, L., Vondrásek, J., and Svatos, A. (2003). Sexual attraction in the silkworm moth. Nature of binding of bombykol in pheromone binding protein–an ab initio study. Chem. Biol. 10, 331–340. doi: 10.1016/S1074-5521(03)00074-7
Konno, Y., Arai, K., and Sekiguchi, K. (1982). (E)-10-hexadecenal, a sex pheromone component of the yellow peach moth, Dichocrocis punctiferalis. Appl. Entomol. Zool. 17, 207–217. doi: 10.1303/aez.17.207
Krieger, J., von Nickisch-Rosenegk, E., Mameli, M., Pelosi, P., and Breer, H. (1996). Binding proteins from the antennae of Bombyx mori. Insect Biochem. Mol. Biol. 26, 297–307. doi: 10.1016/0965-1748(95)00096-8
Kruse, S. W., Zhao, R., Smith, D. P., and Jones, D. N. M. (2003). Structure of a specific alcohol-binding site defined by the odorant binding protein LUSH from Drosophila melanogaster. Nat. Struct. Biol. 10, 694–700. doi: 10.1038/nsb960
Lartigue, A., Gruez, A., Briand, L., Blon, F., Bezirard, V., Walsh, M., et al. (2004). Sulfur single-wavelength anomalous diffraction crystal structure of a pheromone-binding protein from the honeybee Apis mellifera L. J. Biol. Chem. 279, 4459–4464. doi: 10.1074/jbc.M311212200
Lartigue, A., Gruez, A., Spinelli, S., Riviere, S., Brossut, R., Tegoni, M., et al. (2003). The crystal structure of a cockroach pheromone-binding protein suggests a new ligand binding and release mechanism. J. Biol. Chem. 278, 30213–30218. doi: 10.1074/jbc.M304688200
Leal, W. S., Chen, A. M., Ishida, Y., Chiang, V. P., Erickson, M. L., Morgan, T. I., et al. (2005). Kinetics and molecular properties of pheromone binding and release. Proc. Natl. Acad. Sci. U.S.A. 102, 5386–5391. doi: 10.1073/pnas.0501447102
Liu, S. J., Liu, N. Y., He, P., Li, Z. Q., Dong, S. L., and Mu, L. F. (2012). Molecular characterization, expression patterns, and ligand-binding properties of two odorant-binding protein genes from Orthaga achatina (Butler) (Lepidoptera: Pyralidae). Arch. Insect Biochem. Physiol. 80, 123–139. doi: 10.1002/arch.21036
Luo, Z., and Honda, H. (2015). Olfactory and biophysical assessment of the oviposition stimulating potential of host and non-host plants for the yellow peach moth, Conogethes punctiferalis, (Lepidoptera: Crambidae). Appl. Entomol. Zool. 50, 183–189. doi: 10.1007/s13355-014-0320-9
Maibèche-Coisne, M., Longhi, S., Jacquin-Joly, E., Brunel, C., Egloff, M. P., Gastinel, L., et al. (1998). Molecular cloning and bacterial expression of a general odorant-binding protein from the cabbage armyworm Mamestra brassicae. Eur. J. Biochem. 258, 768–774. doi: 10.1046/j.1432-1327.1998.2580768.x
Mao, A., Zhou, J., Bin, M., Zheng, Y., Wang, Y., Li, D., et al. (2016). Sex pheromone recognition and characterization of three pheromone-binding proteins in the legume pod borer, Maruca vitrata Fabricius (Lepidoptera: Crambidae). Sci. Rep. 6:34484. doi: 10.1038/srep34484
Mohanty, S., Zubkov, S., and Gronenborn, A. M. (2004). The solution NMR structure of Antheraea polyphemus PBP provides new insight into pheromone recognition by pheromone-binding proteins. J. Mol. Biol. 337, 443–451. doi: 10.1016/j.jmb.2004.01.009
Sandler, B. H., Nikonova, L., Leal, W. S., and Clardy, J. (2000). Sexual attraction in the silkworm moth: structure of the pheromone-binding-protein-bombykol complex. Chem. Biol. 7, 143–151. doi: 10.1016/S1074-5521(00)00078-8
Sun, M. J., Liu, Y., and Wang, G. R. (2013). Expression patterns and binding properties of three pheromone binding proteins in the diamondback moth, Plutella xylostella. J. Insect Physiol. 59, 46–55. doi: 10.1016/j.jinsphys.2012.10.020
Sun, Y. L., Huang, L. Q., Pelosi, P., and Wang, C. Z. (2013). A lysine at the C-terminus of an odorant-binding protein is involved in binding aldehyde pheromone components in two Helicoverpa species. PLoS ONE 8:e55132. doi: 10.1371/journal.pone.0055132
Tamura, K., Peterson, D., Peterson, N., Stecher, G., Nei, M., and Kumar, S. (2011). MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol. Biol. Evol. 28, 2731–2739. doi: 10.1093/molbev/msr121
Thode, A. B., Kruse, S. W., Nix, J. C., and Jones, D. N. (2008). The role of multiple hydrogen-bonding groups in specific alcohol binding sites in proteins: insights from structural studies of LUSH. J. Mol. Biol. 376, 1360–1376. doi: 10.1016/j.jmb.2007.12.063
Vogt, R. G., Rybczynski, R., and Lerner, M. R. (1991). Molecular cloning and sequencing of general odorant-binding proteins GOBP1 and GOBP2 from the tobacco hawk moth Manduca sexta: comparisons with other insect OBPs and their signal peptides. J. Neurosci. 11, 2972–2984.
Wang, S. Y., Gu, S. H., Han, L., Guo, Y. Y., Zhou, J. J., and Zhang, Y. J. (2013). Specific involvement of two amino acid residues in cis-nerolidol binding to odorant-binding protein 5 AlinOBP5 in the alfalfa plant bug, Adelphocoris lineolatus (Goeze). Insect Mol. Biol. 22, 172–182. doi: 10.1111/imb.12012
Willett, C. S., and Harrison, R. G. (1999). Pheromone binding proteins in the European and Asian corn borers: no protein change associated with pheromone differences. Insect Biochem. Mol. 29, 277–284. doi: 10.1016/S0965-1748(99)00003-X
Xiu, W. M., and Dong, S. L. (2007). Molecular characterization of two pheromone binding proteins and quantitative analysis of their expression in the beet armyworm, Spodoptera exigua Hübner. J. Chem. Ecol. 33, 947–961. doi: 10.1007/s10886-007-9277-2
Xu, X., Xu, W., Rayo, J., Ishida, Y., Leal, W. S., and Ames, J. B. (2010). NMR structure of navel orangeworm moth pheromone-binding protein (AtraPBP1): implications for pH-sensitive pheromone detection. Biochemistry 49, 1469–1476. doi: 10.1021/bi9020132
Yang, M. W., Dong, S. L., and Chen, L. (2009). Electrophysiological and behavioral responses of female beet armyworm Spodoptera exigua (Hübner) to the conspecific female sex pheromone. J. Insect Behav. 22, 153–164. doi: 10.1007/s10905-008-9162-z
Yin, J., Zhuang, X., Wang, Q., Cao, Y., Zhang, S., Xiao, C., et al. (2015). Three amino acid residues of an odorant-binding protein are involved in binding odours in Loxostege sticticalis. Insect Mol. Biol. 24, 528–538. doi: 10.1111/imb.12179
Zhang, T., Gu, S., Wu, K., Zhang, Y., and Guo, Y. (2011). Construction and analysis of cdna libraries from the antennae of male and female cotton bollworms Helicoverpa armigera (Hübner) and expression analysis of putative odorant-binding protein genes. Biochem. Biophys. Res. Commun. 407, 393–399. doi: 10.1016/j.bbrc.2011.03.032
Zhou, J. J., Robertson, G., He, X., Dufour, S., Hooper, A. M., Pickett, J. A., et al. (2009). Characterisation of Bombyx mori odorant-binding proteins reveals that a general odorant-binding protein discriminates between sex pheromone components. J. Mol. Biol. 389, 529–545. doi: 10.1016/j.jmb.2009.04.015
Keywords: pheromone binding proteins, Conogethes punctiferalis, fluorescence competitive binding assays, molecular docking, qPCR
Citation: Ge X, Ahmed T, Zhang T, Wang Z, He K and Bai S (2018) Binding Specificity of Two PBPs in the Yellow Peach Moth Conogethes punctiferalis (Guenée). Front. Physiol. 9:308. doi: 10.3389/fphys.2018.00308
Received: 23 October 2017; Accepted: 14 March 2018;
Published: 03 April 2018.
Edited by:Nicolas Durand, Université Pierre et Marie Curie, France
Reviewed by:Guan-Heng Zhu, University of Kentucky, United States
Loic Briand, UMR6265 Centre des Sciences du Goût et de l'Alimentation (CSGA), France
Herbert Venthur, Universidad de La Frontera, Chile
Copyright © 2018 Ge, Ahmed, Zhang, Wang, He and Bai. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.