Skip to main content

Identification of a neutralizing linear epitope within the VP1 protein of coxsackievirus A10



Coxsackievirus A10 (CV-A10) is a leading cause of hand, foot, and mouth disease (HFMD). It is necessary to identify neutralizing epitopes to investigate and develop an epitope-based vaccine against CV-A10. The viral protein VP1 is the immunodominant capsid protein and contains the critical neutralizing epitope. However, neutralizing epitopes within VP1 protein of CV-A10 have not been well characterized.


Bioinformatics techniques were applied to predict linear epitopes on the CV-A10 VP1 protein. The advanced structural features of epitopes were analyzed by three-dimensional (3D) modeling. The anticipated epitope peptides were synthesized and used to immunize mice as antigens. ELISA and micro-neutralization assay were used to determine the specific IgG antibody and neutralizing antibody titers. The protective efficacy of the epitope peptides in vivo was evaluated using a passive immunization/challenge assay.


Three linear epitopes (EP3, EP4, and EP5) were predicted on CV-A10 VP1, all spatially exposed on the capsid surface, and exhibited adequate immunogenicity. However, only EP4, corresponding to residues 162–176 of VP1, demonstrated potent neutralization against CV-A10. To determine the neutralizing capacity of EP4 further, EP4 double-peptide was synthesized and injected into mice. The mean neutralizing antibody titer of the anti-EP4 double-peptide sera was 1:50.79, which provided 40% protection against lethal infection with CV-A10 in neonatal mice. In addition, sequence and advanced structural analysis revealed that EP4 was highly conserved among representative strains of CV-A10 and localized in the EF loop region of VP1, like EV-A71 SP55 or CV-A16 PEP55.


These data demonstrate that EP4 is a specific linear neutralizing epitope on CV-A10 VP1. Its protective efficacy can be enhanced by increasing its copy number, which will be the foundation for developing a CV-A10 epitope-based vaccine.


Coxsackievirus A10 (CV-A10) is a human enterovirus belonging to the A species of the Enterovirus genus within Picornaviridae family [1]. CV-A10-associated hand, foot, and mouth diseases (HFMD) typically manifest mild and self-limiting symptoms. Rarely, severe clinical manifestations, such as herpangina, acute hemorrhagic conjunctivitis, acute respiratory tract infections, viral myocarditis, and even death, are observed in a few cases [2, 3]. Recently, the increasing incidence of CV-A10 infection has been an alarming trend [4]. In addition, the co-circulation of CV-A10 with other human enteroviruses may increase the likelihood of co-infection and the risk of EV genetic recombination. However, there is no specific medicine or vaccine registered against CV-A10. Several CV-A10 vaccines, including inactivated vaccines, virus-like particles (VLPs), and VLPs derived from non-EV-core viruses, have been developing for years [5].

CV-A10 is a non-enveloped icosahedral particle with a diameter of approximately 24–30 nm. The genome is approximately 7500 nt in length and consists of a single stranded positive sense RNA genome with two open reading frames (ORF1 and ORF2) [6]. ORF1 encodes a polyprotein. Further cleavage of polyprotein generates three precursor proteins. Four structural proteins (VP1, VP2, VP3, and VP4) derived from the precursor protein P1 are assembled to form the virus capsid [7]. P2 and P3 are subdivided into the nonstructural viral proteins, including 2A, 2B, 2C, 3A, 3B, 3C, and 3D [8].

Studies have shown that multiple protruding loops, such as B–C loop (residues 97–105), E–F loop (residues 163–177), G-H loop (residues 208–225), and proximity of C-terminus (residues 253–267) in VP1 of EV-A71 have been identified as the significant antigenic proteins exposed on the viral surface [9,10,11,12]. The neutralizing epitopes in these loop structures have been extensively characterized by neutralizing antisera and monoclonal antibodies. The results indicate that VP1 protein of enterovirus contains potential neutralizing epitopes and is the primary target for diagnostic reagents and vaccine research [13,14,15]. However, neutralizing epitopes in CV-A10 VP1 have not been experimentally confirmed.

The conventional approach for identifying immunodominant linear neutralizing epitopes on the viral capsid is made by screening a peptide library covering the entire protein sequences. These peptides are analyzed with the polyclonal or monoclonal antibody, followed by the determination of their neutralization-inhibitory effects [6, 16], which are generally considered expensive and time consuming. Bioinformatics can be utilized to rapidly predict linear epitopes and present their advanced structural characteristics [17, 18]. When combined with immunological techniques, the efficiency and precision of identifying epitopes can be improved [19]. Current linear epitope prediction includes the single parametric method and the multi-parameter synthetic prediction method utilizing online servers, such as ABCpred, BepiPred, and SVMTrip [20]. The physicochemical properties of epitopes, including secondary structure, hydrophilicity, flexibility, antigenicity, and surface accessibility, can be analyzed by PSIPRED and DNA STAR software. Furthermore, the location of epitopes can be represented visually by three-dimensional (3D) models. These techniques have been extensively applied to investigate the epitopes of influenza virus, hepatitis A virus, SARS-CoV-2, and other pathogens [21,22,23].

