Next Article in Journal
Epstein–Barr Virus and Human Herpesvirus-6 Reactivation in Acute COVID-19 Patients
Next Article in Special Issue
Characterizing Disabilities in Preschool Aged Children with Congenital Zika Virus Syndrome with the ICF Model
Previous Article in Journal
Sacbrood Virus: A Growing Threat to Honeybees and Wild Pollinators
Previous Article in Special Issue
Zika Virus (ZIKV): A New Perspective on the Nanomechanical and Structural Properties
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Selection and Characterization of Single-Stranded DNA Aptamers of Diagnostic Potential against the Whole Zika Virus

by
Liliane Monteiro de Morais
1,
Thiago Santos Chaves
2,
Marco Alberto Medeiros
3,
Kaique Alves Brayner Pereira
4,
Patrícia Barbosa Jurgilas
4,
Sheila Maria Barbosa de Lima
1,
Sotiris Missailidis
3,*,† and
Ana Maria Bispo de Filippis
5,†
1
Laboratório de Tecnologia Virológica—LATEV, Instituto de Tecnologia em Imunobiológicos, Fundação Oswaldo Cruz-Fiocruz, Rio de Janeiro 22725, RJ, Brazil
2
Laboratório de Tecnologia Recombinante—LATER, Instituto de Tecnologia em Imunobiológicos, Fundação Oswaldo Cruz-Fiocruz, Rio de Janeiro 22725, RJ, Brazil
3
Vice-Diretoria de Desenvolvimento Tecnológico—VDTEC, Instituto de Tecnologia em Imunobiológicos, Fundação Oswaldo Cruz-Fiocruz, Rio de Janeiro 22725, RJ, Brazil
4
Laboratório de Macromoléculas—LAMAM, Instituto de Tecnologia em Imunobiológicos, Fundação Oswaldo Cruz-Fiocruz, Rio de Janeiro 22725, RJ, Brazil
5
Laboratório de Flavivirus, Instituto Oswaldo Cruz, Fundação Oswaldo Cruz-Fiocruz, Rio de Janeiro 22725, RJ, Brazil
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Viruses 2022, 14(9), 1867; https://doi.org/10.3390/v14091867
Submission received: 25 July 2022 / Revised: 17 August 2022 / Accepted: 22 August 2022 / Published: 25 August 2022
(This article belongs to the Special Issue Zika Viruses: State-of-the-Art Research in Brazil)

Abstract

:
Zika virus became a major public health problem in early 2015, when cases of Guillain–Barré syndrome and microcephaly were associated with viral infection. Currently, ZIKV is endemic in all tropical areas of the world, and the chance for future Zika epidemics remains very real and accurate diagnosis is crucial. The aim of this work was to select specific ssDNA aptamers that bind to the entire Zika virus and can be used to compose specific diagnostics, without cross-reactivity with other flaviviruses. Zika virus was cultivated in Vero cells and used as a target for aptamer selection. Aptamers specific for the ZIKV were selected using whole-virus SELEX, with counterselection for other flavivirus. Secondary and tertiary structures were evaluated and the molecular anchoring between the aptamers and target were simulated by the HDOCK server. Aptamer interaction was evaluated by ELISA/ELASA and the dissociation constant (Kd) was calculated by thermophoresis. Four ZIKV-specific aptamers were selected. The best two were further characterized and proved to be specific for ZIKV. Aptamers are capable of binding specifically to the ZIKV and differentiate from Dengue virus. The aptamers selected in this work can be used as capture agents in the composition of diagnostic tests to specifically detect ZIKV infection.

1. Introduction

Zika virus (ZIKV) is an arthropod-borne virus (arbovirus) primarily transmitted by the Aedes mosquitoes. ZIKV was discovered in 1947 in primates in the Zika Forest of Uganda during surveillance for yellow fever [1,2]. The first documented ZIKV epidemic occurred in the Yap Island in Micronesia in 2007, followed by outbreaks in the Pacific Islands in 2013–2014 [3,4]. In 2015, ZIKV emerged in the Americas, and spread to 87 countries and territories with autochthonous transmission, according to the WHO, by July 2019 [5,6]. Currently, ZIKV is endemic in all tropical areas of the world, similarly to Dengue virus (DENGV), and the chance for future Zika epidemics remains very real [7].
ZIKV is a small enveloped positive-strand RNA virus belonging to the Flavivirus genus of the Flaviviridae family. The RNA genome encodes three structural (core C, membrane precursor prM, and envelope E) and seven non-structural genes (NS1, NS2A, NS2B, NS3, NS4A, NS4B, and NS5), with untranslated region (UTR) genes flanking the 5′ and 3′ ends [8,9]. High-resolution ZIKV structures identified by cryo-electron microscopy indicate that the overall ZIKV structure is similar to that of other flaviviruses. The flavivirus E protein, the major protein involved in receptor binding and fusion, is formed as a head-to-tail dimer on the surface of viral particles. There are 90 dimers of the E protein on the surfaces of mature viruses [10]. Each monomer has four domains: the stem/transmembrane domain (E-S/E-TM) that anchors the protein into the membrane and domains I, II, and III (ectodomains) that constitute the predominantly surface portion of the protein [10,11].
Although most ZIKV infections are asymptomatic (75–80%), and symptomatic infections are generally mild [12], the Zika virus became a major public health problem in early 2015, when cases of the Guillain–Barré syndrome, microcephaly, and neurological disorders were associated with virus infection [13,14,15]. Due to its association with congenital infections, ZIKV has been considered the newest member of the TORCH pathogens (toxoplasmosis, others, rubella, cytomegalovirus, and herpes) in the Americas [16,17].
Brain changes caused by congenital Zika appear in the second and third trimesters of pregnancy, while other dysmorphic features are often observed in the newborns of mothers infected with the virus [18,19]. As time is a very important factor in the evolution of Zika infection in pregnant women, accurate and early detection becomes of paramount importance for vertical transmission control, and future therapeutic strategies.
There is an increasing need for affordable, rapid, sensitive, and high throughput methods in identifying new or re-emerging infectious agents [20]. The availability of the point of care (POC) tests eliminates diagnostic delays and uncertainties, enabling the timely initiation of proper treatments and the prevention of the spread of diseases [20]. The use of aptamer-based POC tests has been increasingly described in the literature [21,22].
Aptamers are small single-stranded oligonucleotides simultaneously developed by two groups of researchers, Ellington and Szostak, and Tuerk and Gold, in 1990 [23,24]. They are also called chemical antibodies because of their ability to specifically bind to a target, just like antibodies [25]. However, while antibodies recognize and bind to protein epitope sequences, aptamers recognize and bind based on the 3D structure of the target molecule [26]. In addition to specific, high-affinity binding to their targets, aptamers offer several advantages, including in vitro selection, lack of batch-to-batch variation, room temperature transport, and a relatively low cost, making them ideal for use in diagnostic assays [27].
This paper describes the selection and characterization of specific ssDNA aptamers that bind to the whole Zika virus that can be used to compose a specific diagnostic test without cross reactivity with other Flaviviruses.

2. Material and Methods

2.1. Zika Virus Strain

The Zika virus isolate ES2916/2015 (accession number: KX212103.1), from the State of Espírito Santo, Brazil, was used as a target in the selection of aptamers.

2.2. Other Flavivirus Strains

Other circulating Flaviviruses used were yellow fever virus/YFV 17DD (purified vaccine strain); Dengue 1/WestPac Virus; Dengue Virus 2/S16803; Dengue Virus 3/16,562; and Dengue virus 4/TVP36, with accession numbers DQ100292.1; U88535.1; GU289914.1; KU725665.1; and KC963424.1, respectively. All viruses were supplied by the Laboratory of Virological Technology (Bio-Manguinhos/Fiocruz).

2.3. Zika Protein E Gene

