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.


INTRODUCTION
The diagnosis, prevention and treatment of atherosclerosis are important tasks of modern medicine (Byrne et al., 2014;Churov et al., 2019;Solly et al., 2019;Shoeibi, 2020). In most cases, polygenic diseases, including atherosclerosis, develop when the expression of different combinations of candidate genes changes. As a result of a review of various publications, it seems that the approaches used to search for markers for the diagnosis and treatment of atherosclerosis, including miRNA, have not yet solved the problem of diagnosing and treating atherosclerosis Ryu et al., 2020;Sun et al., 2020;You et al., 2020). In human, it is known about 7000 miRNA and more than 20,000 genes, and it is unknown how many miRNA and genes from them participate in the development of atherosclerosis (Byrne et al., 2014;Toba et al., 2014;Lu et al., 2018;Liu et al., 2020;Shi et al., 2020). As a rule, in publications several miRNA and some candidate genes are studied, so the combination of such attempts is large (Wang C. et al., 2020;Wang M. et al., 2020;Zhang et al., 2020). Considering the limitations of existing prediction methods, we can say that the establishment of adequate miRNA markers for diseases is unrealistic using existing approaches in the coming years. This statement is confirmed by the fact that over 20 years of studying miRNA participation in various diseases, diagnostic methods and therapeutic ways of treating diseases with miRNA are not used. Most researches study only the correlation between changes in miRNA concentration and expression of candidate genes, but this approach does not establish specific associations of miRNA and target gene. A number of publications argue that miRNA is the cause of disease without understanding that changes in miRNA concentration can occur through independent expression from intergenic regions or from host gene introns (Madrigal-Matute et al., 2013;Andreou et al., 2015;Feinberg and Moore, 2016;Solly et al., 2019). In all cases, miRNAs appear to modify the expression of target genes and not to cause the disease themselves. It is not taken into account that one miRNA can affect several or even 100s of genes, and one gene may be a target of several miRNAs (Atambayeva et al., 2017;Kondybayeva et al., 2018). In this article, we discuss what should be taken into consideration when reviewing the problem of miRNA interaction with candidate genes and show the need to use a systematic approach in establishing the most probable associations of miRNA and target genes.
It has been found that miRNAs, which are nanoscale regulatory biomolecules, are involved in many biological processes at all stages of the development of atherosclerosis, from early endothelial dysfunction to the erosion and rupture of an unstable atherosclerotic plaque (Cipollone et al., 2011;Madrigal-Matute et al., 2013;Andreou et al., 2015;Maitrias et al., 2015;Feinberg and Moore, 2016;Santovito et al., 2016;Chen et al., 2018;Ren et al., 2018;Wang et al., 2018). In addition, miRNAs are considered novel non-invasive biomarkers of the instability of atherosclerotic plaques and have been associated ischemic disorders (Churov et al., 2019). Their detection in the blood of patients may be a promising direction for the diagnosis of atherosclerosis complications such as ischemic stroke and myocardial infarction (Kanuri et al., 2018;Desjarlais et al., 2019;Moghaddam et al., 2019;Vargas-Alarcon et al., 2019;Velle-Forbord et al., 2019;Wiese et al., 2019;Kondybayeva et al., 2020). When the pathogenic role of a specific miRNA is confirmed, it can be considered a potential therapeutic target (Friedman et al., 2009;Fang et al., 2013;Moghaddam et al., 2019;Kazemi et al., 2020). According to the miRBase database, miRNAs have been found in many human tissues and are able to regulate the expression of more than 60% of all protein-coding genes (Kanuri et al., 2018;Ren et al., 2018;Wiese et al., 2019). Full complementarity between miRNA and mRNA results in degradation (Leidinger et al., 2012;Alagia and Eritja, 2016). However, incomplete complementarity is most often observed, in which case miRNAs suppress translation, generally by binding to the 3 -untranslated regions (3 UTRs) of mRNAs (Ivashchenko et al., 2014a,c;Alagia and Eritja, 2016;Atambayeva et al., 2017).
In addition, miRNAs can bind to other regions of target mRNAs, including the 5 -untranslated regions (5 UTR) and coding sequence (CDS) (Orom et al., 2008;Forman and Coller, 2010;Ivashchenko et al., 2013;Zhou and Rigoutsos, 2014;Niyazova et al., 2015;Kondybayeva et al., 2018;Yurikova et al., 2019). We used a program that effectively determines the quantitative characteristics of the interaction of miRNA with mRNA and allows us to identify fundamentally new properties of the binding of miRNA to mRNA. The aim of this work was to identify associations between miRNAs and mRNA candidate genes of atherosclerosis for use as markers for the diagnosis of this disease.