In the current study, bioinformatics techniques were applied to predict three candidate linear neutralizing epitopes within the VP1 protein of CV-A10. The candidate epitope peptides were subsequently synthesized and conjugated with keyhole limpet hemocyanin (KLH) for the immunization of mice. The levels of specific IgG and neutralizing antibodies were evaluated in the sera of mice. EP4, a highly neutralizing epitope peptide, was identified as a potential neutralizing epitope. Notably, EP4 double-peptide demonstrated enhanced neutralizing protection, which could passively protect neonatal mice from lethal CV-A10 challenges, making it a promising candidate for CV-A10 vaccine.


Cells and virus

Human Rhabdomyoma (RD) cells were cultured in Dulbecco modified eagle medium (DMEM) containing 10% fetal bovine serum (FBS). CV-A10 strain used in this study was the isolate P148/ZS/CHN/2012 (GenBank: MK645898.1).

Linear epitope prediction and secondary structure analysis of CV-A10 VP1 protein

To predict the linear neutralization epitopes of CA-V10 VP1, the amino acid sequence of CV-A10-P148 VP1 protein was obtained from the GenBank database (GenBank ID: QEO24698.1). ABCPreds (, BCPreds ( and SVMTriP ( were used to predict the possible linear neutralizing epitopes. PSIPRED 4.0 and DNA Star Protean modules [24, 25] were utilized to analyze the secondary structure, flexibility, hydrophilicity, surface accessibility, and antigenicity of VP1 protein. According to the epitope prediction and the analysis results of physicochemical properties, a peptide meeting the following conditions was selected as a candidate epitope: (1) it should be located in the high score region predicted by the server; (2) it should avoid α-helical or β-strand structure; (3) it should show good hydrophilicity, flexibility, and antigenicity; (4) loop structures that are exposed to the outside surface are preferentially selected.

Homology modeling and epitope spatial displaying

Homologous sequences were retrieved using the automatic mode of SWISS-MODEL, then the Cryo-EM structure of CV-A10 mature virion (PDB:6IIJ) was selected as a modeling template. The VP1, VP2, VP3, and VP4 amino acid sequences of CV-A10-P148 were submitted to a new modeling project. The PyMOL program was used to visualize the 3D model structure and display the position of the epitope.

Synthetic peptides

All peptides were synthesized by Sener Biotechnology Co., LTD (Hefei, China). An irrelevant peptide (QLINTNGSWHINSTA) from the hepatitis C virus E2 protein was synthesized as a negative control [4]. The purity of all synthetic peptides was > 95%. For immunization, all synthetic peptides were conjugated to KLH.

Mice immunization and serum sample collection

All animal experiments were conducted strictly with experimental protocols approved by the Animal Ethics Committee at Guilin Medical College. BALB/c mice produced in specific pathogen-free (SPF) conditions for this study were provided by SJA Laboratory Animal Co., LTD (Hunan, China). Six-week-old female BALB/C mice were randomly divided into groups with six mice per group. The experimental group, negative control group, and blank group mice were multi-point subcutaneously injected with synthetic peptides (100 µg/dose), irrelevant peptides (100 µg/dose), or PBS, respectively. The first injection of each sample was emulsified with an equal volume of complete Freund's adjuvant (Sigma), and the booster was emulsified with an equal volume of incomplete Freund's adjuvant. Serum samples were harvested two weeks after the last immunization [16].


The peptide-specific IgG antibody in immunized mouse sera was determined by indirect ELISA assay. Briefly, each well in a 96-well plate was coated with 10 ng of a synthetic epitope at 4 °C overnight and then was blocked with 1% bovine serum albumin (BSA) for 1 h at 37 °C. The plates were incubated at 37 °C for 1 h with epitope peptide antisera that diluted by two-fold starting at 1:500 in PBS-Tween-20 (PBST) containing 0.1% BSA and then was incubated with HRP-conjugated goat anti-mouse IgG antibody that diluted at 1:5000 in PBST at 37 °C for 1 h. The plates were washed three times with PBST between each step. After color development, the optical density (OD) was measured on a microplate reader (Bio-Rad) at 450 nm.

Micro-neutralization assay

Referring to previous research [26], the titer of CV-A10 neutralizing antibody in the sera was determined using the cytopathic effect (CPE)-based micro-neutralization assay on RD cells. The neutralization antibody titer was defined as the reciprocal of the highest serum dilutions that could protect 50% of cells from CPE.

In vivo protection analysis

Enterovirus-infected neonatal mouse model is a classic model commonly used in the evaluation of vaccines or antiviral drugs [27, 28]. A standardized CV-A10- infected neonatal mouse model [29] was used to evaluate the protective effects of EP4 double-peptide against a challenge with the mouse adapted CV-A10 strain. Since the immune system of newborn mice is not yet mature enough to elicit immune responses against foreign antigens, protective efficacy was determined by passive immunization/challenge assay [30]. Adult ICR mice were grouped and injected using the same method and dose identical to that used in the immunogenicity experiment. Mice were mated one week after the first booster, and each mouse gave birth to about 10 suckling mice about three weeks after mating. Neonatal mice were challenged intracranially with 50 × LD50 of CV-A10 within 24 h after delivery. They were injected intraperitoneally with pooled anti-EP4 double-peptide sera, anti-irrelevant peptide sera, or anti-PBS sera within 2 h after the challenge. The clinical symptoms and death of mice in each group were monitored for 15 consecutive days.

Sequence alignment and advanced structure analysis

All VP1 protein sequences of different enteroviruses for sequence alignment were retrieved from GenBank. The alignment of VP1 protein sequences was conducted using BioEdit software. The PyMOL program was used for the superposition of VP1 3D models of CV-A10 (PDB: 6IIJ), EV-A71 (PDB: 4aed), and CV-A16 (PDB: 5c4w). The linear neutralizing epitopes were tagged on 3D models of each VP1, respectively.

Statistical analysis

GraphPad Prism 8 software was used for graph generation and statistical analyses. Log-rank test was employed to compare Kaplan Meier survival curves, and the student's two-tailed t-test was utilized to analyze serum-specific IgG antibody titer and neutralizing antibody titer. A p-value of less than 0.05 was considered statistically significant.


Prediction results of the B cell linear epitopes of CV-A10 VP1

B cell linear epitopes are crucial in eliciting humoral immune protection [15]. The virus cannot invade cells if neutralizing antibodies induced by linear neutralizing epitopes are present. Several servers based on neural algorithms and polypeptide feature analysis were implemented to improve the accuracy of epitope prediction. Five candidate linear epitopes on CV-A10 VP1 (Table 1) were selected for further investigation based on the results of ABCpred, BCPred, and SVMTrip servers (Additional file 1: Table S1) and analysis of their secondary structures, flexibility, antigenicity, and hydrophilicity (Additional file 1: Fig. S1).

Table 1 Five possible linear epitopes on CV-A10 VP1 protein

The localization of candidate linear epitopes on viral capsid

The pentamer spatial model of CV-A10 capsid was developed using SWISS-MODEL (Fig. 1). VP4 and five candidate linear epitopes of VP1 were labeled on the model. VP4 protein is located on the inner surface of the capsid. EP1 and EP2 are also located on the inner surface, whereas EP3, EP4, and EP5 are on the outer surface. Neutralizing epitopes, which mediate virus-receptor interactions, were typically located on the external surface of the viral capsid near the vertices; therefore, EP3, EP4, and EP5 were favored as candidate linear epitopes of CV-A10 VP1.

Fig. 1
figure 1

Location of predicted linear epitopes on CV-A10 P1 protein. EP1-EP5 peptide and VP4 protein are colored in red, yellow, blue, green, purple, and orange, respectively. A The tertiary structure is constructed with CV-A10 P1 amino acid sequence. B The top (left) and the bottom (right) of CV-A10 P1 pentamer are shown as a space fill model, where EP3, EP4, and EP5 are exposed on the outer capsid surface, and EP1, EP2, and VP4 proteins are in the inner surface capsid

Analysis of mice IgG responses elicited by synthetic epitopes

For the detection of immunogenicity of the linear epitopes, the three peptides (EP3, EP4, and EP5) were chemically synthesized and conjugated to KLH and administered to mice for immunization. IgG antibody levels in the sera obtained from immunized mice were determined using ELISA with a corresponding synthetic peptide as the capture antigen. The sera of mice immunized with PBS served as blank, whereas the sera of mice immunized with irrelevant peptides served as the negative control. The negative control or blank antisera exhibited no significant reactivity with any peptides. As illustrated in Fig. 2, all antiserum samples against EP3, EP4, or EP5 contained significantly higher IgG levels than the negative control. EP5 demonstrated the highest immunogenicity. The differences in OD values between groups suggest that EP3, EP4, and EP5 may induce varying degrees of immune response in mice.

Fig. 2
figure 2

Peptide-specific IgG antibody levels in the sera obtained from immunized mice. Five adult female BALB/c mice received priming injection at week 0 and boosting injection at weeks 2 and 4 with PBS, irrelevant peptide, EP3, EP4, or EP5, respectively. The end-point dilution titers were tested by peptide ELISA using EP3 (A), EP4 (B), or EP5 (C) peptide as the corresponding antigen. Anti-PBS sera served as blank control, and anti-irrelevant peptide sera served as the negative control (NC). The data are presented as mean ± standard deviation (SD) of at least triplicate wells, and the means of IgG titers were compared using a student's two-tailed t-test. **p < 0.01, ****p < 0.0001

In vitro neutralizing capacity of the candidate linear epitopes antisera

The micro-neutralization assay was utilized to assess the ability of the candidate linear epitope antisera to neutralize CV-A10 infection. As demonstrated in Fig. 3, EP3, EP4, and EP5 antisera neutralized CV-A10 P148 strain with geometric mean titers of 1:4, 1:17.96, and 1:3.56, respectively. In contrast, the control antisera lacked CV-A10 neutralizing activity. There were statistically significant differences between the neutralizing capacity of groups immunized with EP4 and negative control. Since the neutralizing antibody titer elicited by the single epitope peptide was low, mice were immunized with a double-peptide EP4 to determine if EP4 is a dominant neutralizing epitope. Two weeks after three booster injections, sera for the micro-neutralization experiment were collected. As depicted in Fig. 4, sera from EP4 double-peptide group neutralized CV-A10 with titers ranging from 8 to 128, as measured by dilution.

Fig. 3
figure 3

CV-A10-neutralizing antibody titration in anti-peptide sera from immunized mice. Two weeks after the final immunization, anti-peptide sera from individual mice were collected and tested for CV-A10 virus neutralization. Each symbol represents a serum sample, and the horizontal line indicates geometric mean values. According to the statistical analysis, there was a significant difference in neutralizing antibody titers between anti-EP4 sera and negative control sera. The data were analyzed using a student's two-tailed t-test. ***p < 0.001

Fig. 4
figure 4

The neutralization titer of anti-EP4 double-peptide sera. Mice were immunized with EP4 double-peptide and given three booster injections two weeks apart. Sera were obtained two weeks after the final immunization for micro-neutralization experiment. The geometric mean neutralizing antibody titer of anti-EP4 double-peptide sera against CV-A10 reached 1:50.79, and a significant difference was shown between anti-EP4 double-peptide sera and negative control sera. The data were analyzed using a student's two-tailed t-test. **p < 0.01

The in vivo protective efficacy of anti-EP4 double-peptide sera against lethal CV-A10

The protective efficacy of EP4 double-peptide in vivo was evaluated by a lethal dose of a mouse-adapted CV-A10 strain intracranially to ICR neonatal mice [31]. The results demonstrated that the mice treated with anti-PBS sera or anti-irrelevant peptide sera exhibited multiple symptoms, including limb-shake weakness and hind-limb paralysis, and began to die on day 3, with all mice dying by days 7 or 8, respectively (Fig. 5A). The mice treated with anti-EP4 double-peptide sera developed mild or severe clinical symptoms on day 4, and some mice died between day 7 and 10, while the rest gradually regained their health, with a 40% final survival rate on day 15 (Fig. 5A, B).

Fig. 5
figure 5

In vivo protective efficacy of anti-EP4 double-peptide sera against lethal CV-A10 challenge. The neonatal mice were challenged with 50 × LD50 of CV-A10 virus within 24 h after birth (each group = 10). They were then injected with anti-EP4 double-peptide sera, anti-irrelevant peptide sera, or anti-PBS sera. A Survival curve of the neonatal mice. B Clinical scores of the neonatal mice. The severity of clinical symptoms, from mild to severe, was scored in five grades: 0, healthy; 1, reduced mobility; 2, limb weakness; 3, paralysis; 4, death. The Mantel-Cox log-rank test evaluated the survival rates. The clinical scores were compared using Dunn's multiple-comparison test. ***p < 0.001

The sequence and advanced structure analysis of EP4

The CV-A10 VP1 sequences were analyzed by SDT software, and the phylogenetic tree was constructed using MEGA 7.0 software (Additional file 1: Fig. S2). Clustal X2 software was utilized to align the VP1 amino acid sequences of 33 representative strains of CV-A10 (A-F group) with the prototype strains of EV-A71 and CV-A16. The results indicated that EP4 region (residues 162–176 of VP1) was highly conserved among selected CV-A10 isolates. Nevertheless, EP4 was found to be 60% consistent with SP55 and 80% consistent with PEP55 (Fig. 6). EP4, SP55, and PEP55 were all located in EF loop region of VP1 and were highly exposed at the virion surface (Fig. 7A) as demonstrated by the advanced structural analysis by SWISS-MODEL. These three epitopes shared a short peptide segment (WQTATNP) that was identical. EP4 differs from SP55 and PEP55 because it has one additional random curl (Fig. 7B–D). EP4 appears to be a surface-exposed, highly conserved, neutralizing epitope on CV-A10 VP1 protein.

Fig. 6
figure 6

Alignment of CV-A10 VP1 amino acid sequences and the location of CV-A10 EP4. According to the evolutionary analysis of CV-A10 (Additional file 1: Fig. S2), 33 amino acid sequences of CV-A10 VP1 in group A-G were selected and compared with VP1 sequences of EV-A71 and CV-A16 prototype strains. Dots represent residues identical to CV-A10/P148/CHN/2012. EP4 region is boxed with a black line

Fig. 7
figure 7

Advanced structural analysis of CV-A10 EP4, EV-A71 SP55, and CV-A16 PEP55. A Superposition of VP1 of CV-A10 (blue), CV-A16 (yellow) and EV-A71 (light gray). The CV-A10 EP4 is magnified and colored in red. BD represent the secondary structure of CV-A10 EP4, EV-A71 SP55 and CV-A16 PEP, respectively


Neutralizing antibodies are elicited by the neutralizing epitopes. The level of neutralization titer is related to the amount of effective viral antigen, which is suggested as a criterion for evaluating the efficacy of vaccine candidates. Therefore, neutralizing epitopes play a key role in vaccine effectiveness and can serve as biomarkers to monitor vaccine efficacy [32]. The discovery and identification of neutralizing epitopes is critical to vaccine development.

Recently, CV-A10 has become one of the most prevalent pathogens of HFMD [33]. The linear neutralizing epitopes are essential for vaccine development of CV-A10 and its diagnosis [34]. Historically, epitopes were identified and screened by chemo-synthesizing short overlapping peptides based on their sequences. Immunological testing was then used to confirm the antigenic and immunogenic properties of each peptide. This method requires a great deal of time, human resources, and materials, and it may miss some potential neutralizing epitopes despite effectively covering the entire polypeptide chain.

In this study, online servers predicted thirteen epitopes with high scores (Additional file 1: Table S1). Based on the secondary structure analysis of VP1 protein, peptides 15, 11 and 3 with the α-helical or β-strand structure were eliminated (Additional file 1: Fig. S1A). Then, based on hydrophilicity, elasticity, antigenicity, and advanced structure of VP1, five linear neutralizing epitope candidates for CV-A10 VP1 were selected. A 3D model displayed the spatial locations of these five epitopes and demonstrated that EP1 and EP2 were embedded within the capsid, whereas EP3–EP5 were located on the exterior of the capsid. Since neutralizing epitopes are immunoreactive regions on antigens and are responsible for binding to immune cell receptors or free antibodies, they are typically found on the surface of the viral capsid. EP1 and EP2, therefore, do not meet the requirement for neutralization.

Furthermore, Dai's study demonstrated that peptides in the EP1 region are antigenic but lack neutralizing properties [6]. In addition, it is demonstrated that bioinformatics is superior to overlapping peptide mapping for locating the region that corresponds to neutralizing epitope characteristics. EP3–EP5 was synthesized and then administered to mice. ELISA demonstrated that all three candidate epitopes could induce elevated levels of self-specific IgG. EP5 induced the highest IgG titer, possibly because its spatial configuration is more conducive to the recognition by immune cells. In a micro-neutralization assay, EP4 was found to have a higher neutralizing antibody titer than EP3 and EP5. It suggests that EP4 may be a potential neutralizing epitope of CV-A10 VP1. To obtain a high-titer neutralizing antibody against EP4 and demonstrate its protective efficacy in vivo, tandem arrays of EP4 double-peptides conjugated with KLH were synthesized. The geometric mean neutralizing antibody titer against CV-A10 in antiserum from mice immunized with double-peptide EP4 was 1:50.79. To verify the passive protection efficiency of double-peptide EP4 in vivo, an experimental animal model of CV-A10 infection was used [29]. The double-peptide EP4 provided 40% protection against the 50 × LD50 CV-A10 challenge in neonatal mice, whereas the PBS or the irrelevant peptide conferred no protection. After further refinement, EP4 double-peptide could be used to develop a vaccine against CV-A10 epitope based on the above results.

Studies demonstrated that EV-A71 and CV-A16 utilize human scavenger receptor class B, member 2 (hSCARB2), or human P selectin glycoprotein ligand 1 (PSGL-1) as cellular receptors [35]. The locations of epitopes in CV-A6 VP1 are comparable to those in EV-A71, including the BC loop, EF loop, GH loop, and C-terminus [9,10,11,12]. PSGL-1 was considered the cell receptor of CV-A10 [11], so the neutralizing epitopes of CV-A10 VP1 may be located in the same region as EV-A71. Sequence analysis and 3D modeling supported the hypothesis, as mentioned earlier, EP4 was located in EF loop of CV-A10 VP1 [36], the exact location as SP55 of EV-A71 VP1 and PEP55 of CV-A16 VP1. EP4 amino acids were highly conserved in CV-A10 genotype, sharing approximately 60% and 80% sequence identity with EV-A71 SP55 and CV-A16 PEP55, respectively, indicating that EP4 can also be used as a specific target for detecting CV-A10. EP4 contains three amino acids (G164S, A167S, and Q169A) different from SP55 and PEP55 that distinguish it from EV-A71 and CV-A16, suggesting that these three amino acids determine the specificity of EP4. Using a mutagenesis strategy in EP4, the relationship between EP4 and virulence, transmissibility, or pathogenicity of CV-A10 can be further investigated.

In this study, a neutralizing epitope of CV-A10 VP1 was identified. Additional bioinformatic and immunology research is required to determine if neutralizing epitopes exist in BC loop of VP1, GH loop of VP2, and GH loop of VP3. Interestingly, the immune protective range can be expanded by replacing the original epitope with a heterologous neutralizing epitope on the capsid of an EV-A71 VLP vaccine or inactivated vaccine [36, 37], suggesting that EP4 could be used as a candidate epitope for developing the novel polyvalent chimeric vaccine.


On CV-A10 VP1, a highly conserved linear neutralizing epitope (EP4) has been identified. Neutralizing antibodies elicited by EP4 could protect mice from a lethal CV-A10 challenge. It can be used as a guide for designing CV-A10 epitope peptides or enterovirus polyvalent vaccines.

Availability of data and materials

Not applicable.



Coxsackievirus A10


Coxsackievirus A16

EV-A 71:

Enterovirus A71


Hand, foot and mouth disease




Virus-like particles


Open reading frame


Keyhole limpet hemocyanin




Dulbecco’s modified eagle medium


Fetal bovine serum


Specific pathogen-free


Bovine serum albumin




Optical density


Cytopathic effect


Human scavenger receptor class B, member 2


P selectin glycoprotein ligand 1


Institute of Cancer Research


50% Tissue culture infectious dose


  1. Li K, Dong F, Cui B, Cui L, Liu P, Ma C, Zheng H, Wu X, Liang Z. Development of a pseudovirus-based assay for measuring neutralizing antibodies against Coxsackievirus A10. Hum Vaccin Immunother. 2020;16(6):1434–40.

    Article  CAS  PubMed  Google Scholar 

  2. Zhu L, Sun Y, Fan J, Zhu B, Cao L, Gao Q, Zhang Y, Liu H, Rao Z, Wang X. Structures of Coxsackievirus A10 unveil the molecular mechanisms of receptor binding and viral uncoating. Nat Commun. 2018;9(1):4985.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Bian L, Gao F, Mao Q, Sun S, Wu X, Liu S, Yang X, Liang Z. Hand, foot, and mouth disease associated with coxsackievirus A10: more serious than it seems. Expert Rev Anti Infect Ther. 2019;17(4):233–42.

    Article  CAS  PubMed  Google Scholar 

  4. Zhang Z, Dong Z, Li J, Carr MJ, Zhuang D, Wang J, Zhang Y, Ding S, Tong Y, Li D, Shi W. Protective efficacies of formaldehyde-inactivated whole-virus vaccine and antivirals in a murine model of Coxsackievirus A10 infection. J Virol. 2017;91(13):e00333-17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Dai W, Xiong P, Zhang X, Liu Z, Chen J, Zhou Y, Ye X, Zhang C. Recombinant virus-like particle presenting a newly identified coxsackievirus A10 neutralization epitope induces protective immunity in mice. Antiviral Res. 2019;164:139–46.

    Article  CAS  PubMed  Google Scholar 

  6. Guo H, Li Y, Liu G, Jiang Y, Shen S, Bi R, Huang H, Cheng T, Wang C, Wei W. A second open reading frame in human enterovirus determines viral replication in intestinal epithelial cells. Nat Commun. 2019;10(1):4066.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Zhao Y, Zhou D, Ni T, Karia D, Kotecha A, Wang X, Rao Z, Jones EY, Fry EE, Ren J, Stuart DI. Hand-foot-and-mouth disease virus receptor KREMEN1 binds the canyon of Coxsackie Virus A10. Nat Commun. 2020;11(1):38.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Fang CY, Liu CC. Recent development of enterovirus A vaccine candidates for the prevention of hand, foot, and mouth disease. Expert Rev Vaccin. 2018;17(9):819–31.

    Article  CAS  Google Scholar 

  9. Foo DG, Alonso S, Phoon MC, Ramachandran NP, Chow VT, Poh CL. Identification of neutralizing linear epitopes from the VP1 capsid protein of Enterovirus 71 using synthetic peptides. Virus Res. 2007;125(1):61–8.

    Article  CAS  PubMed  Google Scholar 

  10. Foo DG, Ang RX, Alonso S, Chow VT, Quak SH, Poh CL. Identification of immunodominant VP1 linear epitope of enterovirus 71 (EV71) using synthetic peptides for detecting human anti-EV71 IgG antibodies in Western blots. Clin Microbiol Infect. 2008;14(3):286–8.

    Article  CAS  PubMed  Google Scholar 

  11. Liu CC, Chou AH, Lien SP, Lin HY, Liu SJ, Chang JY, Guo MS, Chow YH, Yang WS, Chang KH, Sia C, Chong P. Identification and characterization of a cross-neutralization epitope of Enterovirus 71. Vaccine. 2011;29(26):4362–72.

    Article  CAS  PubMed  Google Scholar 

  12. Kirk K, Poh CL, Fecondo J, Pourianfar H, Shaw J, Grollo L. Cross-reactive neutralizing antibody epitopes against Enterovirus 71 identified by an in silico approach. Vaccine. 2012;30(49):7105–10.

    Article  CAS  PubMed  Google Scholar 

  13. Zhao Q, Syed SF, Zhou EM. Antigenic properties of avian hepatitis E virus capsid protein. Vet Microbiol. 2015;180(1–2):10–4.

    Article  CAS  PubMed  Google Scholar 

  14. Kiener TK, Jia Q, Lim XF, He F, Meng T, Chow VT, Kwang J. Characterization and specificity of the linear epitope of the enterovirus 71 VP2 protein. Virol J. 2012;9:55.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Song Y, Zhang Y, Ji T, Gu X, Yang Q, Zhu S, Xu W, Xu Y, Shi Y, Huang X, Li Q, Deng H, Wang X, Yan D, Yu W, Wang S, Yu D, Xu W. Persistent circulation of Coxsackievirus A6 of genotype D3 in mainland of China between 2008 and 2015. Sci Rep. 2017;7(1):5491.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Banerjee A, Kanwar M, Santra D, Maiti S. Global conserved RBD fraction of SARS-CoV-2 S-protein with T500S mutation in silico significantly blocks ACE2 and rejects viral spike. Transl Med Commun. 2022;7(1):2.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Zhang J, Huang H, Xu L, Lou C, Pan M. Screening and identification of linear B cell epitopes within the nonstructural proteins of enterovirus 71. Viral Immunol. 2019;32(2):84–8.

    Article  CAS  PubMed  Google Scholar 

  18. Liu HB, Yang GF, Liang SJ, Lin J. Bioinformatic analysis of non-VP1 capsid protein of coxsackievirus A6. J Huazhong Univ Sci Technol Med Sci. 2016;36(4):607–13.

    Article  CAS  Google Scholar 

  19. Ahmad I, Ali SS, Zafar B, Hashmi HF, Shah I, Khan S, Suleman M, Khan M, Ullah S, Ali S, Khan J, Ali M, Khan A, Wei DQ. Development of multi-epitope subunit vaccine for protection against the norovirus’ infections based on computational vaccinology. J Biomol Struct Dyn. 2022;40(7):3098–109.

    Article  CAS  PubMed  Google Scholar 

  20. Chen HZ, Tang LL, Yu XL, Zhou J, Chang YF, Wu X. Bioinformatics analysis of epitope-based vaccine design against the novel SARS-CoV-2. Infect Dis Poverty. 2020;9(1):88.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Assmus LM, Guan J, Wu T, Farenc C, Sng XYX, Zareie P, Nguyen A, Nguyen AT, Tscharke DC, Thomas PG, Rossjohn J, Gras S, Croft NP, Purcell AW, La Gruta NL. Overlapping peptides elicit distinct CD8(+) T cell responses following influenza A virus infection. J Immunol. 2020;205(7):1731–42.

    Article  CAS  PubMed  Google Scholar 

  22. Khudyakov YE, Lopareva EN, Jue DL, Fang S, Spelbring J, Krawczynski K, Margolis HS, Fields HA. Antigenic epitopes of the hepatitis A virus polyprotein. Virology. 1999;260(2):260–72.

    Article  CAS  PubMed  Google Scholar 

  23. Olvera A, Noguera-Julian M, Kilpelainen A, Romero-Martin L, Prado JG, Brander C. SARS-CoV-2 consensus-sequence and matching overlapping peptides design for COVID19 immune studies and vaccine development. Vaccines (Basel). 2020;8(3):444.

    Article  CAS  PubMed  Google Scholar 

  24. Buchan DW, Minneci F, Nugent TC, Bryson K, Jones DT. Scalable web services for the PSIPRED Protein Analysis Workbench. Nucl Acids Res. 2013;41(Web Server issue):W349-57.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Plasterer TN. PROTEAN. Protein sequence analysis and prediction. Mol Biotechnol. 2000;16(2):117–25.

    Article  CAS  PubMed  Google Scholar 

  26. Ku Z, Ye X, Huang X, Cai Y, Liu Q, Li Y, Su Z, Huang Z. Neutralizing antibodies induced by recombinant virus-like particles of enterovirus 71 genotype C4 inhibit infection at pre- and post-attachment steps. PLoS One. 2013;8(2):e57601.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Li S, Zhao H, Yang L, Hou W, Xu L, Wu Y, Wang W, Chen C, Wan J, Ye X, Liang Z, Mao Q, Cheng T, Xia N. A neonatal mouse model of coxsackievirus A10 infection for anti-viral evaluation. Antiviral Res. 2017;144:247–55.

    Article  CAS  PubMed  Google Scholar 

  28. Mao Q, Wang Y, Gao R, Shao J, Yao X, Lang S, Wang C, Mao P, Liang Z, Wang J. A neonatal mouse model of coxsackievirus A16 for vaccine evaluation. J Virol. 2012;86(22):11967–76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Liu Q, Dan H, Zhao X, Chen H, Chen Y, Zhang N, Mo Z, Liu H. Construction and characterization of an infectious cDNA clone of coxsackievirus A 10. Virol J. 2019;16(1):98.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Zhang W, Dai W, Zhang C, Zhou Y, Xiong P, Wang S, Ye X, Liu Q, Zhou D, Huang Z. A virus-like particle-based tetravalent vaccine for hand, foot, and mouth disease elicits broad and balanced protective immunity. Emerg Microbes Infect. 2018;7(1):94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Luo Z, Su R, Wang W, Liang Y, Zeng X, Shereen MA, Bashir N, Zhang Q, Zhao L, Wu K, Liu Y, Wu J. EV71 infection induces neurodegeneration via activating TLR7 signaling and IL-6 production. PLOS Pathog. 2019;15(11):e1008142.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Liu CC, Chang HW, Yang G, Chiang JR, Chow YH, Sai IH, Chang JY, Lin SC, Sia C, Hsiao CH, Chou AH, Chong P. Development of a quantitative enzyme linked immunosorbent assay for monitoring the Enterovirus 71 vaccine manufacturing process. J Virol Methods. 2011;176(1–2):60–8.

    Article  CAS  PubMed  Google Scholar 

  33. He SZ, Chen MY, Xu XR, Yan Q, Niu JJ, Wu WH, Su XS, Ge SX, Zhang SY, Xia NS. Epidemics and aetiology of hand, foot and mouth disease in Xiamen, China, from 2008 to 2015. Epidemiol Infect. 2017;145(9):1865–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Gonzalez G, Carr MJ, Kobayashi M, Hanaoka N, Fujimoto T. Enterovirus-associated hand-foot and mouth disease and neurological complications in Japan and the rest of the world. Int J Mol Sci. 2019;20(20):5201.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Baggen J, Thibaut HJ, Strating JRPM, van Kuppeveld FJM. T The life cycle of non-polio enteroviruses and how to target it. Nat Rev Microbiol. 2018;16(6):368–81. (Erratum. In: Nat Rev Microbiol. 2018 May 3).

    Article  CAS  PubMed  Google Scholar 

  36. Meijer A, Van Der Sanden S, Snijders BE, Jaramillo-Gutierrez G, Bont L, Van Der Ent CK, Overduin P, Jenny SL, Jusic E, van der Avoort HG, Smith GJ, Donker GA, Koopmans MP. Emergence and epidemic occurrence of enterovirus 68 respiratory infections in The Netherlands in 2010. Virology. 2012;423(1):49–57.

    Article  CAS  PubMed  Google Scholar 

  37. Ho BC, Yang PC, Yu SL. MicroRNA and pathogenesis of enterovirus infection. Viruses. 2016;8(1):11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references


Not applicable.


This research was supported by grants from the National Natural Science Foundation of China (Nos. 81660280 and 31760262), and the Appropriate Health Technology Promotion Project of Guangxi Health Commission (No. S2019036).

Author information

Authors and Affiliations



QL and HL designed the study and the experiments and revised the manuscript. HZ1, XL, YW and YH performed the experiments. HZ1 contributed to data analysis and the raw draft preparation. HZ3 helped to revise the manuscript. All authors have read and approved the manuscript.

Corresponding authors

Correspondence to Hongbo Liu or Qiliang Liu.

Ethics declarations

Ethics approval and consent to participate

This study was performed in strict accordance with the recommendations in the Guide for the Institutional Animal Care and Use Commission (IACUC). The protocols were approved by the Committee on the Ethics of Animal Experiments of the Guilin Medical University (Number: GLMC201703016).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

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: Table S1

The possible linear epitopes given high scores predicted by three servers. Figure S1: The secondary structure analysis of CV-A10 VP1 protein. (A) The secondary structure prediction of CV-A10 VP1 protein by PSIPRED. α-Helical residues are in pink, β-strand residues are in yellow, putative domain boundaries are indicated in blue, and the locations of the predicted epitopes listed in Table S1 are underlined and numbered. (B) Graphic analysis of secondary structure, flexibility, hydrophilicity, surface accessibility and antigenicity of the CV-A10 VP1 protein using the DNAstar Protean module. Figure S2: Evolutionary analysis of CV-A10 VP1 protein. All available CV-A10 sequences were downloaded from GenBank and saved in “FASTA” format. Sequences with high similarity and no specific separation time and location were excluded, and representative sequences were selected for further analysis. ClustalX2 software was utilized to calibrate and compare the CV-A10 representative sequences. According to the VP1 sequence of CV-A10 prototype strain (Kowalik), the VP1 sequence of each isolate was truncated. The CV-A10 VP1 sequences were analyzed with muscle method by SDT software. The strains with nucleotide sequence consistency greater than 75% in the heat map were classified as the same genotype (A). The phylogenetic tree was inferred with the neighbor-joining (N-J) method and Kimura 2-parameter model by MEGA 7.0 software (B). The EV-A71 and CV-A16 prototype strains were used as outgroups. The reliability of the phylogenetic tree was tested by Bootstrap replication (1000) method. The tree is drawn to scale, with branch lengths in the same units as those of the evolutionary distances used to infer the phylogenetic tree. The evolutionary distances were computed using the Maximum Composite Likelihood method and are in the units of the number of base substitutions per site.

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 The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhu, H., Liu, X., Wu, Y. et al. Identification of a neutralizing linear epitope within the VP1 protein of coxsackievirus A10. Virol J 19, 203 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: