Protein Network Studies on PCOS Biomarkers With S100A8, Druggability Assessment, and RNA Aptamer Designing to Control Its Cyst Migration Effect

The prevalence of polycystic ovary syndrome (PCOS) has been gradually increasing among adult females worldwide. Laparoscopy drilling on ovary is the only available temporary solution with a high incidence of reoccurrence. S100A8 with S100A9 complex is believed to facilitate the cyst migration in PCOS condition. The high evident protein interaction network studies between PCOS biomarkers, cancer invasion markers, and the interactors of S100A8 confirm that this protein has strong interaction with other selective PCOS biomarkers, which may be associative in the immature cyst invasion process. Through the network studies, intensive structural and pathway analysis, S100A8 is identified as a targetable protein. In this research, the non-SELEX in silico method is adapted to construct RNA Library based on the consensus DNA sequence of Glucocorticoid Response Element (GRE) and screened the best nucleotide fragments which are bound within the active sites of the target protein. Selected sequences are joined as a single strand and screened the one which competitively binds with minimal energy. In vitro follow-up of this computational research, the designed RNA aptamer was used to infect the MCF7 cell line through Lipofectamine 2000 mediated delivery to study the anti-cell migration effect. Wound Scratch assay confirms that the synthesized 18-mer oligo has significant inhibition activity toward tumor cell migration at the cellular level.

Introduction

Nucleotide aptamers are successfully explored as better therapeutics to treat diseases and disorders. Time-consuming low-throughput procedures have been in practice to design the aptamers in vitro (Ghavami et al., 2009). Therefore, in silico non-SELEX approach is the better choice to perform the selection of aptamers, which involves the construction of an oligonucleotide library without amplification and binding them with suitable target protein unlike SELEX (Berezovski et al., 2006; Tseng et al., 2011). Designing the RNA aptamer for the validated biomarker helps us to normalize the disease state at the genetic level. Hence, delivering a well-designed aptamer against response elements (REs) can control the strange translation of the target gene. REs are the critical elements involved in the activation of target gene regulation. Inhibiting biomarkers of specific pathophysiological conditions at the molecular level is a better choice to oversee the disease (Strimbu and Tavel, 2010). Target validation is one of the necessary procedures in drug discovery protocol. Since the exact cause of polycystic ovary syndrome (PCOS) is imprecise (Sir-Petermann et al., 2002), it is tedious to identify the best target clinically. Assay on endometrial cell migration is one of the diagnostic tools to identify the complications of this syndrome, and metformin has a proven attenuating effect on the invasion of endometrial cells of diseased women (Tan et al., 2011). Previously, researchers have found that 500 biomarkers are prevalent in PCOS (Dai and Lu, 2012). In this research work, we focused on S100A8 protein which is one of the important biomarkers in PCOS. Protein–protein interaction network (PPIN) is used to identify the associative proteins and its pathways in PCOS. Additionally, druggable properties of S100A8 were studied through pocket analysis. Besides, the aptamer library for specific RE of S100A8 was constructed by a non-SELEX fragment approach. The best aptamer sequence was screened through quality assessments, such as affinity and stability parameters.

Materials and Methods

Network Profile of S100A8 in Polycystic Ovary Syndrome and Enrichment Analysis

Interactors of S100A8 are obtained from BioGRID³.⁵, a dataset repository (Oughtred et al., 2019), and the molecular interaction network was constructed in STRING Database (Szklarczyk et al., 2019). Biomarkers specific to PCOS and cancer cell invasion are retrieved from the recent research articles (Daan et al., 2016; Lu et al., 2017; Gerashchenko et al., 2019) and are used to construct another network. Both the networks were merged to find the first shell interactors of S100A8. Cytoscape 3.7.2 is employed to merge the networks and find proteins which are associated with S100A8. Pathways of S100A8 and its clusters are identified by using ClueGO, a Cytoscape application for clustering the functional network by terms or pathways (Bindea et al., 2009). Molecular functions of Gene Ontology (GO), Reactome Pathway Database (Croft et al., 2011), and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways are used as resources for enrichment analysis. P-value 0.005 is set as a significance to select the clusters of S100A8.

Structure and Druggability Studies on S100A8

Druggability analysis is used to predict the receptiveness and stability of drug target. Physiochemical and geometric properties such as number of pockets, druggable score, and pocket volume determine the efficiency of the target candidate. DoGSite Scorer is used for binding site prediction analysis and druggability assessment, which is based on heavy-atom coordinates employing support vector machines (SVMs) (Volkamer et al., 2012). Pocket volume, lipophilic character, and pocket enclosures were accounted for simple score calculation to suggest the competence of targetability. Three-dimensional structure of the target was retrieved from PDB (ID: 5HLV) and used for the druggable screening.

Glucocorticoid Response Elements for S100A8

REs are the inducers of the receptor and ligand interaction which results in the expression or activation of a particular protein. Since the aptamers are crucial elements in the control of target expression so we decide to design RNA oligomer against specific PCOS targets. Glucocorticoid RE (GRE) (Hsu et al., 2005), hypoxia RE (HRE) (Rees et al., 2001), antioxidant RE (ARE) (Nioi et al., 2003), and interferon gamma (INF-γ) RE (IRE) (Yang et al., 1990) are identified as the influencing REs of the S100A8 gene. Specifically, GREs have the proficiency to inhibit S100A8 through the downregulation of leukocyte transmigration. Glucocorticoids also induce the expression of inhibition factor for macrophage migration, which ultimately downregulates the cyst inflammation. The earlier research report shows that GRE consists of two half-sites with three spacer bases; the consensus pseudo palindromic sequence of GRE is 5′ CAGAACATCATGTTCTGA 3′ (Weikum et al., 2017).

Nucleotide Fragment Library Construction

RNA Composer utilizes the Dot-Bracket format notation of the secondary structure sequence to model the RNA, and the 3D element of modeled RNA was chosen from RNA Frabase (Biesiada et al., 2016). RE is a sequence, which binds with the receptor and plays a crucial role in expression, so the RNA analog library of specific RE was created to mimic the inhibitory action. The consensus sequence was segregated as fragments in such a way that six nucleotides at a stretch were taken per fragment (Figure 4A) for analog library construction. Resulted library sequences were later utilized for binding studies with the target by RNA-Lim method and recognized the various conformations of fragments bound in the active sites of the protein (Hall et al., 2015). Fragments which bound on active sites are selected to design the high précised aptamer model. Diversity in the exhibited conformations of ssRNA–protein complexes was meticulously sampled to construct a fragment library. MC-Fold | MC-Sym pipeline was employed to obtain the secondary and tertiary structures of the constructed aptamer (Parisien and Major, 2008). The proposed mechanism for PCOS control through the aptamer binding on S100A8 is illustrated in Figure 2. Refinement on docking results of their chain-forming poses was done in Molecular Operating Environment (MOE) suit (Ahirwar et al., 2016). The fragment-based approach was adopted for competent docking with S100A8 (Ahirwar et al., 2016); this unusual method has numerous advantages over conventional rigid based docking.

Affinity and Stability Studies of the Designed Aptamer

PatchDock tool is efficiently used to study the binding properties of designed aptamers with target protein (Schneidman-Duhovny et al., 2005). Based on the global binding energy, FireDock is employed to decipher the docked result by flexible refinements rather than the rigidity of protein and also it optimizes the side-chain residues, which minimizes the rigid body conformation of the interactive protein (Mashiach et al., 2008). Previously, it was reported that the stability of RNA will be analyzed by the inverted repeats which form stable hairpin loops (Ahirwar et al., 2016). Oligoanalyzer is an inclusive oligonucleotide scrutinizer employed to check the hairpin loop and stability of the designed aptamer (Owczarzy et al., 2008).

Anti-migration by Scratch Wound Healing Assay on Cell Lines

MCF-7 cells were seeded into a 24-well tissue culture plate containing antibiotic-free minimal essential medium (MEM) and incubated for 24 h at 37°C with 5% CO₂. Sterile microtip was used to a make scratch on the 80% confluence monolayer (Camorani et al., 2014). The culture medium was immediately replaced with fresh medium to remove the dislodged cells. Then, 20 nmol/l of the designed aptamer with the transfecting agent, Lipofectamine 2000 (Invitrogen), was dissolved in dimethyl sulfoxide (DMSO) for timeline studies (Zhou et al., 2008). Cell migration of both sample and control were monitored and compared to study the aptamer effect on cell migration.

Results and Discussion

S100A8 Network With Polycystic Ovary Syndrome Biomarkers

S100A8 interacts with 74 proteins (Supplementary File S1). The high confidential STRING network of S100A8 interactors has 55 nodes with 181 edges. The second network with biomarkers has 98 nodes with 431 edges; among 98, two proteins (SLC35D2, MORC4) are unconnected with the main network. Interactor network and biomarker network details are given in Supplementary Material. Merged network with connected nodes of the interactors and the biomarkers shows 96 nodes with 430 edges (interactions) that are illustrated in Figure 1. Immediate interacting proteins with s100A8 are shown in yellow. By the network analysis, we found that S100A8 directly interact with 10 proteins (Table 1).

FIGURE 1

Figure 1. Proposed mechanism of aptamer inhibition. Aptamer binds to the S100A8 and prevents the intracellular initial complexation with S100A9. Due to aptamer binding, extracellular receptor of advanced glycation end products (RAGE) complexed inflammatory sequences are also prevented.

FIGURE 2

Figure 2. Merged network of Polycystic Ovary Syndrome (PCOS). Combined network of S100A8 integrators with PCOS biomarkers; yellow colored nodes are the first shell interacting proteins with S100A8 that are suspected for cyst migration activity.

TABLE 1

Table 1. Proteins associated with S100A8 in the network.

A total of 246 ontology processes were found within the given significant p-value, among them, S100A8 is found in 88 different processes. Particularly, 10 ovulation and maturation-related GO terms with 14 proteins are identified in the enrichment analysis and are listed in Table 2. Apart from S100A8, RETN and S100A9 are found in both networks and also in enriched biological terms with high distribution. ESR1, GDF9, PDGFRA, and LEP are the other proteins found in a greater number of pathways comparatively in the selected terms. Distribution of proteins is given in the graph of Figure 3A.

TABLE 2

Table 2. Enriched terms of S100A8 in Polycystic Ovary Syndrome (PCOS).

FIGURE 3

Figure 3. Gene ontology enrichment and structural properties of S100A8. (A) Distribution of proteins in polycystic ovary syndrome (PCOS)-related terms in enrichments. (B) Homo dimer crystal structure of S100A8—global symmetric view. (C) High druggable pockets of S100A8—shown in filled purple, sky blue, yellow, and orange colors. (D) Active sites of S100A8.

Target Compatibility Evaluation

S100A8 is involved in seven GO functions which are positively regulated cyst formation and cancer cell migration. Additionally, S100A8 poses 18 pockets, among them, nine are druggable (score > 0.3) and four shows better cutoff scores (Table 3). Targets with low specificity on small-molecule were identified as poorly druggable targets (Barelier et al., 2010). Here we have found four high scored druggable pockets in the selected S100A8 (Figure 3C). Structural features and active sites of S100A8 are shown in Figures 3B,D, respectively. Due to the positive druggable results, it is considered as a notable target to control PCOS. Considerably, calcium-binding protein (S100A8) acts as a ligand for receptor of advanced glycation end products (RAGE) which is involved in many inflammatory and oncogenic pathways. There is evidence that S100A8 has a growth-promoting effect, and it helps cells to acquire cell migration activity through the RAGE binding pathway (Ghavami et al., 2008). S100A8 causes uteroplacental perfusion deficiency which leads to embryo abortion that supports the competence of our target selection (Sir-Petermann et al., 2002). Structural analysis shows that S100A8 has two helix loop helix Ca2+ binding domains known as EF-hands and exists as a complex with S100A9. Calprotectin is present in 1q21 locus of chromosome 1 in humans and has a molecular weight of 10–12 kDa. During tumor development, chromosomal rearrangements take place in the locus of the S100A8 gene and majorly contribute to the cyst formation in PCOS. Also, serum calgranulin (S100A8 and S100A9) levels are higher in women with PCOS than normal women (Dai and Lu, 2012). This evidently shows that binding of S100A8 with RAGE facilitates the p38 mitogen-activated protein (MAP) kinase signaling through calcium phosphorylation which also governs cyst migration.