MATERIALS AND METHODS
The nucleotide sequences of 2565 miRNAs (we name this set as 'old miRNAs') were downloaded from the miRBase database 1 (Release 22.1) and 3707 miRNAs (we name this set as 'new miRNAs') obtained from a report by Londin et al. (2015). Due to this work the number of known miRNAs had more than doubled. The nucleotide sequences of genes were obtained from GenBank 2 . A database of 68 candidate genes including the names of the genes and publication sources was compiled, confirming the associations of these genes with atherosclerosis (Supplementary Table S1). A search for the target genes of miRNAs was performed using the MirTarget program (Ivashchenko et al., 2014b). This program determines the following binding characteristics: the start of the miRNA binding site (BS) of mRNA; the locations of miRNA BS (3 UTR, 5 UTR, CDS); the interaction free energy ( G, kJ/mole); and nucleotide interaction schemes between miRNAs and mRNAs. The ratio of G/ Gm (%) was determined for each BS, where Gm is equal to the free energy binding of miRNA with its full complementary nucleotide sequence. The MirTarget program found hydrogen bonds between adenine (A) and uracil (U), guanine (G) and cytosine (C), G and U, A and C. The distances between A and C were equal 1.04 nanometers, between G and C, and between A and U were equal 1.03 nanometers, between G and U equal to 1.02 nanometers (Leontis et al., 2002). The numbers of hydrogen bonds in the G-C, A-U, G-U, and A-C interactions were found to be 3, 2, 1 and 1, respectively (Kool, 2001;Lemieux and Major, 2002). The MirTarget program determines single miRNA BS in mRNA and miRNA BS which are in clusters (BS arranged with overlapping of nucleotide sequences of the same or several miRNAs) . Predicted by the MirTarget program binding sites in over 30 genes were confirmed experimentally (Yurikova et al., 2019).

