Impact Factor 2020 : 1.545 (@Clarivate Analytics)
  • Users Online: 424
  • Print this page
  • Email this page

Table of Contents
Year : 2020  |  Volume : 10  |  Issue : 1  |  Page : 23-32

Evolution of specific RNA aptamers via SELEX targeting recombinant human CD36 protein: A candidate therapeutic target in severe malaria

Institute for Research in Molecular Medicine (INFORMM), Health Campus, Universiti Sains Malaysia, 16150, Kubang Kerian, Kelantan, Malaysia

Date of Submission15-Jun-2019
Date of Decision20-Jul-2019
Date of Acceptance01-Dec-2019
Date of Web Publication24-Dec-2019

Correspondence Address:
Khairul Mohd Fadzli Mustaffa
Institute for Research in Molecular Medicine (INFORMM), Health Campus, Universiti Sains Malaysia, 16150, Kubang Kerian, Kelantan
Login to access the Email id

Source of Support: This work is supported by an Exploratory Research Grant Scheme (ERGS; 203/CIPPM/6730106) and Higher Institution Centre of Excellent (HICoE; 311/CIPPM/4401005) from the Malaysian, Ministry of Higher Education and UniversitiSains Malaysia Fellowship Scheme, Conflict of Interest: None

DOI: 10.4103/2221-1691.273091

Get Permissions


Objective: To isolate and characterize RNA aptamers that are specific to human CD36 protein using systematic evolution of ligands by exponential enrichment (SELEX) technology to identify candidates for adjunct therapy to reverse the binding of Plasmodium-infected erythrocytes.
Methods: RNA aptamers were isolated using nitrocellulose membrane-based SELEX and binding analysis was screened using an electrophoretic mobility shift assay and enzyme-linked oligonucleotide assay.
Results: Thirteen cycles of nitrocellulose membrane-based SELEX yielded three aptamers (RC60, RC25, RC04) exhibiting high binding against CD36 protein as shown on electrophoretic mobility shift assay. The sequence analysis revealed a G-quadruplex sequence within all the isolated aptamers that might contribute to aptamer binding and thermodynamic stability. The specificity assay further showed that RC60 and RC25 were highly specific to CD36. The competitive inhibition assay demonstrated that RC60 and RC25 shared a similar binding epitope recognized by mAb FA6-152, a specific monoclonal antibody against CD36.
Conclusions: RC60 and RC25 are promising candidates as anti-cytoadherence for severe malaria adjunct therapy.

Keywords: RNA aptamers; SELEX; CD36 protein; Severe malaria; Adjunct therapy; Cytoadherence

How to cite this article:
Nik Kamarudin NA, Sat JN, Mohd Zaidi NF, Mustaffa KM. Evolution of specific RNA aptamers via SELEX targeting recombinant human CD36 protein: A candidate therapeutic target in severe malaria. Asian Pac J Trop Biomed 2020;10:23-32

How to cite this URL:
Nik Kamarudin NA, Sat JN, Mohd Zaidi NF, Mustaffa KM. Evolution of specific RNA aptamers via SELEX targeting recombinant human CD36 protein: A candidate therapeutic target in severe malaria. Asian Pac J Trop Biomed [serial online] 2020 [cited 2022 Jan 28];10:23-32. Available from:

  1. Introduction Top

Malaria is one of the significant parasitic infection in human causing 429 000 deaths worldwide, mostly in the endemic areas especially in Africa region, followed by the South-East Asia region and the Eastern Mediterranean region[1]. Most of the cases are affecting children under five years old. Even though the malaria mortality rate among children is reduced to 35% from 2010 until 2015, but the disease remains a significant killer, claiming the life of 1 child every 2 minutes[2]. Plasmodium falciparum is the most significant agent of malaria in human and responsible for the enormous burden of global mortality and morbidity[3]. A unique characteristic of the mature infected erythrocytes (IEs) to sequestrate the microvascular surface protein has been associated with the development of severe malaria (SM), which leads microcirculatory obstruction, impaired tissue perfusion and inflammatory cells activation[4]. This sequestration mechanism allows the parasite to escape from spleen-dependent killing mechanisms and act as a determinant for adhesion-based complications of malaria infections such as cerebral malaria and pregnancy-associated malaria[5].

The CD36 protein is a multiligand scavenger receptor that facilitates the binding and uptake of a wide variety of particulate ligands such as oxidized low-density lipoproteins, bacteria, β-amyloid plaque, and apoptotic cells by macrophages[6]. The CD36 protein is involved in many biological processes including mediating the sequestration of IEs in microvascular capillaries, thrombosis, anti-angiogenesis, inflammation, phagocytosis and endocytosis[7],[8]. The CD36-mediated cytoadherence of IEs conceivably will contribute to the dysfunction of vital organs such as kidney, lung, and liver by impairing the microcirculatory blood flow in malaria patients[8]. Interestingly, this protein is expressed at a low level in the brain and uninducible by inflammatory cytokines such as tumor necrosis factor alpha (TNF-α). The binding of IEs to CD36 receptor protein commonly caused uncomplicated malaria but it also can indirectly support the cytoadhesion of the IEs to the microvascular brain which can lead to the development of SM through binding to multiple receptors including endothelial protein C receptor[9]. The expressed CD36 on the platelet molecules has been reported to facilitate the platelet clumping and endothelial microvascular obstruction due to sequestered IEs which strongly associated with SM[4]. A recent report using in vivo study of Plasmodium berghei infection in mice showed that sequestration of the parasite is the CD36-dependent type that can promote parasite growth and survival in the host where sequestration occurred at the CIDRα2-6 domain of PfEMP1[10],[11]. Previously, monoclonal antibodies (mAb) FA6-152 and OKM5 have been reported to bind to a specific immune-dominant of 153-183 and 139-184 amino acid of the CD36 protein, respectively to inhibit the cytoadherence of IEs[12]. Even though the successful finding of mAb to inhibit and reverse the cytoadherence, there are some limitations of the use of mAb in immunotherapy including high immunogenicity, expensive cost, batch to batch variation, larger size, thermal instability and chemical modification[13].