TABLE 3

Table 3. Druggability assessment of S100A8 protein.

Construction of RNA Analog Library Using Glucocorticoid Response Element

The fragment-based approach of aptamer docking yielded better interaction with S100A8. By the RNA-Lim method, 18 fragments with the consensus sequence of GRE were constructed and used for binding analysis (Figure 4A). Frag6, Frag9, and Frag10 showed better interaction (Table 4) in the active domain of target with minimal global binding energy. Among the three possible conformations, sequence 1 (Figure 4B) shows better thermal stability and lowest energy than the other two sequences. Optimal structure with a folding simulation at physiological pH shows there are three nucleotides at positions 4–7 that make intramolecular base pairing for loop structure (Figure 4C). Energy minimized aptamers are significantly stable, and the aptamers with a binding energy of ≥-40 are optimal in the therapeutical aspect (Pagano et al., 2008). Oligo fragments selected are by their binding ability on the active sites of the target. The compiled 18-mer binds effectively than the fragments. Stability comparison among the newly constructed aptamer sequences is stated in Table 5. Among the three, Apt1 has high stability with a melting temperature of 41.8°C, and also the simulation studies confirmed that it requires the minimum free energy (–27.93 kcal/mol) for hairpin loop formation. In addition, the Apt1 fragment poses low molecular weight (5,327.4 g/mol) comparatively. Aptamers in practice are available in the range of 15–81 nucleotide length with higher molecular weights (Shigdar et al., 2013), but here, the designed is 18-mer with lower molecular weight so the plasma clearance may be faster.

FIGURE 4

Figure 4. RNA library and structures. (A) Library construction by RNA-Lim method using the consensus sequence of glucocorticoid response element (GRE). (B) 3D structure of potential aptamer candidate. (C) Sequence and physiological structure of aptamer after optimal simulation and folding, and the color of bases indicates their energy levels.

TABLE 4

Table 4. Binding energies of RNA analog fragments with S100A8.

TABLE 5

Table 5. Aptamer stability comparison.

Interaction, Inhibition, and Stability Studies

Primarily, RAGE being the receptor for S100A8 was docked to confirm for its binding ability in the domain, which may affect the binding of the designed aptamer. As a result of protein–protein docking analysis, Arg 114 residue at A domain of RAGE interacting with Gln 44 residue at H domain of S100A8 is found as the most active interaction. To test the comparison of binding interactions, RAGE was docked with GRE, which resulted in a binding energy of –24.38, comparatively higher than its binding with the designed aptamers (–46.33) that is shown in Table 6; this infers that the designed nucleotide aptamer also binds efficiently at the S100A8 binding domain of RAGE (Figures 5B,C). In parallel, interactions of S100A8 with GRE and S100A8 with the designed aptamer were inspected to find the competency, which was found as –22.11 and –45.32 energy levels, respectively. The designed aptamer binds efficiently in the active dimer of the target (Figure 5A).

TABLE 6

Table 6. Docking results.

FIGURE 5

Figure 5. Mechanism of inhibition and docking with RNA aptamer. (A) Aptamer bases binding at the active sites of S100A8. (B) Binding pose of aptamer on receptor of advanced glycation end products (RAGE) protein. (C) Optimal binding pose of aptamer on the S100A8.

Anti-cell Migration Assay on MCF-7 Cell Line

Within 4 h of scratch, development of closure was seen in the control (which does not have aptamer), the wounded area has turned into a normal layer when compared with the initial image of well. In the aptamer well, there is no cell migration observed even after the fourth hour of incubation, it was confirmed in the images of 0 and 4 h of wounded well (Figure 6).

FIGURE 6

Figure 6. Wound scratch assay. Cell migration studies by wound healing assay on MCF7 cell line. Leitz labovert FS inverted microscope was used to view the cell migration. (A) A 4 mm scratch with a sterile tip at 0 h, and the wound is pointed at with an arrowhead in the 10 × magnified illustration. (B) Ten times magnified view of aptamer incubated (4 h) cell line. No cell migration was observed after 4 h of aptamer incubation. (C) Fifty times magnification of scratch in control plate (cell migration seen) after the fourth hour.

Conclusion