RESULTS
The BSs of the miRNAs and mRNAs of the target genes were not uniform along the length of the mRNAs. Both multiple and single BS were identified. BSs could be sequential or overlap with each other. Having overlapping nucleotide sequences in a cluster lead to the compaction of the mRNA sequence that is the target of several miRNAs. When the G and G/ Gm values of miRNA BS are close to each other, it can be assumed that in the presence of equal miRNA concentrations, the miRNAs with a larger number of BS will be more likely to bind to the mRNAs of target genes. When the miRNA-mRNA interaction strength and the degree of their complementarity are similar, the miRNA has the highest concentration upon binding.
Characteristics of miRNA Interactions in the 5 UTRs of the mRNAs of Atherosclerosis Candidate Genes BSs of miRNAs were identified in the 5 UTRs of 14 mRNAs. The GAS6 mRNA contains a large cluster with BS of 21 different miRNAs, half of which have two or three BS, while the rest represent single BS ( Table 1). The cluster size is 39 nt, starting at the 11 nt position and ending at position 49 nt. The total BS length in the cluster is 824 nt, where the degree of compaction is 21. This compaction allowed this number of miRNA BS to be located in a 5 UTR with a length of 153 nt. In an association analysis of seven miRNAs with the GAS6 mRNA the free energy interaction of the miRNAs with the mRNA is more than −130 kJ/mole, which indicates that these associations are promising markers of atherosclerosis.
To prove that the identified clusters are reliable, on the example of GAS6 gene showed that the cluster exists in the mRNA of orthological genes and is conservative in its nucleotide composition (Figure 1). There are mainly new miRNAs in the clusters of BSs in 5 UTR region of mRNA of GAS6, NFE2L2 and SCAP genes. In mRNA of GAS6 gene the cluster includes 20 new miRNA and one old miRNA, in mRNA of NFE2L2 gene the cluster contains 15 BSs of new miRNA; in mRNA of SCAP gene there is a cluster for binding of 13 new and one old miRNA ( Table 1). In the other 11 target genes, there are 28 new and one old miRNAs are linked in 5 UTR region of mRNA (Supplementary Table S2). These data suggest that the use of new miRNAs significantly increased the number of effective gene expression modulators.
Some miRNAs have multiple BSs: ID00061.3p-miR, ID01155.3p-miR, ID01702.3p-miR, ID01804.3p-miR, ID02187.5p-miR, ID00296.3p-miR, ID01155.3p-miR by three, ID01641.3p-miR by four, ID00061.3p-miR, ID00296.3p-miR, ID01702.3p-miR by five BSs. These miRNAs are recommended as associations for the diagnosis of atherosclerosis considering the high free energy interaction with mRNA of candidate genes. ID01935.5p-miR has a complete complementarity interaction with mRNA of NFE2L2 gene, which suggests that ID01935.5p-miR and NFE2L2 association should be recommended as a marker. These miRNA correlations clearly demonstrate the predominant influence of new miRNAs on gene expression with BS in 5 UTR. Consequently, researchers actually do not receive information about the regulation of gene expression by new miRNAs, the effect of which is high when only the old miRNAs are studied.
The following circumstances need to be considered when evaluating the effects of many miRNAs in a cluster of BSs in  mRNA. Primarily, miRNA will bind with higher free energy interaction, so quantitative interaction characteristics are needed. The ratio of miRNA concentrations between alternative miRNAs has a significant importance (miRNAs, which bind in the cluster) and especially in relation to the concentration of mRNA. Obviously, the total concentration of alternative miRNAs should be lower than the mRNA concentration, otherwise the protein will not be synthesized. A decrease of concentration of any miRNA below mRNA concentration will not affect protein synthesis. At the same time, an increase of one miRNA above mRNA concentration can completely suppress gene expression. The mRNA of many genes contains clusters of BSs of several similar miRNAs. In this way, nature has optimized the expression of several genes under common miRNA control. It is necessary to consider that many genes are expressed in cells of different tissues to different degrees, which may affect the dependence of protein synthesis from miRNA. However, a large number of binding miRNA may not allow the target gene expression to increase significantly. Many miRNAs have been found in the blood and serum in combination with AGO proteins and within exosomes, which indicates their capability to circulate freely throughout the body and to reach many organs and tissues. It is necessary to remember that the synthesis of intronic miRNAs depends on the expression of host gene, while the synthesis of other miRNAs is made from transcripts of intergenic regions. The reasons for changes in the expression of such miRNA in case of disease should be known. There are known changes of miRNA expression by several orders. Considering the above information, the reports are perceived cautiously that some researchers have been able to detect marker miRNA in specific diseases.
The above examples of presence of BS clusters in mRNA of some genes for several miRNAs suggest that such genes are expressed under the common control of the miRNA group and, consequently, these genes form a network of interconnected genes controlling key metabolic processes. To confirm the reliability of found miRNA and their BSs, it is possible to use data about the presence of such BSs or their clusters in orthological genes (Figure 1).
With the given examples we show the inadequacy of miRNA BSs establishment by programs based only on a miRNA seed sequence. The whole miRNA nucleotide sequence is important, which is confirmed by the conservation of miRNA nucleotide sequences and corresponding BS during millions of years. The nucleotide sequences of miRNA and their BSs have been conserved in the mRNA genes of animal and plant organisms over 10s of millions of years of evolution (Bari et al., 2013;Ivashchenko et al., 2013Ivashchenko et al., , 2014bYurikova et al., 2019).
ID00061.3p-miR, ID00296.3p-miR, ID00457.3p-miR, ID00522.5p-miR, ID01041.5p-miR, ID01155.3p-miR, ID01641.3p-miR, ID01702.3p-miR, ID01804.3p-miR, ID01873.3p-miR, ID02187.5p-miR and ID03367.5p-miR interact not only with GAS6 mRNA but also with NFE2L2 mRNA, where they form a large cluster together with other miRNAs. The cluster of miRNA BS in the NFE2L2 mRNA is one of the largest clusters of miRNAs associated with atherosclerosis-related genes in which BS are formed in the 5 UTR. The cluster consists of 14 different miRNAs, ranging from the 438 to 482 nt positions, with a size of 45 nt. The total size of all binding sites is 669 nt, resulting in a degree of compaction of 15. Due to compaction, this length of BS is a small fraction of the 5 UTR length of 555 nt. A complete complementarity interaction of ID01935.5p-miR with NFE2L2 mRNA was identified. Eight miRNAs interact with NFE2L2 mRNA with a free energy of more than −130 kJ/mole. These associations are recommended as valuable markers of atherosclerosis. The above examples of associations of miRNAs and target genes show that determining the expression level of one or more miRNA without monitoring target gene expression will not provide adequate data on the specific association of these miRNAs with the disease.
The ADAM10 mRNA contains BS for ID02761.3p-miR with overlapping nucleotide sequences (Supplementary Table S2). The ADCY9 mRNA is characterized by the presence of two clusters, each of which is formed by two single BS. In addition to this gene, clusters were identified in the following mRNAs: IRS2 and PIN1, which contain BS for two miRNAs. The PNPLA3 and mRNAs have a cluster consisting of three single BS. The CXCL12 mRNA contains a cluster of four different miRNAs, of which ID02036.3p-miR, ID01293.5p-miR, ID00417.3p-miR and ID02066.5p-miR have two BS. The cluster starts at 65 nt and ends at 93 nt, with a length of 29 nt and an average G/ Gm of 93%. The sum of the BS lengths of these miRNAs was 170 nt and was six times longer than the cluster length. Due to this compaction, all BS were located in a 5 UTR with a length of 92 nt. The MMP2 mRNAs contained a cluster of six miRNA BS with a total length of 136 nt and a 5 UTR length of 311 nt. The PNPLA3 mRNA contains a cluster of three miRNA BS. Among the clusters in the 5 UTRs of the target genes, one can distinguish the SCAP gene cluster formed by BS from 14 miRNAs from the 99 nt to the 136 nt position, in which both single BS and multiple BS are found, more specifically, seven multiple sites and seven single sites ( Table 1). Of the 14 miRNAs with the exception of ID00756.3p-miR, ID01403.5p-miR, ID01652.3p-miR and ID03151.3p-miR all were in the group of miRNAs binding to the GAS6 mRNA. Furthermore, ID00061.3p-miR, ID00296.3p-miR, and ID01702.3p-miR had five BS each, which significantly increased the likelihood of their binding compared to that of competing miRNAs. The total length of miRNA BS was 735 nt, and the degree of compaction was 19. Without compaction, these sites could not be located in this 255 nt long 5 UTR. This cluster is characterized by a maximum G of −144 kJ/mole for the BS of ID00296.3p-miR and ID01702.3p-miR. Of all the interactions of miRNAs and the 5 UTRs of target mRNAs, five had a free energy of more than −130 kJ/mole. Convincing evidence of the effectiveness of the MirTarget program in determining the characteristics of miRNA binding to the 5 UTR, CDS and 3 UTR of mRNA is given in Figure 2. The diagrams show the interaction of nucleotides along the entire length of the miRNA in the BS of the target mRNA and the characteristics of this interaction.

