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.

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.

In silico Prediction of miRNA Interactions With Candidate Atherosclerosis Gene mRNAs

The involvement of genes and miRNAs in the development of atherosclerosis is a challenging problem discussed in recent publications. It is necessary to establish which miRNAs affect the expression of candidate genes. We used known candidate atherosclerosis genes to predict associations. The quantitative characteristics of interactions of miRNAs with mRNA candidate genes were determined using the program, which identifies the localization of miRNA binding sites in mRNA, the free energy interaction of miRNA with mRNA. In mRNAs of GAS6 and NFE2L2 candidate genes, binding sites of 21 miRNAs and of 15 miRNAs, respectively, were identified. In IRS2 mRNA binding sites of 25 miRNAs were located in a cluster of 41 nt. In ADRB3, CD36, FASLG, FLT1, PLA2G7, and PPARGC1A mRNAs, clusters of miR-466, ID00436.3p-miR, and ID01030.3p-miR BS were identified. The organization of overlapping miRNA binding sites in clusters led to their compaction and caused competition among the miRNAs. The binding of 53 miRNAs to the mRNAs of 14 candidate genes with free energy interactions greater than -130 kJ/mole was determined. The miR-619-5p was fully complementary to ADAM17 and CD36 mRNAs, ID01593.5p-miR to ANGPTL4 mRNA, ID01935.5p-miR to NFE2L2, and miR-5096 to IL18 mRNA. Associations of miRNAs and candidate atherosclerosis genes are proposed for the early diagnosis of this disease.

Prediction of miRNA interaction with mRNA of stroke candidate genes

BACKGROUND

The role of miRNA in tissue affected by stroke is actively studied, but it remains unclear which miRNAs and target genes are involved in the development of stroke.

METHODS

The MirTarget program defines the following features of a miRNA binding to a mRNA: the binding start site, the location of the binding site in mRNA, the free energy of a miRNA binding with a mRNA, and the interaction schemes of miRNA and mRNA.

RESULTS

The interaction of 6565 miRNAs with mRNAs of stroke candidate genes was determined. The association of the mRNAs of stroke candidate genes with miRNAs depends on the level of gene expression. Some highly expressed candidate genes are targets of miR-619-5p and miR-5095, which have binding sites located on overlapping mRNA nucleotide sequences (clusters). miR-619-5p and miR-5095 bind to mRNA of 15 genes. Clusters for the binding of miR-1273f,d,e are in mRNAs of highly expressed genes. The start sites of miR-1273d and miR-1273e binding in all clusters are in sequences with one and ten nucleotides, respectively. The clusters of multiple miR-574-5p and ID00470.5p-miR binding sites and the clusters of the miR-466, ID01030.3p-miR, and ID00436.3p-miR binding sites are in mRNAs of some genes expressed at low levels.

CONCLUSION

The organization of miRNA binding sites into clusters reduces the length of mRNA and creates competition between miRNAs for binding to mRNA of a target gene. The characteristics of miRNA associations with target genes can be used to recommend markers for a diagnosis of stroke.

Prediction of clusters of miRNA binding sites in mRNA candidate genes of breast cancer subtypes

The development of breast cancer (BC) subtypes is controlled by distinct sets of candidate genes, and the expression of these genes is regulated by the binding of their mRNAs with miRNAs. Predicting miRNA associations and target genes is thus essential when studying breast cancer. The MirTarget program identifies the initiation of miRNA binding to mRNA, the localization of miRNA binding sites in mRNA regions, and the free energy from the binding of all miRNA nucleotides with mRNA. Candidate gene mRNAs have clusters (miRNA binding sites with overlapping nucleotide sequences). mRNAs of EPOR, MAZ and NISCH candidate genes of the HER2 subtype have clusters, and there are four clusters in mRNAs of MAZ, BRCA2 and CDK6 genes. Candidate genes of the triple-negative subtype are targets for multiple miRNAs. There are 11 sites in CBL mRNA, five sites in MMP2 mRNA, and RAB5A mRNA contains two clusters in each of the three sites. In SFN mRNA, there are two clusters in three sites, and one cluster in 21 sites. Candidate genes of luminal A and B subtypes are targets for miRNAs: there are 21 sites in FOXA1 mRNA and 15 sites in HMGA2 mRNA. There are clusters of five sites in mRNAs of ITGB1 and SOX4 genes. Clusters of eight sites and 10 sites are identified in mRNAs of SMAD3 and TGFB1 genes, respectively. Organizing miRNA binding sites into clusters reduces the proportion of nucleotide binding sites in mRNAs. This overlapping of miRNA binding sites creates a competition among miRNAs for a binding site. From 6,272 miRNAs studied, only 29 miRNAs from miRBase and 88 novel miRNAs had binding sites in clusters of target gene mRNA in breast cancer. We propose using associations of miRNAs and their target genes as markers in breast cancer subtype diagnosis.