The gene encoding the ectodomain of Zika virus E protein (E-ZIKV) was synthesized from the genetic sequence available from GenBank (https://www.ncbi.nlm.nih.gov/genbank/ (accessed on 23 March 2017)—Accession code YP_009227198.1—Envelope protein E (Zika virus)), by the company Integrated DNA Technology (IDT). A polyhistidine tail was inserted into the C-terminal portion to facilitate the further purification of the protein.

2.4. Production and Purification of Zika Virus

In order to evaluate the replication profile of the Zika viral isolate, growth kinetics were performed in stationary bottles. This experiment was performed in duplicate, and aliquots were collected daily for the evaluation of viral titer (PFU/mL). Vero cells were seeded at a density of 4 × 104 cells/cm2 in 175 cm2 flasks 2 days before infection. Virus was inoculated at two different MOIs (0.02 and 0.002) and allowed to bind to cells for 1 h at 37 °C/5% CO2. Thereafter, the inoculum was removed, the cells were washed with PBS, and VP/SFM medium supplied with 4 mM glutamine and gentamicin 1% (Thermo Fisher Scientific, Rockford, IL, USA) was added. The virus was harvested every 24 h, the first collection point being eight hours post-infection (hpi). The viral titers were determined by plaque-forming unit (PFU) assay on Vero cells. Briefly, 200 µL of 10-fold serial dilutions were incubated for 1 h on Vero cell monolayers in 6-well plates. The virus inoculum was removed, and the cells were overlaid with medium 199 containing 2% Carboxymethyl cellulose (Sigma-Aldrich, Oakville, ON, Canada), 5% FBS, and 1% penicillin/streptomycin. Four days post-infection, cells were fixed with 5% formaldehyde and stained with 0.04% crystal violet (Sigma-Aldrich) to visualize plaques.
The production of ES 2916/2015 Zika virus was performed using a previously selected viral input and time for harvest. The supernatant was collected and clarified using Sartobran P—0.45/0.2 µm filter capsule (Sartorius). After clarification, the viral suspension was submitted to one purification step of liquid chromatography, in Äkta Purifier 10 Chromatograph (GE Healthcare), using ion exchange chromatography (Sartobind Q75 membrane, Sartorius) in 50 mM Tris/300 mM NaCl/sorbitol 8% solution buffer, pH 8.5. The purified viral titer was calculated using the plaque assay.

2.5. Obtaining the Zika Virus Recombinant Protein E (E-ZIKVre)

The E-ZIKVre gene contained in the pUCIDT commercial vector was amplified by PCR with specific oligonucleotides: E-ZIKVre FW (5′-CAT GCC ATG GGC ATT AGG TGC ATA GGC GTT AGC-3′) and E-ZIKVre RV (5′-CCC AAG CTT CTA ATG GTG GTG ATG GTG ATG C-3′) and then cloned into the pET28a expression vector (Novagen, Madison, WI, USA) using restriction endonucleases NcoI (New England Biolabs, Ipswich, MA, USA) and HindIII (New England Biolabs, Ipswich, MA, USA) and the enzyme T4 DNA ligase (Life Technologies, Carlsbad, CA, USA), following the manufacturer’s recommendations. E. coli TOP10 cells were transformed with the pET28a + E-ZIKVre construct by electroporation. Subsequently, the cells were grown on LB agar plates supplemented with 50 µg/mL kanamycin. After overnight incubation at 37 °C, colonies were individually transferred to tubes containing 5 mL of LB medium and 50 μg/mL kanamycin and then incubated for 16 h with shaking at 200 rpm at 37 °C. Plasmid DNA was extracted using the High Pure Plasmid Isolation kit (Roche Life Science, Manheim, Germany). The pET28a + E-ZIKVre purified plasmid was transformed into SHuffle E. coli competent cells (New England Biolabs, Ipswich, MA, USA) using electroporation. A single positive clone was picked to inoculate 10 mL of Terrific Broth (TB) medium, supplemented with 50 μg/mL kanamycin and 1% v/v glucose, and submitted to agitation at 37 °C overnight. Thereafter, 1 mL of bacterial culture was grown in 25 mL of TB medium, containing 50 μg/mL kanamycin and 1% v/v glucose, at 37 °C with 200 RPM shaking until obtaining an optical density (OD) of 0.8 at 600 nm. Before induction, a sample of the culture (1 mL) was collected, centrifuged at 14,000 rpm and stored. For protein expression, the culture was induced with 1 mM isopropyl β-D-1-thiogalactopyranoside (IPTG) at 37 °C for 4 h. Then, the OD600nm of the culture was read. The cells were subsequently pelleted by centrifugation and the supernatants were discarded.

2.6. Analysis of the Expression and Solubility of the Zika Virus Recombinant Protein E (E-ZIKVre)

The pellet collected in the uninduced culture was resuspended in 25 μL of 1× sample buffer—50 mM Tris-HCl, 2% (w/v) sodium dodecyl sulfate (SDS), 0.1% (w/v) blue bromophenol, 10% (v/v) glycerol, 100 mM 2-mercaptoethanol—for every 0.1 OD600nm and boiled at 100 °C for 5 min.
The induced culture pellets were resuspended in 20 mM Tris-HCl/1 mM EDTA/0,1% triton x-100 buffer and lysed by sonication (3 cycles/30 s) on ice. The insoluble fraction was subsequently separated from the soluble cell lysate by centrifugation (14,000 rpm, 10 min, 4 °C). E-ZIKVre was analyzed using SDS-PAGE and Immunoblot.

2.7. E-ZIKVre Protein Solubilization Study

Different formulations of the lysis buffer were used in order to select the best protein solubilization condition. Six buffers containing Tris-HCl, EDTA, triton x-100, and urea at pH 8.0 were tested. The difference between the lysis buffers was only due to the concentration of the chaotropic agent (urea), which ranged from 1 to 6 M. All formulations were added at a rate of 25 μL for each 0.1 OD600nm and subjected to lysis by sonication in 3 cycles of 30 s on ice.

2.8. Purification of E-ZIKVre Protein

The culture conditions in the shake flasks were scaled to a 2 L bioreactor for biomass production in order to obtain biomass for the purification step by metal ion affinity chromatography (IMAC). To express the recombinant E Protein at a large scale, the recombinant strain pET28a + E-ZIKVre was inoculated into E. coli SHuffle cells grown in a 2 L bioreactor (New Brunswick, St Albans, UK). Biomass growth was induced at an OD600nm of 0.8 with 1 mM IPTG at 37 °C for 4 h. Cell culture was centrifuged at 14,000 rpm for 10 min and the pellet resuspended in wash buffer I [20 mM Tris-HCl + 1 mM EDTA + 0.5% triton x-100 + 1 mM phenylmethylsulphonyl fluoride (PMSF), pH 8.0] in the proportion of 10 mL for each 1 g of pellet. The cell suspension was sonicated for two cycles for 15 s with an interval of 30 s on ice and centrifuged at 8000 rpm for 15 min. The pellet was resuspended in washing buffer II (20 mM Tris-HCl + 1 mM EDTA + 0.5% triton x-100 + 1 M Urea + 1 mM PMSF, pH 8.0) and the processes of sonication and centrifugation were repeated. At the end, the soluble fraction was filtered through a 0.22 μm filter and the recombinant protein was purified by IMAC, using the Ni2+ loaded affinity columns of the His Trap HP model (GE Healthcare) in the AKTA Püre equipment (GE Healthcare). Purification took place in steps of increasing imidazole concentration. The purification fractions were analyzed by SDS-PAGE gel and those with the highest concentration were pooled and dialyzed in a step of desalination and refolding of the recombinant protein. The concentration of E-ZIKVre protein was determined by the bicinchoninic acid (BCA) method, using the Pierce BCA Protein Assay kit (Thermo Scientific).

2.9. In Vitro Selection of Zika Virus ssDNA Aptamers Using Conventional SELEX

The Systematic Evolution of Binding by Exponential Enrichment (SELEX) procedure was implemented using the protocol described by Simmons et al., 2012, with modifications. A random single-stranded DNA library was used to select sequences that recognize Zika virus particles. This library contained 1015 variants of 25 random nucleotides flanked by fixed sequences (5′-GGGAGACAAGAATAAACGCTCAA-25n-TTCGACAGGAGGCTCACAACAGGC-3′). Primers annealing to the flanking regions were used during the selection procedure.
The experiment started with an overnight immobilization of the whole Zika virus to the 96-well ELISA plate at 4 °C. One hundred microliters of the aptamer library (1 mM) were added and incubated for 1 h at room temperature (RT). After incubation, target binding sequences were kept, while unbound nucleotide sequences were removed by washing 3 times with PBS-Tween solution, pH 7.4. The elution of bound aptamers was performed using a salt gradient (1.0 M, 1.2 M, and 1.5 M NaCl) and 3.0 M NaSCN. Eluates were desalted using 5 kDa Microcon filters (Millipore, Life Sciences, Burlington, NJ, USA) and amplified by unidirectional PCR using 10 µL 5× reaction buffer, 1.5 U of Taq polymerase (Invitrogen, Waltham, MA, USA), 200 nM primer sense, 2 nM primer antisense, 20 mM dNTPs, and 50 mM MgCl2. Amplification conditions were 10 min at 95 °C; 100 cycles of 90 s at 95 °C; 30 s at 56 °C; 90 s at 72 °C; and a final elongation of 10 min at 72 °C. The SELEX cycle was repeated seven times.
To increase the specificity of aptamers, negative selection was performed with the yellow fever virus (YFV) and the four dengue virus serotypes after the fifth round.
On the last cycle of SELEX, standard PCR of the selected sequences was performed instead, and the quality of the PCR product was verified in 2% agarose gel.

2.10. Cloning and Sequencing of Oligonucleotides Enriched by SELEX

The PCR product was cloned into a pCR2.1 TOPO (Promega) vector. The plasmid was transformed into E. coli TOP10 and after electroporation, the LB medium was added and the material was grown for 1 h at 37 °C, 200 RPM. The culture was then plated onto LB/agar/ampicillin (100 μg/mL) and the material was incubated overnight at 37 °C. The viable clones were then individually seeded from the colonies in 5 mL of LB/ampicillin and incubated overnight at 37 °C, 200 RPM. The plasmids were extracted following the manufacturer protocol of the High Pure Plasmid Isolation Kit (Roche, Manheim, Germany). The extracted plasmids were quantified by Qubit (Invitrogen, Waltham, MA, USA) and amplified using M13 forward and M13 reverse primers. Sequencing was performed using Sanger’s automated methodology using the Big Dye kit, version 3.1, protocol and the 3500 XL equipment (Applied Biosystems, Waltham, MA, USA). Selected ssDNA molecules were aligned by Clustal 1.2 (Version 1.2.4, Institute Pasteur, Paris, France).

2.11. Generation of Two- and Three-Dimensional Structure

The secondary structure of ssDNA was constructed from their nucleotide sequences using the Mfold Web Server [28]. The sequences were selected as linear at a temperature of 25 °C and ionic concentration of 0.1 M of Na+, 0.01 M of Mg2+, computing only fold configurations within 5% from the minimum free energy, and considering a maximum number of 50 folds with no limit to the maximum distance between paired bases. In addition to the predicted secondary structure of ssDNA, the minimum free energy of the predicted structure is provided.
Three-dimensional modeling was performed using the protocol described by JEDDI and SAIZ, with modifications [29]. Briefly, using the RNAalifold Web Server, ssDNA aptamer sequences in FASTA format were converted into ssRNA (Vienna format), containing the coordinates of the three-dimensional structure. The digital platform “RNA Composer” was used to transform structures from ssRNA in Vienna format into ssRNA in PDB format. Equivalent 3D ssRNA models were constructed using the VMD software, and the 3D RNA models were translated into DNA models using the Chimera program. The final 3D ssDNA structures were visualized by the Chimera software.

2.12. Molecular Dynamics Simulations between Zika Virus E Protein and Aptamers

The Zika virus E protein [30], available in the Protein Data Bank (PDB) database [31], and the 3D ssDNA structures of the aptamers, finalized in the Chimera software, were applied to simulate molecular dynamics (molecular docking) using the HDOCK platform.

2.13. Binding Affinity of Selected Aptamers

2.13.1. ELISA

Enzyme immunoassay plates were coated with 50 µM/well of poly-lysine (PL) at 5 μg/mL then stored at 4 °C for 16 h. Aptamers were immobilized on wells and the plate was incubated overnight at 4 °C. To remove unbound aptamers, wells were washed three times with washing buffer (10 mM PBS, 0.05% Tween-20, pH 7.4,). The plate was blocked with 300 μL/well of a blocking agent (10 mM PBS, 1% BSA, pH 7.4) for 1 h at 37 °C; and another cycle of washes followed and 100 μL/well of different concentrations of virus in 10-fold serial dilutions (5,000,000–0 PFU/mL) was added. After incubation for 1 h at 37 °C, the unbound virus was removed by washing three times. One hundred microliters of monoclonal antibody 4G2 in a 10-fold dilution was added and the plate was incubated at 37 °C for another 1 h. Subsequently, the plate was washed again, 100 µL of anti-mouse conjugated with peroxidase at a dilution of 1/30,000 in PBS was added and the plate was incubated for 1 h at 37 °C. The plate was again washed three times and 100 µL of tetramethylbenzidine (TMB) chromogenic solution was added before 20 min of incubation at 37 °C. The color reaction was stopped by adding 100 µL of 1 N HCl and the absorbance was read at 450 nm using the Spectra Max190 reader (Molecular Devices, San Jose, CA, USA) supplied with the Soft Max Pro 6.3 software.

2.13.2. Affinity Binding by Microscale Thermophoresis

Microscale thermophoresis (MST) was performed on a NanoTemper Monolith NT.115 apparatus (Nano Temper Technologies). Briefly, 10 μM Zika virus E protein was labeled using a Monolith Protein Labeling Kit RED-NHS 2nd Generation, with 3-fold excess NHS dye in PBS (pH 7.4). The free dye was removed according to the manufacturer’s instructions, and the protein was eluted in the MST buffer (PBS containing 0.05% Tween 20) and centrifuged at 14,000× g for 10 min. Binding affinity measurements were performed using 20 nM of labeled protein and a serial dilution 1:1 of 16 concentrations of the four aptamers, starting at 50 μM. The protein was maintained in the presence of the aptamers for 15 min prior to transferring to premium Monolith NT.115 capillaries. Experiments were run in duplicate at 76% excitation and a medium MST power at 25 °C. Data were obtained with MO.Control V2.1.1 software (NanoTemper Technologies, Munich, Germany) and recorded data were analyzed with MO.Affinity Analysis V3.0.4 software (NanoTemper Technologies, Munich, Germany). The dissociation constant (Kd) quantifies the equilibrium of the reaction of the labeled molecule A (concentration cA) with its target T (concentration cT) to form the complex AT (concentration cAT) and is defined by the law of mass action as Kd = cAxcT/cAT, where all concentrations are “free” concentrations. During the titration experiments, the concentration of the labeled molecule A is kept constant, and the concentration of added target T is increased.
In order to assess the specificity of aptamers, we used a pool of viral particles containing all four Dengue serotypes. Two aptamers were selected and labeled with the second-generation RED-NHS labeling kit. The reaction was carried out according to the manufacturer’s instructions in buffer 130 mM NaHCO3, 50 mM NaCl, pH 8.3, applying 10 μM of aptamer (in molar ratio to fluorophore ≈ 1:3) at room temperature for 30 min. Unreacted fluorophore was removed using a size exclusion column equilibrated in PBS buffer, 0.05% Tween 20, pH 7.4. Binding affinity analyzes were conducted using labeled aptamers adjusted to 40 nM with PBS buffer, 0.05% Tween 20, 150 mM NaCl, pH 7.4. For measurement, each virus dilution (ZIKV and DENGV pool) was mixed with a volume of labeled aptamer, at a final aptamer concentration of 20 nM and final virus concentrations ranging from 0.00000652 nM to 0.01335 nM. After 15 min of incubation, samples were loaded into standard capillaries. The effect of thermophoresis was measured using the Monolith NT.115 equipment at room temperature of 25 °C. The parameters were set for 100% or 64% excitation (ZIK01/NH2 or ZIK02/NH2 aptamers, respectively) and average LED power. Data from the duplicates were analyzed using the MST signal at 5 s.

3. Results

3.1. Zika Virus Production

Before creating a working virus seed (WVS), a growth kinetics of the virus in Vero cells was performed to evaluate the replication profile of the Zika virus isolate ES2916/2015. The virus was inoculated at the MOIs 0.02 and 0.002 and was collected every 24 h, the first collection point being 8 hpi. Compared to the cell control, ZIKV-infected cells demonstrated no apparent cytopathic effect (CPE) up to 48 hpi for both MOIs. Infected cells demonstrated mild CPE by 72 hpi, moderate CPE by 96 hpi, and extensive CPE by 96–120 hpi (Figure 1A).
Viral kinetics using MOI 0.02 showed faster growth, as expected, and reached a plateau 12 h before the MOI 0.002 curve. The peak of viral production is between 48 and 72 h at MOI 0.02 (6.19 and 6.21 Log10 PFU/mL, respectively) and 72 h, at MOI 0.002 (6.32 Log10 PFU/mL) (Figure 1B).
For working seed bank production, Vero cell cultures were inoculated with a viral input of 0.02 and the virus was collected after 48 h of infection. Following the purification of the ZIKV, a viral titer of 7.44 Log10 PFU/mL was obtained.

3.2. Cloning, Expression, and Purification of E-ZIKVre Recombinant Protein in E. coli SHuffle System

For the expression of the recombinant protein E-ZIKVre, the gene was amplified by PCR. In Figure 2A, it is possible to observe a band of approximately 1247 bp corresponding to the E-ZIKVre gene. This gene was cloned into pET28a vector and E. coli TOP10 cells were transformed with the pET28a+E-ZIKVre construct. Plasmids of possible recombinant clones were treated with restriction endonucleases NcoI and HindIII, where a band corresponding to the E-ZIKVre insert was observed (Figure 2B). Clone 3 was selected to transform E. coli SHuffle cells.
The expression and solubility of the E-ZIKVre protein was evaluated in induced and non-induced cell extracts using SDS-PAGE (Figure 3A) and immunoblotting (Figure 3B). The expected molecular mass (45 kDa) was observed both in the samples referring to the induced cells and in the samples of the insoluble fraction in SDS-PAGE. The immunoblotting assay with anti-his antibody confirmed these results.
Protein solubilization assays were performed and analyzed by SDS-PAGE with different concentrations of chaotropic agent (urea). The condition of 20 mM Tris buffer + 1 mM EDTA + 0.5% triton x-100 + 4 M urea pH 8.0, 3 cycles of 30 s and 30 s intervals on ice was established as ideal to obtain the complete solubilization of the E-ZIKVre protein (Figure 4).
Cell culture conditions were scaled-up to a 2 L bioreactor in order to obtain sufficient biomass for protein purification. The E-ZIKVre protein was purified by IMAC (Figure 5A) and eluted in buffer containing 500 mM imidazole. Fractions collected during elution were analyzed by SDS-PAGE (Figure 5B), with a high degree of homogeneity >85%. After the purification step, the E-ZIKVre protein was dialyzed to remove the chaotropic agent and refold. The quantitation results demonstrated a recovery of E-ZIKVre protein of approximately 37 mg/L.

3.3. ssDNA Aptamers against Zika Virus

SELEX technology was used for the selection of specific DNA aptamers against whole viral particles. A randomized ssDNA library was incubated with Zika virus particles immobilized on ELISA plates and cycles of interaction, and the separation of the bound from unbound species as well as the elution and amplification of the bound species were performed. After the fifth selection round, negative selection was performed with others flaviviruses, the four dengue serotypes and yellow fever virus, to increase the specificity of aptamers. At the end of the SELEX cycles, the PCR-derived molecules were cloned and sequenced. Ninety-four clones containing an aptamer insert were obtained by elution with 3 M NaSCN. The sequencing of 28 of these clones was performed and their homology was evaluated using a Clustal 2.1 alignment (Figure 6). The formation of four aptamers clusters with 100% nucleotide homology was detected: cluster ZK01 (aptamers ZK04 and ZK83); cluster ZK02 (aptamers ZK06, ZK23 and ZK78); cluster ZK03 (aptamers ZK34 and ZK88); and cluster ZK04 (aptamers ZK17 and ZK86). In addition, the homology of 22 nucleotides (GTTTAGAAACCCCACCGTGGCA) was observed in 10 sequences (Figure 6). The four clusters were chosen from the selected pool for further characterization. Clusters ZK01 and ZK02 were analyzed using their complete sequence, whereas the sequences of clusters ZK3 and ZK04 were truncated to their appropriate predetermined size.

3.4. Secondary and Tertiary Structure Prediction of the ssDNA Aptamers

Secondary structures of the ssDNA molecules were predicted using the Mfold Web Server, based on free energy minimization techniques. In mfold, all possible secondary structures are approximated based on Watson–Crick pairing and the most thermodynamically stable structures are selected.
Figure 7A shows the most stable two-dimensional structures of each aptamer group and their respective minimum-fold free energies (ΔG). Clusters ZK02, ZK03, and ZK04 demonstrated a conservation of structure, a hairpin formed by the final 11 nucleotides of the forward primer and the first five nucleotides of the aptamer sequence.
The transformation of secondary sequences into three-dimensional structures consisted of multiple steps and the use of two programs, VMD and chimera. The results of the three-dimensional analysis confirm the structural similarity between aptamers ZK02, ZK03, and ZK04 (Figure 7B).

3.5. Molecular Docking between ssDNA Aptamers and the E-ZIKVre Protein

Protein E is a glycoprotein that involves the entire surface of flavivirus viral particles. Each glycoprotein E monomer forms stable homodimers, and 90 dimers assemble to form the outer envelope of the infectious virus. Protein E is the main target of neutralizing antibodies against the Zika virus and the likely target of the selected aptamers. We verified this assumption using recombinant protein E in thermophoretic analyses, to confirm that the aptamers indeed bind to this protein, verify the binding profile, and calculate the dissociation constants. We also performed the assays with the intact virus to confirm the aptamer-target specificity, using a Dengue virus pool as a control (see next section).
Having confirmed the interaction of the selected aptamers with the E protein, molecular dynamics simulation between the aptamers and ZIKV E protein was performed using the HDOCK server. HDOCK performs global docking to predict binding complexes between two molecules using a hybrid docking strategy.
We analyzed the first 10 predictions and eliminated those that interacted in the unexposed part of the protein, facing the transmembrane region of protein E of ZIKV. Among the models linked to the outer surface of the protein, we selected the one with the lowest docking score. For aptamers ZK01 and ZK02, model 1 was selected. The aptamers ZK03 and ZK04 presented model 3 as the model with the lowest energy that met the pre-established parameters. Figure 8 shows the sites of interaction between each aptamer and the dimeric E protein.
It is observed that all aptamers present binding sites on two adjacent monomers of the same E-dimer (Figure 8 and Table 1).

3.6. ELISA

An in-house ELISA/ELASA was used as a screening tool to assess the performance of the ZK01, ZK02, ZK03, and ZK04 aptamer clusters. Aptamers were tested against the serial dilution of Zika virus at base 10. Optical density data demonstrated that all sequences were only able to recognize Zika virus at a dilution of 5 × 106 PFU/mL (Figure 9).

3.7. Microscale Thermophoresis

Identification and quantification of interactions between ZIKV and aptamers was analyzed by microscale thermophoresis (MST) using a NanoTemper Monolith NT.115 apparatus, aiming to determine the dissociation constants (Kd) and confirm the specificity of the aptamers.
Labeling of the E-ZIKVre protein was performed with reactive dyes using N-Hydroxysuccinimide (NHS)–ester chemistry, which reacts efficiently with the protein’s primary amines to form highly stable dye–protein conjugates. The efficiency of the labeling step was evaluated by preliminary tests available on the equipment, and the intensity of gross fluorescence within the detection limit could be verified. After the preliminary test, the binding check test was performed to assess whether there was a detectable response amplitude between the target (labeled E-ZIKVre protein) and the ligands (aptamers) at the maximum concentration available. In this test, one can select the best assay buffer for the protein, targeting an optimal MST binding assay. With this, it was possible to evaluate the conditions of protein aggregation, adsorption in the capillaries, fluorescence intensity, and homogeneity to obtain a higher signal/noise ratio and establish the assay conditions to carry out the quantification of the affinity through the binding affinity test.
To investigate the binding affinity between E protein and Zika virus aptamers, 20 nM of the labeled protein (target) and serial dilutions (50–1.53 nM) of the four aptamers were used. As shown in Figure 10, the binding affinity of aptamer ZK01 for E-ZIKVre (Kd = 1414.90 ± 0.41 nM) was slightly higher than for the other aptamers with the greater homology (ZK02 Kd = 761.54 ± 0.49 nM; ZK03 K: 919.26 ± 0.68 nM; ZK44 Kd: 820.98 ± 0.34 nM) (Figure 10).
The evaluation of the specificity of the aptamers was verified using a pool of Dengue serotypes in comparison with the Zika virus, as DENGV is well-known to cross-react with Zika in many diagnostic tests. According to Figure 11, the variation in the concentrations of the DENGV pool did not change the normalized fluorescence response, which suggests no interaction in the concentration range evaluated. On the other hand, the variation of ZIKV concentration caused an increase in the response at higher concentrations, which indicates the presence of intermolecular interaction between aptamers and ZIKV (Figure 11). As there was no complete saturation of the curve, due to insufficient viral titters, the dissociation constant (Kd) could not be calculated.

4. Discussion

Aptamers have been extensively used in the study of arbovirus with some very interesting results, both for diagnostic and therapeutic applications, which have been recently reviewed by our group [32]. Particularly in the area of diagnostics, they have been capable of separating between different Dengue serotypes, but have also been applied for ZIKV, tick-borne encephalitis virus (TBEV), Rift Valley fever virus (RVFV) and Chikungunya Virus (CHIKV) with considerable success [32]. In the case of previous Zika studies using aptamers, the authors used the NS1 protein as a target, and in most cases, their selected aptamers are not capable of differentiating between ZIKV and DENGV NS1 proteins, with one exception [32]. In this study, we are presenting for the first time aptamers selected against the whole viral particle, verifying its binding to the whole virus and the E protein specifically, and demonstrating the specificity of that aptamer and its capacity to differentiate between zika and dengue viruses.
Several studies in the literature have reported the expression of structural and non-structural flavivirus proteins using the E. coli prokaryotic system [33,34,35,36]. This system has the advantage of enabling high levels of expression at a relatively low cost. The flavivirus envelope protein is the most abundant on the surface of the viral particle. Despite the system’s limitation in carrying out post-translational modifications, such as glycosylation, immunogenicity assays using flavivirus envelope proteins produced in the prokaryotic system demonstrated the ability of this protein to induce the formation of neutralizing antibodies [37,38], indicating the conservation of antigenic epitopes.
The recombinant protein E-ZIKVre was expressed in the form of inclusion bodies, as shown in other works in the literature [33,37]. Yields of E-ZIKV protein expression in the prokaryotic system of 0.55 mg/mL and 10 mg/L were reported by Han et al. [37] and Liang et al. [39], respectively. In this work, we obtained a final yield after purification of 37 mg/L of culture at a final concentration of 0.76 mg/mL with a purity > 90%, demonstrating the efficiency of the strategy used.
Although different strategies have been adopted for the selection of aptamers against different arboviruses, as shown above [32], their binding affinity can be influenced by the stringency of the selected conditions, such as, for example, reducing the target concentration or modifying the binding and washing conditions [40]. In this work, we used a protocol based on the interruption of interactions formed between the aptamer and the target using buffers with a high salt content and chaotropic agents. In a previous study carried out by Simmons et al. [41], using the same selection conditions, high homology was observed in the sequences generated in both 1.5 M NaCl and 3 M NaSCN elutions, and the performed assays confirmed binding with the target. In our experiments, in addition to the use of the same elution conditions, we introduced negative rounds of selection against the yellow fever virus and the different dengue serotypes, aiming to select not only the highest binding aptamers, but the ones that were only specific to Zika virus and not for other circulating flaviviruses in the country. We used aptamers resulting only from the 3 M NaSCN elution and we achieved the formation of four aptamer clusters and a highly conserved region of 22 nucleotides and a subset of complete sequence identity, indicating a core sequence necessary for binding, with some variability permitted in the formation of the additional binding loop.
The affinity of the specificity of the aptamers for the ZIKV viral particles was studied by ELISA and microscale thermophoresis, confirming the ability of the aptamers to recognize both the entire viral particle, and the E protein in particular. In the ELISA assay, the negative charge on the DNA backbone, due to the phosphates, makes it difficult for aptamers to adhere to polystyrene plates. A practical and inexpensive alternative is the use of polylysine to modify the surface charge of the bottles and plates for cell culture [42] and plates for enzyme assays [43]. We developed an in-house ELISA with polylysine-coated plates that permitted the immobilization of aptamers in order to perform the ELISA assays. Using this method, we demonstrated that the aptamers were capable of recognizing the ZIKV in ELISA assays. The sensitivity of the assay could potentially be improved, as it is understood that the use of lysine permits de deposition of the aptamers in a disorderly manner, not optimizing the exposure of the aptamer binding sites. Alternative methods of immobilization using biotin–streptavidin systems, or gold nanoparticles with thiol-modified aptamers could improve the performance of the system and permit direct use in immunochromatography-based tests. Nevertheless, our current results permit the verification of the aptamers to bind to whole viral particles. This was further confirmed in the MST, where it was possible to identify not only the interaction of the aptamers to ZIKV, but also the inability of the aptamer to interact with DENGV particles at the same viral count.
Finally, it was supposed that the most probable target on the viral particle surface would have been the E protein. This was verified by MST, where the affinity of the aptamer for a recombinant purified ZIKV E protein was calculated. Post-SELEX modifications with functional groups have been described to optimize the detection (fluorescent molecules) or immobilization (biotin, amine) of aptamers [40]. In the development of this project into a viable diagnostic, we aimed to investigate post-SELEX modification to increase aptamer affinity, as well as optimize the assay development format to take advantage of the capability of the aptamer to bind to ZIKV and differentiate from DENGV particles, into developing novel diagnostic strategies for Zika virus infections.

Author Contributions

All authors listed have made a substantial, direct, and intellectual contribution to the work. More specifically, L.M.d.M. is the principal contributor in this work, having personally executed or participated in all aspects of the project. T.S.C. and M.A.M. were responsible for the expression and purification of the E protein, K.A.B.P. and P.B.J. offered support in the analytical techniques and the MST execution, training, and interpretation of the data. S.M.B.d.L. is responsible for the viral research laboratory and contributed to the production of the ZIKV, YFV, and DENGV used in this research. A.M.B.d.F. is the project supervisor who contributed to technical discussions and support on flavivirus, and the revision of the overall project. S.M. is the project supervisor, in addition to conducting the technical discussions and manuscript preparation, was responsible for the training and methodology of the aptamer selection. All authors have read and agreed to the published version of the manuscript.

Funding

Fundação de Amparo à Pesquisa do Estado do Rio de Janeiro: E_18/2015—Research program on Zika, Chikungunya and Dengue in the State of RJ.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We acknowledge Dilson Silva, Henrique Francisco Rocha, Mariana Pierre de Barros Gomes, Noemi Rovaris Gardinali, Ygara da Silva Mendes, and Renata Chagas Bastos, for helpful technical support and discussions. This work was supported by the Fundação de Amparo à Pesquisa do Estado do Rio de Janeiro—FAPERJ.

Conflicts of Interest

The authors declare that they have no conflict of interest related to their research work.

References

  1. Dick, G.W.A. Zika virus (II). Pathogenicity and physical properties. Trans. R. Soc. Trop. Med. Hyg. 1952, 46, 521–534. [Google Scholar] [CrossRef]
  2. Dick, G.W.A.; Kitchen, S.F.; Haddow, A.J. Zika Virus (I). Isolations and serological specificity. Trans. R. Soc. Trop. Med. Hyg. 1952, 46, 509–520. [Google Scholar] [CrossRef]
  3. Duffy, M.R.; Chen, T.H.; Hancock, W.T.; Powers, A.M.; Kool, J.L.; Lanciotti, R.S.; Pretrick, M.; Marfel, M.; Holzbauer, S.; DuBray, C.; et al. Zika Virus Outbreak on Yap Island, Federated States of Micronesia. N. Engl. J. Med. 2009, 360, 2536–2543. [Google Scholar] [CrossRef]
  4. Cao-Lormeau, V.M.; Musso, D. Emerging arboviruses in the Pacific. Lancet 2014, 384, 1571–1572. [Google Scholar] [CrossRef]
  5. World Health Organization. Zika Epidemiology Update. July 2019. Available online: https://www.who.int/emergencies/diseases/zika/zika-epidemiology-update-july-2019.pdf?ua=1 (accessed on 18 July 2019).
  6. Zanluca, C.; de Melo, V.C.A.; Mosimann, A.L.P.; dos Santos, G.I.V.; dos Santos, C.N.D.; Luz, K. First report of autochthonous transmission of Zika virus in Brazil. Mem. Inst. Oswaldo Cruz 2015, 110, 569–572. [Google Scholar] [CrossRef]
  7. Castanha, P.M.S.; Marques, E.T.A. Vaccine development during global epidemics: The Zika experience. Lancet Infect. Dis. 2020, 20, 998–999. [Google Scholar] [CrossRef]
  8. Lindenbach, B.D.; Rice, C.M. Molecular biology of flaviviruses. In Advances in Virus Research; Elsevier: Amsterdam, The Netherlands, 2003; pp. 23–61. [Google Scholar] [CrossRef]
  9. Kuno, G.; Chang, G.J.J. Full-length sequencing and genomic characterization of Bagaza, Kedougou, and Zika viruses. Arch. Virol. 2007, 152, 687–696. [Google Scholar] [CrossRef]
  10. Kostyuchenko, V.A.; Lim, E.X.Y.; Zhang, S.; Fibriansah, G.; Ng, T.S.; Ooi, J.; Shi, J.; Lok, S.-M. Structure of the thermally stable Zika virus. Nature 2016, 533, 425–428. [Google Scholar] [CrossRef] [PubMed]
  11. Sirohi, D.; Chen, Z.; Sun, L.; Klose, T.; Pierson, T.C.; Rossmann, M.G.; Kuhn, R.J. The 3.8 Å resolution cryo-EM structure of Zika virus. Science 2016, 352, 467–470. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Musso, D.; Cao-Lormeau, V.M.; Gubler, D.J. Zika virus: Following the path of dengue and chikungunya? Lancet 2015, 386, 243–244. [Google Scholar] [CrossRef]
  13. Oehler, E.; Watrin, L.; Larre, P.; Leparc-Goffart, I.L.; Lastère, S.; Valour, F.; Baudouin, L.; Mallet, H.P.; Musso, D.; Ghawche, F. Zika virus infection complicated by Guillain-Barre syndrome—Case report, French Polynesia, December 2013. Eurosurveillance 2014, 19, 20720. [Google Scholar] [CrossRef] [PubMed]
  14. Cao-Lormeau, V.M.; Blake, A.; Mons, S.; Lastère, S.; Roche, C.; Vanhomwegen, J.; Dub, T.; Baudouin, L.; Teissier, A.; Larre, P.; et al. Guillain-Barré Syndrome outbreak associated with Zika virus infection in French Polynesia: A case-control study. Lancet 2016, 387, 1531–1539. [Google Scholar] [CrossRef]
  15. Schuler-Faccini, L.; Ribeiro, E.M.; Feitosa, I.M.L.; Horovitz, D.D.; Cavalcanti, D.P.; Pessoa, A.; Doriqui, M.J.R.; Neri, J.I.; Neto, J.M.D.P.; Wanderley, H.Y.; et al. Possible Association Between Zika Virus Infection and Microcephaly—Brazil, 2015. MMWR Morb. Mortal. Wkly. Rep. 2016, 65, 59–62. [Google Scholar] [CrossRef]
  16. Coyne, C.B.; Lazear, H.M. Zika virus—Reigniting the TORCH. Nat. Rev. Microbiol. 2016, 14, 707–715. [Google Scholar] [CrossRef]
  17. Kovacs, A.A.Z. Zika, the Newest TORCH Infectious Disease in the Americas. Clin. Infect. Dis. 2020, 70, 2673–2674. [Google Scholar] [CrossRef]
  18. de Miranda-Filho, D.B.; Martelli, C.M.T.; de Ximenes, R.A.A.; Araújo, T.V.B.; Rocha, M.A.W.; Ramos, R.C.F.; Dhalia, R.; de Oliveira França, R.F.; de Azevedo Marques Júnior, E.T.; Rodrigues, L.C. Initial Description of the Presumed Congenital Zika Syndrome. Am. J. Public Health 2016, 106, 598–600. [Google Scholar] [CrossRef]
  19. Oliveira Melo, A.S.; Malinger, G.; Ximenes, R.; Szejnfeld, P.O.; Alves Sampaio, S.; Bispo de Filippis, A.M. Zika virus intrauterine infection causes fetal brain abnormality and microcephaly: Tip of the iceberg? Ultrasound Obstet. Gynecol. 2016, 47, 6–7. [Google Scholar] [CrossRef]
  20. Li, H.Y.; Jia, W.N.; Li, X.Y.; Zhang, L.; Liu, C.; Wu, J. Advances in detection of infectious agents by aptamer-based technologies. Emerg. Microbes Infect. 2020, 9, 1671–1681. [Google Scholar] [CrossRef]
  21. Ma, X.; Chen, Z.; Zhou, J.; Weng, W.; Zheng, O.; Lin, Z.; Guo, L.; Qiu, B.; Chen, G. Aptamer-based portable biosensor for platelet-derived growth factor-BB (PDGF-BB) with personal glucose meter readout. Biosens. Bioelectron. 2014, 55, 412–416. [Google Scholar] [CrossRef]
  22. Gopinath, S.C.B.; Lakshmipriya, T.; Chen, Y.; Phang, W.M.; Hashim, U. Aptamer-based ‘point-of-care testing’. Biotechnol. Adv. 2016, 34, 198–208. [Google Scholar] [CrossRef]
  23. Ellington, A.D.; Szostak, J.W. In vitro selection of RNA molecules that bind specific ligands. Nature 1990, 346, 818–822. [Google Scholar] [CrossRef]
  24. Tuerk, C.; Gold, L. Systematic evolution of ligands by exponential enrichment: RNA ligands to bacteriophage T4 DNA polymerase. Science 1990, 249, 505–510. [Google Scholar] [CrossRef]
  25. Jayasena, S.D. Aptamers: An Emerging Class of Molecules That Rival Antibodies in Diagnostics. Clin. Chem. 1999, 45, 1628–1650. [Google Scholar] [CrossRef]
  26. Wang, C.H.; Chang, C.P.; Lee, G.B. Integrated microfluidic device using a single universal aptamer to detect multiple types of influenza viruses. Biosens. Bioelectron. 2016, 86, 247–254. [Google Scholar] [CrossRef]
  27. Dunn, M.R.; Jimenez, R.M.; Chaput, J.C. Analysis of aptamer discovery and technology. Nat. Rev. Chem. 2017, 1, 76. [Google Scholar] [CrossRef]
  28. Zuker, M. Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003, 31, 3406–3415. [Google Scholar] [CrossRef]
  29. Jeddi, I.; Saiz, L. Three-dimensional modeling of single stranded DNA hairpins for aptamer-based biosensors. Sci. Rep. 2017, 7, 1178. [Google Scholar] [CrossRef]
  30. Dai, L.; Shi, Y.; Qi, J.; Gao, G.F. 5JHM: Crystal Structure of Zika Virus Envelope Protein. Protein Data Bank. Available online: https://www.rcsb.org/structure/5jhm (accessed on 13 June 2019).
  31. Berman, H.M.; Westbrook, J.; Feng, Z.; Gilliland, G.; Bhat, T.N.; Weissig, H.; Shindyalov, I.N.; Bourne, P.E. The Protein Data Bank. Nucleic Acids Res. 2000, 28, 235–242. [Google Scholar] [CrossRef]
  32. Argondizzo, A.P.C.; Silva, D.; Missailidis, S. Application of Aptamer-based Assays to the Diagnosis of Arboviruses Important for Public Health in Brazil. Int. J. Mol. Sci. 2021, 22, 159. [Google Scholar] [CrossRef]
  33. Dai, L.; Song, J.; Lu, X.; Deng, Y.Q.; Musyoki, A.M.; Cheng, H.; Zhang, Y.; Yuan, Y.; Song, H.; Haywood, J.; et al. Structures of the Zika Virus Envelope Protein and Its Complex with a Flavivirus Broadly Protective Antibody. Cell Host Microbe 2016, 19, 696–704. [Google Scholar] [CrossRef]
  34. Fan, J.; Liu, Y.; Xie, X.; Zhang, B.; Yuan, Z. Inhibition of japanese encephalitis virus infection by flavivirus recombinant e protein domain III. Virol. Sin. 2013, 28, 152–160. [Google Scholar] [CrossRef]
  35. Tripathi, N.K.; Shrivastava, A.; Biswal, K.C.; Rao, P.V.L. Recombinant dengue virus type 3 envelope domain III protein from Escherichia coli. Biotechnol. J. 2011, 6, 604–608. [Google Scholar] [CrossRef]
  36. Ganguly, A.; Malabadi, R.B.; Das, D.; Suresh, M.R.; Sunwoo, H.H. Enhanced Prokaryotic Expression of Dengue Virus Envelope Protein. J. Pharm. Pharm. Sci. 2013, 16, 609–621. [Google Scholar] [CrossRef]
  37. Han, J.F.; Qiu, Y.; Yu, J.Y.; Wang, H.J.; Deng, Y.Q.; Li, X.F.; Zhao, H.; Sun, H.-X.; Qin, C.-F. Immunization with truncated envelope protein of Zika virus induces protective immune response in mice. Sci. Rep. 2017, 7, 10047. [Google Scholar] [CrossRef]
  38. Coller, B.A.G.; Clements, D.E.; Bett, A.J.; Sagar, S.L.; ter Meulen, J.H. The Development of Recombinant Subunit Envelope-Based Vaccines to Protect against Dengue Virus Induced Disease. Vaccine 2011, 29, 7267–7275. [Google Scholar] [CrossRef]
  39. Liang, H.; Yang, R.; Liu, Z.; Li, M.; Liu, H.; Jin, X. Recombinant Zika virus envelope protein elicited protective immunity against Zika virus in immunocompetent mice. PLoS ONE 2018, 13, e0194860. [Google Scholar] [CrossRef]
  40. Stoltenburg, R.; Reinemann, C.; Strehlitz, B. SELEX—A (r)evolutionary method to generate high-affinity nucleic acid ligands. Biomol. Eng. 2007, 24, 381–403. [Google Scholar] [CrossRef]
  41. Simmons, S.C.; McKenzie, E.A.; Harris, L.K.; Aplin, J.D.; Brenchley, P.E.; Velasco-Garcia, M.N.; Missailidis, S. Development of Novel Single-Stranded Nucleic Acid Aptamers against the Pro-Angiogenic and Metastatic Enzyme Heparanase (HPSE1). PLoS ONE 2012, 7, e37938. [Google Scholar] [CrossRef]
  42. Wang, J.; Tian, L.; Chen, N.; Ramakrishna, S.; Mo, X. The cellular response of nerve cells on poly-l-lysine coated PLGA-MWCNTs aligned nanofibers under electrical stimulation. Mater. Sci. Eng. C 2018, 91, 715–726. [Google Scholar] [CrossRef]
  43. Stearns, N.A.; Zhou, S.; Petri, M.; Binder, S.R.; Pisetsky, D.S. The Use of Poly-L-Lysine as a Capture Agent to Enhance the Detection of Antinuclear Antibodies by ELISA. PLoS ONE 2016, 11, e0161818. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Viral kinetics of ES 2916/2015 Zika virus (ZIKV) in VERO cells at MOI 0.02 and MOI 0.002. (A) Cell monolayer images demonstrating the cytopathic effect of ZIKV over time. (B) Viral titers determined by plaque-forming unit (PFU) assay.
Figure 1. Viral kinetics of ES 2916/2015 Zika virus (ZIKV) in VERO cells at MOI 0.02 and MOI 0.002. (A) Cell monolayer images demonstrating the cytopathic effect of ZIKV over time. (B) Viral titers determined by plaque-forming unit (PFU) assay.
Viruses 14 01867 g001
Figure 2. One percent agarose gel. (A) Electrophoresis of the PCR product of the Zika virus recombinant protein E (E-ZIKVre) gene in the commercial vector pUCIDT + E-ZIKVre. Lane 1—Molecular marker 1kb DNA Ladder™ (Life Technologies); line 2–4—PCR product E-ZIKVb; lane 5—negative control. (B) Electrophoresis of the enzymatic cleavage digestion products of the pET28a + E-ZIKVre construct with restriction endonucleases NcoI and HindIII. Lane 1—Marker 1kb DNA Ladder™ (Life Technologies); lane 2—digested clone 1; lane 3—undigested clone 1; lane 4—digested clone 2; lane 5—undigested clone 2; lane 6—digested clone 3; lane 7—undigested clone 3; lane 8—digested clone 4; lane 9—undigested clone 4; lane 10—digested clone 5; lane 11—undigested clone 5; lane 12—digested clone 6; and lane 13—undigested clone 6.
Figure 2. One percent agarose gel. (A) Electrophoresis of the PCR product of the Zika virus recombinant protein E (E-ZIKVre) gene in the commercial vector pUCIDT + E-ZIKVre. Lane 1—Molecular marker 1kb DNA Ladder™ (Life Technologies); line 2–4—PCR product E-ZIKVb; lane 5—negative control. (B) Electrophoresis of the enzymatic cleavage digestion products of the pET28a + E-ZIKVre construct with restriction endonucleases NcoI and HindIII. Lane 1—Marker 1kb DNA Ladder™ (Life Technologies); lane 2—digested clone 1; lane 3—undigested clone 1; lane 4—digested clone 2; lane 5—undigested clone 2; lane 6—digested clone 3; lane 7—undigested clone 3; lane 8—digested clone 4; lane 9—undigested clone 4; lane 10—digested clone 5; lane 11—undigested clone 5; lane 12—digested clone 6; and lane 13—undigested clone 6.
Viruses 14 01867 g002
Figure 3. Evaluation of E-ZIKVre expression and solubility analysis in E. coli SHuffle. (A) SDS-PAGE gel. Line 1—benchmark; line 2—non-induced culture; line 3—induced culture; line 4—soluble fraction; line 5—insoluble fraction. (B) Immunoblot assay. Line 1—dual color; line 2—non-induced culture; line 3—induced culture; line 4—soluble fraction; line 5—insoluble fraction. Primary antibody: Anti-His.
Figure 3. Evaluation of E-ZIKVre expression and solubility analysis in E. coli SHuffle. (A) SDS-PAGE gel. Line 1—benchmark; line 2—non-induced culture; line 3—induced culture; line 4—soluble fraction; line 5—insoluble fraction. (B) Immunoblot assay. Line 1—dual color; line 2—non-induced culture; line 3—induced culture; line 4—soluble fraction; line 5—insoluble fraction. Primary antibody: Anti-His.
Viruses 14 01867 g003
Figure 4. SDS-PAGE of the E-ZIKVre protein solubilization assay expressed in E. coli SHuffle. Line 1—benchmark; line 2—crude extract 4 M urea; line 3—soluble fraction 4 M urea; line 4—insoluble fraction 4M urea; line 5—crude extract 5 M urea; line 6—soluble fraction 5 M urea; line 7—insoluble fraction 5 M urea; line 8—crude extract 6 M urea; line 9—soluble fraction 6 M urea; line 10—insoluble fraction 6 M urea.
Figure 4. SDS-PAGE of the E-ZIKVre protein solubilization assay expressed in E. coli SHuffle. Line 1—benchmark; line 2—crude extract 4 M urea; line 3—soluble fraction 4 M urea; line 4—insoluble fraction 4M urea; line 5—crude extract 5 M urea; line 6—soluble fraction 5 M urea; line 7—insoluble fraction 5 M urea; line 8—crude extract 6 M urea; line 9—soluble fraction 6 M urea; line 10—insoluble fraction 6 M urea.
Viruses 14 01867 g004
Figure 5. Analysis of the E-ZIKVre protein purification step. (A) Chromatogram of purification of protein by IMAC on Ni2+ column. On the vertical axis are the absorbance values in mAU (l = 280 nm) and on the horizontal axis are the volumes in milliliters. The blue peak represents the E-ZIKVre protein. Elution in 20 mM Tris-HCl + 1 mM EDTA + 0.5% triton x-100 + 4M urea (pH 8.0). (B) SDS-PAGE of E-ZIKVre protein purification product by IMAC. Line 1—benchmark; lines 2–5—purified fractions of E-ZIKVre protein.
Figure 5. Analysis of the E-ZIKVre protein purification step. (A) Chromatogram of purification of protein by IMAC on Ni2+ column. On the vertical axis are the absorbance values in mAU (l = 280 nm) and on the horizontal axis are the volumes in milliliters. The blue peak represents the E-ZIKVre protein. Elution in 20 mM Tris-HCl + 1 mM EDTA + 0.5% triton x-100 + 4M urea (pH 8.0). (B) SDS-PAGE of E-ZIKVre protein purification product by IMAC. Line 1—benchmark; lines 2–5—purified fractions of E-ZIKVre protein.
Viruses 14 01867 g005
Figure 6. Alignment of the sequences obtained after selection against the whole Zika virus particle. Red rectangles indicate four clusters of aptamers. Blue square shows nucleotide homology between sequences. The nucleotides in bold were removed.
Figure 6. Alignment of the sequences obtained after selection against the whole Zika virus particle. Red rectangles indicate four clusters of aptamers. Blue square shows nucleotide homology between sequences. The nucleotides in bold were removed.
Viruses 14 01867 g006
Figure 7. Prediction of secondary and tertiary structures of the four selected aptamers for Zika virus. (A) Secondary structure of the aptamers by Mfold web server. Parameters: linear DNA; folding temperature = 25 °C; [Na+] = 0.1; [Mg++] = 0.01; unit = M; ΔG = free energy; the circled areas show the same hairpins structures in 3 aptamers. (B) Computational prediction of the 3D structures of the aptamers using Chimera program. The highlighted areas indicate the same tertiary structure in three aptamers.
Figure 7. Prediction of secondary and tertiary structures of the four selected aptamers for Zika virus. (A) Secondary structure of the aptamers by Mfold web server. Parameters: linear DNA; folding temperature = 25 °C; [Na+] = 0.1; [Mg++] = 0.01; unit = M; ΔG = free energy; the circled areas show the same hairpins structures in 3 aptamers. (B) Computational prediction of the 3D structures of the aptamers using Chimera program. The highlighted areas indicate the same tertiary structure in three aptamers.
Viruses 14 01867 g007
Figure 8. Molecular dynamics simulation of the selected aptamers and dimeric ZIKV E protein by the HDOCK server. Protein E is represented at two angles: the first, showing the aptamer binding on the outside of the protein; the second, with a 90° rotation, detailing the exact location of interaction.
Figure 8. Molecular dynamics simulation of the selected aptamers and dimeric ZIKV E protein by the HDOCK server. Protein E is represented at two angles: the first, showing the aptamer binding on the outside of the protein; the second, with a 90° rotation, detailing the exact location of interaction.
Viruses 14 01867 g008
Figure 9. In-house ELISA: (A) Reactivity of the four aptamer clusters with the ES 2916/2015 Zika virus and negative control (no virus); and (B) Schematic representation of the test.
Figure 9. In-house ELISA: (A) Reactivity of the four aptamer clusters with the ES 2916/2015 Zika virus and negative control (no virus); and (B) Schematic representation of the test.
Viruses 14 01867 g009
Figure 10. Thermophoretic assay to compare affinity between aptamers and E-ZIKVre protein. MST parameters: target concentration: 20 nM; ligand concentration: 50 µM–0.00153 µM; excitation power: 76%; temperature: 25 °C.
Figure 10. Thermophoretic assay to compare affinity between aptamers and E-ZIKVre protein. MST parameters: target concentration: 20 nM; ligand concentration: 50 µM–0.00153 µM; excitation power: 76%; temperature: 25 °C.
Viruses 14 01867 g010
Figure 11. Thermophoretic assay to compare affinity between aptamers and Zika virus (ZIKV) and Dengue virus (DENGV) particles. The figure shows the plot of baseline-corrected normalized fluorescence (ΔFnorm) (‰) versus ligand concentration, in a heating time of five seconds.
Figure 11. Thermophoretic assay to compare affinity between aptamers and Zika virus (ZIKV) and Dengue virus (DENGV) particles. The figure shows the plot of baseline-corrected normalized fluorescence (ΔFnorm) (‰) versus ligand concentration, in a heating time of five seconds.
Viruses 14 01867 g011
Table 1. Docking position of aptamers in the E protein.
Table 1. Docking position of aptamers in the E protein.
ZK01ZK02ZK03ZK04
aaPosition #aaPosition #aaPosition #aaPosition #
SER64AMET68AARG2AER306A
SER66AALA69AVAL46ATHR309A
GLN89ASER70ATHR47AALA310A
LYS118AASP71AGLU62AALA311A
PHE119ASER72ASER64ATHR313A
ALA120AARG73ASER122APHE314A
CYS121ATYR81ALYS123ATHR315A
SER122AASP83ALYS124AILE317A
LYS123AGLY102AARG138AGLN331A
LYS124AASN103AMET140ATYR332A
GLY145AHIS249AGLU162AALA333A
SER146AALA250ALYS166AGLY334A
GLU162ALYS251AASN207ATHR335A
ASN163AARG252AGLU262AASP336A
ARG164AARG2BARG283AGLU367A
GLY182AHIS27BTHR313ASER368A
ASP230AVAL46BPHE314ATHR369A
THR231ATHR47BTHR315AGLU370A
GLY232AMET140BGLN331AASN371A
THR233ASER146BGLU367ASER70B
PRO234AGLU162BGLU370AASP71B
TRP236AASN163BASN371ASER72B
ASP278AARG164BLYS373AARG73B
GLY279ALYS166BGLU62BCYS74B
ARG283AASP278BALA63BTHR76B
LYS301AARG283BSER64BGLN77B
THR309ALYS301BILE65BGLY78B
ALA310AVAL303BSER66BGLU79B
ALA311ATYR305BASP67BTYR81B
THR313AGLN331BMET68BLEU82B
PHE314ATHR335BGLN77BASP83B
THR315APRO338BTYR90BGLY104B
ILE317ATHR366BTRP101B
GLN331AGLU367BGLY102B
ALA333ASER368BASN103B
GLY334ATHR369BLYS118B
THR335AGLU370BSER122B
VAL364AASN371BLYS251B
THR366ASER372BARG252B
GLU367ALYS373BTHR254B
SER368A VAL255B
THR369A VAL256B
GLU370A VAL257B
ASN371A LEU258B
SER372A SER260B
LYS373A
MET374A
MET375A
LYS394A
ARG73B
CYS74B
GLN77B
TYR81B
ASP83B
GLY104B
# A = monomer A; B = monomer B.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Morais, L.M.d.; Chaves, T.S.; Medeiros, M.A.; Pereira, K.A.B.; Jurgilas, P.B.; Barbosa de Lima, S.M.; Missailidis, S.; Bispo de Filippis, A.M. Selection and Characterization of Single-Stranded DNA Aptamers of Diagnostic Potential against the Whole Zika Virus. Viruses 2022, 14, 1867. https://doi.org/10.3390/v14091867

AMA Style

Morais LMd, Chaves TS, Medeiros MA, Pereira KAB, Jurgilas PB, Barbosa de Lima SM, Missailidis S, Bispo de Filippis AM. Selection and Characterization of Single-Stranded DNA Aptamers of Diagnostic Potential against the Whole Zika Virus. Viruses. 2022; 14(9):1867. https://doi.org/10.3390/v14091867

Chicago/Turabian Style

Morais, Liliane Monteiro de, Thiago Santos Chaves, Marco Alberto Medeiros, Kaique Alves Brayner Pereira, Patrícia Barbosa Jurgilas, Sheila Maria Barbosa de Lima, Sotiris Missailidis, and Ana Maria Bispo de Filippis. 2022. "Selection and Characterization of Single-Stranded DNA Aptamers of Diagnostic Potential against the Whole Zika Virus" Viruses 14, no. 9: 1867. https://doi.org/10.3390/v14091867

APA Style

Morais, L. M. d., Chaves, T. S., Medeiros, M. A., Pereira, K. A. B., Jurgilas, P. B., Barbosa de Lima, S. M., Missailidis, S., & Bispo de Filippis, A. M. (2022). Selection and Characterization of Single-Stranded DNA Aptamers of Diagnostic Potential against the Whole Zika Virus. Viruses, 14(9), 1867. https://doi.org/10.3390/v14091867

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop