ORIGINAL RESEARCH article

Front. Bioeng. Biotechnol., 13 May 2020

Sec. Computational Genomics

Volume 8 - 2020 | https://doi.org/10.3389/fbioe.2020.00328

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

  • 1. Centre for Research, Kamaraj College of Engineering and Technology, Madurai, India

  • 2. Department of Biotechnology, Kamaraj College of Engineering and Technology, Madurai, India

Abstract

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 (). 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 (; ). 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 (). Target validation is one of the necessary procedures in drug discovery protocol. Since the exact cause of polycystic ovary syndrome (PCOS) is imprecise (), 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 (). Previously, researchers have found that 500 biomarkers are prevalent in PCOS (). 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 BioGRID3.5, a dataset repository (), and the molecular interaction network was constructed in STRING Database (). Biomarkers specific to PCOS and cancer cell invasion are retrieved from the recent research articles (; ; ) 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 (). Molecular functions of Gene Ontology (GO), Reactome Pathway Database (), 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) (). 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) (), hypoxia RE (HRE) (), antioxidant RE (ARE) (), and interferon gamma (INF-γ) RE (IRE) () 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′ ().

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 (). 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 (). 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 (). 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 (). The fragment-based approach was adopted for competent docking with S100A8 (); 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 (). 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 (). Previously, it was reported that the stability of RNA will be analyzed by the inverted repeats which form stable hairpin loops (). Oligoanalyzer is an inclusive oligonucleotide scrutinizer employed to check the hairpin loop and stability of the designed aptamer ().

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% CO2. Sterile microtip was used to a make scratch on the 80% confluence monolayer (). 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 (). 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 2

TABLE 1

Protein IDNamePolycystic Ovary Syndrome (PCOS) relative function
MMP9Matrix Metallopeptidase 9Proteolytic activity on the extracellular matrix (ECM) and involved in leukocyte migration
RETNResistinPromotes chemotaxis in myeloid cells
CTSSCathepsin SEndo protease involved in the removal of unwanted proteins
S100A9S100 Calcium Binding Protein A9Potent amplifier of inflammation as well as in cancer development and tumor spread
NCF2Neutrophil Cytosolic Factor 2Involved in superoxide generation
MPOMyeloperoxidaseProduce hypochlorous acid and other toxic intermediates which enhance PMN microbicidal activity
ALBAlbuminRegulates blood plasma colloid osmotic pressure and acts as a carrier protein for hormones, fatty acids, metabolites, and exogenous drugs
BECN1Beclin 1Mediates vesicle-trafficking processes, tumorigenesis, neurodegeneration, and apoptosis
TP53Tumor Protein P53Prevents CDK7 kinase activity when associated to CAK complex in response to DNA damage, thus stopping cell cycle progression
GRB2Growth Factor Receptor Bound Protein 2Adapter protein involved in the Ras signaling pathway

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

GO termGroup p-valueNo. of proteinsAssociated proteins found
Ovulation cycle2.11E-139EGFR, ESR1, GDF9, HAS2, LEP, PDGFRA, RETN, S100A8, S100A9
Female sex differentiation6.47E-2112CTNNA1, ESR1, GDF9, ICAM1, LEP, MYC, PDGFRA, RBP4, RETN, S100A8, S100A9, VEGFA
Ovarian follicle development6.47E-217CTNNA1, ESR1, ICAM1, MYC, S100A8, S100A9, VEGFA
Ovulation cycle1.38E-249EGFR, ESR1, GDF9, HAS2, LEP, PDGFRA, RETN, S100A8, S100A9
Ovulation cycle process1.38E-247ESR1, GDF9, LEP, PDGFRA, RETN, S100A8, S100A9
Development of primary female sexual characteristics1.38E-2411CTNNA1, ESR1, GDF9, ICAM1, LEP, MYC, PDGFRA, RETN, S100A8, S100A9, VEGFA
Gonad development1.38E-2411CTNNA1, ESR1, GDF9, ICAM1, LEP, MYC, PDGFRA, RETN, S100A8, S100A9, VEGFA
Ovarian follicle development1.38E-247CTNNA1, ESR1, ICAM1, MYC, S100A8, S100A9, VEGFA
Female gonad development1.38E-2411CTNNA1, ESR1, GDF9, ICAM1, LEP, MYC, PDGFRA, RETN, S100A8, S100A9, VEGFA
Regulation of female gonad development1.38E-244GDF9, RETN, S100A8, S100A9

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

FIGURE 3

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 (). 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 (). S100A8 causes uteroplacental perfusion deficiency which leads to embryo abortion that supports the competence of our target selection (). 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 (). 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

Pocket IDVolume A2Surface A2Drug scoreSimple score
P_02,694.562,620.710.810.61
P_12,652.052,822.450.810.64
P_3209.65195.140.660
P_4180.41168.050.60
P_2252.39481.150.50.14
P_6166.91273.520.370
P_7165.23357.950.370.04
P_5173.21283.310.350
P_8137.56211.730.330
P_9130.58204.550.280
P_10127.43196.550.270
P_11120.91239.320.270
P_14109.89180.880.260
P_12117.99208.390.250
P_13116.19209.360.220
P_15109.89259.050.160
P_16106.85216.050.150
P_17100.78249.250.140

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 (). 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 (), but here, the designed is 18-mer with lower molecular weight so the plasma clearance may be faster.

FIGURE 4

TABLE 4

FragmentsPredicted ΔG (kcal/mol)FragmentsPredicted ΔG (kcal/mol)
Frag 1–16.24Frag 10–30.51
Frag 2–8.55Frag 11–14.34
Frag 3–16.52Frag 12–15.66
Frag 4–51.85Frag 13–15.17
Frag 5–18.46Frag 14–11.21
Frag 6–31.71Frag 15–18.70
Frag 7–16.54Frag 16–10.68
Frag 8–13.17Frag 17–21.76
Frag 9–38.76Frag 18–23.71

Binding energies of RNA analog fragments with S100A8.

TABLE 5

AptamersGC content (%)Tm (°C)Molecular weight (g/mol)ΔGmax (kcal/mol)
Apt1CAUCAUCAUGUUAUGUUC33.341.85327.4–27.93
Apt2AACAUCACAGAAGACAGA38.9375504.7–28.3
Apt3CUGACAACAUCAAUCAUG38.936.65395.5–29.04

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

ProteinTargetBinding siteGlobal energy
RAGES100A8Arg A 114 → Gln H44–25.75
RAGEGREArg B 203 → U14, Arg B 228 → G17–24.38
S100A8GREAsn D 61 → C6, Ala B 1 → U13–22.11
RAGEAptamer• Try B 118 → A13, Arg B 216 → A13–46.33
• Arg B 218 → G10, Asn B 25 → U16
• Gln B 24 → U16
S100A8Aptamer• Lys B 36 → U9, Lys F 48 → U17–45.32
• Ser H 86 → A13, Asp C 32 → G10
• Lys B 18 → G10, Lys B 21 → G10

Docking results.

FIGURE 5

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

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.

