# Cloning and Molecular Characterization of the Recombinant CVB4E2 Immunogenic Viral Protein (rVP1), as a Potential Subunit Protein for Vaccine and Immunodiagnostic Reagent Candidate

Ikbel Hadj Hassine, Jawhar Gharbi, Imene Amara, Ameera Alyami, Reem Subei, Mohammed Almalki, Didier Hober, Manel Ben M'Hadheb

► To cite this version:

Ikbel Hadj Hassine, Jawhar Gharbi, Imene Amara, Ameera Alyami, Reem Subei, et al.. Cloning and Molecular Characterization of the Recombinant CVB4E2 Immunogenic Viral Protein (rVP1), as a Potential Subunit Protein for Vaccine and Immunodiagnostic Reagent Candidate. Microorganisms, 2023, Microorganisms, 11, pp.1192. 10.3390/microorganisms11051192. [https://doi.org/10.3390/microorganisms11051192](https://doi.org/10.3390/microorganisms11051192)

HAL Id: hal-04453486

https://hal.univ-lille.fr/hal-04453486

Submitted on 12 Feb 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Cloning and Molecular Characterization of the Recombinant CVB4E2 Immunogenic Viral Protein (rVP1), as a Potential Subunit Protein for Vaccine and Immunodiagnostic Reagent Candidate

Ikbel Hadj Hassine, Jawhar Gharbi, Imene Amara, Ameera Alyami, Reem Subei, Mohammed Almalki, Didier Hober, Manel Ben M'Hadheb

► To cite this version:

Ikbel Hadj Hassine, Jawhar Gharbi, Imene Amara, Ameera Alyami, Reem Subei, et al.. Cloning and Molecular Characterization of the Recombinant CVB4E2 Immunogenic Viral Protein (rVP1), as a Potential Subunit Protein for Vaccine and Immunodiagnostic Reagent Candidate. Microorganisms, 2023, Microorganisms, 11, pp.1192. 10.3390/microorganisms11051192. [https://doi.org/10.3390/microorganisms11051192](https://doi.org/10.3390/microorganisms11051192)

HAL Id: hal-04453486

https://lilloa.hal.science/hal-04453486v1

Submitted on 12 Feb 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Cloning and Molecular Characterization of the Recombinant CVB4E2 Immunogenic Viral Protein (rVP1), as a Potential Subunit Protein for Vaccine and Immunodiagnostic Reagent Candidate

Ikbel Hadj Hassine<fnref n="1" />, Jawhar Gharbi <fnref n="2" />,*, Imene Amara<fnref n="1" />, Ameera Alyami <fnref n="2" />, Reem Subei <fnref n="2" />, Mohammed Almalki <fnref n="2" />, Didier Hober <fnref n="3" /> and Manel Ben M'hadheb <fnref n="1" />,*

**Abstract:** The aim of the present study was, first, to clone the VP1 gene of the human coxsackievirus B4 strain E2 (CVB4E2) in the prokaryotic pUC19 plasmid expression vector then to compare it with the structural capsid proteins of the same strain using bioinformatic tools. PCR colony amplification followed through a restriction digestion analysis and sequencing process which affirmed the success of the cloning process. SDS-PAGE and Western Blotting were used to characterize the purified recombinant viral protein expressed in bacteria cells. The BLASTN tool revealed that the nucleotide sequence of the recombinant VP1 (rVP1) expressed by pUC19 highly matched the target nucleotide sequence of the diabetogenic CVB4E2 strain. Secondary structure and three-dimension structure prediction suggested that rVP1, such as wild-type VP1, is chiefly composed of random coils and a high percentage of exposed amino acids. Linear B-cell epitope prediction showed that several antigenic epitopes are likely present in rVP1 and CVB4E2 VP1 capsid protein. Additionally, phosphorylation site prediction revealed that both proteins may affect the signal transduction of host cells and can be involved in virus virulence. The present work highlights the usefulness of cloning and bioinformatics characterizations for gene investigation. Furthermore, the collected data are helpful for future experimental research related to the development of immunodiagnostic reagents and subunit vaccines based on the expression of immunogenic viral capsid proteins.

**Keywords:** coxsackievirus B4; recombinant rVP1; molecular characterization; bioinformatics

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Hadj Hassine, I.; Gharbi, J.; Amara, I.; Alyami, A.; Subei, R.; Almalki, M.; Hober, D.; M'hadheb, M.B. Cloning and Molecular Characterization of the Recombinant CVB4E2 Immunogenic Viral Protein (rVP1), as a Potential Subunit Protein for Vaccine and Immunodiagnostic Reagent Candidate. *Microorganisms* **2023**, *11*, 1192. https://doi.org/10.3390/microorganisms11051192

Academic Editor: David J Allen

Received: 14 February 2023

Revised: 19 March 2023

Accepted: 24 March 2023

Published: 1 May 2023

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2023 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

## 1. Introduction

Coxsackievirus B4 (CVB4) belongs to the genus *Enterovirus* of the subfamily of *Picornaviridae*. It has a single-stranded RNA positive-sense genome of approximately 7400 nucleotides long with a single open reading frame (ORF) flanked by untranslated regions (UTRs) at the 5' and 3' ends [1,2]. CVB4 is a non-enveloped small virus with about 30 nm in diameter [3,4]. The genome is translated into a single long polyprotein encoding four structural proteins, VP1 to VP4 [5,6]. The VP1 protein is the most exposed and contains major neutralization epitopes and virulence determinants. This protein has been used by several researchers to study virus evolution, serotype identification, and virulence [4,7,8]. A Diabetogenic variant of CVB4 was isolated in 1979 by Yoon and his colleagues from a 10-year-old boy. The post-mortem collected data and performed studies suggest that the boy's diabetes was virus-induced. It is a virus related to a diabetogenic variant derived from CVB4. Tumors induced in mice by subsequent inoculation are hyperglycemia, beta-cell necrosis, and insulitis. The diabetogenic strain is known as the coxsackievirus B4 E2 strain (CVB4

<footnote n="1">¹ Virology and Antiviral Strategies Research Unit, Institute of Biotechnology, University of Monastir, BP74, Monastir 5000, Tunisia</footnote>

<footnote n="2">² Department of Biological Sciences, College of Science, King Faisal University, P.O. Box 380, Al-Ahsa 31982, Saudi Arabia</footnote>

<footnote n="3">³ Laboratoire de Virologie ULR3610, Université de Lille et CHU de Lille, 59000 Lille, France</footnote>

<footnote>* Correspondence: jagharbi@kfu.edu.sa (J.G.); benmhadhebmanel@yahoo.fr (M.B.M.)</footnote>
<figure><img src="image_2.png" /></figure>

Article

# Cloning and Molecular Characterization of the Recombinant CVB4E2 Immunogenic Viral Protein (rVP1), as a Potential Subunit Protein for Vaccine and Immunodiagnostic Reagent Candidate

Ikbel Hadj Hassine<fnref n="1" />, Jawhar Gharbi <fnref n="2" />,*, Imene Amara<fnref n="1" />, Ameera Alyami <fnref n="2" />, Reem Subei <fnref n="2" />, Mohammed Almalki <fnref n="2" />, Didier Hober <fnref n="3" /> and Manel Ben M'hadheb <fnref n="1" />,*

**Abstract:** The aim of the present study was, first, to clone the VP1 gene of the human coxsackievirus B4 strain E2 (CVB4E2) in the prokaryotic pUC19 plasmid expression vector then to compare it with the structural capsid proteins of the same strain using bioinformatic tools. PCR colony amplification followed through a restriction digestion analysis and sequencing process which affirmed the success of the cloning process. SDS-PAGE and Western Blotting were used to characterize the purified recombinant viral protein expressed in bacteria cells. The BLASTN tool revealed that the nucleotide sequence of the recombinant VP1 (rVP1) expressed by pUC19 highly matched the target nucleotide sequence of the diabetogenic CVB4E2 strain. Secondary structure and three-dimension structure prediction suggested that rVP1, such as wild-type VP1, is chiefly composed of random coils and a high percentage of exposed amino acids. Linear B-cell epitope prediction showed that several antigenic epitopes are likely present in rVP1 and CVB4E2 VP1 capsid protein. Additionally, phosphorylation site prediction revealed that both proteins may affect the signal transduction of host cells and can be involved in virus virulence. The present work highlights the usefulness of cloning and bioinformatics characterizations for gene investigation. Furthermore, the collected data are helpful for future experimental research related to the development of immunodiagnostic reagents and subunit vaccines based on the expression of immunogenic viral capsid proteins.

**Keywords:** coxsackievirus B4; recombinant rVP1; molecular characterization; bioinformatics

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Hadj Hassine, I.; Gharbi, J.; Amara, I.; Alyami, A.; Subei, R.; Almalki, M.; Hober, D.; M'hadheb, M.B. Cloning and Molecular Characterization of the Recombinant CVB4E2 Immunogenic Viral Protein (rVP1), as a Potential Subunit Protein for Vaccine and Immunodiagnostic Reagent Candidate. *Microorganisms* **2023**, *11*, 1192. https://doi.org/10.3390/microorganisms11051192

Academic Editor: David J Allen

Received: 14 February 2023

Revised: 19 March 2023

Accepted: 24 March 2023

Published: 1 May 2023

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2023 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

## 1. Introduction

Coxsackievirus B4 (CVB4) belongs to the genus *Enterovirus* of the subfamily of *Picornaviridae*. It has a single-stranded RNA positive-sense genome of approximately 7400 nucleotides long with a single open reading frame (ORF) flanked by untranslated regions (UTRs) at the 5' and 3' ends [1,2]. CVB4 is a non-enveloped small virus with about 30 nm in diameter [3,4]. The genome is translated into a single long polyprotein encoding four structural proteins, VP1 to VP4 [5,6]. The VP1 protein is the most exposed and contains major neutralization epitopes and virulence determinants. This protein has been used by several researchers to study virus evolution, serotype identification, and virulence [4,7,8]. A Diabetogenic variant of CVB4 was isolated in 1979 by Yoon and his colleagues from a 10-year-old boy. The post-mortem collected data and performed studies suggest that the boy's diabetes was virus-induced. It is a virus related to a diabetogenic variant derived from CVB4. Tumors induced in mice by subsequent inoculation are hyperglycemia, beta-cell necrosis, and insulitis. The diabetogenic strain is known as the coxsackievirus B4 E2 strain (CVB4

<footnote n="1">¹ Virology and Antiviral Strategies Research Unit, Institute of Biotechnology, University of Monastir, BP74, Monastir 5000, Tunisia</footnote>

<footnote n="2">² Department of Biological Sciences, College of Science, King Faisal University, P.O. Box 380, Al-Ahsa 31982, Saudi Arabia</footnote>

<footnote n="3">³ Laboratoire de Virologie ULR3610, Université de Lille et CHU de Lille, 59000 Lille, France</footnote>

<footnote>* Correspondence: jagharbi@kfu.edu.sa (J.G.); benmhadhebmanel@yahoo.fr (M.B.M.)</footnote>
E2) [9]. Polymerase chain reaction (PCR) can amplify small DNA fragments [10]. For this, a primer-based PCR method is crucial to obtain the complete cDNA of the VP1 gene. The amplified gene was subjected to cloning approaches using the pUC19 expression plasmid vector. This small vector can usually carry genes that may benefit the survival of host organisms such as antibiotic resistance [11]. In addition, the pUC19 plasmid is a high copy number plasmid, thereby improving the yield of cloned foreign DNA [12,13]. The sequencing tool is the principal method used in virology to determine and examine viruses. During this study, we employed traditional Sanger sequencing. After that, we performed a bioinformatic analysis to collect useful information for future work.

The present work highlights the usefulness of interest protein cloning and bioin-
formatics for gene investigation. Furthermore, the collected data are helpful for future
experimental research related to the development of immunodiagnostic reagents and the
development of subunit vaccines.

**2. Materials and Methods**

2.1. *Viral Strain, Bacterial Strain, Plasmid and Growth Medium*

The CVB4 E2 diabetogenic strain [9], *Escherichia coli* (E. coli) DH5α strain (F-φ80lacZΔM15 Δ(lacZYA-argF) U169 *recA1 endA1 hsdR17 (r_k-, m_k+) phoA supE44 thi-1 gyrA96 relA1 λ^Invitrogen, Waltham, MD, USA) and pUC19 plasmid (Invitrogen) [14] were used. The viral strain (kindly provided by J. W. Yoon, Julia McFarlane Diabetes Research Centre, Calgary, Alberta, AB, Canada) was propagated in the HeLa cell line from BioWhittaker in Eagle's Minimal Essential Medium (MEM; Gibco BRL, Milano, Italy) supplemented by 10% fetal calf serum (FCS; Gibco BRL), 1% (2 mM) L-glutamine (BioWhittaker, Atlanta, GA, USA), 1% non-essential amino-acids (Gibco BRL), 50 µg/mL streptomycin, 50 IU/mL penicillin (Gibco BRL), and 0.05% (2.5 µg/mL) fungizone (Amphotericin B, Apothecon, Dabhasa, India). The supernatants were collected 3 days after inoculation, clarified at 4000× g for 10 min, divided into aliquots, and finally stored at -80 °C. Virus titer in was stocks were determined on HeLa cells by limiting dilution assay for 50% tissue culture infectious doses (TCID₅₀) using the method of Reed and Muench [15]. The bacterial strain was cultured in Luria-Bertani (LB; Invitrogen) or in LB containing 1.5% agar (Invitrogen). When required, the media was supplemented with ampicillin (Amp; 100 µg/mL).

2.2. *PCR Amplification, Cloning and Sequencing of the CVB4 VP1 Gene*

Total RNA (the genome) was completely extracted using the single-step method based on guanidium isothiocyanate or Tri-Reagent® (Sigma, Taufkirchen, Germany) [16] or the EZ-10 Total RNA Mini-Preps kit (Invitrogen) according to the manufacturer's instructions. The cap-sid VP1 gene of the CVB4 E2 diabetogenic strain was amplified using primers that flanked the VP1 region: CVB4E2-VP1-F: 5'-TTTTATGAAGGGCCAACAGAG-3' (nucleotides 2440-2460, relative to strain CVB4 E2) and CVB4 E2-VP1-R: 5'-CTCGGTGGTGACAGCCTCAA-3' (nucleotides 3251-3270, relative to strain CVB4E2). A pre-denaturation for 3 min at 94 °C was applied then followed by 35 cycles of amplification consisting of 1 min denaturation at 94 °C, 1 min annealing at 54 °C, and 1 min extension at 72 °C. A final extension step was performed at 72 °C for 10 min. Another couple of primers were used to add restriction sites to the VP1 gene: 5'-TATCGAATTCTTTTATGAAGGGCCAACAGAG-3' and 5'-TATCGGATCCCTCGGTGGTGACAGCCTCAA-3'. The underlined sequences represent EcoRI and BamHI restriction sites, respectively. To amplify VP1 containing the restriction sites, an initial denaturing step of 3 min at 94 °C was performed and subjected to 35 cycles of denaturation at 94 °C for 1 min, annealing at 65 °C for 1 min, and an extension at 72 °C for 1 min, followed by a final extension step of 10 min at 72 °C. Then the PCR products were verified on a 2% agarose gel and stained with ethidium bromide. The amplified VP1 band was digested with EcoRI and BamHI (TaKaRa, Kusatsu, Japan), cleaned up using "EZ-10 Spin Column DNA PCR Products Purification" (BIO BASIC, Singapore) and ligated using a T4 DNA ligase (Roche, Basel, Switzerland) with the prokaryotic expression vector pUC19, previously digested with the same restriction enzymes. The competent E. coli DH5α cells were previously prepared by
E2) [9]. Polymerase chain reaction (PCR) can amplify small DNA fragments [10]. For this, a primer-based PCR method is crucial to obtain the complete cDNA of the VP1 gene. The amplified gene was subjected to cloning approaches using the pUC19 expression plasmid vector. This small vector can usually carry genes that may benefit the survival of host organisms such as antibiotic resistance [11]. In addition, the pUC19 plasmid is a high copy number plasmid, thereby improving the yield of cloned foreign DNA [12,13]. The sequencing tool is the principal method used in virology to determine and examine viruses. During this study, we employed traditional Sanger sequencing. After that, we performed a bioinformatic analysis to collect useful information for future work.

