In silico Prediction of miRNA Interactions With Candidate Atherosclerosis Gene mRNAs

Circulating microRNAs in Carotid Atherosclerosis: Complex Interplay and Possible Associations with Atherothrombotic Stroke

Atherosclerosis is a chronic inflammatory disorder which remains the main cause of cardiovascular morbidity and mortality, with carotid atherosclerosis (CA) being a major cause of ischemic stroke. Epigenetic regulation plays a significant role in CA progression and stroke, yet the impact of circulating microRNA expression, associated with atherogenesis, has not been clearly defined. We included 81 patients with moderate-severe CA (mean age 67 ± 7 years, 53% male), 42% of whom had prior ipsilateral ischemic stroke (i.e., were symptomatic). A total of 24 miRs were identified and their plasma expression levels were measured. We observed that several microRNAs were up-regulated in stroke survivors, namely miR-200c-3p (30.6 vs. 29.7, p = 0.047), miR-106b-3p (31.01 vs. 30.25, p = 0.004), and miR-494-5p (39 vs. 33, p < 0.001), while others (miR183-3p [25.5 vs. 28.6, p < 0.001], miR-126-5p [35.6 vs. 37.1, p = 0.03], and miR-216-3p [12.34 vs. 16.2, p < 0.001]) had lower plasma levels in symptomatic patients. In a multivariable logistic regression model for symptomatic CA, the only miRs showing statistical significance were miR-106b-5p, miR-183-3p, miR-216-3p, and miR-494-5p. Cluster analysis demonstrated differential miR expression in CA patients depending on their stroke status. Epigenetic modulation, represented as complex interplay between circulating miRs of different atherogenic potential, may play a significant role in CA development and progression. In our study, we show possible candidates for future research regarding CA and stroke.

New Insight into Mechanisms of Cardiovascular Diseases: An Integrative Analysis Approach to Identify TheranoMiRNAs

MiRNAs regulate both physiological and pathological heart functions. Altered expression of miRNAs is associated with cardiovascular diseases (CVDs), making miRNAs attractive therapeutic strategies for the diagnosis and treatment of heart diseases. A recent publication defined, for the first time, the term theranoMiRNA, meaning the miRNAs that may be used both for diagnosis and treatment. The use of in silico tools may be considered fundamental for these purposes, clarifying several molecular aspects, suggesting future directions for in vivo studies. This study aims to explore different bioinformatic tools in order to clarify miRNA interactions with candidate genes, demonstrating the need to use a computational approach when establishing the most probable associations between miRNAs and target genes. This study focused on the functions of miR-133a-3p, miR-21-5p, miR-499a-5p, miR-1-3p, and miR-126-3p, providing an up-to-date overview, and suggests future lines of research in the identification of theranoMiRNAs related to CVDs. Based on the results of the present study, we elucidated the molecular mechanisms that could be linked between miRNAs and CVDs, confirming that these miRNAs play an active role in the genesis and development of heart damage. Given that CVDs are the leading cause of death in the world, the identification of theranoMiRNAs is crucial, hence the need for a definition of in vivo studies in order to obtain further evidence in this challenging field of research.

Endogenous piRNAs Can Interact with the Omicron Variant of the SARS-CoV-2 Genome

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which caused the COVID-19 pandemic, can still infect populations in many countries around the globe. The Omicron strain is the most mutated variant of SARS-CoV-2. The high transmissibility of the strain and its ability to evade immunity necessitate a priority study of its properties in order to quickly create effective means of preventing its spread. The current research aimed to examine the in silico interaction between PIWI-interacting RNAs (piRNAs) and the SARS-CoV-2 genome (gRNA) to identify endogenous piRNAs and propose synthetic piRNAs with strong antiviral activity for drug development. This study used validated bioinformatic approaches regarding the interaction of more than eight million piRNAs with the SARS-CoV-2 genome. The piRNAs’ binding sites (BSs) in the 5′UTR were located with overlapping nucleotide sequences termed clusters of BSs. Several BSs clusters have been found in the nsp3, nsp7, RNA-dependent RNA polymerase, endoRNAse, S surface glycoprotein, ORF7a, and nucleocapsid. Sixteen synthetic piRNAs that interact with gRNA have been proposed with free binding energy ranging from −170 kJ/mol to −175 kJ/mol, which can be used to create drugs that suppress the reproduction of SARS-CoV-2.

