- Department of Gastroenterology, The First Affiliated Hospital of Nanchang University, Jiangxi Institute of Gastroenterology and Hepatology, Jiangxi Clinical Research Center for Gastroenterology, Nanchang, China
Introduction: Aptamers are valuable for bioassays, but aptamer-target binding is susceptible to reaction conditions. In this study, we combined thermofluorimetric analysis (TFA) and molecular dynamics (MD) simulations to optimize aptamer-target binding, explore underlying mechanisms and select preferred aptamer.
Methods: Alpha-fetoprotein (AFP) aptamer AP273 (as the model) was incubated with AFP under various experimental conditions, and melting curves were measured in a real-time PCR system to select the optimal binding conditions. The intermolecular interactions of AP273-AFP were analysed by MD simulations with these conditions to reveal the underlying mechanisms. A comparative study between AP273 and control aptamer AP-L3-4 was performed to validate the value of combined TFA and MD simulation in selecting preferred aptamers.
Results: The optimal aptamer concentration and buffer system were easily determined from the dF/dT peak characteristics and the melting temperature (Tm) values on the melting curves of related TFA experiments, respectively. A high Tm value was found in TFA experiments performed in buffer systems with low metal ion strength. The molecular docking and MD simulation analyses revealed the underlying mechanisms of the TFA results, i.e., the binding force and stability of AP273 to AFP were affected by the number of binding sites, frequency and distance of hydrogen bonds, and binding free energies; these factors varied in different buffer and metal ion conditions. The comparative study showed that AP273 was superior to the homologous aptamer AP-L3-4.
Conclusion: Combining TFA and MD simulation is efficient for optimizing the reaction conditions, exploring underlying mechanisms, and selecting aptamers in aptamer-target bioassays.
1 Introduction
Nucleic acid aptamers, which are artificial ligands of biomolecules, are selected from synthetic random single-stranded oligonucleotide libraries through systematic evolution of ligands by exponential enrichment (SELEX) (Ellington and Szostak, 1990; Tuerk and Gold, 1990; Wu et al., 2021). Aptamers have similar binding properties to antibodies, but they are more stable and easier to prepare and modify than antibodies. Therefore, aptamers are an ideal tool for biomedical molecular detection and recognition (Zhu et al., 2015; Emrani et al., 2016; Li et al., 2019). However, the binding properties of aptamers to their targets are closely related to their structures (Davydova et al., 2020), and the structures of aptamers are sensitive to the conditions of the reaction system. Hence, the analysis of the three-dimensional structures of aptamers and their influencing factors is important for optimizing their binding properties and for selecting preferred aptamers.
Traditionally, the structural analysis of aptamers mainly depends on nuclear magnetic resonance (NMR) (Longhini et al., 2016; Cai et al., 2018) and X-ray crystallography (Ruigrok et al., 2012; Christina Schmidt et al., 2020), but universal application of these methods is not possible due to expensive instruments and complicated procedures. With the rapid development of artificial intelligence technologies, molecular structure prediction based on bioinformatics and molecular dynamics (MD) simulations is increasingly recognized and applied (Douaki et al., 2022; Poolsup et al., 2023), because structural changes in aptamers (before and after binding to their targets) can be rapidly analyzed without damaging the aptamers, the effects of metal ions in the reaction system on the interaction of aptamers and targets can be analyzed, and information on the stability, binding energy, affinity, and molecular mechanism of aptamer-target binding is provided (Hayashi et al., 2014; Hilder and Hodgkiss, 2017; Cui et al., 2020). Through MD simulations, La Penna et al. (La Penna and Chelli, 2018) found that Mg2+ and Na+ could attach to the aptamer surface to enhance the interaction between the amino acid residues of protein targets and the nucleotides of aptamers; thus, Mg2+ and Na+ facilitated the stability of the binding compounds. Using MD simulations, Vu et al. (Vu et al., 2018) identified the binding site of platelet-derived growth factor-B (PDGF-B) with its specific aptamer, and their binding force was mainly the electrostatic force between the positively charged amino acid of the target and the negatively charged phosphate backbone of the aptamer.
In general, it is necessary to select excellent aptamers from the dozens of aptamers generated by SELEX for subsequent applications (Darmostuk et al., 2015; Vu et al., 2018). Specificity and affinity analyses are traditional methods for selecting aptamers, but these methods are labor-intensive and often inflexible when used in complex assay settings. Thermofluorimetric analysis (TFA) is a method that can be used to evaluate aptamer-target binding and distinguish bound and free aptamers in a convenient and rapid manner. This method is based on the melting curve analysis in a polymerase chain reaction (PCR) system (Hu and Easley, 2017), and it can be used to optimize experimental conditions and select the optimal aptamer in complex reaction systems. TFA has been applied to quantify serum platelet-derived growth factor (PDGF), insulin, and prothrombin based on aptamers (Hu et al., 2015; Kim et al., 2015).
Binding between a target and an aptamer is primarily governed by the mutual adaptation of their conformations (the so-called “lock and key” relationship) (Hayashi et al., 2014), and intermolecular non-covalent bonds such as hydrogen bonds, electrostatic interaction, and van der Waals forces act as auxiliary factors in the binding between aptamer and target (Shigdar, 2019; Morozov et al., 2021), while there are also reports of aptamer and target being bound via covalent bonds (Tivon et al., 2021; Tabuchi et al., 2022). Metal ions can impact the tertiary structure of aptamers and their ability to bind to targets through their charge distribution and potential energy (Moccia et al., 2019). Therefore, it is necessary to optimize the metal ion concentrations in binding buffers to identify aptamers with a wider range of adaptability and high specificity. Conventional methods for optimization are generally based on multifactor, multilevel experiments, which are time-consuming and sometimes yield unsatisfactory results (Cheng et al., 2003; Zhou et al., 2015). Thus, a simple, fast, and effective method is needed to optimize the metal ions in the buffers.
The status of aptamer-target interactions can be simply and rapidly observed through TFA, and aptamer structures and their interactions with targets can be visually analyzed through MD simulation. Thus, we speculated that by combining TFA and MD simulations, the effects of buffers and metal ions on aptamer structures and their interactions with targets can be easily analyzed at both experimental and theoretical levels, which can help to select optimal aptamers and optimize experimental conditions. In the present study, the reported aptamer AP273 against alpha-fetoprotein (AFP) (Dong et al., 2015) and our screened aptamer AP-L3-4 against AFP-L3 were used as models for validation, by which we attempted to provide a simple and feasible method for prioritizing aptamers and optimizing experimental conditions.
2 Materials and methods
2.1 Determination of the optimal aptamer concentration by TFA
Aptamers AP273 (5′-GTG ACG CTC CTA ACG CTG ACT CAG GTG CAG TTC TCG ACT CGG TCT TGA TGT GGG TCC TGT CCG TCC GAA CCA ATC-3′) and AP-L3-4 (5′-ACC GAC CGT GCT GGA CTC TGT CGA AAG GAT ACT GAG TAT TGA GGG GCG TCA GGT GGA AGA GTA TGA GCG AGC GTT GCG-3′) were synthesized (Sangon Biotech (Shanghai) Co., Ltd.) and dissolved in ddH2O to prepare a 100 nM storage solution. AP273 is an ssDNA aptamer against AFP selected using CE-SELEX by Wu’s team (Dong et al., 2015); AP-L3-4 is an ssDNA aptamer against AFP-L3 selected using SELEX in our previous work.
The optimal concentration of aptamer AP273 for binding with AFP was determined by analysing the characteristics of the melting curves generated by different concentrations of aptamer reacting with a fixed concentration of AFP. Gradient concentrations of aptamer working solutions (1.25, 2.5, 5.0, 10, 20, 40, and 80 nM) were prepared by adding the typical buffer used in aptamer screening (HEPES-Na 20 mmol/L, NaCl 120 mmol/L, KCl 4 mmol/L, 2 mmol/L MgCl2, 1 mmol/L CaCl2, pH 7.35). Human recombinant AFP (Nearshore Protein Technology Ltd., Shanghai, China) was dissolved in ddH2O to prepare a solution of 1.45 nM (100 ng/mL). EvaGreen (20×) dye (Biotium, United States of America) was diluted to 8× working solution using ddH2O (EvaGreen is a nucleic acid fluorescent dye commonly used in real-time PCR (Shoute and Loppnow, 2018).
The aptamer working solution was denatured in a metal bath at 95°C for 3 min and then immediately placed in an ice bath for 3 min. Twenty microliters of the denatured solution was added to PCR tubes, followed by the addition of 5 µL of 8× Evagreen working solution and 5 µL of human recombinant AFP solution. This mixture was incubated for 30 min at room temperature and then placed in a StepOnePlus™ real-time PCR system (Applied Biosystems Inc., United States of America) to measure melting curves (from 4°C to 80°C, with a 0.5°C rise every 10 s and fluorescence detection). These experiments were performed in three duplicate tubes. An equal volume of the buffer solution instead of the aptamer working solution was used as a blank control.
The melt region derivative data (dF/dT) of each tube was exported from the PCR system. The average dF/dT value of the three replicate tubes at each temperature point was calculated and then corrected by subtracting the average dF/dT value of the three blank control tubes. The dF/dT values were normalized (0–100). The melting curves were plotted with the temperature as the X-axis and the corrected or normalized dF/dT values as the Y-axis, and the optimal concentration of the aptamer was determined based on the melting temperature (Tm), peak height and peak area.
2.2 Determination of the optimal buffer system by TFA
To determine the optimal buffer for the experiment, the interaction of aptamer AP273 (at the optimal concentration determined above) and AFP (1.45 nM) was performed in three buffers (20 mM HEPES, 10 mM PBS and 20 mM Tris-HCl, pH = 7.35–7.45) with constant metal ion concentrations (Na+ 140 mM, Mg2+ 2 mM, K+ 4 mM, and Ca2+ 1 mM), and then the melting curves were measured via TFA. The experimental procedure and interpretation of the results were the same as described above.
2.3 Determination of optimal metal ion strength by TFA
The optimal metal ion strength for the experiment was also determined via TFA. The concentration points of metal ions were set according to the results reported in the literature as follows (Duan et al., 2012; Kang et al., 2015; Wu et al., 2017; Chen et al., 2020a): 1, 2 and 5 mM for Mg2+; 100, 120 and 140 mM for Na+; 2, 4 and 5 mM for K+; and 1, 2 and 2.5 mM for Ca2+. A 4-factor, 3-level orthogonal test (with a fixed random seed number of 300) was designed to determine the optimal concentration of metal ions. Based on the orthogonal design, nine buffers with different concentrations of metal ions were prepared using the optimal buffer system determined in the previous step (Table 1). The concentrations of aptamers and AFP and the experimental procedures were the same as described above. In subsequent MD analysis, the buffers with the optimal metal ion combination were selected according to the Tm value, peak height and peak area.
2.4 Structure prediction, docking and MD simulations of the aptamers to AFP
Structure prediction and docking of two aptamers and AFP were performed by the following software: 1) The secondary structure of the aptamers was predicted using the online software Mfold (http://www.unafold.org/hybrid2.php) (Zuker, 2003). Aptamer structures with the minimum energy (the lowest ΔG value) were selected. The ct. output formats (dot-bracket notation) were used for the construction of the aptamer 3D structures (Garcia-Recio et al., 2016; Heiat et al., 2016; Subki et al., 2020). 2) The tertiary structure of aptamers was predicted by using online RNA Composer software on the basis of the secondary structure of aptamers (https://rnacomposer.cs.put.poznan.pl) (Popenda et al., 2012; Antczak et al., 2016). 3) The nucleotide change for the aptamer was performed by using DS (Biovia Discovery Studio) software (the uracil was changed to thymine, and the ribose backbone was changed to deoxyribose) (Chang et al., 2020). 4) All-atom (AA) energy minimization simulations of the aptamer were performed in a vacuum by using software NAnoscale Molecular Dynamics (NAMD) (Phillips et al., 2005) with the Visual Molecular Dynamics (VMD) (Humphrey et al., 1996). 5) The 3D model of human AFP was obtained from either the AlphaFoldDB website (https://alphafold.ebi.ac.uk/) using code number AF-J3KMX3-F or the PDBe-KB website (https://www.ebi.ac.uk/pdbe/pdbe-kb/proteins/P02771) with the PDB code 7YIM (UniProt et al., 2023). 6) Molecular docking of the 3D structure of the aptamer and AFP was performed online using ZDOCK 3.0.2 (https://zdock.umassmed.edu/) (Mintseris et al., 2007; Pierce et al., 2011; Pierce et al., 2014). 7) The model with the highest ZDOCK score (Wang et al., 2019) was selected as the initial structure for MD simulations, using Amber 18 software (Lee et al., 2020), in which the aptamer was assigned OL15 nucleic acid force fields (Zgarbova et al., 2015) and the target molecule was assigned ff14SB protein force fields (Tian et al., 2020). The reaction system solvent cassette (Steinberg et al., 2019) was set up according to the metal concentrations screened by the TFA experiment. 8) Docking, MD simulation results, hydrogen bonding, and binding sites were visualized using PyMOL software (The PyMOL Molecular Graphics System, Version 2.0 Schrödinger, LLC) (Lill and Danielson, 2011). The optimal buffer system was determined based on comprehensive analyses of the simulations.
2.5 Binding characteristics of aptamers to AFP under different ionic conditions
With optimal buffering conditions, TFA was performed for both fixed and gradient concentrations of aptamers to verify their binding characteristics to AFP. The experimental procedure and data processing were the same as above, and the binding characteristics of aptamers to AFP under various buffer systems were analyzed according to the Tm value, peak height and peak area.
2.6 Comparisons between different aptamers
In order to verify whether the combination of TFA and MD simulations can distinguish different aptamers and identify the optimal aptamer, the same experiments were performed with the aptamer AP-L3-4. Differences in the binding of aptamers AP273 and AP-L3-4 to AFP were compared, including the binding site, the distance and frequency of hydrogen bonds, binding energy, and binding characteristics for the gradient concentrations of AFP.
3 Results
3.1 Optimal concentration of aptamers determined by TFA
The melt curves of normalized dF/dT value and the heatmaps of peak area for gradient concentrations of aptamers AP273 and AP-L3-4 binding to AFP are shown in Figure 1. The free aptamer peak (peak 1) of AP273 appeared at 2.5 nM and then rose with increasing concentration and peaked at 10 nM, while the bound aptamer peak (peak 2) almost peaked at 2.5 nM and then gradually decreased from 10 nM (Figure 1A). A similar trend was found in aptamer AP-L3-4, but the free aptamer peak topped out and the bound aptamer peak started to decrease earlier (both at 5 nM) (Figure 1B). The heatmaps of the area under the peaks showed that the bound aptamer peak of AP273 was stronger than that of AP-L3-4 at the same aptamer concentration, while the free aptamer peak was the opposite (Figures 1C, D). These findings suggest that TFA can exhibit concentration-dependent changes in aptamer-target interactions and reflect differences between aptamers.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1144347/fchem-11-1144347-HTML/image_m/fchem-11-1144347-g001.jpg)
FIGURE 1. Thermofluorimetric analyses of the gradient concentrations of aptamers binding to AFP. (A) Melting curves of AP273 at 1.25–80 nM. (B) Melting curves of AP-L3-4 at 1.25–80 nM. (C) Heatmap of the area under the peaks of AP273 at 1.25–80 nM. (D) Heatmap of the area under the peaks of AP-L3-4 at 1.25–80 nM.
The results of refinement experiments for both aptamers at concentrations of 1–5 nM are shown in Figure 2. For AP273, the bound aptamer peak (peak 2) appeared at the lowest concentration (1 nM) and was maintained to the maximum concentration (5 nM), and its free aptamer peak (peak 1) appeared at 3 nM and continuously rose (Figure 2A). A similar trend was observed in AP-L3-4, although the free aptamer peak was more pronounced (Figure 2B). The heatmaps of the area under peak showed that AP273 had stronger bound aptamer peaks and weaker free aptamer peaks compared with AP-L3-4 (Figures 2C, D). Considering the need to detect different concentrations of AFP, we selected 4 nM as the optimal concentration for subsequent experiments.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1144347/fchem-11-1144347-HTML/image_m/fchem-11-1144347-g002.jpg)
FIGURE 2. Thermofluorimetric analyses of the refined gradient concentrations of aptamers binding to AFP. (A) Melting curves of AP273 at 1–5 nM. (B) Melting curves of AP-L3-4 at 1–5 nM. (C) Heatmap of the area under the peaks of AP273 at 1–5 nM. (D) Heatmap of the area under the peaks of AP-L3-4 at 1–5 nM.
3.2 Optimal buffer system and ion concentration for aptamer-target interactions determined by TFA
The melting curves of AP273 (4 nM) binding to AFP (1.45 nM) in three buffer systems (HEPES, PBS, Tris-HCl) showed different Tm values (Figure 3A), of which the HEPES buffer system had the largest Tm value and therefore was selected as the optimal buffer system for subsequent experiments. The Tm values of AP273 binding to AFP also varied in HEPES buffer solutions (buffer 1–9) with various metal ion concentrations (according to the orthogonal design) (Figure 3B), with the lowest Tm value in buffer 4 and the highest Tm value in buffer 6. The Tm value was higher in the buffer solution with the lowest ionic strength (buffer 6) than in the buffer solution with the highest ionic strength (buffer 4) (Figure 3C). Higher Tm values indicate more stable binding between aptamers and targets, and thus buffers 5, 6, and 9, which had higher Tm values were selected as buffer solutions for the next molecular dynamics analysis, and buffer 7 with the same concentration of Mg2+ ions, was selected as the control buffer solution.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1144347/fchem-11-1144347-HTML/image_m/fchem-11-1144347-g003.jpg)
FIGURE 3. Thermofluorimetric analyses of AP273 binding to AFP in different buffer and ionic strength conditions. (A) Melting curves in different buffer systems. (B) The melting temperatures in HEPES buffer system with various metal ion concentrations. *p < 0.05, **p < 0.01. (C) Melting curves of high and low ionic strength buffer solutions.
3.3 Secondary structure prediction of AP273 and its free energy
Prediction of the secondary structure is required before the tertiary structure of an aptamer can be determined. In the prediction of secondary structures of aptamer AP273, it was found that the amount and minimum free energy of secondary structure varied with the concentrations of Na+ and Mg2+ ions (Figure 4). Increased concentrations of Mg2+ ions led to fewer predicted secondary structures. The concentration of Na+ ions did not affect the amount of secondary structure, but it impacted the minimum free energy of the secondary structure. Under the same concentration of Mg2+ ions, higher concentrations of Na+ ions led to lower free energy values. The structure with the lowest value of ∆G (the most stable secondary structure) was selected for subsequent prediction of the tertiary structure.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1144347/fchem-11-1144347-HTML/image_m/fchem-11-1144347-g004.jpg)
FIGURE 4. The secondary structure and free energy (∆G) of AP273 predicted at different ion concentrations. The colour figures: secondary structures of AP273. Table: Na+ concentration in the first row, ∆G value in the second row. Red arrow: the secondary structure selected for the prediction of tertiary structure.
3.4 Molecular dynamics simulation of the interaction between AP273 and AFP
The docking models of AP273 and AFP were then subjected to MD simulations in different buffer solutions (Figure 5). The results of root mean square deviation (RMSD) (Figure 5A) and root mean square fluctuation (RMSF) (Figure 5B) showed that the stability of the complex of AP273 and AFP was the highest in buffer 9 and higher than that in buffer 6. The comparisons of AP273, AFP and their complex between buffer 6 and buffer 9 (Figures 5C, D) exhibited that AFP had the highest stability and AP273 had the lowest stability, and that the AP273-AFP complex was more stable in buffer 9 than in buffer 6. These results suggest that the metal ions in buffer systems can affect the stability of aptamer-target complexes.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1144347/fchem-11-1144347-HTML/image_m/fchem-11-1144347-g005.jpg)
FIGURE 5. Molecular dynamics simulations of aptamer AP273 binding to AFP. (A) The root means square deviation (RMSD) of the complex of aptamer AP273 and AFP (CPX) in various buffer solutions. (B) The root mean square fluctuation (RMSF) of residues during MD simulations in various buffer solutions. (C) and (D) The RMSDs of CPX, AFP, and AP273 in buffer 6 and buffer 9.
The free energy of AP273 binding to AFP in various buffer solutions was calculated by the MM/GBSA method (Rastelli et al., 2010; Hou et al., 2011; Genheden and Ryde, 2015; Chen et al., 2020b) and is shown in Table 2. As shown in Table 2, buffer 6 exhibited the largest binding energy but not the best stability (Figure 5A), which is not consistent with the results of TFA experiments, indicating that the binding force and stability of aptamer-target are influenced by other complex environmental factors, such as temperature and pH, etc.
The visualization of the MD simulation results provided qualitative binding information between AP273 and AFP in different buffer solutions (Figure 6), which indicated that the binding sites and hydrogen bonds were different under various buffer solutions and that metal ions were involved in the binding.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1144347/fchem-11-1144347-HTML/image_m/fchem-11-1144347-g006.jpg)
FIGURE 6. Visualization of MD simulation results in various buffer solutions. (A). Binding sites and the hydrogen bonds between AP273 and AFP in four buffer solutions. (B). The interaction of K+ and Mg2+ ions with the complex of aptamer AP273 and AFP in MD simulations.
The hydrogen bonds in AP273 binding to AFP under four buffer solutions were further quantified (Figures 7A, B). The results showed that the frequency of hydrogen bonds was the lowest in buffer 9 and the highest in buffer 6, and that the distance of hydrogen bonds was similar in the four buffer solutions, suggesting that the hydrogen bond frequency rather than hydrogen bond distance plays an important role in the binding between AP273 and AFP. TFA experiments were conducted under buffer 6 and buffer 9, and results showed that buffer 6 was more favorable than buffer 9 for AP273 to bind to AFP in terms of bound aptamer ratio and linear correlation (Figure 7C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1144347/fchem-11-1144347-HTML/image_m/fchem-11-1144347-g007.jpg)
FIGURE 7. Quantitative analyses of hydrogen bonds and linear correlations in the interaction between AP273 and AFP. (A) and (B) The maximum frequency and distance of hydrogen bonds in the four buffer solutions. (C) The linear correlations between AP273 and AFP in buffer 6 and buffer 9.
3.5 Comparison between aptamers AP273 and AP-L3-4
Comparative studies were performed between aptamers AP273 and AP-L3-4 to further verify whether the combination of TFA and MD simulation could select the aptamer with better binding properties. Compared with AP-L3-4, AP273 binding to AFP in MD simulations had more binding sites (Figures 8A, B), a shorter hydrogen bond distance, a higher hydrogen bond frequency, and less free energy (Figure 8C). TFA results showed better linear correlation of AP273 than AP-L3-4 for detecting gradient concentrations for AFP (Figure 8D). These findings demonstrate that the binding property of AP273 to AFP is superior to that of AP-L3-4, suggesting the theoretical and experimental feasibility of combining TFA and MD simulation for preferential aptamer selection.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1144347/fchem-11-1144347-HTML/image_m/fchem-11-1144347-g008.jpg)
FIGURE 8. Comparison of aptamers AP273 and AP-L3-4 binding to AFP in molecular dynamics simulation and thermofluorimetric analysis. (A) and (B) Docking models of AP273 and AP-L3-4 with AFP and their binding sites (green); (C) Distances and frequencies of hydrogen bonds and binding energies of AP273 and AP-L3-4 bound to AFP; (D) The linear correlations of AP273 and AP-L3-4 bound to AFP in thermofluorimetric analysis.
4 Discussion
TFA analyzes the binding of aptamers to their target molecules using melting curves. The free aptamer shows a smaller Tm value than that of the bound aptamer, so two peaks, peak 1 (the free aptamer peak) and peak 2 (the bound aptamer peak), can be observed on the melting curve of TFA (Kim et al., 2015; Mahmoud et al., 2019). In the present study, distinct peak 1 and peak 2 appeared on the melting curves after aptamers AP273 or AP-L3-4 were reacted with AFP, and the Tm value, peak height and peak area were different under various experimental conditions. Using these indicators, the optimal experimental system and the best aptamer can be determined, making TFA useful to optimize the aptamer-target reaction system and to select preferred aptamers. The Tm value of a bound aptamer peak is positively correlated with the stability of aptamer-target binding. In the present study, the Tm values for peak 2 differed across the three buffer systems, with the highest value in the HEPES buffer system, suggesting that the binding of AP273 to AFP is more stable in the HEPES buffer system compared to the other two buffer systems. This may be because HEPES is not prone to form complexes with metal ions in the buffer system and thus exhibits a better buffering effect (Ruzza et al., 2021). Using HEPES as the buffer salt, nine buffer solutions with various metal ion concentrations determined by orthogonal design were formulated for TFA experiments to analyze the effects of metal ions on the binding of AP273 to AFP. The Tm values of peak 2 were maximum in buffer 6 (with the lowest ionic strength) and minimum in buffer 4 (with the highest ionic strength), indicating that the high ionic strength is unfavorable for the binding of AP273 to AFP molecules. This is consistent with previous reports that a buffer solution with relatively low ionic strength is more appropriate for TFA experiments (Hayashi et al., 2014) and that immobilization of aptamers under low ionic strength conditions rather than conventional high ionic strength buffers can greatly improve the performance of the E-AB sensor (Liu et al., 2021).
The metal ions of the buffer system are important in aptamer-target binding. In the secondary structure prediction, we found that the number of secondary structures of AP273 varied with the Mg2+ and Na+ concentrations and that the Mg2+ concentration was proportional to the free energy (absolute value of ΔG) of the aptamer secondary structure. In MD simulations, we found that Mg2+ and K+ contributed to hydrogen bonds at the binding site of the aptamer to the target molecule (Figure 6B). Metal ions are mainly embedded in the helical grooves of aptamer tertiary structures and affect the binding properties of aptamers to the target molecules in different ways (McCluskey et al., 2019). Metal ions can neutralize the negative charge of the phosphate groups of aptamers to affect the amount of charge in aptamer structures, hence affecting the structure of aptamers; in this context, Mg2+ affects the flexibility of aptamers and the stability of aptamer-target compounds (Zhang and Yadavalli, 2010). The aptamer structure is unstable in the binding buffer without Mg2+ (Zhao et al., 2020), and Mg2+ is important for stabilizing the two- and three-dimensional structures of aptamers (Huang et al., 2013). In addition to Mg2+, the Na+, K+ and pH value of a buffer system also have an effect on the binding of aptamers to targets, and the optimization of them can significantly improve the efficiency of aptamer sensors in detecting targets (Hianik et al., 2007; Belleperche and DeRosa, 2018) and significantly increase the sensitivity of aptamer sensors (Yang et al., 2022). The binding of thrombin to its aptamer also depends on pH and electrolytes (Hianik et al., 2007). The pH value may affect the ssDNA conformational changes and the electron transfer between the target protein and the aptamer, which is related to the maintenance of the three-dimensional conformation of the aptamer; and an increase in Na+ concentration leads to a weaker binding of thrombin to the aptamer, possibly due to the shielding effect of Na+ ions on the target. In addition, many studies have confirmed that K+ ions can stabilize the structure of the aptamer, particularly for G-quadruplexes (Santos et al., 2022).
The stability and affinity of aptamer-target complexes involve intermolecular covalent bonds, hydrogen bonds, electrostatic interactions and van der Waals forces, in which hydrogen bonds are one of the strongest non-covalent interactions, and the number and distance of hydrogen bonds contribute to the affinity between aptamers and target molecules (Sabri et al., 2019). In the MD simulations, we found that more hydrogen bonds were formed in the binding of AP273 to AFP in buffer 6 than in buffer 9, indicating a stronger binding affinity between the two molecules in buffer 6. However, the most stable binding of AP273 to AFP was found in buffer 9 rather than buffer 6, where the hydrogen bonds were less frequent but the aptamer-target binding sites were more numerous and dispersed. This may suggest the stability is more related to the binding sites of the aptamer and the target, while the affinity is more related to the hydrogen bonds between the aptamer and the target.
In the comparative study of binding properties between aptamer AP273 and aptamer AP-L3-4, we found that the binding properties of AP273 were superior to those of AP-L3-4 in both TFA and MD simulations, suggesting that the combination of TFA and MD simulation can identify aptamers that exhibit better binding characteristics among homologous aptamers.
Additionally, we found that the stability of AP273 binding to AFP was inconsistent between TFA and MD simulations. In TFA experiments, the most stability was found in buffer 6, followed by buffer 9, buffer 5 and buffer 7, while in MD simulations, the most stability was observed in buffer 9, followed by buffer 6, buffer 7 and buffer 5. This may be due to subtle differences in experimental conditions between MD simulations and TFA. The MD simulations were performed under conditions of constant temperature (25°C), neutral reaction environment, atmospheric pressure, aqueous metal ion solution, and no buffer salt, while the TFA were performed under the HEPES buffer system and constantly increasing temperature. Furthermore, it is important to keep in mind the limitations of molecular simulation methods. The results of MD are microscopic simulations and predictions from large data sets; the prediction of aptamer 3D structure based only on the sequence is still a very unreliable task, which has by no means achieved the maturity of protein structure prediction. There is no guarantee that the prediction of an aptamer structure is close to the correct fold. Certainly, an all-atom energy minimization is not sufficient to assess the stability of the predicted aptamer structure, let alone when done in a vacuum.
5 Conclusion
In this study, the binding of aptamer AP273 to its target AFP was analyzed by TFA and underlying molecular mechanisms were analyzed through MD simulations. On the melting curve of TFA, the free and bound AP273 peaks were clearly visible, and the binding status of AP273 with AFP could be determined based on the peak-related indicators (Tm value, height or area of the peaks, and the ratio of the two peaks). These indicators varied with changes in aptamer-target ratio, buffer system and metal ionic strength, which can be useful for optimizing experimental conditions. The analysis of MD simulation showed that the affinity and stability of AP273 binding to AFP varied at different metal ionic strengths, and the underlying mechanisms were related to the difference in the hydrogen bond frequency and binding sites. In the comparative study of AP273 and the control aptamer AP-L3-4, the two aptamers were different in the peak-related indicators and in hydrogen bonds and binding free energies, which are useful for the preferential selection of aptamers. Although the stability of AP273 binding to AFP was inconsistent between TFA and MD simulation, the combination of the two methods provided dual validation that was more intuitive than mere experimental results and more convincing than mere theoretical analysis, and it will facilitate the translational application of aptamers in bioassays.
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found here: https://alphafold.ebi.ac.uk/; https://www.ebi.ac.uk/pdbe/pdbe-kb/proteins/P02771.
Author contributions
H-LZ: conceptualization, methodology, software, investigation, data interpretation, data curation, writing—original draft. CL: validation, data interpretation. Z-HL: formal analysis. SJ: investigation. DC: investigation. S-SL: investigation. TW: project administration, conceptualization, data interpretation, writing—review and editing, funding acquisition, supervision. K-HZ: conceptualization, data interpretation, resources, writing—review and editing, funding acquisition, supervision.
Funding
This study is supported by the National Natural Science Foundation of China (82160444) and the Science and Technology Project of Jiangxi Province, China (20192BBG70048).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Antczak, M., Popenda, M., Zok, T., Sarzynska, J., Ratajczak, T., Tomczyk, K., et al. (2016). New functionality of RNAComposer: An application to shape the axis of miR160 precursor structure. Acta Biochim. Pol. 63, 737–744. doi:10.18388/abp.2016_1329
Belleperche, M., and DeRosa, M. C. (2018). pH-control in aptamer-based diagnostics, therapeutics, and analytical applications. Pharm. (Basel) 11, 80. doi:10.3390/ph11030080
Cai, S., Yan, J., Xiong, H., Liu, Y., Peng, D., and Liu, Z. (2018). Investigations on the interface of nucleic acid aptamers and binding targets. Analyst 143, 5317–5338. doi:10.1039/c8an01467a
Chang, Y. T., Lin, Y. C., Sun, L., Liao, W. C., Wang, C. C. N., Chou, C. Y., et al. (2020). Wilforine resensitizes multidrug resistant cancer cells via competitive inhibition of P-glycoprotein. Phytomedicine 71, 153239. doi:10.1016/j.phymed.2020.153239
Chen, W., Gao, G., Jin, Y., and Deng, C. (2020). A facile biosensor for Aβ40O based on fluorescence quenching of prussian blue nanoparticles. Talanta 216, 120930. doi:10.1016/j.talanta.2020.120930
Chen, Y., Zheng, Y., Fong, P., Mao, S., and Wang, Q. (2020). The application of the MM/GBSA method in the binding pose prediction of FGFR inhibitors. Phys. Chem. Chem. Phys. 22, 9656–9663. doi:10.1039/d0cp00831a
Cheng, S. A., Fung, W. K., Chan, K. Y., and Shen, P. K. (2003). Optimizing electron spin resonance detection of hydroxyl radical in water. Chemosphere 52, 1797–1805. doi:10.1016/s0045-6535(03)00369-2
Christina Schmidt, M. P., Klussmann, S., Betzel, C., and Betzel, C. (2020). Molecular characterization of a ghrelin-l-aptamer complex. J. Mol. Struct. 1204, 127510. doi:10.1016/j.molstruc.2019.127510
Cui, X., Song, M., Liu, Y., Yuan, Y., Huang, Q., Cao, Y., et al. (2020). Identifying conformational changes of aptamer binding to theophylline: A combined biolayer interferometry, surface-enhanced Raman spectroscopy, and molecular dynamics study. Talanta 217, 121073. doi:10.1016/j.talanta.2020.121073
Darmostuk, M., Rimpelova, S., Gbelcova, H., and Ruml, T. (2015). Current approaches in SELEX: An update to aptamer selection technology. Biotechnol. Adv. 33, 1141–1161. doi:10.1016/j.biotechadv.2015.02.008
Davydova, A. S., Timoshenko, V. V., Lomzov, A. A., Pyshnyi, D. V., Venyaminova, A. G., and Vorobyeva, M. A. (2020). G-quadruplex 2'-F-modified RNA aptamers targeting hemoglobin: Structure studies and colorimetric assays. Anal. Biochem. 611, 113886. doi:10.1016/j.ab.2020.113886
Dong, L., Tan, Q., Ye, W., Liu, D., Chen, H., Hu, H., et al. (2015). Screening and identifying a novel ssDNA aptamer against alpha-fetoprotein using CE-SELEX. Sci. Rep. 5, 15552. doi:10.1038/srep15552
Douaki, A., Garoli, D., Inam, A., Angeli, M. A. C., Cantarella, G., Rocchia, W., et al. (2022). Smart approach for the design of highly selective aptamer-based biosensors. Biosens. (Basel) 12, 574. doi:10.3390/bios12080574
Duan, N., Wu, S., Chen, X., Huang, Y., and Wang, Z. (2012). Selection and identification of a DNA aptamer targeted to Vibrio parahemolyticus. J. Agric. Food Chem. 60, 4034–4038. doi:10.1021/jf300395z
Ellington, A. D., and Szostak, J. W. (1990). In vitro selection of RNA molecules that bind specific ligands. Nature 346, 818–822. doi:10.1038/346818a0
Emrani, A. S., Danesh, N. M., Ramezani, M., Taghdisi, S. M., and Abnous, K. (2016). A novel fluorescent aptasensor based on hairpin structure of complementary strand of aptamer and nanoparticles as a signal amplification approach for ultrasensitive detection of cocaine. Biosens. Bioelectron. 79, 288–293. doi:10.1016/j.bios.2015.12.025
Garcia-Recio, E. M., Pinto-Diez, C., Perez-Morgado, M. I., Garcia-Hernandez, M., Fernandez, G., Martin, M. E., et al. (2016). Characterization of MNK1b DNA aptamers that inhibit proliferation in MDA-mb231 breast cancer cells. Mol. Ther. Nucleic Acids 5, e275. doi:10.1038/mtna.2015.50
Genheden, S., and Ryde, U. (2015). The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities. Expert Opin. Drug Discov. 10, 449–461. doi:10.1517/17460441.2015.1032936
Hayashi, T., Oshima, H., Mashima, T., Nagata, T., Katahira, M., and Kinoshita, M. (2014). Binding of an RNA aptamer and a partial peptide of a prion protein: Crucial importance of water entropy in molecular recognition. Nucleic Acids Res. 42, 6861–6875. doi:10.1093/nar/gku382
Heiat, M., Najafi, A., Ranjbar, R., Latifi, A. M., and Rasaee, M. J. (2016). Computational approach to analyze isolated ssDNA aptamers against angiotensin II. J. Biotechnol. 230, 34–39. doi:10.1016/j.jbiotec.2016.05.021
Hianik, T., Ostatna, V., Sonlajtnerova, M., and Grman, I. (2007). Influence of ionic strength, pH and aptamer configuration for binding affinity to thrombin. Bioelectrochemistry 70, 127–133. doi:10.1016/j.bioelechem.2006.03.012
Hilder, T. A., and Hodgkiss, J. M. (2017). The bound structures of 17β-estradiol-binding aptamers. Chemphyschem 18, 1881–1887. doi:10.1002/cphc.201700363
Hou, T., Wang, J., Li, Y., and Wang, W. (2011). Assessing the performance of the MM/PBSA and MM/GBSA methods. 1. The accuracy of binding free energy calculations based on molecular dynamics simulations. J. Chem. Inf. Model 51, 69–82. doi:10.1021/ci100275a
Hu, J., and Easley, C. J. (2017). Homogeneous assays of second messenger signaling and hormone secretion using thermofluorimetric methods that minimize calibration burden. Anal. Chem. 89, 8517–8523. doi:10.1021/acs.analchem.7b02229
Hu, J., Kim, J., and Easley, C. J. (2015). Quantifying aptamer-protein binding via thermofluorimetric analysis. Anal. Methods 7, 7358–7362. doi:10.1039/c5ay00837a
Huang, R., Xiong, W., Wang, D., Guo, L., Lin, Z., Yu, L., et al. (2013). Label-free aptamer-based partial filling technique for enantioseparation and determination of DL-tryptophan with micellar electrokinetic chromatography. Electrophoresis 34, 254–259. doi:10.1002/elps.201200464
Humphrey, W., Dalke, A., and Schulten, K. (1996). Vmd: Visual molecular dynamics. J. Mol. Graph 14 (33-8), 33–38. doi:10.1016/0263-7855(96)00018-5
Kang, L., Yang, B., Zhang, X., Cui, L., Meng, H., Mei, L., et al. (2015). Enzymatic cleavage and mass amplification strategy for small molecule detection using aptamer-based fluorescence polarization biosensor. Anal. Chim. Acta 879, 91–96. doi:10.1016/j.aca.2015.03.030
Kim, J., Hu, J., Bezerra, A. B., Holtan, M. D., Brooks, J. C., and Easley, C. J. (2015). Protein quantification using controlled DNA melting transitions in bivalent probe assemblies. Anal. Chem. 87, 9576–9579. doi:10.1021/acs.analchem.5b03432
La Penna, G., and Chelli, R. (2018). Structural insights into the osteopontin-aptamer complex by molecular dynamics simulations. Front. Chem. 6, 2. doi:10.3389/fchem.2018.00002
Lee, T. S., Allen, B. K., Giese, T. J., Guo, Z., Li, P., Lin, C., et al. (2020). Alchemical binding free energy calculations in AMBER20: Advances and best practices for drug Discovery. J. Chem. Inf. Model 60, 5595–5623. doi:10.1021/acs.jcim.0c00613
Li, G., Zeng, J., Liu, H., Ding, P., Liang, J., Nie, X., et al. (2019). A fluorometric aptamer nanoprobe for alpha-fetoprotein by exploiting the FRET between 5-carboxyfluorescein and palladium nanoparticles. Mikrochim. Acta 186, 314. doi:10.1007/s00604-019-3403-z
Lill, M. A., and Danielson, M. L. (2011). Computer-aided drug design platform using PyMOL. J. Comput. Aided Mol. Des. 25, 13–19. doi:10.1007/s10822-010-9395-8
Liu, Y., Canoura, J., Alkhamis, O., and Xiao, Y. (2021). Immobilization strategies for enhancing sensitivity of electrochemical aptamer-based sensors. ACS Appl. Mater Interfaces 13, 9491–9499. doi:10.1021/acsami.0c20707
Longhini, A. P., LeBlanc, R. M., Becette, O., Salguero, C., Wunderlich, C. H., Johnson, B. A., et al. (2016). Chemo-enzymatic synthesis of site-specific isotopically labeled nucleotides for use in NMR resonance assignment, dynamics and structural characterizations. Nucleic Acids Res. 44, e52. doi:10.1093/nar/gkv1333
Mahmoud, M., Laufer, S., and Deigner, H. P. (2019). An aptamer based thermofluorimetric assay for ethanolamine. Biochimie 158, 233–237. doi:10.1016/j.biochi.2019.01.014
McCluskey, K., Boudreault, J., St-Pierre, P., Perez-Gonzalez, C., Chauvier, A., Rizzi, A., et al. (2019). Unprecedented tunability of riboswitch structure and regulatory function by sub-millimolar variations in physiological Mg2+. Nucleic Acids Res. 47, 6478–6487. doi:10.1093/nar/gkz316
Mintseris, J., Pierce, B., Wiehe, K., Anderson, R., Chen, R., and Weng, Z. (2007). Integrating statistical pair potentials into protein complex prediction. Proteins 69, 511–520. doi:10.1002/prot.21502
Moccia, F., Platella, C., Musumeci, D., Batool, S., Zumrut, H., Bradshaw, J., et al. (2019). The role of G-quadruplex structures of LIGS-generated aptamers R1.2 and R1.3 in IgM specific recognition. Int. J. Biol. Macromol. 133, 839–849. doi:10.1016/j.ijbiomac.2019.04.141
Morozov, D., Mironov, V., Moryachkov, R. V., Shchugoreva, I. A., Artyushenko, P. V., Zamay, G. S., et al. (2021). The role of SAXS and molecular simulations in 3D structure elucidation of a DNA aptamer against lung cancer. Mol. Ther. Nucleic Acids 25, 316–327. doi:10.1016/j.omtn.2021.07.015
Phillips, J. C., Braun, R., Wang, W., Gumbart, J., Tajkhorshid, E., Villa, E., et al. (2005). Scalable molecular dynamics with NAMD. J. Comput. Chem. 26, 1781–1802. doi:10.1002/jcc.20289
Pierce, B. G., Hourai, Y., and Weng, Z. (2011). Accelerating protein docking in ZDOCK using an advanced 3D convolution library. PLoS One 6, e24657. doi:10.1371/journal.pone.0024657
Pierce, B. G., Wiehe, K., Hwang, H., Kim, B. H., Vreven, T., and Weng, Z. (2014). ZDOCK server: Interactive docking prediction of protein-protein complexes and symmetric multimers. Bioinformatics 30, 1771–1773. doi:10.1093/bioinformatics/btu097
Poolsup, S., Zaripov, E., Huttmann, N., Minic, Z., Artyushenko, P. V., Shchugoreva, I. A., et al. (2023). Discovery of DNA aptamers targeting SARS-CoV-2 nucleocapsid protein and protein-binding epitopes for label-free COVID-19 diagnostics. Mol. Ther. Nucleic Acids 31, 731–743. doi:10.1016/j.omtn.2023.02.010
Popenda, M., Szachniuk, M., Antczak, M., Purzycka, K. J., Lukasiak, P., Bartol, N., et al. (2012). Automated 3D structure composition for large RNAs. Nucleic Acids Res. 40, e112. doi:10.1093/nar/gks339
Rastelli, G., Del Rio, A., Degliesposti, G., and Sgobba, M. (2010). Fast and accurate predictions of binding free energies using MM-PBSA and MM-GBSA. J. Comput. Chem. 31, 797–810. doi:10.1002/jcc.21372
Ruigrok, V. J. B., Levisson, M., Hekelaar, J., Smidt, H., Dijkstra, B. W., and Van der Oost, J. (2012). Characterization of aptamer-protein complexes by X-ray crystallography and alternative approaches. Int. J. Mol. Sci. 13, 10537–10552. doi:10.3390/ijms130810537
Ruzza, P., Honisch, C., Hussain, R., and Siligardi, G. (2021). Free radical generation in far-UV synchrotron radiation circular dichroism assays-protein and buffer composition contribution. Int. J. Mol. Sci. 22, 11325. doi:10.3390/ijms222111325
Sabri, M. Z., Abdul Hamid, A. A., Sayed Hitam, S. M., and Abdul Rahim, M. Z. (2019). In silico screening of aptamers configuration against hepatitis B surface antigen. Adv. Bioinforma. 2019, 6912914. doi:10.1155/2019/6912914
Santos, T., Lopes-Nunes, J., Alexandre, D., Miranda, A., Figueiredo, J., Silva, M. S., et al. (2022). Stabilization of a DNA aptamer by ligand binding. Biochimie 200, 8–18. doi:10.1016/j.biochi.2022.05.002
Shigdar, S. (2019). Aptamer-based diagnostics and therapeutics. Pharm. (Basel) 12, 6. doi:10.3390/ph12010006
Shoute, L. C. T., and Loppnow, G. R. (2018). Characterization of the binding interactions between EvaGreen dye and dsDNA. Phys. Chem. Chem. Phys. 20, 4772–4780. doi:10.1039/c7cp06058k
Steinberg, L., Russo, J., and Frey, J. (2019). A new topological descriptor for water network structure. J. Cheminform 11, 48. doi:10.1186/s13321-019-0369-0
Subki, A., Ho, C. L., Ismail, N. F. N., Zainal Abidin, A. A., and Balia Yusof, Z. N. (2020). Identification and characterisation of thiamine pyrophosphate (TPP) riboswitch in Elaeis guineensis. PLoS One 15, e0235431. doi:10.1371/journal.pone.0235431
Tabuchi, Y., Yang, J., and Taki, M. (2022). Relative nuclease resistance of a DNA aptamer covalently conjugated to a target protein. Int. J. Mol. Sci. 23, 7778. doi:10.3390/ijms23147778
Tian, C., Kasavajhala, K., Belfon, K. A. A., Raguette, L., Huang, H., Migues, A. N., et al. (2020). ff19SB: Amino-Acid-Specific protein backbone parameters trained against quantum mechanics energy surfaces in solution. J. Chem. Theory Comput. 16, 528–552. doi:10.1021/acs.jctc.9b00591
Tivon, Y., Falcone, G., and Deiters, A. (2021). Protein labeling and crosslinking by covalent aptamers. Angew. Chem. Int. Ed. Engl. 60, 16035–16040. doi:10.1002/ange.202101174
Tuerk, C., and Gold, L. (1990). Systematic evolution of ligands by exponential enrichment: RNA ligands to bacteriophage T4 DNA polymerase. Science 249, 505–510. doi:10.1126/science.2200121
UniProt, C., Martin, M. J., Orchard, S., Magrane, M., Ahmad, S., Alpi, E., et al. (2023). UniProt: The universal protein knowledgebase in 2023. Nucleic Acids Res. 51, D523–D531. doi:10.1093/nar/gkac1052
Vu, C. Q., Rotkrua, P., Soontornworajit, B., and Tantirungrotechai, Y. (2018). Effect of PDGF-B aptamer on PDGFRβ/PDGF-B interaction: Molecular dynamics study. J. Mol. Graph Model 82, 145–156. doi:10.1016/j.jmgm.2018.04.012
Wang, Q. L., Cui, H. F., Du, J. F., Lv, Q. Y., and Song, X. (2019). In silico post-SELEX screening and experimental characterizations for acquisition of high affinity DNA aptamers against carcinoembryonic antigen. RSC Adv. 9, 6328–6334. doi:10.1039/c8ra10163a
Wu, H., Liu, R., Kang, X., Liang, C., Lv, L., and Guo, Z. (2017). Fluorometric aptamer assay for ochratoxin A based on the use of single walled carbon nanohorns and exonuclease III-aided amplification. Mikrochim. Acta 185, 27. doi:10.1007/s00604-017-2592-6
Wu, L., Wang, Y., Xu, X., Liu, Y., Lin, B., Zhang, M., et al. (2021). Aptamer-based detection of circulating targets for precision medicine. Chem. Rev. 121, 12035–12105. doi:10.1021/acs.chemrev.0c01140
Yang, F., Li, J., Dong, H., Wang, G., Han, J., Xu, R., et al. (2022). A novel label-free electrochemiluminescence aptasensor using a tetrahedral DNA nanostructure as a scaffold for ultrasensitive detection of organophosphorus pesticides in a luminol-H(2)O(2) system. Analyst 147, 712–721. doi:10.1039/d1an02060a
Zgarbova, M., Sponer, J., Otyepka, M., Cheatham, T. E., Galindo-Murillo, R., and Jurecka, P. (2015). Refinement of the sugar-phosphate backbone torsion beta for AMBER force fields improves the description of Z- and B-dna. J. Chem. Theory Comput. 11, 5723–5736. doi:10.1021/acs.jctc.5b00716
Zhang, X., and Yadavalli, V. K. (2010). Molecular interaction studies of vascular endothelial growth factor with RNA aptamers. Analyst 135, 2014–2021. doi:10.1039/c0an00200c
Zhao, Q., Bai, Y., and Wang, H. (2020). Directing a rational design of aptamer-based fluorescence anisotropy assay for sensitive detection of immunoglobulin E by site-specific binding study. Talanta 217, 121018. doi:10.1016/j.talanta.2020.121018
Zhou, C., Hu, J., Ma, H., Yagoub, A. E., Yu, X., Owusu, J., et al. (2015). Antioxidant peptides from corn gluten meal: Orthogonal design evaluation. Food Chem. 187, 270–278. doi:10.1016/j.foodchem.2015.04.092
Zhu, Y., Cai, Y., Xu, L., Zheng, L., Wang, L., Qi, B., et al. (2015). Building an aptamer/graphene oxide FRET biosensor for one-step detection of bisphenol A. ACS Appl. Mater Interfaces 7, 7492–7496. doi:10.1021/acsami.5b00199
Keywords: aptamers, thermofluorimetric analysis, molecular dynamics simulation, reaction condition optimization, aptamer selection
Citation: Zhang H-L, Lv C, Li Z-H, Jiang S, Cai D, Liu S-S, Wang T and Zhang K-H (2023) Analysis of aptamer-target binding and molecular mechanisms by thermofluorimetric analysis and molecular dynamics simulation. Front. Chem. 11:1144347. doi: 10.3389/fchem.2023.1144347
Received: 14 January 2023; Accepted: 26 April 2023;
Published: 09 May 2023.
Edited by:
Riccardo Nifosì, National Research Council (CNR), ItalyCopyright © 2023 Zhang, Lv, Li, Jiang, Cai, Liu, Wang and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ting Wang, tingwang@ncu.edu.cn; Kun-He Zhang, khzhang@ncu.edu.cn