The present work highlights the usefulness of interest protein cloning and bioin-
formatics for gene investigation. Furthermore, the collected data are helpful for future
experimental research related to the development of immunodiagnostic reagents and the
development of subunit vaccines.

**2. Materials and Methods**

2.1. *Viral Strain, Bacterial Strain, Plasmid and Growth Medium*

The CVB4 E2 diabetogenic strain [9], *Escherichia coli* (E. coli) DH5α strain (F-φ80lacZΔM15 Δ(lacZYA-argF) U169 recA1 endA1 hsdR17 (r_k-, m_k+) phoA supE44 thi-1 gyrA96 relA1 λ^Invitrogen, Waltham, MD, USA) and pUC19 plasmid (Invitrogen) [14] were used. The viral strain (kindly provided by J. W. Yoon, Julia McFarlane Diabetes Research Centre, Calgary, Alberta, AB, Canada) was propagated in the HeLa cell line from BioWhittaker in Eagle's Minimal Essential Medium (MEM; Gibco BRL, Milano, Italy) supplemented by 10% fetal calf serum (FCS; Gibco BRL), 1% (2 mM) L-glutamine (BioWhittaker, Atlanta, GA, USA), 1% non-essential amino-acids (Gibco BRL), 50 µg/mL streptomycin, 50 IU/mL penicillin (Gibco BRL), and 0.05% (2.5 µg/mL) fungizone (Amphotericin B, Apothecon, Dabhasa, India). The supernatants were collected 3 days after inoculation, clarified at 4000× g for 10 min, divided into aliquots, and finally stored at -80 °C. Virus titer in was stocks were determined on HeLa cells by limiting dilution assay for 50% tissue culture infectious doses (TCID₅₀) using the method of Reed and Muench [15]. The bacterial strain was cultured in Luria-Bertani (LB; Invitrogen) or in LB containing 1.5% agar (Invitrogen). When required, the media was supplemented with ampicillin (Amp; 100 µg/mL).

2.2. *PCR Amplification, Cloning and Sequencing of the CVB4 VP1 Gene*

Total RNA (the genome) was completely extracted using the single-step method based on guanidium isothiocyanate or Tri-Reagent® (Sigma, Taufkirchen, Germany) [16] or the EZ-10 Total RNA Mini-Preps kit (Invitrogen) according to the manufacturer's instructions. The cap-sid VP1 gene of the CVB4 E2 diabetogenic strain was amplified using primers that flanked the VP1 region: CVB4E2-VP1-F: 5'-TTTTATGAAGGGCCAACAGAG-3' (nucleotides 2440-2460, relative to strain CVB4 E2) and CVB4 E2-VP1-R: 5'-CTCGGTGGTGACAGCCTCAA-3' (nucleotides 3251-3270, relative to strain CVB4E2). A pre-denaturation for 3 min at 94 °C was applied then followed by 35 cycles of amplification consisting of 1 min denaturation at 94 °C, 1 min annealing at 54 °C, and 1 min extension at 72 °C. A final extension step was performed at 72 °C for 10 min. Another couple of primers were used to add restriction sites to the VP1 gene: 5'-TATCGAATTCTTTTATGAAGGGCCAACAGAG-3' and 5'-TATCGGATCCCTCGGTGGTGACAGCCTCAA-3'. The underlined sequences represent EcoRI and BamHI restriction sites, respectively. To amplify VP1 containing the restriction sites, an initial denaturing step of 3 min at 94 °C was performed and subjected to 35 cycles of denaturation at 94 °C for 1 min, annealing at 65 °C for 1 min, and an extension at 72 °C for 1 min, followed by a final extension step of 10 min at 72 °C. Then the PCR products were verified on a 2% agarose gel and stained with ethidium bromide. The amplified VP1 band was digested with EcoRI and BamHI (TaKaRa, Kusatsu, Japan), cleaned up using "EZ-10 Spin Column DNA PCR Products Purification" (BIO BASIC, Singapore) and ligated using a T4 DNA ligase (Roche, Basel, Switzerland) with the prokaryotic expression vector pUC19, previously digested with the same restriction enzymes. The competent E. coli DH5α cells were previously prepared by
the calcium chloride method as described elsewhere [17]. The recombinant clones resulted from the transformation of the ligation products in the competent cells selected by the PCR colony. For this, each colony was removed from LB plates supplemented with Amp using a sterile tip and dissolved in the appropriate PCR mixture. The insert-specific primers were used. After amplification of 40 cycles, using the same PCR protocol previously applied, the PCR products were visualized by a 1% agarose gel electrophoresis. The recombinant plasmid pUC19-VP1 purification by an "EZ-10 Spin Column Plasmid DNA Kit" (BIO BASIC) was followed by restriction digestion with 15 units of *Eco*RI and *Bam*HI (TaKaRa) for 2 h at 37 °C, and the recombinant VP1 (rVP1) was further sequenced using AB PRISM™ 310 Genetic Analyser (Applied Biosystems, Waltham, MD, USA), according to the Sanger method [18].

## 2.3. Preparation of Cell Lysates for SDS-PAGE

For the preparation of crude cell extracts, 1.5 mL of each grown broth cultures from samples incubated overnight were taken in separate 2 mL Eppendorf tubes and centrifuged at full speed (18,000 × g at 4 °C) for 1 min to obtain a pellet. The supernatant was discarded and the pellet was resuspended in 100 µL of 10% Sodium dodecyl sulphate (SDS) and refrigerated (−20 °C) for 5–10 min to lyse the cells. The same sample incubated overnight post induction was repeatedly frozen and thawed several times after SDS treatment. After thawing the samples were again centrifuged at full speed for 1 min. Fifty µL of the clear supernatant was taken from each tube and placed respectively into four new 2 mL Eppendorf tubes. This was followed by the addition of 50 µL of 2 × Laemmli buffer containing loading dye to each tube and mixed by pipetting. After heating for 10 min at 95 °C, the samples (10 µL each) were loaded in separate wells in 15% SDS-Polyacrylamide gel.

## 2.4. SDS-PAGE and Western Blot Analysis of rVP1 Subunit Protein

Protein samples were separated and analyzed in SDS-PAGE as described by Laemmli (1970). This was performed at room temperature in 15% polyacrylamide gels, using the Mini-Protean tetra cell 4-gel hand casting system (BIO-RAD, China). Ten µL of the sample was loaded in each well along with 5 µL of pre-stained protein ladder (10–250 kD). The gel was run at 25 mA for approximately 2 h, stained with Coomassie Blue staining solution and distained with distain solution (40 mL methanol added to 7 mL glacial acetic acid and 53 mL distilled water to make a final volume of 100 mL), until the gel background was clear. Recombinant CVB4 Viral Protein (rVP1) was resolved on 15% SDS-PAGE and transferred onto nitrocellulose membrane (0.45 µm, SIGMA). After transfer, the membrane was blocked by 5% non-fat dry milk in PBS for 1 h and probed by Mouse anti-enterovirus clone 5-D8/1 (DAKO) diluted 1:500, membranes were thoroughly washed with PBS followed by incubation with the secondary antibody goat anti-mouse Ig-Peroxidase conjugate (SIGMA) with 1:5000 dilution. The membranes were developed with a 4-Chloro-1-Naphthol Solution or by chemiluminescence and colorimetric detection kit (BIO-RAD). All the dilutions were realized with the blocking solution.

## 2.5. Purification of rVP1 Subunit Protein

