Binding Affinity Characterization of Four Antennae-Enriched Odorant-Binding Proteins From Harmonia axyridis (Coleoptera: Coccinellidae)

Harmonia axyridis is an important natural enemy that consumes many agricultural and forestry pests. It relies on a sensitive olfactory system to find prey and mates. Odorant-binding proteins (OBPs) as the first-step of recognizing volatiles, transport odors through sensillum lymph to odorant receptors (ORs). However, little is known about the molecular mechanisms of H. axyridis olfaction. In this study, four H. axyridis antenna specific OBP genes, HaxyOBP3, 5, 12, and 15, were bacterially expressed and the binding features of the four recombinant proteins to 40 substances were investigated using fluorescence competitive binding assays. Three-dimensional structure modeling and molecular docking analysis predicted the binding sites between HaxyOBPs and candidate volatiles. Developmental expression analyses showed that the four HaxyOBP genes displayed a variety of expression patterns at different development stages. The expression levels of HaxyOBP3 and HaxyOBP15 were higher in the adult stage than in the other developmental stages, and HaxyOBP15 was significantly transcriptionally enriched in adult stage. Ligand-binding analysis demonstrated that HaxyOBP3 and HaxyOBP12 only combined with two compounds, β-ionone and p-anisaldehyde. HaxyOBP5 protein displayed binding affinities with methyl salicylate, β-ionone, and p-anisaldehyde (Ki = 18.15, 11.71, and 13.45 μM). HaxyOBP15 protein had a broad binding profile with (E)-β-farnesene, β-ionone, α-ionone, geranyl acetate, nonyl aldehyde, dihydro-β-ionone, and linalyl acetate (Ki = 4.33–31.01 μM), and hydrophobic interactions played a key role in the binding of HaxyOBP15 to these substances according to molecular docking. Taken together, HaxyOBP15 exhibited a broader ligand-binding spectrum and a higher expression in adult stage than HaxyOBP3, 5, and 12, indicating HaxyOBP15 may play a greater role in binding volatiles than other three HaxyOBPs. The results will increase our understanding of the molecular mechanism of H. axyridis olfaction and may also result in new management strategies (attractants/repellents) that increase the biological control efficacy of H. axyridis.

It is helpful for identifying the function to study the OBPs expression patterns (Gong et al., 2014;Tang et al., 2019). OBPs have a variety of functions depending on their distribution (Sun et al., 2012;Xue et al., 2016;Li et al., 2019b). Antennaespecific OBPs play important roles in detecting sex pheromones and plant volatiles (Sun et al., 2014;Liu et al., 2020).
The Asian multicolored ladybird beetle, H. axyridis (Coleoptera: Coccinellidae), as an important natural enemy, can prey on many pests, including aphids, whiteflies, and thrips. Since the early 21 century, this species has been successfully used to control pests of crops (Koch, 2003;Pervez and Omkar, 2006;Wang et al., 2015). Harmonia axyridis is an effective biological control agent, but it can also be a pest in some situations (Soares et al., 2007;Koch and Costamagna, 2016;Ovchinnikov et al., 2019). It may compete with native predators for common food resources, and bring pollution to wine production (Pickering et al., 2004;Katsanis et al., 2013;Grez et al., 2016).
Predators used aphid alarm pheromones and pest-induced volatiles to locate pest (Al Abassi et al., 2000;Hatt et al., 2019), which is an important communication way of pestcrop-natural enemy interactions in agricultural fields. It is necessary for enhancing natural enemies' biological control efficacy to understand their olfactory systems. Therefore, the interaction of H. axyridis with plant volatiles and aphid pheromones may be important for enhancing the effectiveness of H. axyridis as a biological control agent.
We previously identified 19 putative OBPs and characterized their tissue expression patterns by quantitative real-time PCR (qRT-PCR) based on antennae and whole-body transcriptomes of H. axyridis (Qu et al., 2021). HaxyOBP3 (NCBI accession number MT150141), HaxyOBP5 (NCBI accession number MT150143), HaxyOBP12 (NCBI accession number MT150150), and HaxyOBP15 (NCBI accession number MT150153), specifically expressed in adult antennae, may play a more important role in the olfactory perception of H. axyridis. In the present study, these four antennae-specific OBPs were selected for detailed study. The development stage expression profiles of these genes were generated, and their binding characteristics to ligands were also conducted. In addition, protein structures were modeled in three dimensions, and their potential binding sites were studied by molecular docking. The results increase our comprehending of the molecular basis of olfaction of H. axyridis and may help to enhance their biological control effectiveness.

