In-silico engineering of RNA nanoplatforms to promote the diabetic wound healing
BMC Chemistry volume 17, Article number: 52 (2023)
One of the most notable required features of wound healing is the enhancement of angiogenesis, which aids in the acceleration of regeneration. Poor angiogenesis during diabetic wound healing is linked to a shortage of pro-angiogenic or an increase in anti-angiogenic factors. As a result, a potential treatment method is to increase angiogenesis promoters and decrease suppressors. Incorporating microRNAs (miRNAs) and small interfering RNAs (siRNAs), two forms of quite small RNA molecules, is one way to make use of RNA interference. Several different types of antagomirs and siRNAs are now in the works to counteract the negative effects of miRNAs. The purpose of this research is to locate novel antagonists for miRNAs and siRNAs that target multiple genes to promote angiogenesis and wound healing in diabetic ulcers.
In this context, we used gene ontology analysis by exploring across several datasets. Following data analysis, it was processed using a systems biology approach. The feasibility of incorporating the proposed siRNAs and miRNA antagomirs into polymeric bioresponsive nanocarriers for wound delivery was further investigated by means of a molecular dynamics (MD) simulation study. Among the three nanocarriers tested (Poly (lactic-co-glycolic acid) (PLGA), Polyethylenimine (PEI), and Chitosan (CTS), MD simulations show that the integration of PLGA/hsa-mir-422a is the most stable (total energy = -1202.62 KJ/mol, Gyration radius = 2.154 nm, and solvent-accessible surface area = 408.416 nm2). With values of -25.437 KJ/mol, 0.047 nm for the Gyration radius, and 204.563 nm2 for the SASA, the integration of the second siRNA/ Chitosan took the last place. The results of the systems biology and MD simulations show that the suggested RNA may be delivered through bioresponsive nanocarriers to speed up wound healing by boosting angiogenesis.
Approximately 15–25% of type 2 diabetic patients suffer from diabetic foot ulcers . Destruction of the integumentary system (often dermis and epidermis) and subcutaneous tissue occurs in diabetic ulcers in multiple cases, which lead to amputation in 14–24% of affected patients . As a result of diabetic foot ulcers, lower-extremity amputation is a high-risk procedure with a high percentage of five-year mortality (50–59%), higher than many cancers . Due to the massive complications followed by diabetic ulcers, acceleration and improvement in diabetic ulcers healing is a research priority.
Enhancement of angiogenesis, which contributes to wound healing acceleration, is one of the distinguishing characteristics of tissue regeneration. Angiogenesis is regulated by multiple biochemical factors, including hypoxia-inducible-factor-1 (HIF-1), vascular-endothelial-growth-factor (VEGF), transforming-growth-factor (TGF), and fibroblast-growth-factors (FGF), and results in the development of new arteries . Inflammation, proliferation, and remodeling are the three stages of the wound healing process. When skin is damaged, inflammatory cells are driven to the site of the lesion, which is followed by fibroblast growth, and these cells are responsible for generating numerous tissue components . Normal wound healing requires angiogenesis stimulation; however, in diabetic wounds, hyperglycemia as a condition inhibits angiogenesis and transactivates hypoxia-inducible factor (HIF) . A capillary bed is developed in the wound area regarding angiogenic factors, during the proliferative phase of the wound healing. Interestingly, diabetes increases angiogenesis in diabetic retinopathy or nephropathy, although it declines angiogenesis in wound healing . Poor angiogenesis process during diabetic wound healing is associated either with insufficiency in pro-angiogenic or an increase in anti-angiogenic factors . Therefore, improvement of promoters or hallmarks of suppressors of angiogenesis could be used as a therapeutic strategy.
TGF-1 is a versatile cytokine that promotes granulation tissue development and collagen production during wound healing . Polyamines increased the expression of ODC, SSAT, hypoxia-inducible factors-1 (HIF-1), VEGF, and matrix metalloproteinases (MMOs) while decreasing the expression of p27 in Akt1-overexpression cells . The proliferation, differentiation, and migration of multiple cell types are all regulated by FGF2, a growth factor whose expression is decreased in diabetic and pressure ulcer tissues. For instance, wound healing, angiogenesis, bone regeneration and neuroregeneration may all benefit from FGF2’s participation in the formation of new blood vessels and tissue regeneration .
Ribonucleic acid interference (RNAi) is a ubiquitous biological response to double-stranded RNA that offers resistance to endogenous parasites and external harmful nucleic acids and controls protein-coding gene expression. Knocking out or down the genes is associated with numerous biological processes and is regarded as a crucial method in molecular research . The RNAi technology may help accomplish this goal. This naturally occurring technique to sequence-specific gene silencing has the potential to advance experimental biology. It has potential applications in functional genomics, therapeutic intervention, regenerative medicine, and other areas.
RNAi strategy is based on two types of small RNA molecules: microRNA (miRNA) and small interfering RNA (siRNA) . miRNAs are single strands (18–23 nucleotides), small non-coding, and endogenous products which bind to targeted miRNAs and degrade or block the translation of genes . Hence, miRNAs might have significant effects on several pathological and physiological processes. It has been shown that miRNAs upregulate and downregulate those genes associated with angiogenesis, cell proliferation, and collagen synthesis during the wound healing process . Due to the therapeutic potential, the application of anti-miRNAs or antagomirs may inhibit miRNAs’ undesirable effects, including angiogenesis suppression .
The primary function of miRNAs is to regulate gene expression. A miRNA is an RNA that is complementary to a portion of one or more mRNAs. The sequence in the 5′ proximal region between the 2nd and 8th nucleotides is where miRNAs mainly interact with their mRNA targets. It is worth noting that miRNA interacts with expression of target genes by repressing mRNA translation or degrading mRNA through deadenylation from the 3′ end and/or decapping from the 5′ end. According to a recent research, translation suppression regulates 48% of miRNA target genes, mRNA degradation influences 29%, and combined mechanisms control 23% .
Along with miRNAs, siRNAs play a significant role in the RNAi process leading to the knockdown of target genes. The process of gene silencing conducted by siRNAs is more specific than by miRNAs, whereas miRNAs usually affect multiple target genes . Synthetic siRNAs have emerged as prospective therapeutic agents for multiple issues, such as cancer, metabolic disorders, inflammatory disorders, and infectious diseases . Due to the different approaches, the application of siRNAs and miRNAs in pharmaceutical subjects may be considered a combined or parallel therapy.
Multiple therapeutic approaches to siRNAs have been developed and gained momentum in the past years [19, 20]. In the context of diabetic wound healing, siRNAs are designed to minimize undesired effects on healing processes through downregulation of genes such as P53 (increase apoptosis and impair vasculogenesis), matrix metalloproteinase (MMP9) (deteriorate of extracellular matrix), prolyl-hydroxylase domain 2 (PHD-2) (destabilize HIF-1) [21,22,23]. Therefore, utilizing siRNAs could improve and hasten cure for various phases of tissue regeneration and wound healing.
In terms of improving angiogenesis and tissue regeneration in diabetic foot ulcers, multiple miRNAs have been investigated. Diabetes foot ulcer lesions were shown to include extracellular vesicles expressing miR-195-5p and miR-205-5p, which inhibit angiogenesis and tissue regeneration in diabetes foot patients . In both in vivo and in vitro testing, the use of miRNA-497 as a therapy resulted in a reduction of pro-inflammatory cytokines such as IL-1, IL-6, and TNF-. The anti-inflammatory effects of miRNA-497 provide new information on the role it plays in facilitating the speedy recovery of diabetic wounds. As a result of its pro-inflammatory cytokine down-regulation capabilities, miRNA-497 has been indicated as a viable therapy option for diabetic wound healing .
The antiangiogenic microRNA 92a-3p (miR-92a) has been shown to increase angiogenesis in a variety of organ systems, including the heart, the limbs, the blood vessels, and the bones . It has been shown that a synthetic miR-92a inhibitor (MRG-110) can increase ITGA5 expression in vitro in both human vascular endothelial cells and primary human skin fibroblasts, and in mice skin . According to these results, MRG-110 has the potential to improve contextual chronic wound healing as well as acute wound healing .
On the other hand, the basis of biological activity is molecular interactions , which are the product of macromolecular structures . Simulations of molecular dynamics (MD) have grown into a sophisticated approach for determining the links between macromolecular structure and function . Also, physiologically-relevant periods may be compared to the simulation process timeframes. For structural bioinformatics, the knowledge on dynamic macromolecule features is sufficient to transform the paradigm from single-structure studies to conformational ensemble evaluation . MD simulation, capable of modeling the biological macromolecules such as proteins and genes, and their integration with various materials, including polymeric nanocarriers, discusses potential therapeutic approaches for multiple diseases.
To date, multiple antagomirs and siRNAs have been presented to reduce the undesirable effects of miRNAs. To the best of our knowledge, there is not a study that could suggest antagomirs and siRNAs, and discover their delivery via employing bioresponsive nanocarriers at the wound site. Almost previous studies lack the survey of employment possibility of designed and dedicated small molecules in diabetic wound healing. This paper aims to present new antagomirs for miRNAs and siRNAs, which suppress the anti-angiogenesis process by targeting multiple genes and improving the wound-healing process. For this object, the gene ontology algorithms were used by utilizing several databases. Then, a systems biology approach was applied to the whole gathered data. Also, to deliver the designed siRNAs and miRNAs’ antagomirs to the wound site, their integration with three types of various polymeric nanocarriers was considered through an MD simulation study. Figure 1 summarizes the whole study procedure.
Materials and methods
Gene ontology analysis
A systems biology investigation was carried out in this paper, which made use of the regeneration gene database as well as the DisGeNET database (v7.0). All procedures were carried out in conformity with the applicable guidelines and regulations. The regeneration gene database contains genes involved in the regeneration of human tissues and/or organs, while DisGeNET is a discovery platform that contains one of the biggest publicly accessible collections of genes and variations related with human disorders [31, 32].
Using the database for annotation, visualization and integrated discovery (DAVID) database, each of the obtained genes was subjected to an analysis designed to enrich individual gene ontology (GO) terms in order to determine whether or not co-transcriptionally regulated genes share any activities related to wound healing [33, 34]. DAVID distinguishes itself from comparable databases by providing an integrated and enlarged back-end annotation database, innovative modular enrichment techniques, and strong exploratory capability in an integrated data-mining context. One of the most critical elements that directly influences the performance of DAVID’s functional analysis is the quality of massive gene lists obtained from high-throughput biological investigations .
The genes resulted from DAVID entered in another online tool, enriching for GO terms, which is called GeneMANIA  to identify potential partners for angiogenesis, epithelial cell proliferation, regulation of epithelial cell migration, response to fibroblast growth factor, fibroblast migration, epithelial cell migration, and positive regulation of cell migration. A list of genes was loaded into GeneMANIA, and a list of data sets to query was chosen. GeneMANIA then expands the list with functionally related genes or genes that share attributes with the original query genes, displaying a dynamic functional correlation matrix that illustrates the links between the genes and datasets. Another notable feature of GeneMANIA is that it applies weights to data sets depending on their utility for each query .
To confirm the results obtained via the GeneMANIA database, pathway analysis and GO analysis were accomplished using the protein analysis through evolutionary relationships (Panther)  and Enrichr [39, 40] databases. The Panther categorization process is a systematic maintained biological database of gene/protein families and functionally related subfamilies that may be used to categorize and identify gene products’ functions. Panther’s practical application is to correctly estimate the function of uncharacterized genes from any organism based on evolutionary links to genes with known functions . Enrichr, a web-based gadget that provides a variety of graphical summaries of the aggregate functions of gene lists, is very straightforward to use .
Determining appropriate miRNAs
Enrichr database identified miRNAs, which potentially could target provided genes via GeneMANIA and Panther. To confirm the Enrichr database findings, they were further investigated using the miRWalk database , which is an open-source platform with an uncomplicated interface for producing predicted and verified miRNA-binding sites of known genes in humans, mice, rats, dogs, and cows . Verified miRNAs were entered into the mirPath v.3 database to identify the most critical pathways, interfered with by these miRNAs and have a crucial role in wound healing .
To inhibit verified miRNAs, antisense RNAs were designed utilizing the miRNAsong tool . miRNAsong is a web-based application for creating and testing miRNA sponges in a virtual environment. miRNAsong creates miRNA sponge constructs for certain miRNAs and miRNA families or clusters and then tests these constructs to see whether they have the ability to bind to miRNAs in a variety of species . Finally, to identify interactions between selected miRNAs and TFs, the transmiR database was used [45, 46]. It is a significant resource for investigating TF–miRNA regulation and could be used to evaluate a wide range of processes, such as the development of the relationships, expression patterns, and related disorders of miRNAs .
The DAVID database was utilized to identify essential genes in the negative regulation of the apoptotic process. The findings were examined using the GeneMANIA database to identify the most significant genes in anti-apoptotic activities such as apoptotic signaling pathway negative regulation, intrinsic apoptotic signaling pathway negative regulation, and endothelial cell apoptotic process negative regulation. The Enrichr database was utilized to verify the acquired findings.
In this regard, the sequences of NFE2L2 and SIRT1 genes were collected from NCBI nucleotide . To design siRNAs, Eurofins Genomic’s siRNA design tool was used, which is an online gadget providing the use of guidelines initially provided by Tuschl et al. . In addition, it applies the scorings criteria  (including target site, length of siRNA, specificity checking, and nucleotide content of siRNA) developed by Reynolds et al. , Ui-Tei et al. , and Vert et al. . To ensure the safety and prevent off-target effects of the designed siRNAs, the sequences were aligned by the human genome utilizing align sequences nucleotide BLAST .
Molecular dynamics simulations
An MD simulation was performed using the all-atom technique (AA-MD) considering the PLGA, PEI and Chitosan nanocarriers, utilizing the OPLSA force field. Mentioned nanocarriers were designed using the Pepdraw and rcsb online tools and were parameterized with the gmx pdb2 gmx command in GROMAX 2020. Parameters like length, mass, net charge, and hydrophobicity were also provided utilizing Pepdraw tool and molecular dynamics add-ons, such as Avogadro. Polymeric nanocarriers are also designed by Avogadro software, then the parameter is pooled using PolyParGen online tool, while their esp charge were calculated by cp2k software. 18 simulations were performed with a duration of 120 nanoseconds and a time step of 2 femtoseconds. The simulations are then performed at the temperature and pH of 30.91 °C and 6.95, respectively . All the simulation was performed in 4 steps using the coordinate file of the molecules. Initially, the simulation system was placed at an energy level of 100 (kJ/mol) after 100 ps. Next, using the Berendsen thermostat algorithm, the simulation system was placed at 30.91 °C after 1 ns. After the NVT (constant number of atoms, N; constant volume, V; constant temperature, T) stage, using the Parrinello-rahman algorithm and the simulation system equilibrates after 10 ns at a pressure of 1 bar. Finally, the simulation was performed using the LINCS algorithm at 50 ns, the cutoff radius was set to 3 nm and H-bonds were considered for all simulations. The tip3p water model was used for simulation. The simulations are performed with 1080 ti graphics. For the analyzes, the commands gmx gyrate, mmpbsa, gmx sasa were used, and for imaging the simulation system, vmd software was used. For molecular dynamics, GROMACS 2020 was used. Gmx gyrate calculate the gyration radius during the simulations using Eq. (1):
where n is the number of particles, mi, the mass of the particle i, and ri indicate the distance of the particle i from the center of gravity. Rg, helps to understand protein density during the simulation, where for higher Rg, higher density is expected. Similarly, gmx sasa calculates the solvent accessible surface area of the particles in the course of simulations, via the double cubic lattice method .
Results and discussion
Gene ontology analysis and RNAis proposing
All 72 genes that are participated in the wound occurring were gathered from the DisGeNET database21, while all the 26 protein-coding genes that are participated in wound healing were gathered from the regeneration gene database  (Table 1). Analyzing all these 98 genes via the DAVID database , resulted in identifying 15 influential genes possessing significant roles in wound generating and healing (Table S1). The GeneMANIA database reveals that six genes collaborate in angiogenesis and proliferation processes, including AKT1, HIF1A, FGF2, IGF1, TGFB1, and PIK3CA (Fig. 2). Validating obtained results via the Panther database , reveals four genes possessing significant roles in the angiogenesis process, including AKT1, HIF1A, FGF2, and PIK3CA, while the Enrichr database verified the obtained results (Table S2).
Using the Enrichr database, it could be possible to identify miRNAs that target the aforementioned genes (Table S3 and S4). Furthermore, since miRTarBase is an empirically confirmed microRNA-target interactions database , it has about 360.000 miRNA-target interactions (MTIs). The seed region of each miRNA is also used to determine the biological targets of miRNAs using TargetScan. Three selected miRNAs were found by checking them against the miRWalk database (Table 2). According to the mirPath v.3 database, hsa-mir-199a-3p and hsa-mir-422a have the potential to interfere with the hippo signaling pathway . This route is responsible for controlling organ development, tissue regeneration, and stem cell self-renewal. YAP and TAZ, two downstream transcription coactivators with PDZ-binding motifs, govern the Hippo pathway’s essential gene regulation and biological functions .
The antisense RNA molecule is a distinct form of DNA transcript that consists of 19–23 nucleotides and complements miRNA. Antisense RNAs regulate gene expression at numerous levels, including replication, transcription, and translation . As a result, antisense miRNAs were created with the use of the microRNA sponge generator and tester (miRNAsong) program . Furthermore, based on DAVID database findings, transcription factors (TFs) targeted mir-199a-3p, mir-4757-3p, and mir-422a, which play a critical role in wound healing, according to the transmiR database (Figure S1). TFs directly interpret the genome, conducting the initial step of DNA sequence decoding . TFs control turning genes on and off to ensure that they are expressed in the appropriate cell at the correct time and in the proper quantity throughout the cell’s and organism’s life . Based on the findings, three miRNA antisenses, hsa-mir-422a, hsa-mir-199a-3p, and hsa-mir-4757-3p, have been identified as targeting PIK3CA, FGF2, and HIF1A, respectively, while all of them target AKT1.
Wound healing is regulated by a huge proportion of cytokines and biochemicals. As a result, diabetic foot ulcer regeneration necessitates the coordinated activities of a series of diverse cell types . Wounds of multiple kinds, including diabetic foot, chronic venous, and pressure ulcers, all respond favorably to the healing actions of a variety of growth factors . FGFs are involved in a wide variety of developmental abnormalities, including neurodevelopment, keratinocyte organization, angiogenesis, wound healing, and any functional aberration that leads to a diversity of developmental illnesses. Angiogenesis is promoted by FGFs, which are more potent angiogenesis factors than PDGF and VEGF. As a result, granulation tissue may form as a result of angiogenesis as well as fibroblast proliferation .
When the expression of miRNAs linked with a given signaling pathway is reduced, the expression of signaling proteins participating in that system increases noticeably. The subsequent effects include increased cell proliferation and migration, as well as a lowered rate of apoptosis. In diabetes mellitus, reduced cell proliferation and migration, as well as accelerated death, are correlated with increased levels of miRNA proteins that are implicated in the PI3K/mTOR signaling pathway. Targeting miRNAs, which may target genes implicated in these pathways, is likely to promote proliferation and migration . Moreover, Jing et al.  discovered that diabetic mice with low AKT/HIF-1 signaling activity had a substantial reduction in wound healing and angiogenesis.
In terms of successful siRNA design, the DAVID database was used to identify 17 critical genes in the negative regulation of the apoptotic process (Table S5). The GeneMANIA database’s GO analysis revealed the five most influential genes, including AKT1, SCG2, NFE2L2, CXCL12, and SIRT1, that are highly collaborative in the negative regulation of apoptotic processes, including negative regulation of apoptotic signaling pathway, negative regulation of intrinsic apoptotic signaling pathway, and negative regulation of endothelial cell apoptotic process (Fig. 3). The Enrichr database was utilized to confirm the reported findings, which resulted in two genes being blocked by siRNAs. The Enrichr database confirmed that the NFE2L2 and SIRT1 genes play an essential part in the negative regulation of the apoptotic process; hence, these genes are the best candidates for siRNA inhibition. As a result, siRNAs were developed in the Eurofins Genomics database using the siMAX siRNA design tool (Tables S6, S7, and S8) . Furthermore, nucleotide BLAST confirmed three siRNAs that had the least level of off-targeting (Table S9).
Coagulation, inflammation, tissue formation or proliferating, and tissue remodeling or scar formation are the four permanent and temporally overlapping stages of wound healing . During the apoptotic phase, neutrophils are absorbed by macrophages to prevent further inflammation . Apoptosis in immune cells may be the key to ending inflammation and kicking off healing since macrophages are the major phagocytes that help resolve inflammation and promote tissue regeneration . Diabetes reduces the number and activity of macrophages, which in turn reduces the development of lymphatic vessels and hence slows the healing process. NFE2L2 and SIRT1 are two key genes that regulate apoptosis and might be targeted to improve tissue regeneration .
Molecular dynamics (MD) simulation
Based on input data provided in Table 3, MD simulation gives the parameters of length, mass, net charge, and hydrophobicity of each RNAis. The pH and temperature at the diabetic wound area were considered as 6.95 and 30.91 °C, respectively . The key parameters obtained from the software show appropriately the details of the structures of siRNA, miRNA and nanocarriers. Table 3 lists the full names and specifications of each structure.
Table S10 also shows the basic parameters for all six structures of siRNA and miRNA antagomirs, such as the mass and net charge of the molecules, so that they can be used to interpret the output analyses obtained from the simulation of MD. Moreover, Table 4 presents the average amount of Van der Waals and electrostatic energies, total energy, Gyration radius, and solvent-accessible surface area (SASA) for each nanocarrier, miRNAs, and siRNAs.
The van der Waals energies are related to the interactions among atoms, which vary according to the size and mass of each atom. This means that in this analysis, the higher the mass and number of atoms of a molecule, the more van der Waals energy it can have. The calculation of these energies in MD follows the Lennard-Jones equation . Table 4 clearly shows that miRNAs exhibit stronger van der Waals interactions than siRNAs. In general, according to the mass data shown in Table S10, the higher the mass of the molecules, the more van der Waals energy they could have, depending on the type of nanocarrier.
Electrostatic energies are also a function of the net charge of the atoms of each molecule, meaning that molecules that have higher and unlike charges have more electrostatic energy. Molecules that have like charges also possess positive electrostatic energy. In MD, electrostatic energy also follows Coulomb’s law . Therefore, a comparison could be available for these simulations according to the net charge of siRNA and miRNA molecules reported in Table 4. From this comparison, it could be concluded that the net charge of miRNAs is higher than that of siRNAs, and the ratio of their electrostatic energy to the molecular charge of these structures has changed.
Based on the electrostatic energy, Poly (lactic-co-glycolic acid) (PLGA), Polyethylenimine (PEI), and Chitosan nanocarriers were able to have better electrostatic absorption for structures, respectively. The electrostatic adsorption energy is positive for the Chitosan nanocarrier, and this is an unsuitable parameter for the adsorption between the Chitosan nanocarrier and the RNA structures. The total energy is also calculated from the sum of van der Waals and electrostatic energies and is a complete indicator for calculating the amount of adsorption among RNA structures and nanocarriers. Based on the total energy, the PLGA nanocarriers had better adsorption than PEI and Chitosan. Also, miRNAs 1 and 3 had the best adsorption; hence, the highest adsorption energy was dedicated to the adsorption among PLGA and miRNA No. 1, as well as PLGA and miRNA No. 3.
These results are also validated through the analysis of Gyration radius and SASA. The average Gyration radius for the simulations also shows a better aggregation of RNA structures and nanocarriers. The lower value results in more significant aggregation and adsorption. This analysis positively correlates to the results obtained from the total energy analyses.
SASA analysis also shows the surface in contact with the solvent for the structures. Changes in this analysis are also done according to the size of the RNA structures and nanocarriers and the degree of hydrophilicity and hydrophobicity of the molecules. Figure 4 A-F show different stages of simulation for nanocarriers and RNAi structures. It is possible to compare the amount of adsorption in various simulation phases, including input, middle, and output. The analyses were pursued through output phases regarding the superior influence of output mode over input and middle phases in MD simulations. Figure 4G-L demonstrate the output interaction modes of PLGA nanocarriers and RNAis.
The output interaction modes of PEI and Chitosan nanocarriers with RNAi structures are available in Fig. 4M-X. The integration of RNAis and polymeric nanocarriers in these figures investigates that PEI is more integrated with miRNA antagomirs and siRNAs than Chitosan. The comparison of Figs. 4G-L and 5 led to discuss the stability, van der Waals, electrostatic, and total energies in MD simulations. Figure 5 A demonstrates the van der Waals interaction energy of PLGA, PEI, and Chitosan with miRNA antisenses and siRNAs, introducing the integration of PLGA with miRNA antagomir 1 as the optimum delivery system in terms of van der Waals interaction energy. Furthermore, Fig. 5B indicates electrostatic energy, in which the PLGA and PEI are in a similar negative range, while Chitosan is in the positive area, with a severe discrepancy. Finally, the total energy is shown in Fig. 5C, which summarizes all the given analyses in the energy area. Based on this chart, the integration of PLGA is more desirable than other polymeric nanocarriers, while delivering the miRNA antagomirs is more stable than siRNAs to the diabetic ulcer areas.
To confirm the obtained results, analyses of the Gyration radius and SASA are performed and presented in Fig. 5D and E, respectively. The Gyration radius of the PLGA/ miRNA antagomir 1 combination (2.154 nm) is higher than others, while its higher SASA amount (408.416 nm^2) confirms that. Therefore, it would be better to investigate its root-mean-square deviation (RMSD) along the MD simulation (Fig. 5F). RMSD is a prominent stability measure in MD, representing the particle position deviation concerning their reference position . The RMSD results positively correlate with Gyration radius analysis regarding the hydrophilicity performance. Finally, Table S11 summarizes the results according to the MD simulation output analyzes of RNAis and nanocarriers, considering the priority sorted from the highest adsorption to the lowest one.
A systems biology approach was employed in the current work to investigate how some RNAis promote diabetic ulcer regeneration by improving angiogenesis. Moreover, to discover the stable delivery of proposed miRNAs and siRNAs into the wound area, their aggregation with three bioresponsive polymeric nanocarriers (i.e., PLGA, PEI, and Chitosan) were studied through MD simulation. Multiple self-assembly parameters and stable integration were examined, including van der Waals energy, electrostatic energy, total energy, Gyration radius, SASA, and RMSD. Results demonstrate that PLGA, PEI, and Chitosan possess the lowest energy, leading to more stability and safe delivery, respectively. Also, results show that miRNAs prioritize utilize as the therapeutic agents compared with siRNAs. It is proposed that as poor angiogenesis during diabetic wound healing is linked to either a lack of pro-angiogenic factors or an increase in anti-angiogenic factors, improving promoters and decreasing suppressors of the angiogenesis process could accelerate wound healing in diabetic ulcers. The results from the systems biology and MD studies led us to choose utilizing PLGA nanocarrier to deliver the hsa-mir-422a antagomir as the most stable condition to recommend for further development and improvement of angiogenesis applications.
The datasets used and/or analyzed during the current study are available from the corresponding author (Nima Beheshtizadeh) on reasonable request.
Database for annotation, visualization and integrated discovery
Fibroblast growth factors
Prolyl-hydroxylase domain 2
Protein analysis through evolutionary relationships
Small interfering RNA
Solvent-accessible surface area
Transcriptional coactivator with PDZ-binding motif
- TGF- β:
Transforming growth factor β
Vascular endothelial growth factor
Yazdanpanah L, Shahbazian H, Nazari I, Arti HR, Ahmadi F, Mohammadianinejad SE, Cheraghian B, Hesam S. Incidence and Risk Factors of Diabetic Foot Ulcer: A Population-Based Diabetic Foot Cohort (ADFC Study)—Two-Year Follow-Up Study, International Journal of Endocrinology, 2018 (2018) 7631659.
Alavi A, Sibbald RG, Mayer D, Goodman L, Botros M, Armstrong DG, Woo K, Boeni T, Ayello EA, Kirsner RS. Diabetic foot ulcers: part I. Pathophysiology and prevention. J Am Acad Dermatol. 2014;70(1):e1–1. e18.
Armstrong DG, Swerdlow MA, Armstrong AA, Conte MS, Padula WV, Bus SA. Five year mortality and direct costs of care for people with diabetic foot complications are comparable to cancer. J Foot Ankle Res. 2020;13:1–4.
Tonnesen MG, Feng X, Clark RA. Angiogenesis in wound healing, Journal of Investigative Dermatology Symposium Proceedings, Elsevier, 2000, pp. 40–46.
Meng Z, Zhou D, Gao Y, Zeng M, Wang W. miRNA delivery for skin wound healing. Adv Drug Deliv Rev. 2018;129:308–18.
Chaudhary A, Bag S, Banerjee P, Chatterjee J. Wound healing efficacy of Jamun honey in diabetic mice model through reepithelialization, collagen deposition and angiogenesis. J Traditional Complement Med. 2020;10:529–43.
Okonkwo UA, DiPietro LA. Diabetes and wound angiogenesis. Int J Mol Sci. 2017;18:1419.
Martin A, Komada MR, Sane DC. Abnormal angiogenesis in diabetes mellitus. Med Res Rev. 2003;23:117–45.
Galehdari H, Negahdari S, Kesmati M, Rezaie A, Shariati G. Effect of the herbal mixture composed of Aloe Vera, Henna, Adiantum capillus-veneris, and Myrrha on wound healing in streptozotocin-induced diabetic rats. BMC Complement Altern Med. 2016;16:386.
Dai F, Yu W, Song J, Li Q, Wang C, Xie SJA-CD. Extracellular polyamines-induced proliferation and migration of cancer cells by ODC, SSAT, and Akt1-mediated pathway, 28 (2017) 457–64.
Decker CG, Wang Y, Paluck SJ, Shen L, Loo JA, Levine AJ, Miller LS, Maynard HD. Fibroblast growth factor 2 dimer with superagonist in vitro activity improves granulation tissue formation during wound healing. Biomaterials. 2016;81:157–68.
Fakhr E, Zare F, Teimoori-Toolabi L. Precise and efficient siRNA design: a key point in competent gene silencing. Cancer Gene Ther. 2016;23:73–82.
Banerjee J, Sen CK. microRNA and wound healing, microRNA: Medical Evidence, Springer2015, pp. 291–305.
Bhattacharya S, Aggarwal R, Singh VP, Ramachandran S, Datta M. Downregulation of miRNAs during delayed Wound Healing in Diabetes: role of Dicer. Mol Med. 2016;21:847–60.
Veith AP, Henderson K, Spencer A, Sligar AD, Baker AB. Therapeutic strategies for enhancing angiogenesis in wound healing. Adv Drug Deliv Rev. 2019;146:97–125.
Soliman AM, Das S, Abd Ghafar N, Teoh SL. Role of MicroRNA in Proliferation Phase of Wound Healing, 9 (2018).
Lam JK, Chow MY, Zhang Y, Leung SW. siRNA versus miRNA as therapeutics for gene silencing. Mol Therapy-Nucleic Acids. 2015;4:e252.
Potaczek DP, Garn H, Unger SD, Renz H. Antisense molecules: a new class of drugs. J Allergy Clin Immunol. 2016;137:1334–46.
López-Fraga M, Martínez T, Jiménez A. RNA interference technologies and therapeutics: from basic research to products. BioDrugs. 2009;23:305–32.
Alshaer W, Zureigat H, al Karaki A, Al-Kadash A, Gharaibeh L, Hatmal Mm, Aljabali A, Awidi A. siRNA: Mechanism of action, challenges, and therapeutic approaches, European Journal of Pharmacology, 905 (2021) 174178.
Li N, Luo H-C, Ren M, Zhang L-M, Wang W, Pan C-L, Yang L-Q, Lao G-J, Deng J-J. -j. Mai, Efficiency and safety of β-CD-(D3) 7 as siRNA carrier for decreasing matrix metalloproteinase-9 expression and improving wound healing in diabetic rats. ACS Appl Mater Interfaces. 2017;9:17417–26.
Nguyen PD, Tutela JP, Thanik VD, Knobel D, Allen J, Robert J, Chang CC, Levine JP, Warren SM, Saadeh PB. Improved diabetic wound healing through topical silencing of p53 is associated with augmented vasculogenic mediators. Wound Repair and Regeneration. 2010;18:553–9.
Shaabani E, Sharifiaghdam M, Lammens J, De Keersmaecker H, Vervaet C, De Beer T, Motevaseli E, Ghahremani MH, Mansouri P, De Smedt S. Increasing angiogenesis factors in Hypoxic Diabetic Wound Conditions by siRNA delivery: Additive Effect of LbL-Gold Nanocarriers and Desloratadine-Induced lysosomal escape. Int J Mol Sci. 2021;22:9216.
Liu J, Wang J, Fu W, Wang X, Chen H, Wu X, Lao G, Wu Y, Hu M, Yang C, Yan L, Ren M. MiR-195-5p and mir-205-5p in extracellular vesicles isolated from diabetic foot ulcer wound fluid decrease angiogenesis by inhibiting VEGFA expression. Aging. 2021;13:19805–21.
Ban E, Jeong S, Park M, Kwon H, Park J, Song EJ, Kim A. Accelerated wound healing in diabetic mice by miRNA-497 and its anti-inflammatory activity. Biomed Pharmacother. 2020;121:109613.
Gallant-Behm CL, Piper J, Dickinson BA, Dalby CM, Pestano LA, Jackson AL. A synthetic microRNA-92a inhibitor (MRG-110) accelerates angiogenesis and wound healing in diabetic and nondiabetic wounds. Wound Repair Regen. 2018;26:311–23.
Gallant-Behm CL, Piper J, Dickinson BA, Dalby CM, Pestano LA, Jackson AL. A synthetic microRNA-92a inhibitor (MRG-110) accelerates angiogenesis and wound healing in diabetic and nondiabetic wounds. Wound Repair and Regeneration. 2018;26:311–23.
Maximova T, Moffatt R, Ma B, Nussinov R, Shehu A. Principles and overview of sampling methods for modeling Macromolecular structure and Dynamics. PLoS Comput Biol. 2016;12:e1004619.
Ferreira LG, Dos Santos RN, Oliva G, Andricopulo AD. Molecular Docking and Structure-Based Drug Design Strategies, Molecules, 20 (2015).
Hospital A, Goñi JR, Orozco M, Gelpí JL. Molecular dynamics simulations: advances and applications. Adv Appl Bioinform Chem. 2015;8:37–47.
Piñero J, Queralt-Rosinach N, Bravo À, Deu-Pons J, Bauer-Mehren A, Baron M, Sanz F, Furlong LI. DisGeNET: a discovery platform for the dynamical exploration of human diseases and their genes, Database (Oxford), 2015 (2015) bav028.
Piñero J, Ramírez-Anguita JM, Saüch-Pitarch J, Ronzano F, Centeno E, Sanz F, Furlong LI. The DisGeNET knowledge platform for disease genomics: 2019 update. Nucleic Acids Res. 2020;48:D845–d855.
Huang da W, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37:1–13.
Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.
Sherman BT, Huang DW, Tan Q, Guo Y, Bour S, Liu D, Stephens R, Baseler MW, Lane HC, Lempicki RA. DAVID Knowledgebase: a gene-centered database integrating heterogeneous gene annotation resources to facilitate high-throughput gene functional analysis. BMC Bioinformatics. 2007;8:426.
Warde-Farley D, Donaldson SL, Comes O, Zuberi K, Badrawi R, Chao P, Franz M, Grouios C, Kazi F, Lopes CT, Maitland A, Mostafavi S, Montojo J, Shao Q, Wright G, Bader GD, Morris Q. The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. 2010;38:W214–220.
Warde-Farley D, Donaldson SL, Comes O, Zuberi K, Badrawi R, Chao P, Franz M, Grouios C, Kazi F, Lopes CT, Maitland A, Mostafavi S, Montojo J, Shao Q, Wright G, Bader GD, Morris Q. The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. 2010;38:W214–20.
Mi H, Ebert D, Muruganujan A, Mills C, Albou LP, Mushayamaha T, Thomas PD. PANTHER version 16: a revised family classification, tree-based classification tool, enhancer regions and extensive API. Nucleic Acids Res. 2021;49:D394–d403.
Chen EY, Tan CM, Kou Y, Duan Q, Wang Z, Meirelles GV, Clark NR. Ma’ayan, Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinformatics. 2013;14:128.
Kuleshov MV, Jones MR, Rouillard AD, Fernandez NF, Duan Q, Wang Z, Koplev S, Jenkins SL, Jagodnik KM, Lachmann A, McDermott MG, Monteiro CD, Gundersen GW. Ma’ayan, Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 2016;44:W90–97.
Thomas PD, Kejariwal A, Campbell MJ, Mi H, Diemer K, Guo N, Ladunga I, Ulitsky-Lazareva B, Muruganujan A, Rabkin S, Vandergriff JA, Doremieux O. PANTHER: a browsable database of gene products organized by biological function, using curated protein family and subfamily classification. Nucleic Acids Res. 2003;31:334–41.
Sticht C, De La Torre C, Parveen A, Gretz N. miRWalk: an online resource for prediction of microRNA binding sites. PLoS ONE. 2018;13:e0206239.
Vlachos IS, Zagganas K, Paraskevopoulou MD, Georgakilas G, Karagkouni D, Vergoulis T, Dalamagas T, Hatzigeorgiou AG. DIANA-miRPath v3.0: deciphering microRNA function with experimental support. Nucleic Acids Res. 2015;43:W460–466.
Barta T, Peskova L, Hampl A. miRNAsong: a web-based tool for generation and testing of miRNA sponge constructs in silico. Sci Rep. 2016;6:36625.
Tong Z, Cui Q, Wang J, Zhou Y. TransmiR v2.0: an updated transcription factor-microRNA regulation database. Nucleic Acids Res. 2019;47:D253–d258.
Wang J, Lu M, Qiu C, Cui Q. TransmiR: a transcription factor-microRNA regulation database. Nucleic Acids Res. 2010;38:D119–122.
Wang J, Lu M, Qiu C, Cui Q. TransmiR: a transcription factor–microRNA regulation database. Nucleic Acids Res. 2010;38:D119–22.
Tuschl T, Zamore PD, Lehmann R, Bartel DP, Sharp PA. Targeted mRNA degradation by double-stranded RNA in vitro. Genes Dev. 1999;13:3191–7.
Reynolds A, Leake D, Boese Q, Scaringe S, Marshall WS, Khvorova A. Rational siRNA design for RNA interference. Nat Biotechnol. 2004;22:326–30.
Ui-Tei K, Naito Y, Takahashi F, Haraguchi T, Ohki-Hamazaki H, Juni A, Ueda R, Saigo K. Guidelines for the selection of highly effective siRNA sequences for mammalian and chick RNA interference. Nucleic Acids Res. 2004;32:936–48.
Vert J-P, Foveau N, Lajaunie C, Vandenbrouck Y. An accurate and interpretable model for siRNA efficacy prediction. BMC Bioinformatics. 2006;7:520.
Gethin G, O’Connor GM, Abedin J, Newell J, Flynn L, Watterson D, O’Loughlin A. Monitoring of pH and temperature of neuropathic diabetic and nondiabetic foot ulcers for 12 weeks: an observational study. Wound Repair Regen. 2018;26:251–6.
Eisenhaber F, Lijnzaad P, Argos P, Sander C, Scharf M. The double cubic lattice method: efficient approaches to numerical integration of surface area and volume and to dot surface contouring of molecular assemblies. J Comput Chem. 1995;16:273–84.
Zhao M, Rotgans B, Wang T, Cummins SF. REGene: a literature-based knowledgebase of animal regeneration that bridge tissue regeneration and cancer. Sci Rep. 2016;6:23167.
Huang HY, Lin YC, Li J, Huang KY, Shrestha S, Hong HC, Tang Y, Chen YG, Jin CN, Yu Y, Xu JT, Li YM, Cai XX, Zhou ZY, Chen XH, Pei YY, Hu L, Su JJ, Cui SD, Wang F, Xie YY, Ding SY, Luo MF, Chou CH, Chang NW, Chen KW, Cheng YH, Wan XH, Hsu WL, Lee TY, Wei FX, Huang HD. miRTarBase 2020: updates to the experimentally validated microRNA-target interaction database. Nucleic Acids Res. 2020;48:D148–d154.
Zhao B, Tumaneng K, Guan K-L. The Hippo pathway in organ size control, tissue regeneration and stem cell self-renewal. Nat Cell Biol. 2011;13:877–83.
Lee M-J, Byun MR, Furutani-Seiki M, Hong J-H. Jung, YAP and TAZ regulate skin Wound Healing. J Invest Dermatology. 2014;134:518–25.
Xu J-z, Zhang J-l, Zhang W-g. Antisense RNA: the new favorite in genetic research. J Zhejiang University-SCIENCE B. 2018;19:739–49.
Lambert SA, Jolma A, Campitelli LF, Das PK, Yin Y, Albu M, Chen X, Taipale J, Hughes TR, Weirauch MT. Hum Transcription Factors Cell. 2018;172:650–65.
Park JW, Hwang SR, Yoon I-S. Advanced Growth factor Delivery Systems in Wound Management and skin regeneration. Molecules; 2017. p. 22.
Liu Y, Liu Y, Deng J, Li W, Nie X. Fibroblast growth factor in Diabetic Foot Ulcer: progress and therapeutic prospects. Front Endocrinol. 2021;12:1348.
Jere SW, Houreld NN, Abrahamse H. Role of the PI3K/AKT (mTOR and GSK3β) signalling pathway and photobiomodulation in diabetic wound healing. Cytokine Growth Factor Rev. 2019;50:52–9.
Jing L, Li S, Li Q. Akt/hypoxia-inducible factor-1α signaling deficiency compromises skin wound healing in a type 1 diabetes mouse model. Exp Ther Med. 2015;9:2141–6.
Demidova-Rice TN, Hamblin MR, Herman IM. Acute and impaired wound healing: pathophysiology and current methods for drug delivery, part 1: normal and chronic wounds: biology, causes, and approaches to care. Adv Skin Wound Care. 2012;25:304–14.
Brubaker AL, Schneider DF, Kovacs EJ. Neutrophils and natural killer T cells as negative regulators of wound healing. Expert Rev Dermatology. 2011;6:5–8.
Li L, Yan B, Shi Y-Q, Zhang W-Q, Wen Z-L. Live imaging reveals differing roles of macrophages and neutrophils during zebrafish tail fin regeneration *. J Biol Chem. 2012;287:25353–60.
Arya AK, Tripathi R, Kumar S, Tripathi K. Recent advances on the association of apoptosis in chronic non healing diabetic wound. World J Diabetes. 2014;5:756–62.
Kataoka Y, Yamada Y. Van der Waals type equation of state for Lennard-Jones fluid and the fluctuation of the potential energy by molecular dynamics simulations. Mol Simul - MOL SIMULAT. 2011;38:1–6.
Zhou H-X, Pang X. Electrostatic Interactions in Protein Structure, Folding, Binding, and Condensation, Chemical reviews, 118 (2018) 1691–1741.
Kaushik AC, Sahi S. Molecular modeling and molecular dynamics simulation-based structural analysis of GPR3. Netw Model Anal Health Inf Bioinf. 2017;6:9.
This study has been funded and supported by Tehran university of medical sciences (TUMS); Grant no. 54069. This fund was helpful in analysis and interpretation of data.
Ethics approval and consent to participate
This article does not contain any studies with human participants or animals performed by any of the authors.
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
About this article
Cite this article
Beheshtizadeh, N., Salimi, A., Golmohammadi, M. et al. In-silico engineering of RNA nanoplatforms to promote the diabetic wound healing. BMC Chemistry 17, 52 (2023). https://doi.org/10.1186/s13065-023-00969-4
- RNA interferences
- Diabetic foot ulcers
- Wound healing
- Tissue engineering