Statements

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.

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.

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.

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

  • 1

    AhirwarR.SmitaN.ShikhaA.SrinivasanR.SouvikM.PradipN. (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. 10.1038/srep21285

  • 2

    BarelierS.PonsJ.GehringK.LancelinJ.-M.KrimmI. (2010). Ligand specificity in fragment-based drug design.J. Med. Chem.5352565266. 10.1021/jm100496j

  • 3

    BerezovskiM. V.MusheevM. U.DrabovichA. P.JitkovaJ. V.KrylovS. N. (2006). Non-SELEX: selection of aptamers without intermediate amplification of candidate oligonucleotides.Nat. Protoc.1:1359. 10.1038/nprot.2006.200

  • 4

    BiesiadaM.Pachulska-WieczorekK.AdamiakR. W.PurzyckaK. J. (2016). RNAComposer and RNA 3D structure prediction for nanotechnology.Methods103120127. 10.1016/j.ymeth.2016.03.010

  • 5

    BindeaG.BernhardM.HubertH.PornpimolC.MarieT.AmosK.et al (2009). ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks.Bioinformatics2510911093. 10.1093/bioinformatics/btp101

  • 6

    CamoraniS.EspositoC. L.RienzoA.CatuognoS.IaboniM.CondorelliG.et al (2014). Inhibition of receptor signaling and of glioblastoma-derived tumor growth by a novel PDGFRbeta aptamer.Mol. Ther.22828841. 10.1038/mt.2013.300

  • 7

    CroftD.O’KellyG.WuG.HawR.GillespieM.MatthewsL.et al (2011). Reactome: a database of reactions, pathways and biological processes.Nucleic Acids Res.39D691D697. 10.1093/nar/gkq1018

  • 8

    DaanN. M.KosterM. P.de WildeM. A.DalmeijerG. W.EveleinA. M.FauserB. C.et al (2016). Biomarker profiles in women with PCOS and PCOS Offspring; a pilot study.PLoS One11:e0165033. 10.1371/journal.pone.0165033

  • 9

    DaiG.LuG. (2012). Different protein expression patterns associated with polycystic ovary syndrome in human follicular fluid during controlled ovarian hyperstimulation.Reprod Fertil. Dev.24893904. 10.1071/RD11201

  • 10

    GerashchenkoT. S.NikitaM. N.NadezhdaV. K.SofiaY. Z.MarinaV. Z.NadezhdaV. (2019). Markers of cancer cell invasion: are they good enough?J. Clin. Med.8:1092. 10.3390/jcm8081092

  • 11

    GhavamiS.RashediI.DattiloB. M.EshraghiM.ChazinW. J.HashemiM.et al (2008). S100A8/A9 at low concentration promotes tumor cell growth via RAGE ligation and MAP kinase-dependent pathway.J. Leukoc. Biol.8314841492. 10.1189/jlb.0607397

  • 12

    GhavamiS.SethC.MohammadH.MehdiE.WalterJ. C.AndrewJ. H.et al (2009). S100A8/A9: a janus-faced molecule in cancer therapy and tumorgenesis.Eur. J. Pharmacol.6257383. 10.1016/j.ejphar.2009.08.044

  • 13

    HallD.LiS.YamashitaK.AzumaR.CarverJ. A.StandleyD. M. (2015). RNA-LIM: a novel procedure for analyzing protein/single-stranded RNA propensity data with concomitant estimation of interface structure.Anal. Biochem.4725261. 10.1016/j.ab.2014.11.004

  • 14

    HsuK.PasseyR. J.EndohY.RahimiF.YoussefP. P.YenT.et al (2005). Regulation of S100A8 by glucocorticoids.J. Immunol.17423182326. 10.4049/jimmunol.174.4.2318

  • 15

    LuC.XiaoqinL.LinW.NingJ.JunY.XiaoboZ.et al (2017). Integrated analyses for genetic markers of polycystic ovary syndrome with 9 case-control studies of gene expression profiles.Oncotarget831703180. 10.18632/oncotarget.13881

  • 16

    MashiachE.Schneidman-DuhovnyD.AndrusierN.NussinovR.WolfsonH. J. (2008). FireDock: a web server for fast interaction refinement in molecular docking.Nucleic Acids Res.36W229W232. 10.1093/nar/gkn186

  • 17

    NioiP.McMahonM.ItohK.YamamotoM.HayesJ. 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), 337348. 10.1042/bj20030754

  • 18

    OughtredR.ChrisS.Bobby-JoeB.JenniferR.LorrieB.ChristieC.et al (2019). The BioGRID interaction database: 2019 update.Nucleic Acids Res.47D529D541. 10.1093/nar/gky1079

  • 19

    OwczarzyR.TataurovA. V.WuY.MantheyJ. A.McQuistenK. A.AlmabraziH. G.et al (2008). IDT SciTools: a suite for analysis and design of nucleic acid oligomers.Nucleic Acids Res.36W163W169. 10.1093/nar/gkn198

  • 20

    PaganoB.MartinoL.RandazzoA.GiancolaC. (2008). Stability and binding properties of a modified thrombin binding aptamer.Biophys. J.94562569. 10.1529/biophysj.107.117382

  • 21

    ParisienM.MajorF. (2008). The MC-Fold and MC-Sym pipeline infers RNA structure from sequence data.Nature4525155. 10.1038/nature06684

  • 22

    ReesB. B.BowmanJ. A.SchulteP. M. (2001). Structure and sequence conservation of a putative hypoxia response element in the lactate dehydrogenase-B gene of Fundulus.Biol. Bull.200247251. 10.2307/1543505

  • 23

    Schneidman-DuhovnyD.InbarY.NussinovR.WolfsonH. J. (2005). PatchDock and SymmDock: servers for rigid and symmetric docking.Nucleic Acids Res.33W363W367.

  • 24

    ShigdarS.QiaoL.ZhouS.-F.XiangD.WangT.LiY.et al (2013). RNA aptamers targeting cancer stem cell marker CD133.Cancer Lett.3308495. 10.1016/j.canlet.2012.11.032

  • 25

    Sir-PetermannT.MaliqueoM.AngelB.LaraH. E.Pérez-BravoF.RecabarrenS. E.et al (2002). Maternal serum androgens in pregnant women with polycystic ovarian syndrome: possible implications in prenatal androgenization.Hum. Reprod.1725732579. 10.1093/humrep/17.10.2573

  • 26

    StrimbuK.TavelJ. A. (2010). What are biomarkers?Curr. Opin. HIV AIDS5463466. 10.1097/COH.0b013e32833ed177

  • 27

    SzklarczykD.AnnikaL. G.DavidL.AlexanderJ.StefanW.JaimeH.et al (2019). STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets.Nucleic Acids Res.47D607D613. 10.1093/nar/gky1131

  • 28

    TanB. K.RaghuA.JingC.HendrikL.LouisJ. S. C.HarpalS. R.et al (2011). Metformin treatment exerts antiinvasive and antimetastatic effects in human endometrial carcinoma cells.J. Clin. Endocrinol. Metab.96808816. 10.1210/jc.2010-1803

  • 29

    TsengC. Y.AshrafuzzamanM.ManeJ. Y.KaptyJ.MercerJ. R.TuszynskiJ. A. (2011). Entropic fragment-based approach to aptamer design.Chem. Biol. Drug Des.78113. 10.1111/j.1747-0285.2011.01125.x

  • 30

    VolkamerA.KuhnD.RippmannF.RareyM. (2012). DoGSiteScorer: a web server for automatic binding site prediction, analysis and druggability assessment.Bioinformatics2820742075. 10.1093/bioinformatics/bts310

  • 31

    WeikumE. R.KnueselM. T.OrtlundE. A.YamamotoK. R. (2017). Glucocorticoid receptor control of transcription: precision and plasticity via allostery.Nat. Rev. Mol. Cell Biol.18159174. 10.1038/nrm.2016.152

  • 32

    YangZ.SugawaraM.PonathP. D.WessendorfL.BanerjiJ.LiY.et al (1990). Interferon gamma response region in the promoter of the human DPA gene.Proc. Natl. Acad. Sci. U.S.A.8792269230. 10.1073/pnas.87.23.9226

  • 33

    ZhouJ.LiH.LiS.ZaiaJ.RossiJ. J. (2008). Novel dual inhibitory function aptamer-siRNA delivery system for HIV-1 therapy.Mol. Ther.1614811489. 10.1038/mt.2008.92

Summary

Keywords

network analysis, druggability, RNA aptamer, lim method, pcos targets, protein network

Citation

Manibalan S, Shobana A, Kiruthika M, Achary A, Swathi M, Venkatalakshmi R, Thirukumaran K, Suhasini K and Roopathy S (2020) Protein Network Studies on PCOS Biomarkers With S100A8, Druggability Assessment, and RNA Aptamer Designing to Control Its Cyst Migration Effect. Front. Bioeng. Biotechnol. 8:328. doi: 10.3389/fbioe.2020.00328

Received

04 October 2019

Accepted

25 March 2020

Published

13 May 2020

Volume

8 - 2020

Edited by

Quan Zou, University of Electronic Science and Technology of China, China

Reviewed by

Ankush Sharma, University of Oslo, Norway; Qing Li, Huntsman Cancer Institute, University of Utah, United States

Updates

Copyright

*Correspondence: Subramaniyan Manibalan,

This article was submitted to Bioinformatics and Computational Biology, a section of the journal Frontiers in Bioengineering and Biotechnology

Disclaimer

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics