- Case Report
- Open access
- Published:
Screening and identification of cyprinid herpesvirus 2 (CyHV-2) ORF55-interacting proteins by phage display
Virology Journal volume 20, Article number: 66 (2023)
Abstract
Background
Cyprinid herpesvirus 2 (CyHV-2) is a pathogenic fish virus belonging to family Alloherpesviridae. The CyHV-2 gene encoding thymidine kinase (TK) is an important virulence-associated factor. Therefore, we aimed to investigate the biological function of open reading frame 55 (ORF55) in viral replication.
Methods
Purified CyHV-2 ORF55 protein was obtained by prokaryotic expression, and the interacting peptide was screened out using phage display. Host interacting proteins were then predicted and validated.
Results
ORF55 was efficiently expressed in the prokaryotic expression system. Protein and peptide interaction prediction and dot-blot overlay assay confirmed that peptides identified by phage display could interact with the ORF55 protein. Comparing the peptides to the National Center for Biotechnology Information database revealed four potential interacting proteins. Reverse transcription quantitative PCR results demonstrated high expression of an actin-binding Rho-activating protein in the latter stages of virus-infected cells, and molecular docking, cell transfection and coimmunoprecipitation experiments confirmed that it interacted with the ORF55 protein.
Conclusion
During viral infection, the ORF55 protein exerts its biological function through interactions with host proteins. The specific mechanisms remain to be further explored.
Introduction
Herpesviral haematopoietic necrosis (HVHN), a highly lethal disease of crucian carp (Carassius carassius), is caused by cyprinid herpesvirus 2 (CyHV-2). This disease was first detected in goldfish (Carassius auratus) cultured in Aichi and Nara prefectures in Japan in 1992, with a mortality rate of 100% [1]. Since 2009, CyHV-2 has been detected in Gibel Carp (Carassius auratus gibelio) cultured in some areas of Jiangsu Province in China, and the virus rapidly became endemic throughout the main carp farming areas in China [2]. CyHV-2 has since dramatically spread across the world [3,4,5,6]. CyHV-2 is a member of the Alloherpesviridae family (genus Cyprinivirus), which also comprises cyprinid herpesvirus 3 (CyHV-3) and carp pox virus (CyHV-1) [7]. CyHV-2 is a DNA virus with a vesicular membrane and a large linear double-stranded DNA genome. The virions are oval, with a diameter of approximately 170 − 200 nm [8].
Developing the molecular epidemiological characterisation of CyHV-2 is an important basis for research on the infection mechanism, and for prevention and control of the virus. Some genes of CyHV-2 have been used for research on virus detection, preparation of antibody, and genotype distinction and diversity analysis [9,10,11]. Additionally, the CyHV-2 open reading frame 55 (ORF55) gene was predicted to encode thymidine kinase (TK), and enzyme participating in nucleic acid metabolism [12]. In previous studies, the function of TK was investigated in herpes simplex virus 1 (HSV-1) [13,14,15], duck plague virus (DPV) [16] and Epstein-Barr (EB) virus [17]. The TK gene has since been studied in Pseudorabies virus (PRV) [18]. The TK gene was not essential for the growth of PRV, but it was essential for viral proliferation, and it was also the most important gene determining PRV critical virulence. Additionally, deletion of the TK gene could significantly decrease the virulence of PRV, while also greatly reducing the replication, transmission and lethal ability of the virus in nerve tissues, but immunogenicity was not affected [19].
Thus, TK plays an important role during viral infection. Although the experimental evidence mentioned above linked the TK gene with virulence, its overall specific mechanism of action remains elusive. From previous experimental results, we deduced that ORF55 also has a virulence effect, and is an important virulence gene in CyHV-2. Therefore, further exploration of the ORF55 gene is important.
Phage display technology (PDT) is a powerful approach first reported in 1985 [20]. PDT is a well-established method for selecting peptides and proteins with novel binding functions [21]. A foreign protein is first expressed in fusion with one of the phage capsid proteins by cloning the foreign DNA sequence in frame with the gene encoding the chosen capsid protein [22]. PDT allows a direct link and identification between a large number of random peptides and DNA-coding sequences. It is used extensively in medicine, biology and chemistry [23, 24]. However, how the CyHV-2 ORF55 protein performs its biological functions in goldfish cells remains poorly understood. Screening and identifying peptides or proteins that may interact with the ORF55 protein could help to elucidate its mechanism of action.
In this study, we cloned the ORF55 gene into a prokaryotic expression vector using seamless cloning, and used PDT to screen interacting partners. Bioinformatics prediction identified four potential interacting proteins. This study lays a foundation for further exploration of the biological functions of the ORF55 protein during viral infection.
Materials and methods
Bacterial strains, cloning vectors, cells and virus
Escherichia coli DH5α and E. coli BL21 (DE3) strains were obtained from TransGen Biotech (Beijing, China). The prokaryotic expression vector pMAL-c5X was purchased from New England Biolabs (NEB; Beijing, China). The CyHV-2 YC-01 strain (GenBank accession No. MN593216; 275,367 bp) was isolated and purified from Crucian Carp showing obvious symptoms of the early stages gill haemorrhage in the laboratory at Shanghai Ocean University, and preserved at − 80 °C. A cell line from the caudal fin of C. auratus gibelio (GiCF) [25] was established in previous experiments and cultured in M199 (Genomcell, Zhejiang, China) medium supplemented with 12% fetal bovine serum (FBS; Every Green, Zhejiang, China) at 27 °C. This cell line could be stably infected with CyHV-2 and used for in vitro virus replication. HEK-293 T cells were cultured in Dulbecco’s modified Eagle’s medium with high glucose containing 10% FBS in a temperature-controlled incubator at 37 °C with 5% CO2.
Sequence-based prediction
The National Center for Biotechnology Information (NCBI) database of the National Library of Medicine was used to obtain the sequence of the TK gene from the CyHV-2 YC-01 strain. The ExPASy ProtParam tool (https://web.expasy.org/protparam/), a server that calculates theoretical protein parameters including molecular weight, amino acids, isoelectric point (pI) and instability index, was used assess the protein physical and chemical properties [26]. In order to explore the functions of the ORF55 protein, we applied different in silico approaches using Cell-PLoc (http://www.csbio.sjtu.edu.cn/bioinf/Cell-PLoc/), PSORTb (https://www.psort.org/psortb/), lncLocator (http://www.csbio.sjtu.edu.cn/bioinf/lncLocator/), iLoc-LncRNA (http://lin-group.cn/server/iLoc-LncRNA/home.php) and Euk-mPLOC (http://www.csbio.sjtu.edu.cn/bioinf/euk-multi-2/) to predict localisation and other characteristics [27,28,29,30,31]. The secondary structure of the ORF55 protein sequence was predicted by various web servers including SOPMA (https://npsa-pbil.ibcp.fr/cgi-bin/npsa_automat.pl?page=npsa_sopma.html), NovoPro (https://novopro.cn/tools/secondary-structure-prediction.html), GOR (http://gor.bb.iastate.edu/) and PSSpred (https://zhanggroup.org/PSSpred/) [32,33,34,35]. Also, we downloaded the TK gene sequences of other species and performed evolutionary tree analysis using MEGA version 7.0. For homology analysis as well as sequence alignment, TK genes of all published CyHV-2 strains and representative CyHV-1 and CyHV-3 strains were selected for comparison. Protein tertiary structure was predicted using I-TASSER (https://zhanggroup.org/I-TASSER/).
Seamless cloning
The ORF55 gene was amplified using PCR with sense and antisense primers 5′-GCGATATCGTCGACGGATCCATGGCTTTTCTGGAGTTGGTGCTGG-3′ and 5′-CCTGCAGGGAATTCGGATCCTCATACAAAAGAACAAGGGGCATCC-3′. Restrictions sites for BamHI were inserted in both primes. Thermal cycling included an initial denaturation step at 98 °C for 5 min, followed by 34 cycles of denaturation at 94 °C for 30 s, annealing at 58 °C for 30 s, and extension at 72 °C for 30 s, and a final extension at 72 °C for 5 min. The subsequent ligation reaction was carried out at 50 °C according to the pEASY-Basic Seamless Cloning and Assembly Kit (TransGen Biotech), after which ligation products were transformed into E. coli DH5α. Positive colonies were screened and plasmid was isolated and sent to Sangon (Shanghai, China) for DNA sequencing.
Protein purification
Transformed E. coli BL-21 (DE3) cells were grown in Luria Bertani (LB) broth media containing 100 µg/mL ampicillin, and isopropylthiogalactoside (IPTG) was added to induce protein expression. Cells were collected by centrifugation and resuspended in 40 mL of 30 mM Tris–HCl, 20% sucrose, 1 mM EDTA, and mixed with 40 mL MgSO4. Following complete ultrasonic lysis and clarification by centrifugation, the bacterial supernatant was incubated with Amylose Resin (Sangon) on ice and transferred to an affinity chromatography column. After column washing, protein was eluted in elution buffer (20 mM Tris–HCl, 0.2 mol/L NaCl and 1 mM EDTA) containing 10 mM maltose.
Western blotting
The bacterial solution before and after induction, as well as the purified protein, were separated by 12% SDS-PAGE at 80 V for 30 min followed by 120 V for 60 min. Gels were stained with Coomassie Brilliant Blue. For western blotting, purified protein was subjected to electroblotting on a polyvinylidene difluoride (PVDF) membrane, which was subsequently blocked with 5% nonfat milk. The membrane was incubated with anti-MBP-Tag monoclonal antibody as primary antibody (Proteintech, USA) then with goat anti-mouse IgG (H + L)-horseradish peroxidase (HRP) conjugate (Bioworld Technology, USA) as secondary antibody. Results were detected using an enhanced chemiluminescent (ECL) detection system.
Phage display
The ORF55 protein concentration was determined from a standard curve prepared using bovine serum albumin (Sangon) and a Bicinchoninic Acid (BCA) Protein Quantification Kit (Beijing Solarbio Science & Technology Co. Ltd., Beijing, China). Protein was seeded in a 12-well plate at a concentration of 100 μg/mL, with three replicates for each experimental group. The phage library (NEB) was bound to the ORF55 protein according to the instruction manual of the Random Ph.D.-12 Phage Display Peptide Library Kit. Following three rounds of panning, the plate with plaques was sent to Genewiz Biotechnology Co. Ltd. (Tianjin, China) for sequencing and identification. The plaque sequencing were compared with NCBI crucian carp and goldfish databases.
Interaction validation
As described above, we obtained polypeptides that may interact with the ORF55 protein through PDT. The molecular docking tool from Home for researchers was used to predict interactions between the ORF55 protein and polypeptides, and to confirm their specific interaction sites. Dot blot experiments were also performed to verify interactions. Serially diluted peptides were blotted onto a PVDF membrane, soaked in ORF55 protein solution, and successively incubated with primary and secondary antibodies as described above for western blotting.
Protein prediction
Sequences of the four potential interacting proteins were downloaded from NCBI databases, and their sequences in other hosts were obtained for alignment, tertiary structure prediction, evolutionary tree analysis, and prediction of other physicochemical properties.
Preliminary in vitro validation
GiCF cells were infected with the supernatant of CyHV-2 YC-01 and sampled at various timepoints to observe whether the virus solution could alter expression of the above genes. RNA was extracted from sampled cells using TRIzol reagent and reverse-transcribed to cDNA as a template for gene amplification. Gene expression was measured using PrimeScript RT Master Mix (Takara, Beijing, China) from amplified cDNA. A real-time system instrument was used to quantify gene expression levels. Reactions contained 6.5 µL of TB Green Premix Ex Taq II (2 × ; Takara), 0.5 µL forward/reverse primers, and 500 ng cDNA. Reactions were repeated in triplicate.
Subcellular localisation analysis
Based on the above results, two proteins most likely to interact with the ORF55 protein were selected for further validation. Eukaryotic expression plasmids were constructed using seamless cloning technology as described above. Two luciferase reporter genes were selected as vectors. Unlike prokaryotic expression, a free endotoxin plasmid extraction kit was required for amplification. Recombinant plasmids were co-transfected into HEK-293Â T cells. After 40Â h, cell nuclei were stained using DAPI, and photographed using a confocal microscope.
Coimmunoprecipitation analysis
Samples of cells co-transfected as described above were collected to determine whether proteins were successfully expressed. Cells were lysed using NP40 lysate, and green fluorescent protein (GFP) fusion-tagged proteins and their interacting partners were selected by rabbit anti-GFP immunomagnetic beads. Western blotting was performed to probe interactions between proteins using anti-RFP Tag mouse monoclonal antibody as primary antibody as described above.
Statistical analysis
Statistical analysis was performed using GraphPad Prism (https://www.graphpad.com/scientific-software/prism/) and IBM SPASS statistics (SPASS 19.0). Results are presented as mean ± standard error (SE). Statistical significance was calculated using one-way analysis of variance (ANOVA). Values were considered significant (*) for p-values of 0.01 to 0.05, very significant (**) for p-values of 0.001 to 0.01, and extremely significant (***) for p-values 0 to 0.001.
Results
Physiochemical parameters
The physicochemical properties of the CyHV-2 YC-01 ORF55 protein were predicted using the Protparam web tool (Table 1). The predicted pI of ORF55 was < 7, suggesting it is acidic. The instability index value indicated that it is an unstable protein in CyHV-2 and CyHV-3, but stable in CyHV-1. The calculated aliphatic index indicated that it was thermally over a wide temperature range. ORF55 was classified as polar based on the grand average hydropathy (GRAVY) value. Knowing where the ORF55 protein is localised is critical to predicting its functions. All different algorithm-based servers predicted that the ORF55 protein is in the cytoplasm (Table 2). Almost all servers predicted that the ORF55 protein is rich in the alpha-helix secondary structure, which is likely important for protein structure and function (Table 3).
Bioinformatics analysis of the ORF55 protein
A phylogenetic tree of TKs from different viruses was constructed using the neighbour-joining method in MEGA 7.0 software, revealing the evolutionary journey and relatedness of TK genes (Fig. 1A). Homology analysis in MegAlign 7.1 showed that all CyHV-2 strain types shared 100% sequence identity and > 50% identity with CyHV-1 and CyHV-3 (Additional file 2: Fig. 1G). All CyHV-2 strain types shared 98.6% similarity (Fig. 1B). Among them, ST-J1 and YC-01 strains were completely identical, with nine more bases (GGACGAGGA) than other strains. The results indicate that the ORF55 gene is highly conserved between CyHV-2 strains. The protein tertiary structure analysis (Fig. 1C) showed that there was a Pfam: TK domain located at 2–176, which laid a foundation for subsequent in-depth analysis of the reciprocal binding sites.
Validation of protein expression
PCR results showed that the ORF55 gene was 645 bp in length, and 904 bp after connecting the prokaryotic expression vector, as expected (Fig. 2A). Expression of the ORF55 recombinant plasmid resulted in protein expression in E. coli BL-21 (DE3) following induction with 0.5 mM IPTG (Fig. 2B). SDS-PAGE results revealed a single band corresponding to MBP-tagged recombinant ORF55 protein, with molecular weights of 45 kDa for MBP and 70 kDa for ORF55 (Fig. 2C). Western blotting displayed a single band for purified recombinant ORF55 protein (Fig. 2D).
Screening of interacting polypeptides
A standard curve was established using BSA and the ORF55 protein concentration was measured at 562 nm using a spectrophotometer (Additional file 2: Fig. 2). During the three rounds of phage screening, the screening conditions were optimised by gradually reducing the amount of coated protein, increasing the concentration of Tween-20 and washing with detergent, and reducing the incubation duration with phage to obtain positive clones with high affinity (Table 4A). To determine phage enrichment for specific binding to the ORF55 protein in each round of screening requires a quantifiable indicator, in addition to the visual counting of phage plaques obtained in each round of phage panning (Fig. 3), hence phage titres after each round of screening were measured. The results showed that the correlation between the number of selection rounds and the titres of phages bound to ORF55 protein was positive (Table 4B). The recovery yield of phages was increased approximately 84-fold after the third round of panning compared with the first round. This indicates an obvious enrichment in the specific binding capability of phages.
Prediction of interacting proteins
After the third round of screening, 50 plaque samples were randomly selected and sent to Genewiz for DNA sequencing. Based on the resulting data, nucleotide sequences were translated into amino acid sequences using DNAMAN (version 6). The results showed that the polypeptide sequence LSPGANSHVSRH had the highest homology, with 38 repetitions (Table 4C), suggesting it binds strongly to the ORF55 protein.
Validation of polypeptide-protein interactions
Through analysis of the phage display results described above, we obtained a polypeptide sequence of 12 amino acids. In this experiment, we used the dot blot method for validation. When the concentration of the ORF55 protein was > 0.5 mg/mL, the dot blot results showed a significant band, with the MBP protein tag alone serving as a negative control (Fig. 4A).
In order to further understand the role of the ICG (ligand) in nanostructure formation, ORF55 was used as a macromolecule and the polypeptide was used as a ligand for in silico analysis of the interaction between them during NP formation. The docking results revealed active binding between the ORF55 protein and the polypeptide. Specifically, binding between the two was due to interactions between a lysine residue of the ORF55 protein and a sulphonate group of the polypeptide (Fig. 4B). It was found that the docked ligand is close to the lysine and lies within hydrogen bonding distance of the sulphonate group. The distance from the lysine residue of PLL to the adjacent oxygen atom of ICG is in the range of 1.8 Å to 2.2 Å. Docking studies revealed that the docked ICG forms two hydrogen bonds with the lysine, with a bond length < 2.2 Å, suggesting possible aggregation of the ICG molecule within PLL/salt aggregates.
Bioinformatics analysis of potential interacting proteins
Amino acid sequences were analysed by alignment to the NCBI crucian carp and goldfish databases, and the ORF55 protein was predicted to interact with four proteins (Table 4D). They were gastrula zinc finger protein XlCGF8.2DB-like (ZFP, GenBank accession NO. XP_026094447.1), WD repeat-containing protein 7 isoform X1 (WDR7, XP_026052722.1), actin-binding Rho-activating protein (ABRA, XP_026140919.1) and zinc finger protein 516-like Isoform X1 (ZFP516, XP_026087732.1). A brief analysis of the respective physical and chemical parameters of the above proteins was performed (Additional file 1: Table S1). The Swiss-model (https://swissmodel.expasy.org/) was used to predict the tertiary structures (Fig. 5) and the evolutionary tree (Additional file 2: Fig. 3G) of the four proteins. The results showed that all four shared at least 70% sequence similarity between homologs in different species, indicating high conservation between their respective gene families. The results of conserved domain analysis showed that all four proteins had at least one domain or superfamily structure (Additional file 1: Table S2) potentially playing a specific role in the process of viral infection.
RT-qPCR verification
GiCF cells were infected with the supernatant of CyHV-2 YC-01 and sampled at different timepoints (0 h, 12 h, 24 h, 36 h and 48 h). A significant CPE was observed under the microscope ~ 48 h after viral infection, and all cells had died ~ 144 h after infection (Fig. 6A). According to the nucleic acid sequences published by NCBI, RT-qPCR primers were designed (Table 5). The results showed that the relative expression of ZFP, ABRA and ZFP516 was upregulated in the latter stages of viral infection, without a significant change in WDR7 expression at any of the timepoints (Fig. 6B − E). Melting curves of RT-qPCR experiments are shown in Additional file 2: Fig. 4. Additionally, viral titres were tested, and the results showed a significant increase with increasing time, as predicted (Fig. 6F).
Subcellular localisation
Through the results of the relative expression of four genes described above, we inferred that ABRA might interact with the ORF55 protein, and attempted to validate this. Both pEGFP-N1 and pDsRed-Express-N1 were chosen to connect ORF55 and ABRA genes, using primers listed in Table 6. Both were cotransfected into HEK-293 T cells, and RFP empty vector served as a control. After 40 h, cells were transferred to a confocal microscope, observed, and photographed. The results showed that the ORF55 protein colocalised with ABRA in the nuclei of 293 T cells, but it did not colocalise with RFP, indicating that the RFP tag and ORF55 protein were located in different organelles, respectively. (Fig. 7).
Coimmunoprecipitation
Cell samples after successful transfection as described above were analysed by western blotting. The experiment was divided into two parts; the first set of samples contained RFP and GFP-ORF55 protein, and the second set of samples contained RFP-ABRA and GFP-ORF55 protein. Both were recognised by GFP antibodies in the input (IP) group, while the RFP antibody only recognised 26 kDa and 43 kDa proteins, respectively. The IP group confirmed that the GFP magnetic beads had successfully adsorbed the GFP-ORF 55 protein. The final coimmunoprecipitation results revealed a clear band at 43 kDa in the second set of samples, implying an interaction between the GFP-ORF55 protein and the RFP-ABRA protein, but no interaction with the RFP protein (Fig. 8A).
Furthermore, we used molecular docking to verify the interaction between the ORF55 protein and ABRA, and to make specific binding site predictions. ZDOCK Score values and their best pose interactions were calculated (Table 7). PROK2 forms hydrogen bonds with amino acids including ARG167-ASP109 and ARG167-ASP112. Comprehensive analysis revealed that proteins PROK2 and CSF3 formed stable protein docking models (Fig. 8B).
Discussion
Outbreaks of CyHV-2 with a mortality rate of > 90% and no effective treatment measures have been reported in many countries and regions. There are three virus types in the Alloherpesviridae family, among which CyHV-1 is the least harmful, hence it is rarely reported. CyHV-3, also known as Koi herpesvirus (KHV) [36], causes a large number of koi deaths worldwide, hence it has received international attention. There have been some studies on virulence genes of KHV [37], but research on CyHV-2 protein is limited, and the main focus has been on the identification and functions of major membrane proteins and capsid proteins, as well as immunological detection methods [38]. The function and virulence mechanism of the CyHV-2 virulence remain unknown.
In the field of virus research, the E. coli prokaryotic expression system has many advantages including relatively low cost and simple operation, hence it is widely employed, as demonstrated for Duck Tembusu virus (DTMUV) non-structural protein 2A (NS2A) [39]. The CyHV-2 protein received some attention [40]. Therefore, we wished to set up a prokaryotic expression system and purify recombinant ORF55 in the present work. According to the preliminary analysis of physicochemical parameters, the ORF55 protein was unstable and not highly hydrophilic (Table 1), indicating that the final purified protein may not be highly soluble and/or easily degraded. Based on these analyses, the pMAL-c5X vector was selected for gene cloning because MBP is known to be highly soluble [41] and can enhance the biological activity of fused target proteins [42]. We successfully cloned the ORF55 gene into the pMAL-c5X vector, and expressed and purified the MBP-ORF55 fusion protein.
We then explored the reference lists of relevant studies on potential interacting proteins (Table 4D). Zinc finger proteins (ZFPs) are a class of proteins that contains short, stable, self-folding ‘finger’ structures that bind zinc ions. Studies have shown that zinc-finger proteins are natural host endogenous immune factors that regulate expression of cytokines such as interleukins and interferons, and directly regulate lymphocyte growth and differentiation [43]. The human specificity protein 1 (SP1) transcription factor belongs to the C2H2-type zinc finger family, and has been widely studied in HeLa cells (human cervical cancer cells) [44]. This protein regulates gene expression by interacting with other proteins, acting as a negative or positive regulator of gene expression [45]. In addition, a class of zinc-finger antiviral proteins has been widely studied in land animal viruses, including Small Ruminant Morbillivirus (SRMV) [46], Hepatitis E virus (HEV) [47] and Human immunodeficiency virus [48]. Members of this class of proteins can have specific inhibitory effects on viral replication.
Actin-binding Rho-activating Protein (ABRA, also known as Stars) is a muscle-specific actin-binding protein, a novel and evolutionarily conserved actin-binding protein (ABP) that binds to I-Band in sarcomeres and actin filaments in transfected cells to activate Rho signalling events [49]. ABRA may be involved in regulating actin cytoskeleton functions and actin polymerisation in various cells. Actin is a globular multifunctional protein that forms microfilaments, present in essentially all eukaryotic cells, and is the most abundant cytoskeletal protein in eukaryotic cells. Actin plays a regulatory role in viral infection, and is closely related to virus adsorption and entry, virus replication, virus assembly, virus release, and virus transmission. Studies have shown that actins including Profilins [50], Cofilin [51, 52] and Drebrin [53] can play a key role in viral infection [54]. Rho family proteins were the first cloned proteins in the Ras superfamily; they have GTP hydrolysis activity and play an important role in the regulation of cytoskeletal recombination. Rho is highly expressed in various malignancies and is closely related to tumorigenesis, invasion and metastasis [55].
We also investigated WD repeat-containing protein 7 (WDR7). F-box and WD repeat domain containing 7 (FBXW7) is a member of the F-box protein family that functions as the substrate recognition component of the SCF E3 ubiquitin ligase. FBXW7 is a critical tumour suppressor and is able to control proteasome-mediated degradation of oncoproteins [56]. Herein, the WDR7 gene did not show an obvious change in expression, indicating that it might not be altered during viral infection of cells.
In this experiment, PDT was used to explore a polypeptides interacting with the ORF55 protein, and predict potential interacting host proteins. Putative peptides can then be tested in cells to assess inhibitory or promoting effects on the virus, and further explore its mechanism. Through molecular docking and in vitro experiments, we determined that ABRA could interact with the ORF55 protein. In future work, the mechanism of interaction will be further investigated. For example, RNA interference (RNAi) and overexpression can be used to determine the nature of interactions, and signalling pathways can be studied by transcriptome analysis. For the other three identified interacting proteins, other methods can employed to probe interactions with the ORF55 protein. This screening results provide a basis for exploring the specific biological functions of the ORF55 protein in the mechanism of virus infection.
Conclusion
In summary, we successfully constructed a prokaryotic expression system for the CyHV-2 YC-01 TK gene, and applied the resulting purified ORF55 protein in phage display experiments. A polypeptide sequence interacting with the ORF55 protein was successfully screened, and four potential interacting proteins were predicted using the NCBI database. Following validation by RT-qPCR and coimmunoprecipitation experiments, we determined that the ORF55 protein could interact with ABRA. The findings lay a foundation for the further investigation of the infection mechanism of ORF55.
Availability of data and materials
Not applicable.
References
Jung SJ, Miyazaki T. Herpesviral haematopoietic necrosis of goldfish, Carassius auratus (L.). J Fish Diseases. 1995;18(3):211–20.
Wu T, Ding Z, Ren M, An L, Xiao Z, Liu P, Gu W, Meng Q, Wang W. The histo- and ultra-pathological studies on a fatal disease of Prussian carp (Carassius gibelio) in mainland China associated with cyprinid herpesvirus 2 (CyHV-2). Aquaculture. 2013;412–413:8–13. https://doi.org/10.1016/j.aquaculture.2013.07.004.
Goodwin AE, Merry GE, Sadler J. Detection of the herpesviral hematopoietic necrosis disease agent (Cyprinid herpesvirus 2) in moribund and healthy goldfish: validation of a quantitative PCR diagnostic method. Dis Aquat Org. 2006;69:137–43. https://doi.org/10.3354/dao069137.
Jeffery KR, Bateman K, Bayley A, Feist SW, Hulland J, Longshaw C, Stone D, Woolford G, Way K. Isolation of a cyprinid herpesvirus 2 from goldfish, Carassius auratus (L.), in the UK. J Fish Dis. 2007;30:649–56. https://doi.org/10.1111/j.1365-2761.2007.00847.x.
Boitard PM, Baud M, Labrut S, de Boisseson C, Jamin M, Bigarre L. First detection of cyprinid herpesvirus 2 (CyHV-2) in goldfish (Carassius auratus) in France. J Fish Dis. 2016;39:673–80. https://doi.org/10.1111/jfd.12400.
Giovannini S, Bergmann SM, Keeling C, Lany C, Schütze H, Schmidt-Posthaus H. Herpesviral hematopoietic necrosis in goldfish in Switzerland: early lesions in clinically normal goldfish (Carassius auratus). Vet Pathol. 2016;53:847–52. https://doi.org/10.1177/0300985815614974.
Tang R, Lu L, Wang B, Yu J, Wang H. Identification of the immediate-early genes of cyprinid herpesvirus 2. Viruses-Basel. 2020. https://doi.org/10.3390/v12090994.
Xu J, Zeng L, Zhang H, Zhou Y, Ma J, Fan Y. Cyprinid herpesvirus 2 infection emerged in cultured gibel carp, Carassius auratus gibelio in China. Vet Microbiol. 2013;166:138–44. https://doi.org/10.1016/j.vetmic.2013.05.025.
Marsh M, Helenius A. Virus entry: open sesame. Cell. 2006;124:729–40. https://doi.org/10.1016/j.cell.2006.02.007.
Kong SY, Jiang YS, Wang Q, Lu JF, Xu D, Lu LQ. Detection methods of Cyprinid herpesvirus 2 infection in silver crucian carp (Carassius auratus gibelio) via a pORF72 monoclonal antibody. J Fish Dis. 2017;40:1791–8. https://doi.org/10.1111/jfd.12648.
Zhou Y, Jiang N, Ma J, Fan Y, Zhang L, Xu J, Zeng L. Protective immunity in gibel carp, Carassius gibelio of the truncated proteins of cyprinid herpesvirus 2 expressed in Pichia pastoris. Fish Shellfish Immunol. 2015;47:1024–31. https://doi.org/10.1016/j.fsi.2015.11.012.
Davison AJ, Kurobe T, Gatherer D, Cunningham C, Korf I, Fukuda H, Hedrick RP, Waltzek TB. Comparative genomics of carp herpesviruses. J Virol. 2013;87:2908–22. https://doi.org/10.1128/jvi.03206-12.
Graessmann A, Sandberg G, Guhl E, Graessmann M. Methylation of single sites within the herpes simplex virus tk coding region and the simian virus 40 T-antigen intron causes gene inactivation. Mol Cell Biol. 1994;14:2004–10. https://doi.org/10.1128/mcb.14.3.2004.
Hay KA, Gaydos A, Tenser RB. The role of herpes simplex thymidine kinase expression in neurovirulence and latency in newborn vs adult mice. J Neuroimmunol. 1995;61(1):41–52. https://doi.org/10.1016/0165-5728(95)00071-9.
Yi QY, Bai ZS, Cai B, Chen N, Chen LS, Yuan T, Mao JH. HSV-TK/GCV can induce cytotoxicity of retinoblastoma cells through autophagy inhibition by activating MAPK/ERK. Oncol Rep. 2018;40:682–92. https://doi.org/10.3892/or.2018.6454.
Wen Y, Cheng A, Wang M, Ge H, Shen C, Liu S, Xiang J, Jia R, Zhu D, Chen X, et al. A Thymidine Kinase recombinant protein-based ELISA for detecting antibodies to duck plague virus. Virol J. 2010. https://doi.org/10.1186/1743-422x-7-77.
Gill MB, Kutok JL, Fingeroth JD. Epstein-Barr virus thymidine kinase is a centrosomal resident precisely localized to the periphery of centrioles. J Virol. 2007;81:6523–35. https://doi.org/10.1128/jvi.00147-07.
Zhang L, Ruan K, Sang G, Xu Z, Tong W, Yu H, Shan T, Gao F, Li L, Kong N, et al. Tk-deleted pseudorabies virus retains high pathogenicity in rats. J Vet Res. 2021;65:401–5. https://doi.org/10.2478/jvetres-2021-0056.
Zhao Y, Wang LQ, Zheng HH, Yang YR, Liu F, Zheng LL, Jin Y, Chen HY. Construction and immunogenicity of a gE/gI/TK-deleted PRV based on porcine pseudorabies virus variant. Mol Cell Probes. 2020;53:101605. https://doi.org/10.1016/j.mcp.2020.101605.
Smith GP. Filamentous fusion phage: novel expression vectors that display cloned antigens on the virion surface. Science. 1985;228:1315–7. https://doi.org/10.1126/science.4001944.
Dias RL, Fasan R, Moehle K, Renard A, Obrecht D, Robinson JA. Protein ligand design: from phage display to synthetic protein epitope mimetics in human antibody Fc-binding peptidomimetics. J Am Chem Soc. 2006;128:2726–32. https://doi.org/10.1021/ja057513w.
Mandrup OA, Friis NA, Lykkemark S, Just J, Kristensen P. A novel heavy domain antibody library with functionally optimized complementarity determining regions. PLoS ONE. 2013;8:e76834. https://doi.org/10.1371/journal.pone.0076834.
Lee KJ, Lee JH, Chung HK, Choi J, Park J, Park SS, Ju EJ, Park J, Shin SH, Park HJ, et al. Novel peptides functionally targeting in vivo human lung cancer discovered by in vivo peptide displayed phage screening. Amino Acids. 2015;47:281–9. https://doi.org/10.1007/s00726-014-1852-6.
Wang G, Yin P, Wang J, Ma P, Wang Y, Cai Y, Qi H, Liu A. Specific heptapeptide screened from pIII phage display library for sensitive enzyme-linked chemiluminescence immunoassay of vascular endothelial growth factor. Sensors Actuators B-Chem. 2021. https://doi.org/10.1016/j.snb.2021.129555.
Lu J, Xu D, Lu L. A novel cell line established from caudal fin tissue of Carassius auratus gibelio is susceptible to cyprinid herpesvirus 2 infection with the induction of apoptosis. Virus Res. 2018;258:19–27. https://doi.org/10.1016/j.virusres.2018.09.010.
Wilkins MR, Gasteiger E, Bairoch A, Sanchez JC, Williams KL, Appel RD, Hochstrasser DF. Protein identification and analysis tools in the ExPASy server. Methods Mol Biol. 1999;112:531–52. https://doi.org/10.1385/1-59259-584-7:531.
Chou KC, Shen HB. Cell-PLoc: a package of Web servers for predicting subcellular localization of proteins in various organisms. Nat Protoc. 2008;3:153–62. https://doi.org/10.1038/nprot.2007.494.
Yu NY, Wagner JR, Laird MR, Melli G, Rey S, Lo R, Dao P, Sahinalp SC, Ester M, Foster LJ, et al. PSORTb 3.0: improved protein subcellular localization prediction with refined localization subcategories and predictive capabilities for all prokaryotes. Bioinformatics. 2010;26:1608–15. https://doi.org/10.1093/bioinformatics/btq249.
Cao Z, Pan X, Yang Y, Huang Y, Shen HB. The lncLocator: a subcellular localization predictor for long non-coding RNAs based on a stacked ensemble classifier. Bioinformatics. 2018;34:2185–94. https://doi.org/10.1093/bioinformatics/bty085.
Su ZD, Huang Y, Zhang ZY, Zhao YW, Wang D, Chen W, Chou KC, Lin H. iLoc-lncRNA: predict the subcellular location of lncRNAs by incorporating octamer composition into general PseKNC. Bioinformatics. 2018;34:4196–204. https://doi.org/10.1093/bioinformatics/bty508.
Chou KC, Shen HB. A new method for predicting the subcellular localization of eukaryotic proteins with both single and multiple sites: Euk-mPLoc 20. PLoS ONE. 2010;5:e9931. https://doi.org/10.1371/journal.pone.0009931.
Geourjon C, Deléage G. SOPMA: significant improvements in protein secondary structure prediction by consensus prediction from multiple alignments. Comput Appl Biosci. 1995;11:681–4. https://doi.org/10.1093/bioinformatics/11.6.681.
Jones DT. Protein secondary structure prediction based on position-specific scoring matrices. J Mol Biol. 1999;292:195–202. https://doi.org/10.1006/jmbi.1999.3091.
Sen TZ, Jernigan RL, Garnier J, Kloczkowski A. GOR V server for protein secondary structure prediction. Bioinformatics. 2005;21:2787–8. https://doi.org/10.1093/bioinformatics/bti408.
Yan R, Xu D, Yang J, Walker S, Zhang Y. A comparative assessment and analysis of 20 representative sequence alignment methods for protein structure prediction. Sci Rep. 2013;3:2619. https://doi.org/10.1038/srep02619.
Hedrick RP, Gilad O, Yun S, Spangenberg JV, Marty GD, Nordhausen RW, Kebus MJ, Bercovier H, Eldar A. A herpesvirus associated with mass mortality of Juvenile and Adult Koi, a strain of common carp. J Aquat Anim Health. 2000;12:44–57. https://doi.org/10.1577/1548-8667(2000)012%3c0044:Ahawmm%3e2.0.Co;2.
Boutier M, Gao Y, Vancsok C, Suárez NM, Davison AJ, Vanderplasschen A. Identification of an essential virulence gene of cyprinid herpesvirus 3. Antiviral Res. 2017;145:60–9. https://doi.org/10.1016/j.antiviral.2017.07.002.
Shen Z, Jiang Y, Lu J, Sano M, Xu D, Lu L. Application of a monoclonal antibody specific for the ORF92 capsid protein of Cyprinid herpesvirus 2. J Virol Methods. 2018;261:22–7. https://doi.org/10.1016/j.jviromet.2018.07.012.
Jiang B, Zhang W, Liu P, Yang C, Wang M, Jia R, Zhu D, Liu M, Yang Q, Wu Y, et al. The prokaryotic expression, polyclonal antibody preparation, and subcellular localization of the transmembrane protein NS2A of the duck Tembusu virus. Acta Virol. 2020;64:380–5. https://doi.org/10.4149/av_2020_308.
Haq WY, Kang SK, Lee SB, Kang HC, Choi YJ, Lee CN, Kim YS. High-level soluble expression of bioactive porcine myostatin propeptide in E coli. Appl Microbiol Biotechnol. 2013;97:8517–27. https://doi.org/10.1007/s00253-013-5134-0.
Waugh DS. The remarkable solubility-enhancing power of Escherichia coli maltose-binding protein. Postepy Biochem. 2016;62:377–82.
Lee SB, Park SK, Kim YS. Maltose binding protein-fusion enhances the bioactivity of truncated forms of pig myostatin propeptide produced in E coli. PLoS ONE. 2017;12:e0174956. https://doi.org/10.1371/journal.pone.0174956.
Laity JH, Lee BM, Wright PE. Zinc finger proteins: new insights into structural and functional diversity. Curr Opin Struct Biol. 2001;11:39–46. https://doi.org/10.1016/s0959-440x(00)00167-6.
Li L, He S, Sun JM, Davie JR. Gene regulation by Sp1 and Sp3. Biochem Cell Biol. 2004;82:460–71. https://doi.org/10.1139/o04-045.
Liu H, He J, Koh SP, Zhong Y, Liu Z, Wang Z, Zhang Y, Li Z, Tam BT, Lin P, et al. Reprogrammed marrow adipocytes contribute to myeloma-induced bone disease. Sci Transl Med. 2019. https://doi.org/10.1126/scitranslmed.aau9087.
Tang A, Tang J, Miao Q, Zhu J, Guo H, Liu C, Meng C, Li C, Chen Z, Liu G. Zinc finger antiviral protein (ZAP) inhibits small ruminant morbillivirus replication in vitro. Vet Microbiol. 2021;260:109163. https://doi.org/10.1016/j.vetmic.2021.109163.
Yu W, Ji H, Long F, Chen S, He Q, Xia Y, Cong C, Yang C, Wei D, Huang F. Inhibition of hepatitis E virus replication by zinc-finger antiviral Protein synergizes with IFN-β. J Viral Hepat. 2021;28:1219–29. https://doi.org/10.1111/jvh.13522.
Ficarelli M, Antzin-Anduetza I, Hugh-White R, Firth AE, Sertkaya H, Wilson H, Neil SJD, Schulz R, Swanson CM. CpG dinucleotides inhibit HIV-1 replication through zinc finger antiviral protein (ZAP)-dependent and -independent mechanisms. J Virol. 2020. https://doi.org/10.1128/jvi.01337-19.
Arai A, Spencer JA, Olson EN. STARS, a striated muscle activator of Rho signaling and serum response factor-dependent transcription. J Biol Chem. 2002;277:24453–9. https://doi.org/10.1074/jbc.M202216200.
Ohta K, Matsumoto Y, Nishio M. Profilin2 is required for filamentous actin formation induced by human parainfluenza virus type 2. Virology. 2019;533:108–14. https://doi.org/10.1016/j.virol.2019.05.013.
Lv X, Li Z, Guan J, Hu S, Zhang J, Lan Y, Zhao K, Lu H, Song D, He H, et al. Porcine hemagglutinating encephalomyelitis virus activation of the integrin α5β1-FAK-cofilin pathway causes cytoskeletal rearrangement to promote its invasion of N2a cells. J Virol. 2019. https://doi.org/10.1128/jvi.01736-18.
Zheng K, Xiang Y, Wang X, Wang Q, Zhong M, Wang S, Wang X, Fan J, Kitazato K, Wang Y Epidermal growth factor receptor-PI3K signaling controls cofilin activity to facilitate herpes simplex virus 1 entry into neuronal cells. mBio 2014, 5, e00958–00913, .
Li B, Ding S, Feng N, Mooney N, Ooi YS, Ren L, Diep J, Kelly MR, Yasukawa LL, Patton JT, et al. Drebrin restricts rotavirus entry by inhibiting dynamin-mediated endocytosis. Proc Natl Acad Sci U S A. 2017;114:E3642-e3651. https://doi.org/10.1073/pnas.1619266114.
Borovac J, Bosch M, Okamoto K. Regulation of actin dynamics during structural plasticity of dendritic spines: Signaling messengers and actin-binding proteins. Mol Cell Neurosci. 2018;91:122–30. https://doi.org/10.1016/j.mcn.2018.07.001.
Malliri A, Collard JG. Role of Rho-family proteins in cell adhesion and cancer. Curr Opin Cell Biol. 2003;15:583–9. https://doi.org/10.1016/s0955-0674(03)00098-x.
Yeh CH, Bellon M, Nicot C. FBXW7: a critical tumor suppressor of human cancers. Mol Cancer. 2018;17:115. https://doi.org/10.1186/s12943-018-0857-2.
Acknowledgements
This work was supported by grants from the Key Research and Development Plan of Jiangsu Province (No. BE2021369) and the earmarked fund for CARS (CARS-45-16).
Author information
Authors and Affiliations
Contributions
Conceptualization, H.W. and L. L.; Methodology, M.Q.; Software, M.Q. and S. X.; Validation, M.Q, S. X. and J.W.; Formal Analysis, M.Q.; Investigation, M.Q.; Data Curation, M.Q. and Y. Y; Writing – Original Draft Preparation, M.Q. and H.W.; Writing – Review & Editing, H.W.; Supervision, H.W.; Project Administration, H.W.; Funding Acquisition, H.W. and F. Y. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent to publication
Not applicable.
Competing interests
The authors declare that they have no competing of interest.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1.
Analysis of physical and chemical parameters and conserved domain of interactive protein.
Additional file 2.
Bioinformatic analysis and melting curves.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Qian, M., Xiao, S., Yang, Y. et al. Screening and identification of cyprinid herpesvirus 2 (CyHV-2) ORF55-interacting proteins by phage display. Virol J 20, 66 (2023). https://doi.org/10.1186/s12985-023-02026-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12985-023-02026-x