Insect Samples
Harmonia axyridis was obtained from Beijing Kuoye Tianyuan Biological Technology Co., Ltd., rearing in a growth chamber of the Beijing Academy of Agriculture and Forestry Sciences with the temperature of 23 ± 1°C, 16:8 h (L:D) photoperiod and 70% relative humidity. The adults and larvae were fed with aphid Aphis craccivora Koch (Qu et al., 2018). To determine the transcript levels of HaxyOBP3, 5, 12, and 15 under various developmental stages (eggs, first, second, third, and fourth instar, pupae, and male and female adults), samples were collected and stored at −80°C. Three biological replicates were conducted.

Specific Expression of OBP Genes
The TRIzol reagent (Invitrogen, Carlsbad, CA, United States) was used to extract total RNA samples based on the manufacturer's instructions. The first-strand cDNA was synthesized using the PrimeScript™ RT reagent Kit (TAKARA, Japan) following the provided protocol. The development stage expression pattern of HaxyOBPs was assessed by qRT-PCR. qRT-PCR was performed on ABI PRISM 7500 (Applied Biosystems, United States). The reaction consisted of 10 μl SYBR Premix Ex TaqTM II (TaKaRa, Japan), 1 μl of each primer (10 μmol L −1 ), 2 μl cDNA, 0.4 μl Rox Reference Dye II (Takara, Japan), and 5.6 μl nuclease free water. The reaction conditions were 95°C for 30 s, followed by 40 cycles of 95°C for 5 s and 60°C for 34 s. Primers of HaxyOBPs were based on Qu et al. (2021). EF1A and RPS13 genes were used as housekeeping genes (Qu et al., 2018). All samples were tested in three biological replicates. The 2 -△△CT method was used for relative quantification (Schmittgen and Livak, 2008). The differences in the transcript levels of HaxyOBPs in different developmental stages were compared by One-way ANOVA (SPSS 19.0, Chicago, IL, United States), followed by Tukey's test. Heat map illustrating the log 2 transformation of HaxyOBPs mRNA expression levels in different developmental stages.

Expression and Purification of Recombinant OBPs
The DNA sequences that encode the HaxyOBP3, 5, 12, and 15 proteins were chemically synthesized and cloned into pET30a (+) by GenScript (Nanjing, China; Wang et al., 2020b). The positive plasmid was then transformed into BL21 (DE3) cells for the expression of recombinant proteins, and proteins induced with 0.5 mmol/L isopropyl β-D-1-thiogalactopyranoside (IPTG) for 4 h at 37°C (HaxyOBP3) or 16 h at 15°C (HaxyOBP5, 12, and 15). HaxyOBP5 was expressed in the supernatant. HaxyOBP3, 12, and 15 were mainly found in inclusion bodies. Inclusion bodies were denatured by 8 M urea. Recombinant proteins of HaxyOBP3, 12, and 15 were dissolved and refolded based on the reported methods (Zhang et al., 2012).
Protein purification was performed with His-Tag Purification Resin column (Genscript Biology Company, Nanjing, Jiangsu, China) and purified by gradient imidazole buffer (20, 50, 100, 250, and 500 mmol·L −1 ). The purity and size of proteins were detected by SDS-PAGE, and concentrations of proteins were measured with bicinchoninic acid (BCA) Protein Assay Kit (ThermoFisher Scientific-Life Technologies, Carlsbad, CA, United States).

Competitive Fluorescence Binding Assay
A Cary Eclipse Fluorescence Spectrophotometer (Agilent Technologies, United States) was used to determine the results of the binding assay. N-phenyl-1-naphthylamine (1-NPN) for HaxyOBP15 and 4,4′-Dianilino-1,1′-binaphthyl-5,5′-disulfonic acid dipotassium salt (bis-ANS) for HaxyOBP3, 5, and 12 were chosen as the fluorescent probe. The excitation wavelength was 337 nm of 1-NPN and 295 nm of bis-ANS, and the emission spectrum was recorded between 350 and 500 nm for 1-NPN and between 300 and 550 nm for bis-ANS. The recombinant proteins prepared in Tris-HCl (50 mM, pH 7.4) was titrated with aliquots of 1 mM 1-NPN or bis-ANS to final concentrations ranging from 2 to 16 μM to measure the binding affinity. To further measure the binding affinity of ligands to HaxyOBPs, proteins and fluorescent probe at 2 μM were titrated with aliquots of 1 mM odorants. The binding constant (K 1-NPN/bis-ANS ) of 1-NPN or bis-ANS to HaxyOBPs was calculated by GraphPad