From the network analysis, S100A8 is identified as a targetable protein to control PCOS. The druggable property of the target was validated by topological measures. S100A8 acts as a ligand for RAGE to promote cell migration in cancers and PCOS conditions. GRE inhibits S100A8 by competitive binding at the minimal level through a feedback mechanism. Additionally, S100A9 and resistin were also found along with S100A8 as associative proteins. We adopted a computational method to develop an RNA aptamer and designed 18 oligos based on the consensus sequences of GRE, which binds to both RAGE and S100A8. In addition to the computational studies, the cell line studies proved the anti-migration activity of the designed aptamer at minimal dose delivery with Lipofectamine 2000. The newly designed 18mer effectively stopped the cancer cell migration through dual action, and it is identified as a potential therapeutic to control PCOS and cancers.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, to any qualified researcher.

Author Contributions

All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.

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.

Acknowledgments

We thank the scientific committee of Kamaraj College of Engineering and Technology for their encouragement with best project award in TECHNOVISION’17. Our sincere thanks to Tamil Nadu State Council for Science and Technology (TNSCST), Directorate of Technical Education, Chennai, for their partial financial support through student project scheme. I dedicate this manuscript as a gift to my student A. Shobana for her recent birthday. This manuscript has been released as a pre-print at bioRxiv (https://www.biorxiv.org/content/10.1101/603357v1), doi: https://doi.org/10.1101/603357 entitled “Identification of target candidate in polycystic ovarian syndrome and in vitro evaluation of therapeutic activity of the designed RNA aptamer,” and this is a revised version.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbioe.2020.00328/full#supplementary-material

References

Ahirwar, R., Smita, N., Shikha, A., Srinivasan, R., Souvik, M., and Pradip, N. (2016). In silico selection of an aptamer to estrogen receptor alpha using computational docking employing estrogen response elements as aptamer-alike molecules. Sci. Rep. 6:21285. doi: 10.1038/srep21285

PubMed Abstract | CrossRef Full Text | Google Scholar

Barelier, S., Pons, J., Gehring, K., Lancelin, J.-M., and Krimm, I. (2010). Ligand specificity in fragment-based drug design. J. Med. Chem. 53, 5256–5266. doi: 10.1021/jm100496j

PubMed Abstract | CrossRef Full Text | Google Scholar

Berezovski, M. V., Musheev, M. U., Drabovich, A. P., Jitkova, J. V., and Krylov, S. N. (2006). Non-SELEX: selection of aptamers without intermediate amplification of candidate oligonucleotides. Nat. Protoc. 1:1359. doi: 10.1038/nprot.2006.200

PubMed Abstract | CrossRef Full Text | Google Scholar

Biesiada, M., Pachulska-Wieczorek, K., Adamiak, R. W., and Purzycka, K. J. (2016). RNAComposer and RNA 3D structure prediction for nanotechnology. Methods 103, 120–127. doi: 10.1016/j.ymeth.2016.03.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Bindea, G., Bernhard, M., Hubert, H., Pornpimol, C., Marie, T., Amos, K., et al. (2009). ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics 25, 1091–1093. doi: 10.1093/bioinformatics/btp101

PubMed Abstract | CrossRef Full Text | Google Scholar

Camorani, S., Esposito, C. L., Rienzo, A., Catuogno, S., Iaboni, M., Condorelli, G., et al. (2014). Inhibition of receptor signaling and of glioblastoma-derived tumor growth by a novel PDGFRbeta aptamer. Mol. Ther. 22, 828–841. doi: 10.1038/mt.2013.300

PubMed Abstract | CrossRef Full Text | Google Scholar

Croft, D., O’Kelly, G., Wu, G., Haw, R., Gillespie, M., Matthews, L., et al. (2011). Reactome: a database of reactions, pathways and biological processes. Nucleic Acids Res. 39, D691–D697. doi: 10.1093/nar/gkq1018

PubMed Abstract | CrossRef Full Text | Google Scholar

Daan, N. M., Koster, M. P., de Wilde, M. A., Dalmeijer, G. W., Evelein, A. M., Fauser, B. C., et al. (2016). Biomarker profiles in women with PCOS and PCOS Offspring; a pilot study. PLoS One 11:e0165033. doi: 10.1371/journal.pone.0165033

PubMed Abstract | CrossRef Full Text | Google Scholar

Dai, G., and Lu, G. (2012). Different protein expression patterns associated with polycystic ovary syndrome in human follicular fluid during controlled ovarian hyperstimulation. Reprod Fertil. Dev. 24, 893–904. doi: 10.1071/RD11201

PubMed Abstract | CrossRef Full Text | Google Scholar

Gerashchenko, T. S., Nikita, M. N., Nadezhda, V. K., Sofia, Y. Z., Marina, V. Z., and Nadezhda, V. (2019). Markers of cancer cell invasion: are they good enough? J. Clin. Med. 8:1092. doi: 10.3390/jcm8081092

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghavami, S., Rashedi, I., Dattilo, B. M., Eshraghi, M., Chazin, W. J., Hashemi, M., et al. (2008). S100A8/A9 at low concentration promotes tumor cell growth via RAGE ligation and MAP kinase-dependent pathway. J. Leukoc. Biol. 83, 1484–1492. doi: 10.1189/jlb.0607397

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghavami, S., Seth, C., Mohammad, H., Mehdi, E., Walter, J. C., Andrew, J. H., et al. (2009). S100A8/A9: a janus-faced molecule in cancer therapy and tumorgenesis. Eur. J. Pharmacol. 625, 73–83. doi: 10.1016/j.ejphar.2009.08.044

PubMed Abstract | CrossRef Full Text | Google Scholar

Hall, D., Li, S., Yamashita, K., Azuma, R., Carver, J. A., and Standley, D. M. (2015). RNA-LIM: a novel procedure for analyzing protein/single-stranded RNA propensity data with concomitant estimation of interface structure. Anal. Biochem. 472, 52–61. doi: 10.1016/j.ab.2014.11.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Hsu, K., Passey, R. J., Endoh, Y., Rahimi, F., Youssef, P. P., Yen, T., et al. (2005). Regulation of S100A8 by glucocorticoids. J. Immunol. 174, 2318–2326. doi: 10.4049/jimmunol.174.4.2318

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, C., Xiaoqin, L., Lin, W., Ning, J., Jun, Y., Xiaobo, Z., et al. (2017). Integrated analyses for genetic markers of polycystic ovary syndrome with 9 case-control studies of gene expression profiles. Oncotarget 8, 3170–3180. doi: 10.18632/oncotarget.13881

PubMed Abstract | CrossRef Full Text | Google Scholar

Mashiach, E., Schneidman-Duhovny, D., Andrusier, N., Nussinov, R., and Wolfson, H. J. (2008). FireDock: a web server for fast interaction refinement in molecular docking. Nucleic Acids Res. 36, W229–W232. doi: 10.1093/nar/gkn186

PubMed Abstract | CrossRef Full Text | Google Scholar

Nioi, P., McMahon, M., Itoh, K., Yamamoto, M., and Hayes, J. D. (2003). Identification of a novel Nrf2-regulated antioxidant response element (ARE) in the mouse NAD(P)H:quinone oxidoreductase 1 gene: reassessment of the ARE consensus sequence. Biochem. J. 374(Pt 2), 337–348. doi: 10.1042/bj20030754

PubMed Abstract | CrossRef Full Text | Google Scholar

Oughtred, R., Chris, S., Bobby-Joe, B., Jennifer, R., Lorrie, B., Christie, C., et al. (2019). The BioGRID interaction database: 2019 update. Nucleic Acids Res. 47, D529–D541. doi: 10.1093/nar/gky1079

PubMed Abstract | CrossRef Full Text | Google Scholar

Owczarzy, R., Tataurov, A. V., Wu, Y., Manthey, J. A., McQuisten, K. A., Almabrazi, H. G., et al. (2008). IDT SciTools: a suite for analysis and design of nucleic acid oligomers. Nucleic Acids Res. 36, W163–W169. doi: 10.1093/nar/gkn198

PubMed Abstract | CrossRef Full Text | Google Scholar

Pagano, B., Martino, L., Randazzo, A., and Giancola, C. (2008). Stability and binding properties of a modified thrombin binding aptamer. Biophys. J. 94, 562–569. doi: 10.1529/biophysj.107.117382

PubMed Abstract | CrossRef Full Text | Google Scholar

Parisien, M., and Major, F. (2008). The MC-Fold and MC-Sym pipeline infers RNA structure from sequence data. Nature 452, 51–55. doi: 10.1038/nature06684

PubMed Abstract | CrossRef Full Text | Google Scholar

Rees, B. B., Bowman, J. A., and Schulte, P. M. (2001). Structure and sequence conservation of a putative hypoxia response element in the lactate dehydrogenase-B gene of Fundulus. Biol. Bull. 200, 247–251. doi: 10.2307/1543505

PubMed Abstract | CrossRef Full Text | Google Scholar

Schneidman-Duhovny, D., Inbar, Y., Nussinov, R., and Wolfson, H. J. (2005). PatchDock and SymmDock: servers for rigid and symmetric docking. Nucleic Acids Res. 33, W363–W367.

PubMed Abstract | Google Scholar

Shigdar, S., Qiao, L., Zhou, S.-F., Xiang, D., Wang, T., Li, Y., et al. (2013). RNA aptamers targeting cancer stem cell marker CD133. Cancer Lett. 330, 84–95. doi: 10.1016/j.canlet.2012.11.032

PubMed Abstract | CrossRef Full Text | Google Scholar

Sir-Petermann, T., Maliqueo, M., Angel, B., Lara, H. E., Pérez-Bravo, F., Recabarren, S. E., et al. (2002). Maternal serum androgens in pregnant women with polycystic ovarian syndrome: possible implications in prenatal androgenization. Hum. Reprod. 17, 2573–2579. doi: 10.1093/humrep/17.10.2573

PubMed Abstract | CrossRef Full Text | Google Scholar

Strimbu, K., and Tavel, J. A. (2010). What are biomarkers? Curr. Opin. HIV AIDS 5, 463–466. doi: 10.1097/COH.0b013e32833ed177

PubMed Abstract | CrossRef Full Text | Google Scholar

Szklarczyk, D., Annika, L. G., David, L., Alexander, J., Stefan, W., Jaime, H., et al. (2019). STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 47, D607–D613. doi: 10.1093/nar/gky1131

PubMed Abstract | CrossRef Full Text | Google Scholar

Tan, B. K., Raghu, A., Jing, C., Hendrik, L., Louis, J. S. C., Harpal, S. R., et al. (2011). Metformin treatment exerts antiinvasive and antimetastatic effects in human endometrial carcinoma cells. J. Clin. Endocrinol. Metab. 96, 808–816. doi: 10.1210/jc.2010-1803

PubMed Abstract | CrossRef Full Text | Google Scholar

Tseng, C. Y., Ashrafuzzaman, M., Mane, J. Y., Kapty, J., Mercer, J. R., and Tuszynski, J. A. (2011). Entropic fragment-based approach to aptamer design. Chem. Biol. Drug Des. 78, 1–13. doi: 10.1111/j.1747-0285.2011.01125.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Volkamer, A., Kuhn, D., Rippmann, F., and Rarey, M. (2012). DoGSiteScorer: a web server for automatic binding site prediction, analysis and druggability assessment. Bioinformatics 28, 2074–2075. doi: 10.1093/bioinformatics/bts310

PubMed Abstract | CrossRef Full Text | Google Scholar

Weikum, E. R., Knuesel, M. T., Ortlund, E. A., and Yamamoto, K. R. (2017). Glucocorticoid receptor control of transcription: precision and plasticity via allostery. Nat. Rev. Mol. Cell Biol. 18, 159–174. doi: 10.1038/nrm.2016.152

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Z., Sugawara, M., Ponath, P. D., Wessendorf, L., Banerji, J., Li, Y., et al. (1990). Interferon gamma response region in the promoter of the human DPA gene. Proc. Natl. Acad. Sci. U.S.A. 87, 9226–9230. doi: 10.1073/pnas.87.23.9226

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, J., Li, H., Li, S., Zaia, J., and Rossi, J. J. (2008). Novel dual inhibitory function aptamer-siRNA delivery system for HIV-1 therapy. Mol. Ther. 16, 1481–1489. doi: 10.1038/mt.2008.92

PubMed Abstract | CrossRef Full Text | Google Scholar