Characteristics of miRNA Interactions With the CDS of mRNAs of Atherosclerosis Candidate Genes
Most genes contained more one BS for one miRNA in the CDS of mRNAs. There were some genes with several BS or multiple BS where clusters could form. It is worth noting that IRS2 is characterized by the presence of three clusters and a large number of BS compared to other genes. This suggests that this gene is more susceptible to regulation by miRNA ( Table 2). The first cluster was composed of 25 different miRNAs, 13 of which had two or more BS (ID00061.3p-miR, ID00457.3p-miR, ID00756.3p-miR, ID01155.3p-miR, ID01702.3p-miR, ID01804.3p-miR, ID01873.3p-miR, ID01879.5p-miR ID02064.5p-miR, ID02187.5p-miR, ID03229.5p-miR, ID03367.5p-miR, miR-3960), and the rest were composed of single BS. The cluster was 41 nt long, starting at 2586 nt and ending at 2626 nt. The total length of the BS in the cluster was 1114 nt, while the degree of cluster compaction was 28. The 2-nd cluster was formed by single BS of seven miRNAs. The 3-day cluster consisted of the BS of nine different miRNAs. ID00296.3p-miR and ID01702.3p-miR had three binding sites. In general, the cluster was 41 nt long, extending from position 4304 to 4344 nt. The total length of all the BS of the three clusters was 1607 nt, which is 40% of the total CDS length of 4017 nt. Due to compaction, the total length of all clusters is 110 nt, that is, only 2.7%. Thirteen miRNAs interacted with the IRS2 mRNA with a free energy of more than −130 kJ/mole, which gives reason to recommend these interactions as diagnostic markers of atherosclerosis. Each of the three clusters in the MR of the orthologous IRS2 genes encoded different highly homologous oligopeptides (Supplementary Table S3).
Clusters consisting of single miRNA BS were found in the following genes (Supplementary Table S4): two BS in ACE with an average G of −128 kJ/mole; five BS in ADRB3 for five miRNAs with an average G of −115 kJ/mole and BS in FASLG for ID00061.3p-miR, ID00296.3p-miR, ID01641.3p-miR, and ID01702.3p-miR. SIRT1 contained two sites for ID03332.3p-miR and single BS for ID00278.3p-miR and ID00811.3p-miR, with an average G of −129 kJ/mole.
TBC1D10B gene is involved in the functioning of cardiovascular endothelial cells (Li et al., 2017) and it is a target of miR-762 with fully complementary binding (Supplementary Table S4). This site encodes the APAPAPAPAPAPA oligopeptide, with the flanking oligopeptides AWVPGSAQTS and VTGSTVVVLTL (Supplementary Table S5).
MRNA of CDKN1C gene (Rodriguez et al., 2007) contains three clusters consisting of both single and multiple BS. The first cluster consists of two BS for miR-762 and ID03129.3p-miR and has a size of 36 nt, extending from 738 to 753 nt. The second cluster is composed of 17 miR-762 BS, 3 ID00099.3p-miR BS, and 1 ID02682.5p-miR BS. ID00036.3p-miR, ID01075.3p-miR and ID00411.5p-miR cover a region of 31 nt, extending from the 888 nt position to 918 nt with an average G value of −125 kJ/mole. The mRNA of the candidate gene CDKN1C has 17 sequentially located binding sites for miR-762 that encode the oligopeptide (AP)18, which indicates a strong   Table S6). These animals can be used to study the role of miR-762 in the development of atherosclerosis. A similar type of organization of miRNA binding multiple sites has been found in the mRNAs of many genes involved in the development of cardiovascular and other diseases (Kondybayeva et al., 2018;Aisina et al., 2019;Yurikova et al., 2019).

dependence of gene expression on this miRNA (Supplementary
KLF2 is characterized by the presence of two clusters consisting of the BS for 14 and 2 different miRNAs. The first of these clusters is formed by single BS starting from the 264 to 296 nt position with a size of 33 nt ( Table 2). The total length of all BS is 305 nt. The degree of compaction of the binding sites in the cluster is nine. The second cluster of this gene is formed only by two single BS. In the CDS of KLF2 mRNA BS were detected for ID00061.3p-miR, ID00457.3p-miR, ID01155.3p-miR, ID01702.3p-miR, ID01804.3p-miR, ID02187.5p-miR, ID03367.5p-miR which have BS in the 5 UTR of GAS6 mRNA. This indicates an overall control of these genes expression by these miRNAs. The maximum G/ Gm value (100%) among these genes was detected at the ID01593.5p-miR BS, which interact with the ANGPTL4 mRNA.
PDE4D formed two clusters. The first cluster contained the BS for ID00061.3p-miR, ID01155.3p-miR, ID03064.3p-miR, miR-3960 and ID01702.3p-miR, which had three BS with a free energy value of −134 kJ/mole. ID01641.3p-miR, miR-3960 and ID01702.3p-miR BS were found in the second cluster. The size of the first cluster was 34 nt, extending from position 336 to 369 nt. The second cluster consisted of five multiple BS and eight single BS; the cluster size was 49 nt, extending from 391 to 439 nt. The total length of miRNA BS was 491 nt, and the degree of compaction was 10 (Supplementary Table S4).
There were 39 BS for new miRNAs and two BS for old miRNA in IRS2 mRNA. Two clusters of binding sites in the KLF2 mRNA interacted with 14 new miRNAs and one old miRNA. In addition, there were 15 BS for new miRNA and two BS for miR-3960 in the PDE4D mRNA. A total in CDS mRNA of 40 target genes BS were identified for 148 new miRNAs and 32 old miRNAs.

Characteristics of miRNA Interactions in the 3 UTRs of mRNAs of Atherosclerosis Candidate Genes
The 40 genes involved in the development of atherosclerosis were characterized by BS in the 3 UTRs. A feature of miRNA interactions with atherosclerosis candidate genes is the presence of BS for miR-466, ID00436.3p-miR, and ID01030.3p-miR in .3p-miR. Another feature of candidate genes is the presence of clusters of binding sites for miR-5095 and miR-619-5p in their mRNA. There are two clusters of BS of miR-5095 and miR-619-5p located 12 nucleotides apart in IL18 mRNA. In these clusters, the start sites of miR-5095 and miR-619-5p BS are located six nucleotides apart, which is probably due to their common origin. In the mRNAs of CD36 and ADAM17 completely complementary BS with miR-619-5p and BS between the mRNA of IL18 with miR-5096 were revealed. The BS of miR-1273f and miR-1273e form a cluster in the mRNAs of IGF1 gene (Supplementary Table S7).
The CD36 mRNA cluster, consisting of BS to two different miRNAs, miR-619-5p and miR-5585-3p, extended from 4168 to 4197 nt. The CD36 mRNA cluster consisted of both single BS and multiple BS. In general, the cluster consisted of five ID01030.3p-miR, five miR-466, six ID00436.3p-miR and one ID01727.5p-miR sequences. The size was 38 nt, extending from 3529 to 3566 nt, with an average G value of −106 kJ/mole. Associations of miRNA and target genes having BS in 3 UTR mRNA are very different from those in 5 UTR and CDS. In the 3 UTR mRNA there are no clusters of BS over four miRNAs. The number of old and new miRNAs is comparable. The total number of binding sites of new miRNA is 57 and 63 old miRNA per mRNA of 40 target genes. This is probably due to the fact that the GC content of old miRNAs is comparable to that of 3 UTRs, while in 5 UTR and CDS it is about 10% higher. In the new miRNA GC content is comparable to that in 5 UTR and CDS, so the number of BS for new miRNA in these regions of mRNA is much larger.
The choice of miRNA-gene associations is difficult to offer as miRNAs such as miR-619-5p, miR-5095-, miR-5096-and miR-5585-3p, which have targets in many genes and most likely play the role of stabilizers of protein expression. Similar to miR-466, ID00436.3p-miR and ID01030.3p-miR, which have multiple BS in the mRNA of many genes, which also gives them the role of stabilizing the expression of their targets. A significant increase in the concentration of any these miRNAs will lead to numerous metabolic disorders and to different diseases consequently. If about 100-200 key miRNAs are identified with miRNA-chips, the range of candidate genes will be significantly reduced because with bioinformatic approaches it could be done easily and quickly.

DISCUSSION
Since it would be advisable to have lower costs for the diagnosis of the disease, it is desirable to minimize the number of associations of miRNA and target genes. However, in the case of polygenic diseases, such methods could be difficult to develop because it is not known which gene and which miRNA cause a particular patient's disease. Therefore, now it is necessary to select a list of priority genes for diagnosis and then check which miRNA could influence protein expression. With the MirTarget program, this is easily to install and, if alternative miRNAs are involved, it is also simple to identify miRNAs competitiveness. The results of the present research demonstrate the effectiveness of this approach.
It was found that the BS of some miRNAs formed clusters in the 5 UTR and CDS sequences of mRNA, resulting in competition between these miRNAs for binding to mRNA and, accordingly, for suppressing the expression of the target gene. The quantitative characteristics of miRNA BS with mRNA make it possible to predict which miRNAs can more efficiently bind to mRNA at equal miRNA concentrations. At different concentrations of these miRNAs, kinetic equations can be used to predict their suppressive effect on the expression of the target gene. The revealed associations of several miRNAs that bind to mRNAs of different candidate genes make it possible to predict the effect of these miRNAs on the corresponding genes, which can be expressed to different degrees. For example, GAS6, NFE2L2, and IRS2 are targets of several identical miRNAs and bind with these miRNAs depending on their concentrations in the cell. In cells of different tissues, the ratio of miRNA and target gene concentrations will be different, and the effect of the miRNAs on these target genes will also be different. This brief discussion of the different interactions of miRNAs and candidate genes under different circumstances shows the complexity of these networks and begins to paint a picture of the interactions of miRNA groups with groups of target genes. Without this assumption, it is difficult to predict or evaluate the involvement of a single miRNA or a single gene in the development of a disease. From the identified interactions, it is necessary to choose those that control more genes and miRNAs and reflect the development of atherosclerosis. A number of mRNAs and miRNAs strongly interacting with each other are proposed in this study. For example, PDE4D and ID00296.3p-miR; SCAP and ID00296.3p-miR and ID01702.3p-miR; ADCY9 and ID00296.3p-miR; IRS2 and ID01702.3p-miR, ID03064.3p-miR, ID01804.3p-miR, and ID00296.3p-miR; NFE2L2 and ID01935.5p-miR, ID00296.3p-miR, ID01702.3p-miR, and ID01804.3p-miR. In addition, miRNAs that bind with mRNAs of candidate genes with complete complementarity are of interest. Within addition to this condition, it is necessary to consider miRNAs interacting with a large number of genes, since such miRNAs can simultaneously control the expression of many genes. Genes that are regulated by a large number of miRNAs can also be considered as effective markers. Unfortunately, our knowledge is still insufficient to determine which candidate genes play a key role in the development of the disease, and it is not known which miRNAs are the most likely signals for the development of atherosclerosis. However, one thing is obvious -it is necessary to simultaneously study miRNA and gene interactions in order to improve the likelihood of obtaining a valid result. From the established interactions, it is first necessary to consider the interactions between candidate genes with multiple miRNAs that form clusters of BS and miRNA interactions with high free energy. For example, IRS2 is the target of 41 miRNAs with BS located in the CDS, of which ID00061.3p-miR, ID00296.3p-miR, ID00457.3p-miR, ID01041.5p-miR, ID01155.3p-miR, ID01641.3p-miR, ID01702.3p-miR, ID01778.3p-miR, ID01804.3p-miR, ID01895.5p-miR, ID02052.5p-miR, ID02064.5p-miR, ID02950.3p-miR, and ID03064.3p-miR have free energy interactions above −130 kJ/mole. The mRNA of GAS6 interacts with 21 miRNAs, of which those with ID00296.3p-miR, ID01041.5p-miR, ID01106.5p-miR, ID01641.3p-miR, ID01702.3p-miR, ID01804.3p-miR, ID02084.3p-miR, and ID02294.5p-miR are more efficient. Associations of this type were observed for NFE2L2, ADCY9, PNPLA3, SCAP, KLF2 and PDE4D, which bind to many miRNAs with high free energy. Another type of association is the interaction of one or more miRNA with many genes. For example, miR-466, ID00436.3p-miR and ID01030.3p-miR bind with the mRNAs of ADRB3, CD36, FASLG, FLT1, PLA2G7, and PPARGC1A. In addition, miR-619-5p targets several candidate atherosclerosis genes: CD36, IL18, ADAM17, ADAM33, APOL1, BRCA1, F11R, ITGA2, LDLR, MTHFR, PNPLA3, SOAT1, and TNFSF10. The mRNA of ANGPTL4 is fully complementary with ID01593.5p-miR; the mRNAs of ADAM17 and CD36 with miR-619-5p; the mRNA of IL18 with miR-5096; and the mRNA of NFE2L2 with ID01935.5p-miR.
Here, miRNA interactions with candidate atherosclerosis genes have been established, which consist of one gene and several miRNAs, one gene and one miRNA, one miRNA and several genes, and two or more miRNAs with two or more candidate genes. The revealed cluster organization of miRNA BS in mRNA candidate genes contributes to a more accurate diagnosis of the participation of competing miRNAs in the development of atherosclerosis.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/ Supplementary Material.

AUTHOR CONTRIBUTIONS
AI and DM conceived of the study, drafted the manuscript, and gave final approval of the version to be published. SL, DA, AP, and AR conceived of the study and drafted the manuscript. All authors made substantial contributions to acquisition of data, to interpretation and modification of the data, were involved in subsequent rounds of revisions, and read and approved the final manuscript.

FUNDING
The work was carried out with the financial support of the Ministry of Education and Science of the Republic of Kazakhstan within the framework of the grant no. AP05132460. We thank the European Union Horizon 2020 research and innovation program (grant agreement no. 645648) 'Muscle Stress Relief.'