Three-Dimensional Modeling and Molecular Docking
Three-dimensional structure of HaxyOBP12 and HaxyOBP15, more than 30% homology with the OBP templates in the Protein Database, 1 was modeled by Program MODELLER (Martí-Renom et al., 2000;Fiser et al., 2010), while HaxyOBP3 and HaxyOBP5 that had less than 30% homology were also generated using a deep residual neural network trRosetta (https://yanglab.nankai.edu. cn/trRosetta; Yang et al., 2020). Three methods, including Verify_3D, Procheck, and ERRAT were used to assess the final 3D model of HaxyOBPs protein (Laskowski et al., 1996;Webb and Sali, 2016). AutoDock Vina (version 1.1.2) was selected to analyze the binding mode between HaxyOBPs protein and compounds with the default parameters (Morris et al., 2009). The top ranked binding mode was evaluated according to the Vina docking score, and visually analyzed by PyMOL (version 1.9.0; http://www.pymol.org/).

Developmental Stage Expression of HaxyOBPs
qRT-PCR was used to determine the expression levels of HaxyOBP3, 5, 12, and 15 in different developmental stages (Figure 1). HaxyOBP3 and HaxyOBP15 were both highly expressed in adults, and HaxyOBP15 had a significantly higher expression level in this stage. Transcripts of HaxyOBP5 were especially abundant in the first instar. In addition, HaxyOBP12 showed similar relative transcript levels in all developmental stages.

Expression and Purification of HaxyOBPs
The recombinant proteins of HaxyOBP3, 5, 12, and 15 were successfully expressed in the E. coli system induced by IPTG. HaxyOBP5 was mainly detected in the supernatant, while HaxyOBP3, 12, and 15 were present in inclusion bodies (Supplementary Figure S1). Therefore, 8 mol/L urea was used to extract the protein of HaxyOBP3, 12, and 15 before the purification. Renaturation, dialysis, and ultrafiltration were then used to obtain the purified target proteins of HaxyOBP3, 12, and 15. The OBPs were purified by nickel affinity chromatography. SDS-PAGE analysis revealed the final purified proteins as a single band, a molecular weight of about 15 kDa, consistent with the predicted molecular mass (Figure 2).

Binding Characteristics of HaxyOBPs
To determine the binding spectra of four HaxyOBPs recombinant proteins, fluorescence competitive binding assay was conducted. The binding characteristics of HaxyOBP3, 5, and 12 with fluorescent probe bis-ANS were detected by molecular fluorescence spectrometry, and the dissociation constants (K i value) of HaxyOBP3, HaxyOBP5, and HaxyOBP12 were 3.07 ± 057, 2.23 ± 0.36, and 1.84 ± 0.10 μM, respectively. Using the same method, we detected the binding characteristics of HaxyOBP15 with 1-NPN, and the dissociation constant was 5.07 ± 0.31 μM (Figure 3).

Homology Modeling and Molecular Docking
Sequence alignments showed that HaxyOBP12 and HaxyOBP15 share 44 and 31% amino acid identities with 6JPM and 4Z45, respectively. Sequence alignments showed that HaxyOBP3 and HaxyOBP5 share 29.27 and 28.93% amino acid identities, respectively, with the templates 1C3Y and 6QQ4, less than 30.00% ( Table 2). The low identity may decrease the accuracy of the predicted model. So, we used the other method, trRosetta, to predict the model of HaxyOBP3 and HaxyOBP5. The trRosetta can predict the protein more accurately for the low identity sequence. The models predicted by Homology modeling were named Mod-HaxyOBP12 and Mod-HaxyOBP15. The models predicted by trRosetta were named trR-HaxyOBP3 and trR-HaxyOBP5.
The protein structures of HaxyOBP3, 5, 12, and 15 were composed of six typical α-helices, forming a hydrophobic binding cavity, which are the important features of insect OBPs (Figure 5).

DISCUSSION
Odorant-binding proteins are the front-line environmental odorant sensors, playing an essential role in insect behavior (Pelosi et al., 2006). Temporal and spatial expression patterns of OBPs in insects are interrelated with their specific physiological functions (Song et al., 2014;Sun et al., 2016;Zhang et al., 2017b). The transcripts of HaxyOBP3, 5, 12, and 15 were mainly restricted to adult antennae in our previous study (Qu et al., 2021), implying a role of these proteins in olfactory chemoreception. Clarifying the expression characteristics of insects' OBPs at different developmental stages can also help to understand their functions in olfactory recognition (Ju et al., 2014). In this study, qRT-PCR indicated that HaxyOBP3, 5, 12, and 15 had different transcript levels during the different developmental stages of H. axyridis. HaxyOBP5 and HaxyOBP12 were abundant in the larval stage, indicating their connection to larval biological characteristics of H. axyridis. However, HaxyOBP3 and HaxyOBP15 were both highly expressed in adult stage, and the expression level of HaxyOBP15 was significantly higher in this stage, indicating they might be involved in adult-specific behaviors.
The results of the fluorescence binding assay also showed that HaxyOBP15 had a broader ligand-binding affinity, and it could bind seven substances including (E)-β-Farnesene, β-ionone, α-ionone, geranyl acetate, dihydro-β-ionone, nonyl aldehyde, and linalyl acetate, comparing with HaxyOBP3, 5, 12. These results were consistent with the fact that HaxyOBP15 gene showed significantly higher expression level in adult stage, indicating that HaxyOBP15 played a key role in olfactory communication of adult H. axyridis. Harmonia axyridis is an important natural enemy in many crops (Koch, 2003;Pervez and Omkar, 2006). Plant volatiles and sex pheromone are essential signal chemicals in pest-crop-natural enemy interactions Frontiers in Physiology | www.frontiersin.org 8 March 2022 | Volume 13 | Article 829766 FIGURE 6 | Binding pattern of HaxyOBP3 with p-anisaldehyde and β-ionone. . Compared with larvae, adults of H. axyridis have a wider range of activities due to their ability to fly, and they need to recognize more odorants, so as to detect food sources or find suitable oviposition sites. Homology modeling and molecular docking were used to further study the specific binding characteristics of OBPs. Classic OBPs usually has six α-helical domains, and fold together to form a compact pocket for combing odors (Pelosi et al., 2014). The present study showed that the predicted 3D structures of the four HaxyOBPs were consistent with those of classic OBPs, having six α-helical domains. Ligands are usually bound in a hydrophobic cavity of insect OBPs (Wogulis et al., 2006;Northey et al., 2016). In this study, molecular docking results showed that HaxyOBP15 broadly bound with more substances, suggesting that HaxyOBP15 may have adapted to binding to substances with different shapes and sizes. Some residues of HaxyOBP15 may specifically interact with functional groups of substances. For instance, HaxyOBP15 possess a key amino acid residue, H53, which appears to be involved in the recognition of a broad range of substances.
Among the 40 candidate compounds, nine compounds, including β-ionone, α-ionone, dihydro-β-ionone, geranyl acetate, nonanal, linalyl acetate, EBF, p-anisaldehyde, and methyl salicylate, bound to four HaxyOBPs according to fluorescence binding assay and molecular docking. β-ionone as a fragrance compound, existing in the flowers and fruits of many plants (Beekwilder et al., 2014;Fu et al., 2019;Guarino et al., 2021) and having a strong repellent effect on flea beetles, spider mites, and whiteflies (Caceres et al., 2016), could bind with the four HaxyOBPs proteins. Moreover, HaxyOBP15 could also bind with α-ionone and dihydro-β-ionone, the analogs of β-ionone. For HaxyOBP15, hydrophobic interactions played a key role in the binding of HaxyOBP15 to three substances based on molecular docking. Three substances have similar chemical structures, and also commonly exist in plant volatiles, playing an important role in interactions between plants and insects (Li et al., 2019a). For example, Dihydro-β-ionone is attractive to the crucifer flea beetle (Caceres et al., 2016). The bouquet of Philodendron adamantinum is mainly composed of dihydro-β-ionone, which can attract the beetle Erioscelis emarginata and promote the pollination process (Pereira et al., 2014). The α-ionone was widely used as a male attractant for Bractocera latifrons (Nishida et al., 2009).
HaxyOBP15 could also bind geranyl acetate, nonanal, and linalyl acetate, and H53 was the key residue between HaxyOBP15 and these three substances by molecular docking. Geranyl acetate is similar to (E)-ß-farnesene (EBF) in structure, but the polarity and hydrophilicity of two compounds are different. Geranyl acetate is an ester, and EBF is a hydrocarbon. Previous studies reported that geranyl acetate is a strong activator of OR5 of aphids and also has binding affinity to OBP3 and OBP7, demonstrating features shared by several other behaviorally active repellents (Zhang et al., 2017a). Linalyl acetate is a monoterpene ester, which can be isolated from essential oils of Chrysactinia mexicana, Lavandula angustifolia, and Thymus leptophyllus (Jan et al., 2016).
In summary, the expression levels of the four HaxyOBPs genes showed different in development stages, and HaxyOBP15 was significantly higher expressed in the adult of H. axyridis based on the results of qRT-PCR. Ligand binding assays and molecular docking demonstrated HaxyOBP15 exhibited high specificity for more substances, comparing with HaxyOBP3, 5, and 12, suggesting HaxyOBP15 may play the most prominent role in the olfactory chemoreception of H. axyridis. These results can provide insight into the mechanism of olfactory communication of H. axyridis and enhance the biological control effectiveness of H. axyridis.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.

AUTHOR CONTRIBUTIONS
CQ and CL conceived and designed the experiments. CQ performed the experiments. CQ and F-qL analyzed the data. Z-kY and CQ provided the homology modeling and molecular docking and wrote the initial manuscript. CQ and H-pZ prepared the figures. CL, X-lY, and SW edited and reviewed the manuscript. All authors accepted the final version of the manuscript.