piRNA and miRNA Can Suppress the Expression of Multiple Sclerosis Candidate Genes

Multiple sclerosis (MS) is a common inflammatory demyelinating disease with a high mortality rate. MS is caused by many candidate genes whose specific involvement has yet to be established. The aim of our study was to identify endogenous miRNAs and piRNAs involved in the regulation of MS candidate gene expression using bioinformatic methods. A program was used to quantify the interaction of miRNA and piRNA nucleotides with mRNA of the target genes. We used 7310 miRNAs from three databases and 40,000 piRNAs. The mRNAs of the candidate genes revealed miRNA binding sites (BSs), which were located separately or formed clusters of BSs with overlapping nucleotide sequences. The miRNAs from the studied databases were generally bound to mRNAs in different combinations, but miRNAs from only one database were bound to the mRNAs of some genes. For the first time, a direct interaction between the complete sequence of piRNA nucleotides and the nucleotides of their mRNA BSs of target genes was shown. One to several clusters of BSs of miRNA and piRNA were identified in the mRNA of ADAM17, AHI1, CD226, EOMES, EVI5, IL12B, IL2RA, KIF21B, MGAT5, MLANA, SOX8, TNFRSF1A, and ZBTB46 MS candidate genes. These piRNAs form the expression regulation system of the MS candidate genes to coordinate the synthesis of their proteins. Based on these findings, associations of miRNAs, piRNAs, and candidate genes for MS diagnosis are recommended.

piRNAs may regulate expression of candidate genes of esophageal adenocarcinoma

Elucidation of ways to regulate the expression of candidate cancer genes will contribute to the development of methods for cancer diagnosis and therapy. The aim of the present study was to show the role of piRNAs as efficient regulators of mRNA translation of esophageal adenocarcinoma (EAC) candidate genes. We used bioinformatic methods to study the interaction characteristics of up to 200 thousand piRNAs with mRNAs of 38 candidate EAC genes. The piRNAs capable of binding to mRNA of AR, BTG3, CD55, ERBB3, FKBP5, FOXP1, LEP, SEPP1, SMAD4, and TP53 genes with high free energy by the formation of hydrogen bonds between canonical (G-C, A-U) and noncanonical (G-U, A-C) piRNA and mRNA nucleotide pairs were revealed. The organization of piRNA binding sites (BSs) in the mRNA of candidate genes was found to overlap nucleotide sequences to form clusters. Clusters of piRNA BSs were detected in the 5'-untranslated region, coding domain sequence, and 3'-untranslated region of mRNA. Due to the formation of piRNA binding site clusters, compaction of BSs occurs and competition between piRNAs for binding to mRNA of candidate EAC genes occurs. Associations of piRNA and candidate genes were selected for use as markers for the diagnosis of EAC.

In silico and in vitro analysis of microRNAs with therapeutic potential in atherosclerosis

Atherosclerosis is a chronic inflammatory disease in which aberrant lipid metabolism plays a key role. MicroRNAs (miRNAs), micro-coordinators of gene expression, have been recently proposed as novel clinical biomarkers and potential therapeutic tools for a broad spectrum of diseases. This study aimed to identify miRNAs with therapeutic potential in atherosclerosis. Bioinformatic databases, including experimentally validated and computational prediction tools as well as a novel combination method, were used to identify miRNAs that are able to simultaneously inhibit key genes related to the pathogenesis of atherosclerosis. Further validation of genes and miRNAs was conducted using the STRING online tool, KEGG pathway analysis and DIANA-miRPath. The inhibitory effects of the identified miRNAs in HepG2 and Huh7 cells were verified by real-time PCR. The MTT assay was utilized to evaluate cell cytotoxicity effects of miRNAs. Atherosclerotic drug-targeted genes were selected as key genes. Strong interactions between genes were confirmed using STRING. These genes were shown to be integral to critical pathological processes involved in atherosclerosis. A novel combined method of validated and predicted tools for the identification of effective miRNAs was defined as the combination score (C-Score). Bioinformatic analysis showed that hsa-miR-124-3p and hsa-miR-16-5p possessed the best C-Score (0.68 and 0.62, respectively). KEGG and DIANA-miRPath analysis showed that selected genes and identified miRNAs were involved in atherosclerosis-related pathways. Compared with the controls in both HepG2 and Huh7 cell lines, miR-124 significantly reduced the expression of CETP, PCSK9, MTTP, and APOB, and miR-16 significantly reduced the expression of APOCIII, CETP, HMGCR, PCSK9, MTTP, and APOB, respectively. The cytotoxicity assay showed that miR-124 reduced cell viability, especially after 72 h; however, miR-16 did not show any significant cytotoxicity in either cell line. Our findings indicate that hsa-miR-124 and miR-16 have potential for use as therapeutic candidates in the treatment of atherosclerosis.

Profile analysis and functional modeling identify circular RNAs in nonalcoholic fatty liver disease as regulators of hepatic lipid metabolism

Nonalcoholic fatty liver disease (NAFLD) is the leading cause of chronic liver disease, associated with an outcome of hepatic fibrosis/cirrhosis and hepatocellular carcinoma. However, limited exploration of the underlying mechanisms hinders its prevention and treatment. To investigate the mechanisms of epigenetic regulation in NAFLD, the expression profile of circular RNA (circRNA) of rodents in which NAFLD was induced by a high-fat, high-cholesterol (HFHC) diet was studied. Modeling of the circRNA-microRNA (miRNA) -mRNA regulatory network revealed the functional characteristics of NAFLD-specific circRNAs. The targets and effects in the liver of such NAFLD-specific circRNAs were further assessed. Our results uncovered that the downregulation of 28 annotated circRNAs characterizes HFHC diet-induced NAFLD. Among the downregulated circRNAs, long intergenic non-protein coding RNA, P53 induced transcript (LNCPINT) -derived circRNAs (circ_0001452, circ_0001453, and circ_0001454) targeted both miR-466i-3p and miR-669c-3p. Their deficiency in NAFLD abrogated the circRNA-based inhibitory effect on both miRNAs, which further inactivated the AMPK signaling pathway via AMPK-α1 suppression. Inhibition of the AMPK signaling pathway promotes hepatic steatosis, depending on the transcriptional and translational upregulation of lipogenic genes, such as those encoding sterol regulatory element-binding protein 1 (SREBP1) and fatty acid synthase (FASN) in hepatocytes. The levels of LNCPINT-derived circRNAs displayed a negative association with hepatic triglyceride (TG) concentration. These findings suggest that loss of LNCPINT-derived circRNAs may underlie NAFLD via miR-466i-3p- and miR-669c-3p-dependent inactivation of the AMPK signaling pathway.

Bioinformatics Analysis of the Interaction of miRNAs and piRNAs with Human mRNA Genes Having di- and Trinucleotide Repeats

The variability of nucleotide repeats is considered one of the causes of diseases, but their biological function is not understood. In recent years, the interaction of miRNAs and piRNAs with the mRNAs of genes responsible for developing neurodegenerative and oncological diseases and diabetes have been actively studied. We explored candidate genes with nucleotide repeats to predict associations with miRNAs and piRNAs. The parameters of miRNAs and piRNA binding sites with mRNAs of human genes having nucleotide repeats were determined using the MirTarget program. This program defines the start of the initiation of miRNA and piRNA binding to mRNAs, the localization of miRNA and piRNA binding sites in the 5′-untranslated region (5′UTR), coding sequence (CDS) and 3′-untranslated region (3′UTR); the free energy of binding; and the schemes of nucleotide interactions of miRNAs and piRNAs with mRNAs. The characteristics of miRNAs and piRNA binding sites with mRNAs of 73 human genes were determined. The 5′UTR, 3′UTR and CDS of the mRNAs of genes are involved in the development of neurodegenerative, oncological and diabetes diseases with GU, AC dinucleotide and CCG, CAG, GCC, CGG, CGC trinucleotide repeats. The associations of miRNAs, piRNAs and candidate target genes could be recommended for developing methods for diagnosing diseases, including neurodegenerative diseases, oncological diseases and diabetes.

Identification of Bovine miRNAs with the Potential to Affect Human Gene Expression