In past 20 years, Ellington and Szostak had discovered new small molecules that can be used for therapy and diagnostic which mimic mAb called aptamer[13]. The aptamer is a synthetic oligonucleotide in the form of single-stranded DNA (ssDNA) or RNA, which could bind to a variety of targets at high specificity and affinity. Aptamers bind to their target via an “induced fit” mechanism, similar to antibodies[14]. Even though both DNA and RNA aptamer can form complex structures, the RNA aptamer can form a more diverse tertiary structure as RNA contains 2’-OH group on their ribose sugar that contributes to tighter and specific binding to the target[15]. The modification of RNA aptamer at 2’OH with 2’ fluoro, 2’amino, and 2’O-methoxy has been found to greatly increase resistance to nuclease environment within the bloodstream. Systematic evolution of ligands by exponential enrichment (SELEX) technology is a method used to isolate aptamers, with high affinity and specificity to the target, and having a dissociation constant in the low nanomolar to picomolar level[14],[15]. Aptamers have more advantages than antibodies regarding production, stability, binding affinity, small size, and less immunogenicity[16]. Previous aptamer development for malaria treatment is focused on disruption of rosette formation and heme metabolism[13]. Unfortunately, an aptamer targeting the host endothelial surface receptors to block or reverse the binding of IEs has not yet been explored.

The present work describes the use of a SELEX strategy to isolate RNA aptamers against the recombinant (rhCD36). The isolated potential RNA aptamer candidates could be further studied on their potential to inhibit and reverse the IEs binding to human endothelial surface CD36 protein.

  2. Materials and methods Top

2.1. Preparation of RNA library

The aptamer library (5’-GCGAATTCTGGGGCGATATATCC(N40) CCAAATAGCCGAATTCGCACG-3’), forward primer (5’-GCTA ATACGACTCACTATAGGCGAATTCTGGGGCGATATATCC-3’) and reverse primer (5’-CGTGCGAATTCGGCTATTTGG-3’) were synthesized by the Integrated DNA Technologies (United Kingdom) and PAGE purified. The T7 promoter region (underlined) and N40 denotes 40 consecutive A, G, C or T nucleotides at 1:1:1:1 ratio, thus providing for approximately 1015 random RNA library. Amplification of ssDNA aptamer library was done by mixed 1×Green GoTaq® Flexi Buffer, 1.5 mM MgCl2, 0.2 mM dNTPs, 0.5 μM forward primer, 0.5 μM reverse primer, 10 nM oligonucleotide and 5 units GoTaq® DNA Polymerase (Promega, United States) in 100 μL total reaction. Then, the reaction mixture was amplified using the following condition; pre-denaturation at 94 °C for 5 min, 5 cycles of denaturation at 94 °C for 30 s, annealing at 61°C for 30 s, elongation at 72°C for 30 s and post-elongation at 72°C for 7 min using Veriti® Therma Cycler (Thermo Fisher Scientific, US). The PCR product was ethanol precipitated and followed by transcription using AmpliScribe™ T7-Flash™ Transcription Kit (Epicentre Biotechnology, USA) according to the manufacturer’s protocol. Briefly, the transcription reaction was performed in 20 μL containing 1× of AmpliScribe™ T7-Flash™ Reaction Buffer, 9 mM of each ATP, CTP, GTP, and UTP, 20 units of RiboGuard RNase Inhibitor (40 U/μL), 1 μg of dsDNA sample and 20 units of AmpliScribe™ T7-Flash™ Enzyme Solution (10 U/μL). The reaction mixture was mixed and incubated at 37° C for 16 h (overnight) using AccuBlock Digital Dry Bath (Labnet International. Inc, Edison). On the following day, the reaction tube was added with 10 units of RNase-free DNase I and incubated again at 37°C for 20 min to remove the dsDNA that existed in the reaction mixture. The reaction mixture was purified by using phenol: chloroform (24:1) and precipitated by using ethanol.

2.2. Selection of RNA aptamer

SELEX experiment was performed using nitrocellulose membrane filter immobilization methods as described by Hall and colleagues[17]. Briefly, the total RNA pool was resuspended in 100 μL of 1°SELEX buffer (SB) (137 mM of NaCl, 2.7 mM of KCl, 10 mM of Na2HPO4 and 1.7 mM of KH2PO4, pH 7.4) and heated at 95°C for 2 min, followed by cold down at room temperature for 10 min to allow RNA refolded into the secondary structure. The RNA pool was subjected for negative selection by filtering through pre-equilibrated (in 1×SB) 0.45 μm nitrocellulose membrane filter (Millipore, USA) using ‘pop-top’ filter holder (Whatman, United Kingdom) to remove RNA/filter bound and collect the filtrate (non-filter binder). The negative selection was performed in every SELEX cycle. During first SELEX cycle, 700 nM of the filtrate was mixed 20 nM of yeast tRNA (Life Technologies, USA) and 650 nM rhCD36 protein (R&D system, United Kingdom) in 100 μL of 1xSB. Then, the reaction was incubated at room temperature for 30 min on MACSmix tube rotator (MiltenyiBiotec, Germany). The binding reaction was filtered through a pre-wet 0.45 μm nitrocellulose membrane filter. After washed with 2× volume of 1× SB, the nitrocellulose membrane filter was collected and transferred to a new 1.5 mL tube. The RNAs/protein complex that bound to nitrocellulose membrane filter was recovered by heating in 400 μL elution buffer containing 7 M Urea at 95°C for 5 min. The bound RNA was transferred into a new centrifuge tube followed by ethanol precipitation. The stringency of selection was done by reducing the protein concentration and the incubation time and increasing the yeast tRNA concentration. Addition of nucleic acid competitors such as yeast tRNA or spermidine has the advantage to reduce nonspecific binding of RNA library[18]. The SELEX experiments were performed until significant binding affinity towards the target molecule was observed. In this study, the amount of the nucleic acid retained as protein/RNA complexes on the surface of the nitrocellulose membrane filter were collected, reverse transcribed, amplified and gel electrophoresed before quantified using the ImageJ software[19].

2.3. Amplification of bound aptamer by RT-PCR

The recovered RNA was reverse-transcribed using reverse transcriptase (Promega, United State). The reaction mixture was conducted in 20 μL containing l×AMV reverse transcriptase, 0.5 μM of reverse primer, 1.0 mM of dNTPs and 20 units of AMV reverse transcriptase. The reaction mixture was firstly heated at 95°C for 2 min following cold down at room temperature for 10 min before dNTPs and enzyme were added. The reaction mixture was incubated at 42°C for 1 h. The first cDNA was amplified using PCR. The PCR cycle was optimized until the band of the right size appeared. The PCR product was confirmed by electrophoresis of 10 μL of the PCR mixture on 2% (w/v) agarose gel and visualized using GeneFlash image analyzer (Syngene, USA). The PCR product was ethanol precipitated and directly used for the in vitro transcription using the AmpliScribe T7-Flash Transcription Kit. The isolated RNA was used for the next cycle of SELEX selection.

2.4. Cloning and sequencing

After 13 rounds of selection, the RT-PCR product was cloned using the TOPO TA Cloning Kit (Invitrogen, USA). There were 66 colonies that were randomly selected from positive individual colonies by using blue and white screening and plasmids were extracted using DNA-spin Plasmid DNA Extraction Kit (iNtRON Biotechnology, Korea). The plasmid was sent for sequencing using M13 forward primer (5’-GTAAAACGACGGCCAGT-3’) and M13 reverse primer (5’-CAGGAAACAGCTATGAC-3’) (First BASE Laboratories Sdn Bhd, Malaysia).

2.5. Prediction of secondary structure

The sequencing results were analyzed and aligned using ClustaW Multiple Alignment provided by MEGA 6 software (Version 6.06). The RNAs were clustered based on sequence alignment and phylogenetic analysis (Neighbor-Joining Tree) provided by MEGA 6 software. The online mfold software (http://unafold.rna. was used to predict the secondary structures of RNA aptamers. The secondary structure is essential to predict the binding region of the aptamer with the target and its functionality. Therefore, the analysis was conducted using free online mfold software[20]. The predicted secondary structures that have the lowest Gibbs free energy were selected because it will form the most stable secondary structure predicted by the software[20]. The aptamer sequences were further analyzed their motif consensus using online MEME (Multiple Em for Motif Elicitation) software at http:// (Version 4.12.0). The E-value generated by the software with the value less than or equal to 0.05 indicates that the motif identified is significant[21]. The online QGRS Mapper ( was used to determine the G-quadruplex sequence of isolated RNA aptamer[22].

2.6. Electrophoretic mobility shift assay (EMSA)

The plasmids of the selected aptamer cluster were amplified by using specific primers as described above. The products were gel purified and used for in vitro transcription. Binding activity of selected aptamers and target protein-ligand was analyzed using fluorescence-based EMSA. This method uses two fluorescent dyes for detection, which are PeqGreen nucleic acid stain for RNA and SYPRO Ruby protein gel stain for protein. The binding assay was prepared by mixing 1 μM of selected RNA aptamer clusters (initially denatured at 95°C for 2 min in 1xSB before cooling to room temperature for 10 min) with 0.5 μM of rhCD36. The binding reactions were then incubated at room temperature for 20 min. Following the addition of 4 μL of 6xloading dye, the RNA and protein complexes were separated from the free RNA using 5% (v/v) native PAGE. After electrophoresis at 140 V for 40 min in 0.5×Tris- Borate-EDTA (TBE) buffer, the gel was stained with a 1×PeqGreen nucleic acid stain in 0.5×TBE buffer for 10 min. After that, the gel was stained with SYPRO Ruby Protein stain for 3 h followed by destaining using destaining solution [10% (v/v) 2-Isopropanol, 10% (v/v) glacial acetic acid] for 1 h. The gel was visualized using a standard 300 nm UV transilluminator for both stains.

2.7. Semi-quantitative RT-PCR (sqRT-PCR) assay

The semi-quantitative RT-PCR assay was carried out to further evaluate the binding ability of isolated RNA aptamer to the target protein. This assay was used to quantify the bound RNA to the target by measuring the intensity of reverse transcribed and amplified product from the gel electrophoresis using ImageJ software[19]. The binding reactions were prepared at a fixed concentration of RNA and rhCD36 protein for all selected RNA aptamers. RNA was diluted to a final concentration of 1 μM and initially denatured at 95°C for 2 min in 1×SB before cooling down to room temperature for 10 min. Subsequently, 0.5 μM of target proteins were added to corresponding tubes, and the reaction mixture was incubated at room temperature for 20 min. The RNA/protein complexes were separated using a nitrocellulose membrane filter. The eluted (bound) and filtrate (unbound) RNA/s were reverse transcribed and amplified. The amplification was performed at the 13th PCR cycle to avoid over-amplification of PCR products. The PCR product was confirmed by electrophoresis of 10 μL of the PCR mixture on 2% (w/v) agarose gel and visualized using an image analyzer. The band’s intensity of amplified products was measured using ImageJ software (Version 1.46r).

2.8. Determination of dissociation constant value

The modified agarose-based EMSA assay was used to determine dissociation constant (Kd) as described by Ream et al[23]. Briefly, the series of rhCD36 protein concentrations were used from 0, 25, 50, 100, 250 and 500 nM and fixed concentration of RNA (100 nM). The complexes were resolved on 1% (w/v) agarose gel containing 1×Diamond Gel Stain (Promega, United States) (in 0.5×TBE buffer), ran 150 volts for 10 min. The gel was visualized using a standard 300 nm UV transilluminator. The signal intensity of bound and free RNA was measured using the ImageJ program and calculated Equation 1. The equilibrium dissociation (Kd) constant was calculated using nonlinear, curve-fitting algorithm [Y = Bmax*X / (Kd + X)] (GraphPad Prism 7.0 Software). Equation 1: Percentage of fraction bound aptamer.

2.9. The specificity of isolated RNA aptamers

The RNA aptamer with a polyA tail was synthesized using a TGCGAATTCGGCTATTTGG-3’). Enzyme-linked oligonucleotide assays (ELONA) based binding assay was conducted in order to determine the binding specificity of rhCD36-isolated RNA aptamer cluster to bovine serum albumin (BSA) protein and ICAM-1 protein. Briefly, 1 μg/mL rhCD36, recombinant human intercellular adhesion molecule 1 (rhICAM-1) (R&D system, United Kingdom) and BSA (Sigma-Aldrich, United States) (diluted in bicarbonate/carbonate coating buffer, pH 9.6) were coated on 96 wells microplate at 37°C for 1 h. The plate was then washed with ELISA washing buffer (1xPBS, pH 7.4 containing 0.02% (v/v) Tween-20) three times using ELISA plate washer. Afterward, the plate was blocked with 2% (w/v) non-fat milk (Oxoid, United Kingdom) diluted in ELISA washing buffer. After that, the mixture of 100 nM of polyA tail RNA aptamer and 100 nM of biotinylated polydT (5’-/BiosG/ TTTTTTTTTTTTTTTTTTTT3’) (previously diluted in 1xSB denatured at 95°C for 2 min and reannealed at room temperature for 10 min) was added and the plate was incubated at 37°C for 1 h. The bound RNA was detected using HRP-conjugated streptavidin (Thermo Fisher Scientific, United States) (diluted 1:10 000 in 1× PBS, pH 7.4). After washing five times, the color was developed by adding TMB substrate and reaction was stopped using 2N HCl. The absorbance was read at 450 nm using SpectraMax M5 (Molecular Device, United States).

2.10. Competitive inhibition assay

Ten μg/mL of avidin (Thermo Fisher Scientific, United States) (diluted in bicarbonate/carbonate coating buffer, pH 7.4) was coated on 96 wells microplate at 4°C overnight. The plate was then washed with ELISA washing buffer (1×PBS, pH 7.4 containing 0.02% (v/ v) Tween-20) three times using ELISA plate washer. The plate was then blocked with 2% (w/v) non-fat milk diluted in ELISA washing buffer. Subsequently, the mixture of 100 nM of polyA tail RNA aptamer and 100 nM of biotinylated polydT (5’-/BiosG/ TTTTTTTTTTTTTTTTTTTT3 ' ) (previously diluted in 1×SB, denatured at 95°C for 2 min and reannealed at room temperature for 10 min) was added and the plate was incubated at 37°C for 1 h. After washing, the mixture of 1 μg/mL of rhCD36 and 1 μg/mL of mAb FA6-152 (Abcam, United Kingdom) was added and the plate was incubated at 37°C for 1 h. As the control, the polyA tail RNA aptamer was allowed to bind with rhCD36 in the absence of mAb FA6-152. The bound rhCD36 was detected using HRP-conjugated anti-human IgG (Thermo Fisher Scientific, United States) (diluted 1:100 000 in 1×PBS, pH 7.4). After washing five times, the color was developed by adding TMB substrate and reaction was stopped using 2N HCl. The absorbance was read at 450 nm using SpectraMax M5 (Molecular Device, United States).

2.11 Statistical analysis

Data analysis was done using Two-way ANOVA and unpaired t-test using GraphPad Prism (Version 7.0). Error bar indicates standard deviation(±SD) of the triplicate experiments and P-value less than 0.05 was considered statistically significant.

  3. Results Top

3.1. Selection of RNA aptamers against human rhCD36

Isolation of RNA aptamers specifically bind to recombinant human endothelial surface receptor CD36 protein was conducted using nitrocellulose membrane-based SELEX. PCR initially prepared the random library and followed by in vitro transcription. The stringency of selection was increased by reducing the concentration of rhCD36 protein and incubation time while increasing the concentration of yeast tRNA. Introducing the BSA protein after five rounds of the SELEX cycle might increase the chances to isolate high affinity and specificity of RNA aptamer to rhCD36 protein. In our study, the counter SELEX with BSA protein was conducted at cycle 8th, 10th, and 12th. [Table 1] shows the SELEX condition in this study and the optimized number of PCR cycles.
Table 1: Condition for SELEX experiment of RNA aptamers against rhCD36 and optimized number of PCR cycles.

Click here to view

3.2. Evaluation of SELEX enrichment

SELEX enrichment of RNA pool binding with rhCD36 was evaluated using sqRT-PCR at 1st, 3rd, 7th, 9th, 11th, and 13th SELEX cycle. [Figure 1] shows the band intensity of amplified cDNA of collected RNA from selected SELEX cycle. The equal concentration of RNA from selected SELEX cycles and rhCD36 protein were used in the binding assay. The complexes were separated using nitrocellulose and reverse transcribed the bound RNA. Amplification of cDNA was limited to thirteen PCR cycle, and PCR products were loaded into the gel followed by ImageJ quantification [Supplementary Figure 1] [Additional file 1]. The increasing recovery of bound RNA was observed starting at 3rd SELEX cycle and reached a maximum at 13th SELEX cycle.
Figure 1: The SELEX enrichment of RNA pool targeting protein using semiquantitative PCR. The image represents the percentage of RNA that bound to rhCD36 from selected SELEX cycles. The ImageJ quantification shows that the RNA pool from 13th SELEX cycle demonstrated the highest percentage of bound RNA to rhCD36 protein.

Click here to view

3.3. Sequence and structure analysis of aptamer

After the 13th SELEX cycle, the amplified cDNA was cloned and sequenced. Sixty-six clones of rhCD36 aptamer candidates were successfully sequenced, aligned and classified into several groups based on their sequence frequency and homology. From [Table 2], eleven clusters of aptamer sequence have been identified targeting rhCD36. Cluster 1 (15.15%), C1_RC60 showed the highest frequency of occurrence, which was 13.64%. However, C1_RC17 that had a single frequency of occurrence was grouped similar to C1_RC60 as it only showed a single mutation at position (43C>U) within the randomized region. There were four sequences of cluster 2 (7.58%) (C2_RC15, C2_RC04, C2_RC71, and C2_ RC51) with single mutation at position (48U>C for C2_RC04 and 27A>G for C2_RC71) compared with C2_RC15, both present within randomized region. However, C2_RC15 showed a single base deletion at position 39delG. In cluster 3 (6.06%), there was one mutation site within randomized region identified in C3_RC57 that occurred at position 50U>C compared with C3_RC25. Other isolated clusters were C4_RC26 (6.06%), C5_RC16 (4.55%), C6_RC46 (4.55%), C7_RC14 (3.03%), C8_RC11 (3.03%), C9_ RC52 (3.03%), C10_RC67 (3.03%) and C11_RC55 (3.03%). The aptamer clusters were also confirmed using Neighbor-Joining (NJ) tree [Supplementary Figure 2] [Additional file 2]. [Figure 2] shows the significant motif sequence found by MEME software was located within a randomized region (green) instead of two constant regions (red and light blue). The motif sequence of rhCD36 aptamer clusters was discovered at ‘[U/G/A][U/G]GG[A/U/G]GG[A/U/G]GG[A/U][G/ U]G[A/U/G]’. Besides, the sequences of isolated RNA aptamers were analyzed using online QGRS Mapper in order to determine the location of the G-quadruplex sequence. We also discovered that all the isolated RNA aptamers possess G-quadruplex sequence within the randomized region [Supplementary Figure 3] [Additional file 3].
Figure 2: The sequence analysis of isolated aptamers by online MEME software. (A) Significant motif sequence was found within the randomized region (green) instead of constant sequence (red and light blue). The purple and yellow rectangle show other motifs found by software with insignificant E-value. (B) The thickness of a letter indicates the probability of base insertions between aptamer clusters.

Click here to view
{Figure 1
Table 2: The random sequence of isolated RNA aptamer clusters targeted rhCD36 protein according to the frequency number and degree of sequence similarities.

Click here to view

3.4. Secondary structure prediction

The predicted secondary structure of 6 selected RNA aptamers represented for cluster 1-6 as showed in [Figure 3]. Most of the aptamer clusters had a hairpin loop structure connected to the multi-branch or core loop. Also, there were also several clusters having bulge loop, and interior loop connected to core loop structure. The secondary structure analysis showed most of the motif sequences were located within a single-stranded of hairpin loop structure and double-stranded of the hairpin stem. Therefore, it can be hypothesized that the discovered motif sequence might play a role in the aptamer-target interactions as most of the motifs were located between hairpin loop structure and base pairs region within the randomized sequence.
Figure 3: Secondary structure prediction of RNA aptamers that bind to rhCD36 protein using the online mfold software. The arrows indicate the discovered motif sequence.

Click here to view

3.5. Binding assay of isolated RNA aptamer

The binding assay was performed to confirm whether isolated RNA aptamers can bind to rhCD36 proteins. In this study, EMSA and sqRT-PCR were performed to confirm the aptamer binding. The interaction between isolated RNA aptamer with the target rhCD36 protein was observed on EMSA [Figure 4] and [Supplementary Figure 4] [Additional file 4]. Unfortunately, some of the complexes retained in the well, and this might be due to complexes aggregation. Statistical analysis of sqRT-PCR assay using unpaired t-test showed significant binding of isolated RNA aptamer compared with the RNA library (initial library) in nine aptamer candidates for targeted rhCD36 [Supplementary Figure 5] [Additional file 5]. The result obtained from the sqRT-PCR assay was in concordance with the EMSA-based binding assay. Therefore, three aptamer candidates (RC60, RC04, and RC25) targeted rhCD36 proteins were selected to represent cluster 1 to cluster 3 and further evaluated for their binding affinity as they showed higher binding to rhCD36 compared to other aptamer clusters.
Figure 4: Electrophoretic mobility shift assay (EMSA) of selected RNA aptamers targeting rhCD36.

Click here to view

3.6. Dissociation constant (Kd) value of isolated aptamers

The Kd value was determined using a modified agarose-based EMSA assay. The concentration of RNA aptamer was fixed while the concentration of protein was varied (25 to 500 nM). Thus, the estimated Kd values for RC04, RC60 and RC25 were (17.49±4.16) nM, (21.90±2.40) nM, and (23.07±3.49) nM, respectively [Supplementary Figure 6] [Additional file 6], [Supplementary Figure 7] [Additional file 7], [Supplementary Figure 8] [Additional file 8]. The results demonstrated that the binding affinity of RC04 had the highest binding affinity followed by RC60 and RC25.

3.7. The specificity of isolated RNA aptamer

In this study, ELONA was used to rapidly assess the relative binding specificity of several isolated RNA aptamers to the rhCD36, BSA and rhICAM-1. Previously, RC04 showed the highest binding affinity compared to RC60 and RC25. However, in the binding specificity analysis, RC04 candidate showed cross-reactivity to BSA and rhICAM-1 protein while RC60 and RC25 showed significantly no cross-reactivity to BSA and rhICAM-1 protein [Figure 5]. Negative binding assay of RC60 and RC25 to nitrocellulose membrane filter also showed no binding [Supplementary Figure 9] [Additional file 9].
Figure 5: Binding specificity of isolated RNA aptamers to rhCD36. The binding specificity of isolated RNA aptamers was estimated by comparing the binding between rhCD36, BSA and rhICAM-1 protein. Data analysis was done using Two-way ANOVA (Multiple comparisons test) where ***P < 0.001. Error bar indicates ±SD of the triplicate experiment

Click here to view

3.8. The binding orientation of isolated RNA aptamer to rhCD36

Aptamer competitive binding assay was designed in order to determine the binding epitope of RNA aptamer to rhCD36 protein. Significant reduction of absorbance was observed for RC60, RC04, and RC25 in the presence of mAb FA6-152 [Figure 6]. The reduction of absorbance might be caused by the competitive inhibition of mAb FA6-152 at the same epitope site of RNA aptamer binding to rhCD36. Other RNA aptamers (RC26, RC16, and RC46) showed a reduction on the signal intensity with no significant difference. The most convincing hypothesis is that the immobilized RNA aptamer increases the captured protein when there are more free epitopes compared to protein immobilization.
Figure 6: Enzyme-linked oligonucleotide assays (ELONA) for binding site confirmation. The polyA tail RNA aptamer (100 nM) was immobilized onto ELISA microplate. The competitive inhibition of aptamer binding was measured by quantifying the absorbance of bound rhCD36 in the presence of mAb FA6-152. The data obtained were compared with control (in the absence of mAb FA6-152). Data analysis was done using Two-way ANOVA (Multiple comparisons test) where *P < 0.05 and **P < 0.01. Error bar indicates ±SD of the triplicate experiment.

Click here to view

  4. Discussion Top

Previous studies have reported that there is an association between cytoadherence of IEs to the host endothelial surface receptor protein and parasite burden in the blood vessel. The parasite surface protein (PfEMP1) was reported to bind to several hosts endothelial surface protein including CD36 which linked to uncomplicated malaria, but it also can indirectly support the cytoadhesion of the IEs to the microvascular brain which leads to the development of severe malaria. Thus, this study was proposed to isolate and characterize a group of novel RNA aptamers that bind to CD36 protein as a potential malaria adjunct therapy. Nitrocellulose membrane-based SELEX was performed to isolate the RNA aptamers that specifically bind to recombinant human endothelial surface receptor CD36 protein. In each round of SELEX cycle, the amounts of specific binders are present at low copy number compared to a non-specific binder which then requires optimum PCR method to increase amplification efficiency to avoid amplification bias and minimize the loss of potential high-affinity binders. Over-amplification of PCR product hypothetically tends to cause DNA mispairing as oligonucleotide template containing an abundance of different sequences of the randomized region that can anneal with each other, leading to more round of SELEX cycle[24],[25]. Therefore, to avoid PCR over-amplification throughout this SELEX experiment, the number of PCR cycles was optimized in every SELEX cycle in the increment of two PCR cycles until the expected band size appears[26]. Then, the scale-up of the DNA library for the next SELEX cycle was carried out using an optimized number of PCR cycle, followed by ethanol precipitation before subjected to in vitro transcription. The SELEX enrichment analysis showed that the increased amount of specific RNA bound to rhCD36 protein starting from round three and it reached maximum binding at around thirteen. The results obtained were in concordance with another study when they found increasing RNA or DNA recovery after 3 to 5 rounds of the SELEX[27].

The sequence analysis revealed a G-quadruplex sequence between 11 and 12 bases was found within all the isolated RNA aptamers. As predicted, the presence of potassium ion (K +) at 2.7 mM in the selection buffer has increased the chances of G-quadruplex formation that plays an important role in aptamer binding[28]. The aptamer that contains G-quadruplex was reported to be more thermodynamically and chemically stable, which increases stability within the nuclease environment of the bloodstream[29]. Some studies reported that the contact of the aptamer and the target occurs in the single-stranded region, but it also can happen in the base-paired region[30]. Nevertheless, further study is needed to confirm this hypothesis using RNase/DNase footprinting assay and in silico Aptamer Docking[31],[32].

There were three aptamer candidates (RC60, RC04, and RC25) that show a significant binding with the rhCD36 protein. The RC04 aptamer shows the highest binding affinity compared to RC60 and RC25 even though it shows only one sequence frequency in the cluster. This finding correlates with another study which reported that the aptamer sequence frequency is a poor predictor of aptamer affinity. Thus, evaluation of the cycle-to-cycle enrichment of aptamer is a better predictor of binding compared to frequency determination in the final pool of the SELEX cycle[33]. Assuming the added target proteins were active, and the stoichiometric of the RNA aptamer and protein formation was at a similar ratio, the Kd was estimated[34].

The binding specificity analysis of RC60, RC04 and RC25 showed that RC60 and RC25 demonstrated significantly no cross-reactivity. Meanwhile, RC04 showed cross-reactivity against BSA and rhICAM-1 protein. The ICAM-1 protein is an immunoglobulin (Ig) superfamily member that is expressed on the surface of endothelial membrane cells and leucocytes encompassing five extracellular Ig-like domain (D1-D5), a transmembrane domain and a cytoplasmic domain[35]. There is a relationship of high cytoadherence and cerebral malaria that associated with adhesion of IE to ICAM-1 protein[36],[37]. Furthermore, the binding orientation analysis showed that RC60, RC04, and RC25 have the same binding epitope with mAb FA6-152. This monoclonal antibody was known to be able to inhibit the binding of ox-LDL and recognition of apoptotic neutrophils by residues 155-183[38]. Moreover, several studies have reported that mAb FA6-152 also efficiently inhibited the binding of infected erythrocyte to CD36[12].

In conclusion, this study successfully isolated and characterized two potential RNA aptamers (RC60 and RC25) that have affinity against rhCD36. These isolated RNA aptamer candidates could be further investigated on its future use as anti-cytoadherence for severe malaria adjunct therapy.

Conflict of interest statement

The authors declare that there is no conflict of interest.


This work is supported by an Exploratory Research Grant Scheme (ERGS; 203/0ΡΡΜ/6730106) and Higher Institution Centre of Excellent (HICoE; 311/ΏΡΡΜ/4401005) from the Malaysian, Ministry of Higher Education and UniversitiSains Malaysia Fellowship Scheme.

Authors’ contributions

NAANK, JNAS, NFMZ and KMFM performed the experiment. NAANK and KMFM performed the data and statistical analysis, manuscript preparation and editing. Both NAANK and KMFM authors contributed to the final version of the manuscript. KMFM designed and supervised the project. All authors approved the final version for submission.

  References Top

World Health Organization. World malaria report 2016. Geneva Switzerland; 2016, p.186.  Back to cited text no. 1
World Health Organization. Fact sheet: World malaria report 2016. [Online] Available from: [Accessed on 10th September 2017].  Back to cited text no. 2
Cowman AF, Healer J, Marapana D, Marsh K. Malaria: Biology and disease. Cell 2016; 167(3): 610-624.  Back to cited text no. 3
Wassmer SC, Grau GE. Severe malaria: What’s new on the pathogenesis front? Int JParasitol 2017; 47(2-3): 145-152.  Back to cited text no. 4
Deroost K, Pham TT, Opdenakker G, Van den Steen PE. The immunological balance between host and parasite in malaria. FEMS Microbiol Rev 2016; 40(2): 208-257.  Back to cited text no. 5
Prabhu Das MR, Baldwin CL, Bollyky PL, Bowdish DME, Drickamer K, Febbraio M, et al. A consensus definitive classification of scavenger receptors and their roles in health and disease. J Immunol 2017; 198(10): 3775-3789.  Back to cited text no. 6
Gong QW, Liao MF, Liu L, Xiong XY, Zhang Q, Zhong Q, et al. CD36 gene polymorphisms are associated with intracerebral hemorrhage susceptibility in a hanchinese population. Biomed Res Int 2017; 2017: 7.  Back to cited text no. 7
Gowda DC, Wu XZ. Parasite recognition and signaling mechanisms in innate immune responses to malaria. Front Immunol 2018; 9: 3006.  Back to cited text no. 8
Tuikue Ndam N, Moussiliou A, Lavstsen T, Kamaliddin C, Jensen ATR, Mama A, et al. Parasites causing cerebral falciparum malaria bindmultiple endothelial receptors and express EPCR and ICAM-1-binding PfEMP1. J Infect Dis 2017; 215(12): 1918-1925.  Back to cited text no. 9
[Lagassé HA, Anidi IU, Craig JM, Limjunyawong N, Poupore AK, Mitzner W, et al. Recruited monocytes modulate malaria-induced lung injury through CD36-mediated clearance of sequestered infected erythrocytes. JLeukoc Biol 2016; 99(5): 659-671.  Back to cited text no. 10
Hsieh Turner U Bolla JR, Robinson CV, Lavstsen T, Higgins MK. The structural basis for CD36 binding by the malaria parasite. Nat Commun 2016; 7: 12837.  Back to cited text no. 11
Mustaffa KMF, Storm J, Whittaker M, Szestak T, Craig AG. In vitro inhibition and reversal of Plasmodium falciparum cytoadherence to endothelium by monoclonal antibodies to ICAM-1 and CD36. Malar J 2017; 16(1): 279.  Back to cited text no. 12
Nik Kamarudin NAA, Mohammed NA, Mustaffa KMF. Aptamer technology: Adjunct therapy for malaria. Biomedicines 2017; 5(1): E1.  Back to cited text no. 13
Zhou G, Wilson G, Hebbard Duan W, Liddle C, George J, et al. Aptamers: A promising chemical antibody for cancer therapy. Oncotarget 2016; 7(12): 13446-13463.  Back to cited text no. 14
Sun HG, Zhu X, Lu PY, Rosato RR, Tan W, Zu ΥΕ. Oligonucleotide aptamers: New tools for targeted cancer therapy. Mol Ther Nucleic Acids 2014; 3: e182. Doi: 10.1038/mtna.2014.32.  Back to cited text no. 15
Zhou JH, Rossi J. Aptamers as targeted therapeutics: Current potential and challenges. Nat Rev Drug Discov 2017; 16(6): 440.  Back to cited text no. 16
Hall B, Arshad S, Seo K, Bowman C, Corley M, Jhaveri SD, et al. In vitro selection of RNA aptamers to a protein target by filter immobilization. Curr Protoc Nucleic Acid Chem 2010; 89(1). Doi: 10.1002/0471142700. nc0903s40.  Back to cited text no. 17
Sevostyanova A, Groisman EA. An RNA motif advances transcription by preventing Rho-dependent termination. Proc Natl Acad Sci USA 2015; 112(50): E6835-E6843.  Back to cited text no. 18
Kwon HM, Lee KH, Han BW, Han MR, Kim DH, Kim DE. An RNA aptamer that specifically binds to the glycosylated hemagglutinin of avian influenza virus and suppresses viral infection in cells. PLoS One 2014; 9(5): e97574. Doi: 10.1371/journal.pone.0097574.  Back to cited text no. 19
Zuker M. Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res 2003; 31(13): 3406-3415.  Back to cited text no. 20
Bailey TL, Boden M, Buske FA, Frith M, Grant CE, Clementi L, et al. MEME SUITE: Tools for motif discovery and searching. Nucleic Acids Res 2009; 37: W202-W208.  Back to cited text no. 21
Kikin O, D’Antonio Ε, Bagga PS. QGRS Mapper: A web-based server for predicting G-quadruplexes in nucleotide sequences. Nucleic Acids Res 2006; 34: W676-W682.  Back to cited text no. 22
Ream JA, Lewis LK, Lewis KA. Rapid agarose gel electrophoretic mobility shift assay for quantitating protein: RNA interactions. Anal Biochem 2016; 511: 36-41.  Back to cited text no. 23
Liu LJ, Chen Y, Wang W, Chen C, Gao MH, Zhang XQ, et al. Screening and identification of aptamers against pulmonary surfactant protein A. Chinese J Anal Chem 2013; 41(11): 1659-1663.  Back to cited text no. 24
Tolle F, Wilke J, Wengel J, Mayer G. By-product formation in repetitive PCR amplification of DNA libraries during SELEX. PLoS One 2014; 9(12): e114693. Doi: 10.1371/journal.pone.0114693.  Back to cited text no. 25
Marimuthu C. Generation and characterization of RNA aptamer againtrHuEPO-α by SELEX technology. Advanced Medical and Dental Institute: Universiti Sains Malaysia; 2013.  Back to cited text no. 26
Zhang YJ, You YD, Xia ZW, Han XY, Tian YP, Zhou ND. Graphene oxide-based selection and identification of ofloxacin-specific single- stranded DNA aptamers. RSC Adv 2016; 6(101): 99540-99545.  Back to cited text no. 27
Szameit K, Berg K, Kruspe S, Valentini E, Magbanua E, Kwiatkowski M, et al. Structure and target interaction of a G-quadruplex RNA-aptamer. RNA Biol 2016; 13(10): 973-987.  Back to cited text no. 28
Kwok CK, Merrick CJ. G-quadruplexes: Prediction, characterization, and biological application. Trends Biotechnol 2017; 35(10): 997-1013.  Back to cited text no. 29
Choi SJ, Ban C. Crystal structure of a DNA aptamer bound to PÆDH elucidates novel single-stranded DNA structural elements for folding and recognition. Sci Rep 2016; 6: 34998.  Back to cited text no. 30
Ruscito A, DeRosa MC. Small-molecule binding aptamers: Selection strategies, characterization, and applications. Front Chem 2016; 4: 14.  Back to cited text no. 31
Ahirwar R, Nahar S, Aggarwal S, Ramachandran S, Maiti S, Nahar P. In silico selection of an aptamer to estrogen receptor alpha using computational docking employing estrogen response elements as aptamer-alike molecules. Sci Rep 2016; 6: 21285.  Back to cited text no. 32
Hoinka J, Berezhnoy A, Sauna ZE, Gilboa E, Przytycka TM. AptaCluster - A method to cluster HT-SELEX aptamer pools and lessons from its application. Res Comput Mol Biol 2014; 8394: 115-128. Doi:10.1007/978- 3-319-05269-4_9.  Back to cited text no. 33
Daniel C, Roupioz Y, Gasparutto D, Livache T, Buhot A. Solution- phase vssurface-phase aptamer-protein affinity from a label-free kinetic biosensor. PLoS One 2013; 8(9): e75419. Doi:10.1371/journal. pone.0075419.  Back to cited text no. 34
Bengtsson A, Joergensen L, Rask TS, Olsen RW, Andersen MA, Turner L, et al. A novel domain cassette identifies Plasmodium falciparum PfEMP1 proteins binding ICAM-1 and is a target of cross-reactive, adhesion-inhibitory antibodies. J Immunol 2013; 190(1): 240-249.  Back to cited text no. 35
Ochola LB, Siddondo BR, Ocholla H, Nkya S, Kimani EN, Williams TN, et al. Specific receptor usage in Plasmodium falciparum cytoadherence is associated with disease outcome. PLoS One 2011; 6(3): e14741. Doi: 10.1371/journal.pone.0014741.  Back to cited text no. 36
Grau GE, Craig AG. Cerebral malaria pathogenesis: Revisiting parasite and host contributions. Future Microbiol 2012; 7(2): 291-302.  Back to cited text no. 37
Kar NS, Ashraf MZ, Valiyaveettil M, Podrez EA. Mapping and characterization of the binding site for specific oxidized phospholipids and oxidized low density lipoprotein of scavenger receptor CD36. J Biol Chem 2008; 283(13): 8765-8771.  Back to cited text no. 38


  [Figure 1], [Figure 2], [Figure 3], [Figure 4], [Figure 5], [Figure 6]

  [Table 1], [Table 2]

This article has been cited by
1 Aptamers isolated against mosquito-borne pathogens
Tholasi Nadhan Navien,Tzi Shien Yeoh,Andrew Anna,Thean-Hock Tang,Marimuthu Citartan
World Journal of Microbiology and Biotechnology. 2021; 37(8)
[Pubmed] | [DOI]


    Similar in PUBMED
   Search Pubmed for
   Search in Google Scholar for
 Related articles
    Access Statistics
    Email Alert *
    Add to My List *
* Registration required (free)  

  2. Materials and...
  In this article
1. Introduction
3. Results
4. Discussion
Article Figures
Article Tables

 Article Access Statistics
    PDF Downloaded326    
    Comments [Add]    
    Cited by others 1    

Recommend this journal