The pellet obtained from the recombinant bacterial culture was treated with lysis CellyticB buffer (SIGMA), as per the protocol. Then batch purification of recombinant rVP1 subunit protein from the *E. coli* DH5α strain was carried out under native conditions using Ni-NTA affinity chromatography (QIAGEN, Hilden, Germany). The supernatant was carefully removed and 10–15 µL was checked for the presence of the concerned protein by SDS-PAGE analysis. The supernatant was subjected to purification using Ni-NTA affinity chromatography. One mL of 50% Ni-NTA slurry was added to 4 mL cleared lysate and mixed gently by shaking it at 4 °C for 60 min. The lysate-Ni-NTA mixture was loaded into a column with the bottom outlet capped. The bottom cap was removed and column flow-through was collected. The flow-through was saved for the SDS-PAGE analysis. The column was washed twice with 4 mL wash buffer and the collected wash fractions were also analyzed using SDS-PAGE. One mL elution buffer was used to elute the protein four
the calcium chloride method as described elsewhere [17]. The recombinant clones resulted from the transformation of the ligation products in the competent cells selected by the PCR colony. For this, each colony was removed from LB plates supplemented with Amp using a sterile tip and dissolved in the appropriate PCR mixture. The insert-specific primers were used. After amplification of 40 cycles, using the same PCR protocol previously applied, the PCR products were visualized by a 1% agarose gel electrophoresis. The recombinant plasmid pUC19-VP1 purification by an "EZ-10 Spin Column Plasmid DNA Kit" (BIO BASIC) was followed by restriction digestion with 15 units of *Eco*RI and *Bam*HI (TaKaRa) for 2 h at 37 °C, and the recombinant VP1 (rVP1) was further sequenced using AB PRISM™ 310 Genetic Analyser (Applied Biosystems, Waltham, MD, USA), according to the Sanger method [18].

## 2.3. Preparation of Cell Lysates for SDS-PAGE

For the preparation of crude cell extracts, 1.5 mL of each grown broth cultures from samples incubated overnight were taken in separate 2 mL Eppendorf tubes and centrifuged at full speed (18,000 × g at 4 °C) for 1 min to obtain a pellet. The supernatant was discarded and the pellet was resuspended in 100 µL of 10% Sodium dodecyl sulphate (SDS) and refrigerated (−20 °C) for 5–10 min to lyse the cells. The same sample incubated overnight post induction was repeatedly frozen and thawed several times after SDS treatment. After thawing the samples were again centrifuged at full speed for 1 min. Fifty µL of the clear supernatant was taken from each tube and placed respectively into four new 2 mL Eppendorf tubes. This was followed by the addition of 50 µL of 2 × Laemmli buffer containing loading dye to each tube and mixed by pipetting. After heating for 10 min at 95 °C, the samples (10 µL each) were loaded in separate wells in 15% SDS-Polyacrylamide gel.

## 2.4. SDS-PAGE and Western Blot Analysis of rVP1 Subunit Protein

Protein samples were separated and analyzed in SDS-PAGE as described by Laemmli (1970). This was performed at room temperature in 15% polyacrylamide gels, using the Mini-Protean tetra cell 4-gel hand casting system (BIO-RAD, China). Ten µL of the sample was loaded in each well along with 5 µL of pre-stained protein ladder (10–250 kD). The gel was run at 25 mA for approximately 2 h, stained with Coomassie Blue staining solution and distained with distain solution (40 mL methanol added to 7 mL glacial acetic acid and 53 mL distilled water to make a final volume of 100 mL), until the gel background was clear. Recombinant CVB4 Viral Protein (rVP1) was resolved on 15% SDS-PAGE and transferred onto nitrocellulose membrane (0.45 µm, SIGMA). After transfer, the membrane was blocked by 5% non-fat dry milk in PBS for 1 h and probed by Mouse anti-enterovirus clone 5-D8/1 (DAKO) diluted 1:500, membranes were thoroughly washed with PBS followed by incubation with the secondary antibody goat anti-mouse Ig-Peroxidase conjugate (SIGMA) with 1:5000 dilution. The membranes were developed with a 4-Chloro-1-Naphthol Solution or by chemiluminescence and colorimetric detection kit (BIO-RAD). All the dilutions were realized with the blocking solution.

## 2.5. Purification of rVP1 Subunit Protein

The pellet obtained from the recombinant bacterial culture was treated with lysis CellYTicB buffer (SIGMA), as per the protocol. Then batch purification of recombinant rVP1 subunit protein from the *E. coli* DH5α strain was carried out under native conditions using Ni-NTA affinity chromatography (QIAGEN, Hilden, Germany). The supernatant was carefully removed and 10–15 µL was checked for the presence of the concerned protein by SDS-PAGE analysis. The supernatant was subjected to purification using Ni-NTA affinity chromatography. One mL of 50% Ni-NTA slurry was added to 4 mL cleared lysate and mixed gently by shaking it at 4 °C for 60 min. The lysate-Ni-NTA mixture was loaded into a column with the bottom outlet capped. The bottom cap was removed and column flow-through was collected. The flow-through was saved for the SDS-PAGE analysis. The column was washed twice with 4 mL wash buffer and the collected wash fractions were also analyzed using SDS-PAGE. One mL elution buffer was used to elute the protein four
times. The purity of the recombinant protein rVP1 in different fractions was checked by electrophoresis on 15% SDS-PAGE gel after staining with Coomassie brilliant blue stain.

## 2.6. *Bioinformatic Analysis*

For sequences homology and conserved domains analysis, NCBI BLASTN (https://blast.ncbi.nlm.nih.gov/Blast.cgi?PAGE_TYPE=BlastSearch, accessed on 15 January 2023) and Conserved Domain Database (CDD) (https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi, accessed on 15 January 2023) [19] were used. For further proteomic analysis, the nucleotide sequence was converted to an amino-acids (AA) format by the EMBOSS Seqret tool (https://www.ebi.ac.uk/Tools/sfc/emboss_seqret/, accessed on 15 January 2023). The ProtParam tool (https://web.expasy.org/protparam/, accessed on 15 January 2023) was utilized to analyze the physicochemical properties of the rVP1 and CVB4 VP1-VP4. The properties include the relative molecular weights, the atomic composition of VP proteins' AA, molar extinction coefficient, theoretical isoelectric point (Pi), hydrophilic parameters, instability index (II), and so on [20]. SignalP-4.0 server (http://www.cbs.dtu.dk/services/SignalP-4.1/, accessed on 15 January 2023) was used for predicting the signal peptides [21]; NetPhos3.1 server (http://www.cbs.dtu.dk/services/NetPhos/, accessed on 15 January 2023) was utilized to determine phosphorylation sites [22,23]. InterProScan (https://www.ebi.ac.uk/interpro/, accessed on 15 January 2023), and SMART (http://smart.embl-heidelberg.de/, accessed on 15 January 2023) were employed for the analysis of the structural and functional domains [24,25]. The secondary structure of each VP protein and the solvent accessibility of its AA residues were predicted by SOPMA (https://npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=/NPSA/npsa_sopma.html, accessed on 15 January 2023), and Predictprotein (https://www Predictprotein (http://smart.embl-heidelberg.de/, accessed on 15 January 2023) were employed for the analysis of the structural and functional domains [24,25]. The secondary structure of each VP protein and the solvent accessibility of its AA residues were predicted by SOPMA (https://npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=/NPSA/npsa_sopma.html, accessed on 15 January 2023), and Predictprotein (https://www predictprotein.org/, accessed on 15 January 2023) websites. Secondary structure prediction is useful for 3D structure prediction [26]. Protein modeling was determined by the automatic mode of SWISS-MODEL (https://swissmodel.expasy.org/, accessed on 15 January 2023), and model structures were visualized using RasMol 2.7.5 Software (http://www.rasmol.org/software/RasMol_2.7.5/, accessed on 15 January 2023). BepiPred-2.0 was used to predict β-cell epitopes of the proteins (http://www.cbs.dtu.dk/services/BepiPred/, accessed on 15 January 2023).

# 3. Results

## 3.1. *PCR Amplification, Cloning, and Sequencing of the rVP1 Gene*

The amplification of the VP1 gene, with or without EcoRI/BamHI restriction sites, was carried out using the diabetogenic CVB4E2 strain. No specific band was amplified from the mock (Figure 1A, Lane 1). Whereas, a target fragment of approximately 850 bp was amplified and then purified from the positive sample (Figure 1A, Lanes 2 and 3, respectively). Before proceeding to the cloning step, a purification process is essential to remove organic and inorganic materials and undesired chemical compounds. No specific bands were observed from the mock (Figure 1B, Lanes 1 and 3). After purification of VP1 with or without restriction sites, the same bands are observed as in Figure 1 (Figure 1B, Lanes 2 and 4, respectively). To clone the DNA fragment into prokaryotic expression vector pUC19 to yield pUC19-VP1, first, pUC19 was first extracted and purified from E. coli transformed colonies, thus giving the expected fragment-size (Figure 2A, Lane 2). The purification was followed by restriction enzyme digestion using EcoRI and BamHI restriction enzymes, giving a fragment with a different size (Figure 2A, Lane 1) compared to the original (Figure 2A, Lane 2). To end the cloning by a ligation reaction between the flanked and purified VP1, the digested pUC19 was also purified (Figure 2B, Lane 1). To confirm the cloning process, a PCR colony was amplified from the transformed colonies (Figure 3). As expected, the amplification of negative control (Figure 3A, Lane 1) and non-transformed colonies (Figure 3A, Lanes 2 to 6) did not show the target fragment. A fragment of approximately 850 bp was found predicting the presence of the cloned DNA fragment in the recombinant colony (Figure 3A, Lane 7). After that, the DNA from the
times. The purity of the recombinant protein rVP1 in different fractions was checked by electrophoresis on 15% SDS-PAGE gel after staining with Coomassie brilliant blue stain.

## 2.6. Bioinformatic Analysis

For sequences homology and conserved domains analysis, NCBI BLASTN (https://blast.ncbi.nlm.nih.gov/Blast.cgi?PAGE_TYPE=BlastSearch, accessed on 15 January 2023) and Conserved Domain Database (CDD) (https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi, accessed on 15 January 2023) [19] were used. For further proteomic analysis, the nucleotide sequence was converted to an amino-acids (AA) format by the EMBOSS Seqret tool (https://www.ebi.ac.uk/Tools/sfc/emboss_seqret/, accessed on 15 January 2023). The ProtParam tool (https://web.expasy.org/protparam/, accessed on 15 January 2023) was utilized to analyze the physicochemical properties of the rVP1 and CVB4 VP1-VP4. The properties include the relative molecular weights, the atomic composition of VP proteins' AA, molar extinction coefficient, theoretical isoelectric point (Pi), hydrophilic parameters, instability index (II), and so on [20]. SignalP-4.0 server (http://www.cbs.dtu.dk/services/SignalP-4.1/, accessed on 15 January 2023) was used for predicting the signal peptides [21]; NetPhos3.1 server (http://www.cbs.dtu.dk/services/NetPhos/, accessed on 15 January 2023) was utilized to determine phosphorylation sites [22,23]. InterProScan (https://www.ebi.ac.uk/interpro/, accessed on 15 January 2023), and SMART (http://smart.embl-heidelberg.de/, accessed on 15 January 2023) were employed for the analysis of the structural and functional domains [24,25]. The secondary structure of each VP protein and the solvent accessibility of its AA residues were predicted by SOPMA (https://npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=/NPSA/npsa_sopma.html, accessed on 15 January 2023), and Predictprotein (https://www Predictprotein (http://smart.embl-heidelberg.de/, accessed on 15 January 2023) were employed for the analysis of the structural and functional domains [24,25]. The secondary structure of each VP protein and the solvent accessibility of its AA residues were predicted by SOPMA (https://npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=/NPSA/npsa_sopma.html, accessed on 15 January 2023), and Predictprotein (https://www predictprotein.org/, accessed on 15 January 2023) websites. Secondary structure prediction is useful for 3D structure prediction [26]. Protein modeling was determined by the automatic mode of SWISS-MODEL (https://swissmodel.expasy.org/, accessed on 15 January 2023), and model structures were visualized using RasMol 2.7.5 Software (http://www.rasmol.org/software/RasMol_2.7.5/, accessed on 15 January 2023). BepiPred-2.0 was used to predict β-cell epitopes of the proteins (http://www.cbs.dtu.dk/services/BepiPred/, accessed on 15 January 2023).

# 3. Results

## 3.1. PCR Amplification, Cloning, and Sequencing of the rVP1 Gene

The amplification of the VP1 gene, with or without EcoRI/BamHI restriction sites, was carried out using the diabetogenic CVB4E2 strain. No specific band was amplified from the mock (Figure 1A, Lane 1). Whereas, a target fragment of approximately 850 bp was amplified and then purified from the positive sample (Figure 1A, Lanes 2 and 3, respectively). Before proceeding to the cloning step, a purification process is essential to remove organic and inorganic materials and undesired chemical compounds. No specific bands were observed from the mock (Figure 1B, Lanes 1 and 3). After purification of VP1 with or without restriction sites, the same bands are observed as in Figure 1 (Figure 1B, Lanes 2 and 4, respectively). To clone the DNA fragment into prokaryotic expression vector pUC19 to yield pUC19-VP1, first, pUC19 was first extracted and purified from E. coli transformed colonies, thus giving the expected fragment-size (Figure 2A, Lane 2). The purification was followed by restriction enzyme digestion using EcoRI and BamHI restriction enzymes, giving a fragment with a different size (Figure 2A, Lane 1) compared to the original (Figure 2A, Lane 2). To end the cloning by a ligation reaction between the flanked and purified VP1, the digested pUC19 was also purified (Figure 2B, Lane 1). To confirm the cloning process, a PCR colony was amplified from the transformed colonies (Figure 3). As expected, the amplification of negative control (Figure 3A, Lane 1) and non-transformed colonies (Figure 3A, Lanes 2 to 6) did not show the target fragment. A fragment of approximately 850 bp was found predicting the presence of the cloned DNA fragment in the recombinant colony (Figure 3A, Lane 7). After that, the DNA from the
same colony was extracted, attesting to the presence of pUC19-VP1 (Figure 3B, Lane 1). Furthermore, a restriction analysis showed two bands, each corresponding to pUC19 and the cloned VP1 (Figure 3C, Lane 1). DNA sequencing confirmed the success of the cloning process. The nucleotide sequence has been submitted in GenBank under the Accession Nb. MW193767.

<figure><img src="image_5.png" /></figure>

**Figure 1.** PCR amplification and Purification of the VP1 gene. (A) Lane M: GeneRuler 100 bp DNA Ladder (ThermoFisher, Waltham, MD, USA); Lane 1: negative control; Lane 2: PCR amplification product of the VP1 gene; Lane 3: PCR amplification product of VP1 with primers containing *Eco*RI and *Bam*HI recognition sites. (B) Lane M: GeneRuler 100 bp DNA Ladder (ThermoFisher); Lanes 1 and 3: negative control; Lanes 2 and 4: Purification products of VP1 and VP1 flanked by the restriction enzyme sites, respectively.

<figure><img src="image_6.png" /></figure>

**Figure 2.** Preparation of pUC19 for the cloning process. (A) Lane M: Lambda DNA/HindIII Marker (ThermoFisher); Lane 1: Digestion product of purified pUC19; Lane 2: Purified pUC19 extracted from transformed bacterial cells. (B) Lane M: Lambda DNA/HindIII Marker (ThermoFisher); Lane 1: Purification product of the digested pUC19.
same colony was extracted, attesting to the presence of pUC19-VP1 (Figure 3B, Lane 1). Furthermore, a restriction analysis showed two bands, each corresponding to pUC19 and the cloned VP1 (Figure 3C, Lane 1). DNA sequencing confirmed the success of the cloning process. The nucleotide sequence has been submitted in GenBank under the Accession Nb. MW193767.

<figure><img src="image_5.png" /></figure>

**Figure 1.** PCR amplification and Purification of the VP1 gene. (A) Lane M: GeneRuler 100 bp DNA Ladder (ThermoFisher, Waltham, MD, USA); Lane 1: negative control; Lane 2: PCR amplification product of the VP1 gene; Lane 3: PCR amplification product of VP1 with primers containing *Eco*RI and *Bam*HI recognition sites. (B) Lane M: GeneRuler 100 bp DNA Ladder (ThermoFisher); Lanes 1 and 3: negative control; Lanes 2 and 4: Purification products of VP1 and VP1 flanked by the restriction enzyme sites, respectively.

<figure><img src="image_6.png" /></figure>

**Figure 2.** Preparation of pUC19 for the cloning process. (A) Lane M: Lambda DNA/HindIII Marker (ThermoFisher); Lane 1: Digestion product of purified pUC19; Lane 2: Purified pUC19 extracted from transformed bacterial cells. (B) Lane M: Lambda DNA/HindIII Marker (ThermoFisher); Lane 1: Purification product of the digested pUC19.
<figure><img src="image_7.png" /><figcaption>**Figure 3.** Confirmation of the cloning process. (A) Lane M: GeneRuler 100 bp DNA Ladder (ThermoFisher); Lane 1: Negative control; Lanes 2 to 6: PCR colony of Non-transformed colonies; Lane 7: PCR colony of a transformed colony. (B) Lane M: Lambda DNA/HindIII Marker (ThermoFisher); Lane 1: The extracted DNA of the same colony confirmed the cloning of VP1 in pUC19. (C) Lane M: Lambda DNA/HindIII Marker (ThermoFisher); Lane 1: The restriction analysis showing two bands corresponding to pUC19 and VP1.</figcaption></figure>

### 3.2. SDS-PAGE Analysis, Western Blotting and Purification of the rVP1 Subunit Protein

The molecular mass of recombinant rVP1 subunit protein, obtained through the separation of different bacterial proteins from cell lysates through using a 15% SDS-PAGE technique was estimated to be approximately 30 kD, which is very close to the expected molecular weight. Recombinant *E. coli* DH5α host cells expressing rVP1 subunit protein, subjected to overnight incubation as well as repeated freeze-thawing post-induction, lead to increased yield of recombinant VP1 subunit viral protein (Figure 4A).

<figure><img src="image_8.png" /><figcaption>**Figure 4.** Molecular characterization of the recombinant CVB4 E2 rVP1 subunit viral protein. (A) SDS-PAGE analysis of the recombinant rVP1 protein. MW: Protein Weight Marker; 1: Control DH5α cells; 2–3: Transformed DH5α cells expressing VP1 gene. (B) Purification of recombinant rVP1 subunit protein using Ni-NTA affinity chromatography. 1: Clarified Bacteria cells expressing VP1; 2–3: Wash fractions; 4–6: Elution fractions. (C) Western Blotting analysis of the recombinant rVP1 protein. 1: Negative control; 2–3: Non-transformed bacteria cell lysates; 4: Transformed DH5α cell lysate expressing VP1 gene.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 3.** Confirmation of the cloning process. (A) Lane M: GeneRuler 100 bp DNA Ladder (ThermoFisher); Lane 1: Negative control; Lanes 2 to 6: PCR colony of Non-transformed colonies; Lane 7: PCR colony of a transformed colony. (B) Lane M: Lambda DNA/HindIII Marker (ThermoFisher); Lane 1: The extracted DNA of the same colony confirmed the cloning of VP1 in pUC19. (C) Lane M: Lambda DNA/HindIII Marker (ThermoFisher); Lane 1: The restriction analysis showing two bands corresponding to pUC19 and VP1.</figcaption></figure>

### 3.2. SDS-PAGE Analysis, Western Blotting and Purification of the rVP1 Subunit Protein

The molecular mass of recombinant rVP1 subunit protein, obtained through the separation of different bacterial proteins from cell lysates through using a 15% SDS-PAGE technique was estimated to be approximately 30 kD, which is very close to the expected molecular weight. Recombinant *E. coli* DH5α host cells expressing rVP1 subunit protein, subjected to overnight incubation as well as repeated freeze-thawing post-induction, lead to increased yield of recombinant VP1 subunit viral protein (Figure 4A).

<figure><img src="image_8.png" /><figcaption>**Figure 4.** Molecular characterization of the recombinant CVB4 E2 rVP1 subunit viral protein. (A) SDS-PAGE analysis of the recombinant rVP1 protein. MW: Protein Weight Marker; 1: Control DH5α cells; 2–3: Transformed DH5α cells expressing VP1 gene. (B) Purification of recombinant rVP1 subunit protein using Ni-NTA affinity chromatography. 1: Clarified Bacteria cells expressing VP1; 2–3: Wash fractions; 4–6: Elution fractions. (C) Western Blotting analysis of the recombinant rVP1 protein. 1: Negative control; 2–3: Non-transformed bacteria cell lysates; 4: Transformed DH5α cell lysate expressing VP1 gene.</figcaption></figure>
After batch purification of the recombinant rVP1 recombinant protein of CVB4E2 from *E. coli* DH5α cells using Ni-NTA affinity chromatography, a relatively purified fraction of the recombinant protein rVP1 was obtained in elution fractions 4, 5, and 6, as confirmed by the 15% SDS-PAGE analysis of different fractions (Figure 4B).

The exact identity of the expressed protein was established by subjecting the recombinant *E. coli* DH5α cells to Western transfer analysis, along with negative control (non-transformed *E. coli* DH5α). Mouse anti-enterovirus clone 5-D8/1 (DAKO) conjugated with anti-mouse Ig-Peroxidase (SIGMA) directed against VP1 viral protein picked up a band at approximately 30 kD regions from the lanes loaded with cell lysate from the host cells transformed to express the recombinant VP1 viral protein. However, no such signal was produced in lanes from non-transformed cell lysates (Figure 4C).

### 3.3. Bioinformatic Analysis of the Sequence and Structure of the rVP1 Subunit Protein

#### 3.3.1. Bioinformatic Analysis of the Nucleotide Sequence of rVP1

The NCBI BLASTN tool showed that the nucleotide sequence of the rVP1 highly matched the target nucleotide sequence of the human coxsackievirus B4 strain E2 variant (GenBank Accession Nb. AF311939). The nucleotide sequences shared 96% identity with each other. Moreover, CDD analysis indicated that the rVP1 belonged to the rhv-like superfamily and contained a conserved domain pfam00073 (Figure 5A). As expected, the gene expression leads to a coxsackievirus capsid protein.

<figure><img src="image_9.png" /><figcaption>**Figure 5.** Prediction of the conserved domain, signal peptides sequence, and phosphorylation sites of the rVP1. The conserved domain, Signal peptide sequence, and phosphorylation site of rVP1 were analyzed using NCBI Conserved Domains search tool (A), SignalP-4.0 Server (B), and NetPhos 2.0 program (C), respectively.</figcaption></figure>
After batch purification of the recombinant rVP1 recombinant protein of CVB4E2 from *E. coli* DH5α cells using Ni-NTA affinity chromatography, a relatively purified fraction of the recombinant protein rVP1 was obtained in elution fractions 4, 5, and 6, as confirmed by the 15% SDS-PAGE analysis of different fractions (Figure 4B).

The exact identity of the expressed protein was established by subjecting the recombinant *E. coli* DH5α cells to Western transfer analysis, along with negative control (non-transformed *E. coli* DH5α). Mouse anti-enterovirus clone 5-D8/1 (DAKO) conjugated with anti-mouse Ig-Peroxidase (SIGMA) directed against VP1 viral protein picked up a band at approximately 30 kD regions from the lanes loaded with cell lysate from the host cells transformed to express the recombinant VP1 viral protein. However, no such signal was produced in lanes from non-transformed cell lysates (Figure 4C).

### 3.3. Bioinformatic Analysis of the Sequence and Structure of the rVP1 Subunit Protein

#### 3.3.1. Bioinformatic Analysis of the Nucleotide Sequence of rVP1

The NCBI BLASTN tool showed that the nucleotide sequence of the rVP1 highly matched the target nucleotide sequence of the human coxsackievirus B4 strain E2 variant (GenBank Accession Nb. AF311939). The nucleotide sequences shared 96% identity with each other. Moreover, CDD analysis indicated that the rVP1 belonged to the rhv-like superfamily and contained a conserved domain pfam00073 (Figure 5A). As expected, the gene expression leads to a coxsackievirus capsid protein.

<figure><img src="image_9.png" /><figcaption>**Figure 5.** Prediction of the conserved domain, signal peptides sequence, and phosphorylation sites of the rVP1. The conserved domain, Signal peptide sequence, and phosphorylation site of rVP1 were analyzed using NCBI Conserved Domains search tool (A), SignalP-4.0 Server (B), and NetPhos 2.0 program (C), respectively.</figcaption></figure>
### 3.3.2. Bioinformatic Analysis of the Polypeptide Sequence of rVP1

Protparam tool indicated that rVP1 expressed by the pUC19 expressing vector and VP1 protein of the Human Coxsackievirus B4E2 strain have many identical or close physicochemical properties compared to the other proteins. All the proteins are hydrophilic with a positive polarity, nevertheless, have some differences. Unlike the others, VP4 is unstable. PIs of VP3 and VP4 are less than 7 and significantly different from that of rVP1 (7.8), VP1 (8.81), and VP2 (8.83). rVP1 and VP1 have similar amino acid composition, atomic composition, and molecular weight. The analysis of CVB4 rVP1 and VP1-4 by applying different tools revealed that all the capsid proteins did not include potential signal peptide sequences (Figure 5B). However, all the proteins harbored phosphorylation sites. rVP1 includes 35 potential phosphorylation sites including 21 serine, 9 threonine, and 4 tyrosine; VP1 contains 23 potential phosphorylation sites, including 14 serine, 5 threonine, and 4 tyrosine (Figure 5C).

VP2-4 have a fewer number of phosphorylation sites. The subcellular location predicted by TargetP demonstrated that all the proteins could be at any extra-mitochondrial location. The prediction of the secondary structure and the solvent accessibility of AA residues are presented in Table 1. SOPMA and Predictprotein showed/revealed that random coils were the principal conformation of the capsid proteins, a higher ratio of β-folding, and lower ratios of β-turns. PredictProtein demonstrated that most of the AA residues in the VP4 protein were intermediate to solvent (83.82%). Percentages of exposed and buried AA residues in the other VP proteins were equivalent. rVP1 and VP1 have 47.65% and 45.42% of exposed AA, respectively. The protein structure homology-modeling was determined by employing the SWISS-MODEL server. The selection of templates depends on the similarity and coverage of AA sequences. Based on evaluation parameters, high standards models were applied as the tertiary structure of the VP proteins (Table 2). Model structures were observed by using RasMol 2.7.5 Software (Figure 6). The prediction of Linear B cell epitopes of rVP1 and VP1-VP4 capsid proteins was determined by the BepiPred server (Table 3). Among the potential linear B cell epitopes, VP4 included 2 epitopes with over 10 AA residues, respectively located at AA 5-19 and AA 29-54; VP3 contained 3 epitopes, respectively situated at AA 7-39, AA 55-67 and AA 139-150; VP2 included 5 epitopes respectively situated at AA 5-27, 38-60, 66-92, 130-165, 223-244; VP1 included 7 epitopes with over 10 AA residues, respectively located at 5-31, 33-50, 82-96, 127-138, 152-165, 199-209, 251-280; rVP1 6 epitopes 18-63, 83-93, 130-140, 155-168, 202-213, 258-273 (Table 3). The model structures were observed by using RasMol 2.7.5 Software (Figure 7).

<table><caption>**Table 1.** Prediction of secondary structures and solvent accessibility for the recombinant rVP1 and the CVB4E2 capsid proteins.</caption><thead><tr><th></th><th>rVP1</th><th>VP1</th><th>VP2</th><th>VP3</th><th>VP4</th></tr></thead><tbody><tr><td></td><td></td><td></td><td colspan="3">SOPMA</td></tr><tr><td>α-helix</td><td>23.47%</td><td>20.07%</td><td>17.34%</td><td>17.15%</td><td>17.56%</td></tr><tr><td>β-folding</td><td>25.63%</td><td>23.24%</td><td>26.61%</td><td>29.29%</td><td>22.06%</td></tr><tr><td>β-turns</td><td>3.25%</td><td>2.11%</td><td>7.66%</td><td>5.44%</td><td>4.41%</td></tr><tr><td>Random coils</td><td>47.65%</td><td>54.58%</td><td>48.39%</td><td>48.12%</td><td>55.88%</td></tr><tr><td></td><td></td><td></td><td colspan="3">Predictprotein</td></tr><tr><td>α-helix</td><td>4.69%</td><td>8.45%</td><td>2.82%</td><td>8.37%</td><td>13.24%</td></tr><tr><td>β-folding</td><td>23.47%</td><td>26.41%</td><td>24.60%</td><td>23.43%</td><td>16.18%</td></tr><tr><td>Random coils</td><td>71.84%</td><td>65.14%</td><td>72.58%</td><td>68.20%</td><td>70.59%</td></tr><tr><td>Exposed AA</td><td>47.65%</td><td>45.42%</td><td>43.95%</td><td>42.68%</td><td>16.18%</td></tr><tr><td>Buried AA</td><td>7.22%</td><td>7.04%</td><td>6.05%</td><td>11.30%</td><td>-</td></tr><tr><td>Intermediate AA</td><td>45.13%</td><td>47.54%</td><td>50%</td><td>46.03%</td><td>83.82%</td></tr></tbody></table>
### 3.3.2. Bioinformatic Analysis of the Polypeptide Sequence of rVP1

Protparam tool indicated that rVP1 expressed by the pUC19 expressing vector and VP1 protein of the Human Coxsackievirus B4E2 strain have many identical or close physicochemical properties compared to the other proteins. All the proteins are hydrophilic with a positive polarity, nevertheless, have some differences. Unlike the others, VP4 is unstable. PIs of VP3 and VP4 are less than 7 and significantly different from that of rVP1 (7.8), VP1 (8.81), and VP2 (8.83). rVP1 and VP1 have similar amino acid composition, atomic composition, and molecular weight. The analysis of CVB4 rVP1 and VP1-4 by applying different tools revealed that all the capsid proteins did not include potential signal peptide sequences (Figure 5B). However, all the proteins harbored phosphorylation sites. rVP1 includes 35 potential phosphorylation sites including 21 serine, 9 threonine, and 4 tyrosine; VP1 contains 23 potential phosphorylation sites, including 14 serine, 5 threonine, and 4 tyrosine (Figure 5C).

VP2-4 have a fewer number of phosphorylation sites. The subcellular location predicted by TargetP demonstrated that all the proteins could be at any extra-mitochondrial location. The prediction of the secondary structure and the solvent accessibility of AA residues are presented in Table 1. SOPMA and Predictprotein showed/revealed that random coils were the principal conformation of the capsid proteins, a higher ratio of β-folding, and lower ratios of β-turns. PredictProtein demonstrated that most of the AA residues in the VP4 protein were intermediate to solvent (83.82%). Percentages of exposed and buried AA residues in the other VP proteins were equivalent. rVP1 and VP1 have 47.65% and 45.42% of exposed AA, respectively. The protein structure homology-modeling was determined by employing the SWISS-MODEL server. The selection of templates depends on the similarity and coverage of AA sequences. Based on evaluation parameters, high standards models were applied as the tertiary structure of the VP proteins (Table 2). Model structures were observed by using RasMol 2.7.5 Software (Figure 6). The prediction of Linear B cell epitopes of rVP1 and VP1-VP4 capsid proteins was determined by the BepiPred server (Table 3). Among the potential linear B cell epitopes, VP4 included 2 epitopes with over 10 AA residues, respectively located at AA 5-19 and AA 29-54; VP3 contained 3 epitopes, respectively situated at AA 7-39, AA 55-67 and AA 139-150; VP2 included 5 epitopes respectively situated at AA 5-27, 38-60, 66-92, 130-165, 223-244; VP1 included 7 epitopes with over 10 AA residues, respectively located at 5-31, 33-50, 82-96, 127-138, 152-165, 199-209, 251-280; rVP1 6 epitopes 18-63, 83-93, 130-140, 155-168, 202-213, 258-273 (Table 3). The model structures were observed by using RasMol 2.7.5 Software (Figure 7).

<table><caption>**Table 1.** Prediction of secondary structures and solvent accessibility for the recombinant rVP1 and the CVB4E2 capsid proteins.</caption><thead><tr><th></th><th>rVP1</th><th>VP1</th><th>VP2</th><th>VP3</th><th>VP4</th></tr></thead><tbody><tr><td></td><td></td><td></td><td colspan="3">SOPMA</td></tr><tr><td>α-helix</td><td>23.47%</td><td>20.07%</td><td>17.34%</td><td>17.15%</td><td>17.56%</td></tr><tr><td>β-folding</td><td>25.63%</td><td>23.24%</td><td>26.61%</td><td>29.29%</td><td>22.06%</td></tr><tr><td>β-turns</td><td>3.25%</td><td>2.11%</td><td>7.66%</td><td>5.44%</td><td>4.41%</td></tr><tr><td>Random coils</td><td>47.65%</td><td>54.58%</td><td>48.39%</td><td>48.12%</td><td>55.88%</td></tr><tr><td></td><td></td><td></td><td colspan="3">Predictprotein</td></tr><tr><td>α-helix</td><td>4.69%</td><td>8.45%</td><td>2.82%</td><td>8.37%</td><td>13.24%</td></tr><tr><td>β-folding</td><td>23.47%</td><td>26.41%</td><td>24.60%</td><td>23.43%</td><td>16.18%</td></tr><tr><td>Random coils</td><td>71.84%</td><td>65.14%</td><td>72.58%</td><td>68.20%</td><td>70.59%</td></tr><tr><td>Exposed AA</td><td>47.65%</td><td>45.42%</td><td>43.95%</td><td>42.68%</td><td>16.18%</td></tr><tr><td>Buried AA</td><td>7.22%</td><td>7.04%</td><td>6.05%</td><td>11.30%</td><td>-</td></tr><tr><td>Intermediate AA</td><td>45.13%</td><td>47.54%</td><td>50%</td><td>46.03%</td><td>83.82%</td></tr></tbody></table>
<table><caption>**Table 2.** Prediction of the sequences homologous to the recombinant rVP1 and the CVB4E2 capsid proteins.</caption><thead><tr><th></th><th>Template</th><th>GMQE</th><th>QMEAN</th></tr></thead><tbody><tr><td>rVP1</td><td>1h8t.1A</td><td>0.85</td><td>-3.29</td></tr><tr><td>VP1</td><td>4GB3</td><td>0.78</td><td>-3.06</td></tr><tr><td>VP2</td><td>1mqt.1.B</td><td>0.57</td><td>-6.83</td></tr><tr><td>VP3</td><td>1oop.1.C</td><td>0.93</td><td>-1.90</td></tr><tr><td>VP4</td><td>1pov.1.A</td><td>0.57</td><td>-5.03</td></tr></tbody></table>

<table><caption>**Table 3.** Potential B cell epitope sites on the recombinant rVP1 and CVB4E2 capsid proteins.</caption><thead><tr><th>Capsid</th><th>N°</th><th>Start</th><th>End</th><th>Peptide</th><th>Length</th></tr></thead><tbody><tr><td rowspan="7">rVP1</td><td>1</td><td>7</td><td>7</td><td>E</td><td>1</td></tr><tr><td>2</td><td>18</td><td>63</td><td>AATLCPGPSNSEQNPALTAVESGHSSQVDPSDSMQTRHVHNYHSRS</td><td>46</td></tr><tr><td>3</td><td>83</td><td>93</td><td>SSAESNNLKRY</td><td>11</td></tr><tr><td>4</td><td>30</td><td>40</td><td>QEMSSASNSDV</td><td>11</td></tr><tr><td>5</td><td>155</td><td>168</td><td>PVPTSVNDYVWQTS</td><td>13</td></tr><tr><td>6</td><td>202</td><td>213</td><td>WSNFSRDGIYGY</td><td>12</td></tr><tr><td>7</td><td>230</td><td>234</td><td>SSPGG</td><td>5</td></tr><tr><td>8</td><td>258</td><td>273</td><td>LCQYKKAKNGNFDVEA</td><td>16</td></tr><tr><td rowspan="9">VP1</td><td>1</td><td>5</td><td>31</td><td>ESVERAMGRVADTIARGPSNSEQIPAL</td><td>27</td></tr><tr><td>2</td><td>33</td><td>50</td><td>AVETGHTSQVDPSDTMQT</td><td>18</td></tr><tr><td>3</td><td>53</td><td>53</td><td>V</td><td>1</td></tr><tr><td>4</td><td>55</td><td>61</td><td>NYHSRSE</td><td>7</td></tr><tr><td>5</td><td>82</td><td>96</td><td>AESNNLKRYAEWVIN</td><td>15</td></tr><tr><td>6</td><td>127</td><td>138</td><td>QEMSTATNSVVP</td><td>12</td></tr><tr><td>7</td><td>152</td><td>165</td><td>PVPTSVNDYVWQTS</td><td>14</td></tr><tr><td>8</td><td>176</td><td>176</td><td>N</td><td>1</td></tr><tr><td>9</td><td>199</td><td>209</td><td>WSNFSRDGIYG</td><td>11</td></tr><tr><td>10</td><td>227</td><td>231</td><td>SSPGGL</td><td>6</td></tr><tr><td>11</td><td>251</td><td>280</td><td>PPRRLCQYKKAKNVNFDVEAVTTERASLVTT</td><td>31</td></tr><tr><td rowspan="7">VP2</td><td>1</td><td>5</td><td>27</td><td>EECGYSDRVRITSITLGNSTITTQE</td><td>13</td></tr><tr><td>2</td><td>38</td><td>60</td><td>WPDYLSDEEATAEDQPTQPDVAT</td><td>23</td></tr><tr><td>3</td><td>66</td><td>92</td><td>LKSVKWRCSQRGGGSSQMRCQRWVSL</td><td>27</td></tr><tr><td>4</td><td>130</td><td>165</td><td>WGALTRKMRPHMVIYAEGKQQNNLNKMQSQVRLLCK</td><td>36</td></tr><tr><td>5</td><td>213</td><td>215</td><td>SHY</td><td>3</td></tr><tr><td>6</td><td>223</td><td>244</td><td>RERPPTFPSLQLLQALNTMACD</td><td>22</td></tr><tr><td>1</td><td>7</td><td>39</td><td>TPGSTQFLTSDDFQSPSAMPQFDVNPEMNIPGR</td><td>33</td></tr><tr><td rowspan="5">VP3</td><td>2</td><td>55</td><td>67</td><td>INNLQANLKTMEA</td><td>13</td></tr><tr><td>3</td><td>72</td><td>80</td><td>VRSTDEMGQ</td><td>9</td></tr><tr><td>4</td><td>91</td><td>94</td><td>SSVL</td><td>4</td></tr><tr><td>5</td><td>139</td><td>150</td><td>GAGAPDSRKNAM</td><td>12</td></tr><tr><td>6</td><td>202</td><td>206</td><td>PAEAQ</td><td>5</td></tr><tr><td>7</td><td>228</td><td>236</td><td>DTQFIKQDT</td><td>9</td></tr><tr><td rowspan="3">VP4</td><td>1</td><td>5</td><td>19</td><td>STQKTGAHETSLSAT</td><td>15</td></tr><tr><td>2</td><td>29</td><td>54</td><td>INYYKDAASNSANRQDFTQDPSKFTE</td><td>26</td></tr><tr><td>3</td><td>61</td><td>62</td><td>IK</td><td>2</td></tr></tbody></table>
<table><caption>**Table 2.** Prediction of the sequences homologous to the recombinant rVP1 and the CVB4E2 capsid proteins.</caption><thead><tr><th></th><th>Template</th><th>GMQE</th><th>QMEAN</th></tr></thead><tbody><tr><td>rVP1</td><td>1h8t.1A</td><td>0.85</td><td>-3.29</td></tr><tr><td>VP1</td><td>4GB3</td><td>0.78</td><td>-3.06</td></tr><tr><td>VP2</td><td>1mqt.1.B</td><td>0.57</td><td>-6.83</td></tr><tr><td>VP3</td><td>1oop.1.C</td><td>0.93</td><td>-1.90</td></tr><tr><td>VP4</td><td>1pov.1.A</td><td>0.57</td><td>-5.03</td></tr></tbody></table>

<table><caption>**Table 3.** Potential B cell epitope sites on the recombinant rVP1 and CVB4E2 capsid proteins.</caption><thead><tr><th>Capsid</th><th>N°</th><th>Start</th><th>End</th><th>Peptide</th><th>Length</th></tr></thead><tbody><tr><td rowspan="7">rVP1</td><td>1</td><td>7</td><td>7</td><td>E</td><td>1</td></tr><tr><td>2</td><td>18</td><td>63</td><td>AATLCPGPSNSEQNPALTAVESGHSSQVDPSDSMQTRHVHNYHSRS</td><td>46</td></tr><tr><td>3</td><td>83</td><td>93</td><td>SSAESNNLKRY</td><td>11</td></tr><tr><td>4</td><td>30</td><td>40</td><td>QEMSSASNSDV</td><td>11</td></tr><tr><td>5</td><td>155</td><td>168</td><td>PVPTSVNDYVWQTS</td><td>13</td></tr><tr><td>6</td><td>202</td><td>213</td><td>WSNFSRDGIYGY</td><td>12</td></tr><tr><td>7</td><td>230</td><td>234</td><td>SSPGG</td><td>5</td></tr><tr><td>8</td><td>258</td><td>273</td><td>LCQYKKAKNGNFDEA</td><td>16</td></tr><tr><td rowspan="10">VP1</td><td>1</td><td>5</td><td>31</td><td>ESVERAMGRVADTIARGPSNSEQIPAL</td><td>27</td></tr><tr><td>2</td><td>33</td><td>50</td><td>AVETGHTSQVDPSDTMQT</td><td>18</td></tr><tr><td>3</td><td>53</td><td>53</td><td>V</td><td>1</td></tr><tr><td>4</td><td>55</td><td>61</td><td>NYHSRSE</td><td>7</td></tr><tr><td>5</td><td>82</td><td>96</td><td>AESNNLKRYAEWVIN</td><td>15</td></tr><tr><td>6</td><td>127</td><td>138</td><td>QEMSTATNSVVP</td><td>12</td></tr><tr><td>7</td><td>152</td><td>165</td><td>PVPTSVNDYVWQTS</td><td>14</td></tr><tr><td>8</td><td>176</td><td>176</td><td>N</td><td>1</td></tr><tr><td>9</td><td>199</td><td>209</td><td>WSNFSRDGIYG</td><td>11</td></tr><tr><td>10</td><td>227</td><td>231</td><td>SSPGGL</td><td>6</td></tr><tr><td>11</td><td>251</td><td>280</td><td>PPRRLCQYKKAKNVNFDEAVTTERASLVTT</td><td>31</td></tr><tr><td rowspan="7">VP2</td><td>1</td><td>5</td><td>27</td><td>EECGYSDRVRITSITLGNSTITTQE</td><td>13</td></tr><tr><td>2</td><td>38</td><td>60</td><td>WPDYLSDEEATAEDQPTQPDVAT</td><td>23</td></tr><tr><td>3</td><td>66</td><td>92</td><td>LKSVKWRCSQRGGGSSQMRCQRWVSL</td><td>27</td></tr><tr><td>4</td><td>130</td><td>165</td><td>WGALTRKMRPHMVIYAEGKQQNNLNKMQSQVRLLCK</td><td>36</td></tr><tr><td>5</td><td>213</td><td>215</td><td>SHY</td><td>3</td></tr><tr><td>6</td><td>223</td><td>244</td><td>RERPPTFPSLQLLQALNTMACD</td><td>22</td></tr><tr><td>7</td><td>7</td><td>39</td><td>TPGSTQFLTSDDFQSPSAMPQFDVNPEMNIPGR</td><td>33</td></tr><tr><td rowspan="5">VP3</td><td>2</td><td>55</td><td>67</td><td>INNLQANLKTMEA</td><td>13</td></tr><tr><td>3</td><td>72</td><td>80</td><td>VRSTDEMGQ</td><td>9</td></tr><tr><td>4</td><td>91</td><td>94</td><td>SSVL</td><td>4</td></tr><tr><td>5</td><td>139</td><td>150</td><td>GAGAPDSRKNAM</td><td>12</td></tr><tr><td>6</td><td>202</td><td>206</td><td>PAEAQ</td><td>5</td></tr><tr><td>7</td><td>228</td><td>236</td><td>DTQFIKQDT</td><td>9</td></tr><tr><td rowspan="3">VP4</td><td>1</td><td>5</td><td>19</td><td>STQKTGAHETSLSAT</td><td>15</td></tr><tr><td>2</td><td>29</td><td>54</td><td>INYYKDAASNSANRQDFTQDPSKFTE</td><td>26</td></tr><tr><td>3</td><td>61</td><td>62</td><td>IK</td><td>2</td></tr></tbody></table>
<figure><img src="image_10.png" /><figcaption>**Figure 6.** Prediction of the tertiary structure of rVP1 and CVB4E2 capsid proteins. The colored ribbons show the tertiary structure of the viral proteins. The red ones show the α helices, the yellows are the β-folding, and the white ones correspond to other residues.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Figure 7.** The possible B-epitope regions of the proteins and their presentation in tertiary structures. RasMol generated the tertiary ribbon and molecular surface representations. B cell epitopes are shown in red, blue, green, gold, purple, orange, and light green based on the order: AA 18–63, 83–93, 130–140, 155–168, 202–213, 258–273 of rVP1; AA 5–31, 33–50, 82–96, 127–138, 152–165, 199–209, 251–280 of VP1; AA 5–27, 38–60, 66–92, 130–165, 223–244 of VP2; AA 7–39, AA 55–67 and AA 139–150 of VP3; and AA 5–19 and AA 29–54 of VP4.</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>**Figure 6.** Prediction of the tertiary structure of rVP1 and CVB4E2 capsid proteins. The colored ribbons show the tertiary structure of the viral proteins. The red ones show the α helices, the yellows are the β-folding, and the white ones correspond to other residues.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Figure 7.** The possible B-epitope regions of the proteins and their presentation in tertiary structures. RasMol generated the tertiary ribbon and molecular surface representations. B cell epitopes are shown in red, blue, green, gold, purple, orange, and light green based on the order: AA 18–63, 83–93, 130–140, 155–168, 202–213, 258–273 of rVP1; AA 5–31, 33–50, 82–96, 127–138, 152–165, 199–209, 251–280 of VP1; AA 5–27, 38–60, 66–92, 130–165, 223–244 of VP2; AA 7–39, AA 55–67 and AA 139–150 of VP3; and AA 5–19 and AA 29–54 of VP4.</figcaption></figure>
**4. Discussion**

The aim of this study was first to clone, express, purify and characterize the VP1 gene of the CVB4E2 strain in the pUC19 plasmid and then to compare it with the structural capsid proteins of the same strain by bioinformatic tools. PCR colony amplification followed by restriction digestion analysis and sequencing process affirmed the success of the cloning process. After preparing the bacterial lysate, the expression of CVB4E2 recombinant protein was confirmed by separating different proteins in the bacterial cell lysate using the SDS-PAGE technique. The molecular mass of recombinant VP1 subunit protein, obtained by separating different bacterial proteins through a 15% SDS-PAGE technique was estimated to be approximately 30 kD, which is very close to the expected molecular weight. Bands intensity in the PAGE gel is usually proportional to the amount of protein [27]. The exact identity of the expressed protein was established by subjecting the recombinant *E. coli* DH5α cells expressing the target protein to Western Blot transfer analysis. The recombinant rVP1 subunit protein was identified by using the anti-enterovirus monoclonal antibody mAb 5-D8/1 conjugated to the anti-mouse peroxydase. This mAb picked up a band of approximately 30 kD in the elution fractions from the bacteria cell lysates. However, no such signal was produced in the non-transformed cell lysates. Purification of the recombinant rVP1 subunit protein from *E. coli* DH5α was carried out using Ni-NTA affinity chromatography. Upon purification of the recombinant rVP1 subunit protein, the purified fraction was obtained in the elution fractions 4, 5 and 6 as confirmed by 15% SDS-PAGE analysis of different fractions.

The encoded protein of rVP1 is a member of the rhv_like superfamily with a pfam00073 conserved domain. rVP1 is a Coxsackievirus capsid protein of the *Picornaviridae* family. According to the negative value of GRAVY, all the viral proteins are hydrophilic. Furthermore, they have positive polarity. However, the theoretical isoelectric point, molecular extinction coefficient, and instability index were relatively different. The theoretical isoelectric point is an important parameter to identify and purify proteins. The values of this parameter revealed that rVP1, VP1, and VP2 are basic. On the contrary, the other viral proteins are acidic. The molecular extinction coefficient is known as the absorbance of a protein at 280 nm. It depends on the number of aromatic residues, Tryptophan (Trp), tyrosine (Tyr), and cystine (Cys). The molar extinction coefficient cannot be determined by Protparam when the AA sequence does not contain aromatic residues. That is the case with the VP4 protein. If the instability index is higher than 40, it indicates that the protein is unstable, and if it is lower than 40, the protein is stable [27,28]. Unlike the other proteins, VP4 is stable. Although rVP1 and VP1 had similar molecular weights, the other proteins have significantly different molecular weights. Phosphorylation site prediction indicated that rVP1 and VP1 have more phosphorylation sites than the other viral proteins. Consequently, they are all able to affect the signal transduction of host cells and can be involved in virus virulence. The secondary structure prediction of rVP1 and VP of CV-B4 E2 principally includes random coils, also known as loops. Furthermore, unlike VP4, all the other proteins had a higher percentage of exposed AA. For several biomedical applications such as vaccine development, the determination of antigenic determinant is crucial. B-cell epitopes are hydrophilic, and have more random coils and fewer secondary structures. Furthermore, epitopes have a significantly higher ratio of exposed AA and a lesser percentage of buried AA [29]. Prediction of B-cell epitopes by applying Bepipred revealed that all the compared proteins harbored some antigenic determinants (Table 3, Figure 7). According to this, rVP1 and VP1 have the highest number of significant epitopes with similar locations. However, VP4 has a lesser number of neutralizing epitopes.

The major immunogenic protein VP1 of the CVB4E2 strain was successfully expressed in *E. coli* DH5α cells using expression vector pUC19. As a widely used method, the construction of a pUC19 vector for protein expression and purification in *E. coli* is fast, inexpensive and scalable. Ni-NTA affinity chromatography proved a rapid and excellent technique to purify the rVP1 subunit protein expressed in *E. coli*. Bioinformatic characterization proved a well-studied configuration of the expressed protein, which preserves most characteristics of
**4. Discussion**

The aim of this study was first to clone, express, purify and characterize the VP1 gene of the CVB4E2 strain in the pUC19 plasmid and then to compare it with the structural capsid proteins of the same strain by bioinformatic tools. PCR colony amplification followed by restriction digestion analysis and sequencing process affirmed the success of the cloning process. After preparing the bacterial lysate, the expression of CVB4E2 recombinant protein was confirmed by separating different proteins in the bacterial cell lysate using the SDS-PAGE technique. The molecular mass of recombinant VP1 subunit protein, obtained by separating different bacterial proteins through a 15% SDS-PAGE technique was estimated to be approximately 30 kD, which is very close to the expected molecular weight. Bands intensity in the PAGE gel is usually proportional to the amount of protein [27]. The exact identity of the expressed protein was established by subjecting the recombinant *E. coli* DH5α cells expressing the target protein to Western Blot transfer analysis. The recombinant rVP1 subunit protein was identified by using the anti-enterovirus monoclonal antibody mAb 5-D8/1 conjugated to the anti-mouse peroxydase. This mAb picked up a band of approximately 30 kD in the elution fractions from the bacteria cell lysates. However, no such signal was produced in the non-transformed cell lysates. Purification of the recombinant rVP1 subunit protein from *E. coli* DH5α was carried out using Ni-NTA affinity chromatography. Upon purification of the recombinant rVP1 subunit protein, the purified fraction was obtained in the elution fractions 4, 5 and 6 as confirmed by 15% SDS-PAGE analysis of different fractions.

The encoded protein of rVP1 is a member of the rhv_like superfamily with a pfam00073 conserved domain. rVP1 is a Coxsackievirus capsid protein of the *Picornaviridae* family. According to the negative value of GRAVY, all the viral proteins are hydrophilic. Furthermore, they have positive polarity. However, the theoretical isoelectric point, molecular extinction coefficient, and instability index were relatively different. The theoretical isoelectric point is an important parameter to identify and purify proteins. The values of this parameter revealed that rVP1, VP1, and VP2 are basic. On the contrary, the other viral proteins are acidic. The molecular extinction coefficient is known as the absorbance of a protein at 280 nm. It depends on the number of aromatic residues, Tryptophan (Trp), tyrosine (Tyr), and cystine (Cys). The molar extinction coefficient cannot be determined by Protparam when the AA sequence does not contain aromatic residues. That is the case with the VP4 protein. If the instability index is higher than 40, it indicates that the protein is unstable, and if it is lower than 40, the protein is stable [27,28]. Unlike the other proteins, VP4 is stable. Although rVP1 and VP1 had similar molecular weights, the other proteins have significantly different molecular weights. Phosphorylation site prediction indicated that rVP1 and VP1 have more phosphorylation sites than the other viral proteins. Consequently, they are all able to affect the signal transduction of host cells and can be involved in virus virulence. The secondary structure prediction of rVP1 and VP of CV-B4 E2 principally includes random coils, also known as loops. Furthermore, unlike VP4, all the other proteins had a higher percentage of exposed AA. For several biomedical applications such as vaccine development, the determination of antigenic determinant is crucial. B-cell epitopes are hydrophilic, and have more random coils and fewer secondary structures. Furthermore, epitopes have a significantly higher ratio of exposed AA and a lesser percentage of buried AA [29]. Prediction of B-cell epitopes by applying Bepipred revealed that all the compared proteins harbored some antigenic determinants (Table 3, Figure 7). According to this, rVP1 and VP1 have the highest number of significant epitopes with similar locations. However, VP4 has a lesser number of neutralizing epitopes.

The major immunogenic protein VP1 of the CVB4E2 strain was successfully expressed in *E. coli* DH5α cells using expression vector pUC19. As a widely used method, the construction of a pUC19 vector for protein expression and purification in *E. coli* is fast, inexpensive and scalable. Ni-NTA affinity chromatography proved a rapid and excellent technique to purify the rVP1 subunit protein expressed in *E. coli*. Bioinformatic characterization proved a well-studied configuration of the expressed protein, which preserves most characteristics of
the wild-type CVB4 virus proteins. Results of this preliminary study support the potential strategy to produce related subunit vaccines, and immunodiagnostic reagents based on the recombinant viral protein.

**Author Contributions:** Conceptualization, I.H.H., D.H. and M.B.M.; methodology, I.H.H., I.A., A.A., R.S. and M.A.; writing, I.H.H.; funding acquisition, M.B.M. and J.G. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was supported by the Deanship of Scientific Research, Vice Presidency for Graduate Studies and Scientific Research, King Faisal University, Saudi Arabia [GRANT Nb. 2400] and the Tunisian Ministry of Higher Education and Scientific Research, Tunisia [GRANT UR17/ES30].

**Acknowledgments:** This work was supported by the Deanship of Scientific Research, Vice Presidency for Graduate Studies and Scientific Research, King Faisal University, Saudi Arabia [GRANT Nb. 2400] and the Tunisian Ministry of Higher Education and Scientific Research, Tunisia [GRANT UR17/ES30].

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Tapparel, C.; Siegrist, F.; Petty, T.J.; Kaiser, L. Picornavirus and enterovirus diversity with associated human diseases. *Infect. Genet. Evol.* **2013**, *14*, 282–293. [CrossRef]

2. Zell, R.; Delwart, E.; Gorbalenya, A.E.; Hovi, T.; King, A.M.Q.; Knowles, N.J.; Lindberg, A.M.; Pallansch, M.A.; Palmenberg, A.C.; Reuter, G.; et al. ICTV virus taxonomy profile: Picornaviridae. *J. Gen. Virol.* **2017**, *98*, 2421–2422. [CrossRef]

3. Solomon, T.; Lewthwaite, P.; DaCardosa, M.J.; McMinn, P.; Ooi, M.H. Virology, epidemiology, pathogenesis, and control of enterovirus 71. *Lancet Infect. Dis.* **2010**, *10*, 778–790. [CrossRef] [PubMed]

4. Hassine, I.H.; Gharbi, J.; Hamrita, B.; Almalki, M.A.; Rodríguez, J.F.; Ben M'hadheb, M. Characterization of Coxsackievirus B4 virus-like particles VLP produced by the recombinant baculovirus-insect cell system expressing the major capsid protein. *Mol. Biol. Rep.* **2020**, *47*, 2835–2843. [CrossRef]

5. Ben M'hadheb, M.; Souii, A.; Harrabi, M.; Jrad-Battikh, N.; Gharbi, J. In Vitro-reduced translation efficiency of coxsackievirus B3 Sabin3-like strain is correlated to impaired binding of cellular initiation factors to viral IRES RNA. *Curr. Microbiol.* **2015**, *70*, 756–761. [CrossRef]

6. Gharbi, J.; Hiar, R.; Ben M'hadheb, M.; Jaïdane, H.; Bouslama, L.; N'saïbia, S. Aouni, M. Nucleotide sequences of IRES domains IV and V of natural ECHO virus type 11 isolates with different replicative capacity phenotypes. *Virus Genes* **2006**, *32*, 269–276. [CrossRef]

7. Acharya, R.; Fry, E.; Stuart, D.; Fox, G.; Rowlands, D.; Brown, F. The three-dimensional structure of foot-and-mouth disease virus at 2.9 Å resolution. *Nature* **1989**, *337*, 709–716. [CrossRef] [PubMed]

8. Oberste, M.S.; Maher, K.; Kilpatrick, D.R.; Flemister, M.R.; Brown, B.A.; Pallansch, M.A. Typing of human enteroviruses by partial sequencing of VP1. *J. Clin. Microbiol.* **1999**, *37*, 1288–1293. [CrossRef] [PubMed]

9. Yoon, J.W.; Austin, M.; Onodera, T.; Notkins, A.L. Isolation of a virus from the pancreas of a child with diabetic ketoacidosis. *N. Engl. J. Med.* **1979**, *300*, 1173–1179. [CrossRef]

10. Hoover, D.M.; Lubkowski, J. DNAWorks: An automated method for designing oligonucleotides for PCR-based gene synthesis. *Nucleic. Acids. Res.* **2002**, *30*, 1–7. [CrossRef]

11. Niu, X.N.; Wei, Z.Q.; Zou, H.F.; Xie, G.G.; Wu, F.; Li, K.J.; Jiang, W.; Tang, J.L.; He, Y.Q. Complete sequence and detailed analysis of the first indigenous plasmid from Xanthomonas oryzae pv. Oryzicola. *BMC Microbiol.* **2015**, *15*, 1–15. [CrossRef] [PubMed]

12. Messing, J.; Vieira, J. A new pair of M13 vectors for selecting either DNA strand of double-digest restriction fragments. *Gene* **1982**, *19*, 269–276. [CrossRef]

13. Messing, J. New M13 vectors for cloning. *Methodes Enzymol.* **1983**, *101*, 20–78.

14. Yanisch-perron, C.; Vieira, J.; Messing, J. Improved M13 phage cloning vectors and host strains: Nucleotide sequences of the M13mp18 and pUC19 vectors. *Gene* **1985**, *33*, 103–119. [CrossRef]

15. Reed, L.J.; Muench, H. A simple method of estimating fifty percent endpoints. *Am. J. Epidemiol.* **1983**, *27*, 493–497. [CrossRef]

16. Chomczynski, P.; Sacchi, N. Single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction. *Anal. Biochem.* **1987**, *162*, 156–159. [CrossRef] [PubMed]

17. Chung, C.T.; Niemela, S.L.; Miller, R.H. One-step preparation of competent *Escherichia coli*: Transformation and storage of bacterial cells in the same solution. *Proc. Natl. Acad. Sci. USA* **1989**, *86*, 2172–2175. [CrossRef] [PubMed]

18. Sanger, F.; Nicklen, S.; Coulson, A.R. DNA sequencing with chain-terminating inhibitors. *Proc. Natl. Acad. Sci. USA* **1977**, *74*, 5463–5467. [CrossRef]

19. Marchler-Bauer, A.; Bryant, S.H. CD-Search: Protein domain annotations on the fly. *Nucleic Acids Res.* **2004**, *32*, W327–W331. [CrossRef]
the wild-type CVB4 virus proteins. Results of this preliminary study support the potential strategy to produce related subunit vaccines, and immunodiagnostic reagents based on the recombinant viral protein.

**Author Contributions:** Conceptualization, I.H.H., D.H. and M.B.M.; methodology, I.H.H., I.A., A.A., R.S. and M.A.; writing, I.H.H.; funding acquisition, M.B.M. and J.G. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was supported by the Deanship of Scientific Research, Vice Presidency for Graduate Studies and Scientific Research, King Faisal University, Saudi Arabia [GRANT Nb. 2400] and the Tunisian Ministry of Higher Education and Scientific Research, Tunisia [GRANT UR17/ES30].

**Acknowledgments:** This work was supported by the Deanship of Scientific Research, Vice Presidency for Graduate Studies and Scientific Research, King Faisal University, Saudi Arabia [GRANT Nb. 2400] and the Tunisian Ministry of Higher Education and Scientific Research, Tunisia [GRANT UR17/ES30].

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Tapparel, C.; Siegrist, F.; Petty, T.J.; Kaiser, L. Picornavirus and enterovirus diversity with associated human diseases. *Infect. Genet. Evol.* **2013**, *14*, 282–293. [CrossRef]

2. Zell, R.; Delwart, E.; Gorbalenya, A.E.; Hovi, T.; King, A.M.Q.; Knowles, N.J.; Lindberg, A.M.; Pallansch, M.A.; Palmenberg, A.C.; Reuter, G.; et al. ICTV virus taxonomy profile: Picornaviridae. *J. Gen. Virol.* **2017**, *98*, 2421–2422. [CrossRef]

3. Solomon, T.; Lewthwaite, P.; DaCardosa, M.J.; McMinn, P.; Ooi, M.H. Virology, epidemiology, pathogenesis, and control of enterovirus 71. *Lancet Infect. Dis.* **2010**, *10*, 778–790. [CrossRef] [PubMed]

4. Hassine, I.H.; Gharbi, J.; Hamrita, B.; Almalki, M.A.; Rodríguez, J.F.; Ben M'hadheb, M. Characterization of Coxsackievirus B4 virus-like particles VLP produced by the recombinant baculovirus-insect cell system expressing the major capsid protein. *Mol. Biol. Rep.* **2020**, *47*, 2835–2843. [CrossRef]

5. Ben M'hadheb, M.; Souii, A.; Harrabi, M.; Jrad-Battikh, N.; Gharbi, J. In Vitro-reduced translation efficiency of coxsackievirus B3 Sabin3-like strain is correlated to impaired binding of cellular initiation factors to viral IRES RNA. *Curr. Microbiol.* **2015**, *70*, 756–761. [CrossRef]

6. Gharbi, J.; Hiar, R.; Ben M'hadheb, M.; Jaïdane, H.; Bouslama, L.; N'saïbia, S. Aouni, M. Nucleotide sequences of IRES domains IV and V of natural ECHO virus type 11 isolates with different replicative capacity phenotypes. *Virus Genes* **2006**, *32*, 269–276. [CrossRef]

7. Acharya, R.; Fry, E.; Stuart, D.; Fox, G.; Rowlands, D.; Brown, F. The three-dimensional structure of foot-and-mouth disease virus at 2.9 Å resolution. *Nature* **1989**, *337*, 709–716. [CrossRef] [PubMed]

8. Oberste, M.S.; Maher, K.; Kilpatrick, D.R.; Flemister, M.R.; Brown, B.A.; Pallansch, M.A. Typing of human enteroviruses by partial sequencing of VP1. *J. Clin. Microbiol.* **1999**, *37*, 1288–1293. [CrossRef] [PubMed]

9. Yoon, J.W.; Austin, M.; Onodera, T.; Notkins, A.L. Isolation of a virus from the pancreas of a child with diabetic ketoacidosis. *N. Engl. J. Med.* **1979**, *300*, 1173–1179. [CrossRef]

10. Hoover, D.M.; Lubkowski, J. DNAWorks: An automated method for designing oligonucleotides for PCR-based gene synthesis. *Nucleic. Acids. Res.* **2002**, *30*, 1–7. [CrossRef]

11. Niu, X.N.; Wei, Z.Q.; Zou, H.F.; Xie, G.G.; Wu, F.; Li, K.J.; Jiang, W.; Tang, J.L.; He, Y.Q. Complete sequence and detailed analysis of the first indigenous plasmid from Xanthomonas oryzae pv. Oryzicola. *BMC Microbiol.* **2015**, *15*, 1–15. [CrossRef] [PubMed]

12. Messing, J.; Vieira, J. A new pair of M13 vectors for selecting either DNA strand of double-digest restriction fragments. *Gene* **1982**, *19*, 269–276. [CrossRef]

13. Messing, J. New M13 vectors for cloning. *Methodes Enzymol.* **1983**, *101*, 20–78.

14. Yanisch-perron, C.; Vieira, J.; Messing, J. Improved M13 phage cloning vectors and host strains: Nucleotide sequences of the M13mp18 and pUC19 vectors. *Gene* **1985**, *33*, 103–119. [CrossRef]

15. Reed, L.J.; Muench, H. A simple method of estimating fifty percent endpoints. *Am. J. Epidemiol.* **1983**, *27*, 493–497. [CrossRef]

16. Chomczynski, P.; Sacchi, N. Single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction. *Anal. Biochem.* **1987**, *162*, 156–159. [CrossRef] [PubMed]

17. Chung, C.T.; Niemela, S.L.; Miller, R.H. One-step preparation of competent *Escherichia coli*: Transformation and storage of bacterial cells in the same solution. *Proc. Natl. Acad. Sci. USA* **1989**, *86*, 2172–2175. [CrossRef] [PubMed]

18. Sanger, F.; Nicklen, S.; Coulson, A.R. DNA sequencing with chain-terminating inhibitors. *Proc. Natl. Acad. Sci. USA* **1977**, *74*, 5463–5467. [CrossRef]

19. Marchler-Bauer, A.; Bryant, S.H. CD-Search: Protein domain annotations on the fly. *Nucleic Acids Res.* **2004**, *32*, W327–W331. [CrossRef]
20. Gasteiger, E.; Hoogland, C.; Gattiker, A.; Duvaud, S.; Wilkins, M.R.; Appel, R.D.; Bairoch, A. Protein Identification and Analysis Tools on the ExPASy server. In *The proteomics Protocols Handbook*; John, M.W., Ed.; Humana Press: Totowa, NJ, USA, 2005; pp. 571–607.

21. Petersen, T.N.; Brunak, S.; Von Heijne, G.; Nielsen, H. Signal P 4.0: Discriminating signal peptides from transmembrane regions. *Nat. Methods* **2011**, *8*, 785–786. [CrossRef]

22. Blom, N.S.; Gammeltoft, S.; Brunak, S. Sequence and structure-based prediction of eukaryotic protein phosphorylation sites. *J. Mol. Biol.* **1999**, *294*, 1351–1362. [CrossRef] [PubMed]

23. Blom, N.; Sicheritz-Pontén, T.; Gupta, R.; Gammeltoft, S.; Brunak, S. Prediction of post-translational glycosylation and phosphorylation of proteins from the amino acid sequence. *Proteomics* **2004**, *4*, 1633–1649. [CrossRef]

24. Letunic, I.; Bork, P. 20 years of the SMART protein domain annotation resource. *Nucleic Acids Res.* **2018**, *46*, D493–D496. [CrossRef] [PubMed]

25. Letunic, I.; Doerks, T.; Bork, P. SMART: Recent updates, new developments and status in 2015. *Nucleic Acids Res.* **2015**, *43*, D257–D260. [CrossRef]

26. Selbig, J.; Mevissen, T.; Lengauer, T. Decision tree-based formation of consensus protein secondary structure prediction. *Bioinformatics* **1999**, *15*, 1039–1046. [CrossRef]

27. Idicula-thomas, S.; Balaji, P.V. Understanding the relationship between the primary structure of proteins and its propensity to be soluble on overexpression in Escherichia coli. *Protein Sci.* **2005**, *14*, 582–592. [CrossRef] [PubMed]

28. Bradford, M.M. 1976. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. *Anal. Biochem.* **1976**, *72*, 248–254. [CrossRef]

29. Dennis, M.K.; Field, A.S.; Burai, R.; Ramesh, C.; Whitney, K.; Bologa, C.G.; Oprea, T.I.; Yamaguchi, Y.; Hayashi, S.; Sklar, L.A.; et al. Identification of a GPER/GPR30 antagonist with improved estrogen receptor counterselectivity. *J. Steroid. Biochem. Mol. Biol.* **2011**, *127*, 358–366. [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
20. Gasteiger, E.; Hoogland, C.; Gattiker, A.; Duvaud, S.; Wilkins, M.R.; Appel, R.D.; Bairoch, A. Protein Identification and Analysis Tools on the ExPASy server. In *The proteomics Protocols Handbook*; John, M.W., Ed.; Humana Press: Totowa, NJ, USA, 2005; pp. 571–607.

21. Petersen, T.N.; Brunak, S.; Von Heijne, G.; Nielsen, H. Signal P 4.0: Discriminating signal peptides from transmembrane regions. *Nat. Methods* **2011**, *8*, 785–786. [CrossRef]

22. Blom, N.S.; Gammeltoft, S.; Brunak, S. Sequence and structure-based prediction of eukaryotic protein phosphorylation sites. *J. Mol. Biol.* **1999**, *294*, 1351–1362. [CrossRef] [PubMed]

23. Blom, N.; Sicheritz-Pontén, T.; Gupta, R.; Gammeltoft, S.; Brunak, S. Prediction of post-translational glycosylation and phosphorylation of proteins from the amino acid sequence. *Proteomics* **2004**, *4*, 1633–1649. [CrossRef]

24. Letunic, I.; Bork, P. 20 years of the SMART protein domain annotation resource. *Nucleic Acids Res.* **2018**, *46*, D493–D496. [CrossRef] [PubMed]

25. Letunic, I.; Doerks, T.; Bork, P. SMART: Recent updates, new developments and status in 2015. *Nucleic Acids Res.* **2015**, *43*, D257–D260. [CrossRef]

26. Selbig, J.; Mevissen, T.; Lengauer, T. Decision tree-based formation of consensus protein secondary structure prediction. *Bioinformatics* **1999**, *15*, 1039–1046. [CrossRef]

27. Idicula-thomas, S.; Balaji, P.V. Understanding the relationship between the primary structure of proteins and its propensity to be soluble on overexpression in Escherichia coli. *Protein Sci.* **2005**, *14*, 582–592. [CrossRef] [PubMed]

28. Bradford, M.M. 1976. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. *Anal. Biochem.* **1976**, *72*, 248–254. [CrossRef]

29. Dennis, M.K.; Field, A.S.; Burai, R.; Ramesh, C.; Whitney, K.; Bologa, C.G.; Oprea, T.I.; Yamaguchi, Y.; Hayashi, S.; Sklar, L.A.; et al. Identification of a GPER/GPR30 antagonist with improved estrogen receptor counterselectivity. *J. Steroid. Biochem. Mol. Biol.* **2011**, *127*, 358–366. [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.