Milk and other products from large mammals have emerged during human evolution as an important source of nutrition. Recently, it has been recognized that exogenous miRNAs (mRNA inhibited RNA) contained in milk and other tissues of the mammalian body can enter the human body, which in turn have the ability to potentially regulate human metabolism by affecting gene expression. We studied for exogenous miRNAs from Bos taurus that are potentially contain miRNAs from milk and that could act postprandially as regulators of human gene expression. The interaction of 17,508 human genes with 1025 bta-miRNAs, including 245 raw milk miRNAs was studied. The milk bta-miR-151-5p, bta-miR-151-3p, bta-miRNA-320 each have 11 BSs (binding sites), and bta-miRNA-345-5p, bta-miRNA-614, bta-miRNA-1296b and bta-miRNA-149 has 12, 14, 15 and 26 BSs, respectively. The bta-miR-574-5p from cow’s milk had 209 human genes in mRNAs from one to 25 repeating BSs. We found 15 bta-miRNAs that have 100% complementarity to the mRNA of 13 human target genes. Another 12 miRNAs have BSs in the mRNA of 19 human genes with 98% complementarity. The bta-miR-11975, bta-miR-11976, and bta-miR-2885 BSs are located with the overlap of nucleotide sequences in the mRNA of human genes. Nucleotide sequences of BSs of these miRNAs in 5′UTR mRNA of human genes consisted of GCC repeats with a total length of 18 nucleotides (nt) in 18 genes, 21 nt in 11 genes, 24 nt in 14 genes, and 27–48 nt in nine genes. Nucleotide sequences of BSs of bta-miR-11975, bta-miR-11976, and bta-miR-2885 in CDS mRNA of human genes consisted of GCC repeats with a total length of 18 nt in 33 genes, 21 nt in 13 genes, 24 nt in nine genes, and 27–36 nt in 11 genes. These BSs encoded polyA or polyP peptides. In only one case, the polyR (SLC24A3 gene) was encoded. The possibility of regulating the expression of human genes by exogenous bovine miRNAs is discussed.

Evolutionary Changes in the Interaction of miRNA With mRNA of Candidate Genes for Parkinson's Disease

Parkinson’s disease (PD) exhibits the second-highest rate of mortality among neurodegenerative diseases. PD is difficult to diagnose and treat due to its polygenic nature. In recent years, numerous studies have established a correlation between this disease and miRNA expression; however, it remains necessary to determine the quantitative characteristics of the interactions between miRNAs and their target genes. In this study, using novel bioinformatics approaches, the quantitative characteristics of the interactions between miRNAs and the mRNAs of candidate PD genes were established. Of the 6,756 miRNAs studied, more than one hundred efficiently bound to mRNA of 61 candidate PD genes. The miRNA binding sites (BS) were located in the 5′-untranslated region (5′UTR), coding sequence (CDS) and 3′-untranslated region (3′UTR) of the mRNAs. In the mRNAs of many genes, the locations of miRNA BS with overlapping nucleotide sequences (clusters) were identified. Such clusters substantially reduced the proportion of nucleotide sequences of miRNA BS in the 5′UTRs, CDSs, and 3′UTRs. The organization of miRNA BS into clusters leads to competition among miRNAs to bind mRNAs. Differences in the binding characteristics of miRNAs to the mRNAs of genes expressed at different rates were identified. Single miRNA BS, polysites for the binding for one miRNA, and multiple BS for two or more miRNAs in one mRNA were identified. Evolutionary changes in the BS of miRNAs and their clusters in 5′UTRs, CDSs and 3′UTRs of mRNA of orthologous candidate PD genes were established. Based on the quantitative characteristics of the interactions between miRNAs and mRNAs candidate PD genes, several associations recommended as markers for the diagnosis of PD.

Predicting Associations of miRNAs and Candidate Gastric Cancer Genes for Nanomedicine

Nanoscale miRNAs regulate the synthesis of most human proteins involved in differentiation, proliferation, cell cycle, apoptosis, and other processes associated with the growth and the development of an organism. miRNAs also play a number of important roles in the development of gastric cancer. In this work, we studied the quantitative characteristics of miRNA interactions with 69 candidate gastric cancer genes using bioinformatics approaches. To this end, the MirTarget program was used, which determines the characteristics of miRNA binding to mRNA in the 5′UTR, CDS, and 3′UTR. Associations of miRNAs with alternative target genes and associations of genes with alternative miRNAs were established. The cluster organization of miRNA binding sites (BSs) in mRNA was revealed, leading to the emergence of miRNA competition for binding to the mRNA of a target gene. Groups of target genes with clusters of overlapping BSs include miR-5095, miR-619-5p, miR-1273 family, miR-466, ID01030.3p-miR, ID00436.3p-miR, miR-574-5p, and ID00470.5p-miR. In the defined associations of target genes and miRNAs, miRNA BSs are organized into clusters of multiple BSs, which facilitate the design and the development of a system of chips that can be used to control the state of miRNA and target genes associations in gastric cancer.