Next Article in Journal
Effect of High Glucose-Induced Oxidative Stress on Paraoxonase 2 Expression and Activity in Caco-2 Cells
Next Article in Special Issue
Intragenic MicroRNAs Autoregulate Their Host Genes in Both Direct and Indirect Ways—A Cross-Species Analysis
Previous Article in Journal
Rare Human Diseases: Model Organisms in Deciphering the Molecular Basis of Primary Ciliary Dyskinesia
Previous Article in Special Issue
Identification of Transcriptional Markers and microRNA–mRNA Regulatory Networks in Colon Cancer by Integrative Analysis of mRNA and microRNA Expression Profiles in Colon Tumor Stroma
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Relationship Between the miRNA Sequence and Disease May be Revealed by Focusing on Hydrogen Bonding Sites in RNA–RNA Interactions

by
Tatsunori Osone
1,* and
Naohiro Yoshida
1,2
1
School of Materials and Chemical Technology, Tokyo Institute of Technology, Yokohama 226-8502, Japan
2
Earth-Life Science Institute, Tokyo Institute of Technology, Tokyo 152-8550, Japan
*
Author to whom correspondence should be addressed.
Cells 2019, 8(12), 1615; https://doi.org/10.3390/cells8121615
Submission received: 30 October 2019 / Revised: 26 November 2019 / Accepted: 9 December 2019 / Published: 11 December 2019
(This article belongs to the Special Issue microRNA Bioinformatics)

Abstract

:
MicroRNAs are important genes in biological processes. Although the function of microRNAs has been elucidated, the relationship between the sequence and the disease is not sufficiently clear. It is important to clarify the relationship between the sequence and the disease because it is possible to clarify the meaning of the microRNA genetic code consisting of four nucleobases. Since seed theory is based on sequences, its development can be expected to reveal the meaning of microRNA sequences. However, this method has many false positives and false negatives. On the other hand, disease-related microRNA searches using network analysis are not based on sequences, so it is difficult to clarify the relationship between sequences and diseases. Therefore, RNA–RNA interactions which are caused by hydrogen bonding were focused on. As a result, it was clarified that sequences and diseases were highly correlated by calculating the electric field in microRNA which is considered as the torus. It was also suggested that four diseases with different major classifications can be distinguished. Conventionally, RNA was interpreted as a one-dimensional array of four nucleobases, but a new approach to RNA from this study can be expected to provide a new perspective on RNA-RNA interactions.

1. Introduction

According to the World Health Organization, the leading causes of death worldwide are stroke (ST), chronic obstructive pulmonary disease (COPD), Alzheimer’s disease (AD), and tuberculosis (TB) [1]. ST has been shown to have a higher incidence in developing countries than in developed countries. Of the various types of ST, 85% are ischemic. It has also been suggested that both genetic and environmental factors are involved in pathophysiology [2]. COPD is caused by an inflammatory response in the lungs to harmful gas or particles. COPD is progressive and often coexists with other diseases. The most common risk factor is smoking [3]. AD, the most common cause of dementia, is a progressive neurodegenerative disease that develops due to neuronal death and loss or atrophy of synaptic functions [4]. TB is an infectious disease caused by Mycobacterium tuberculosis that causes death beyond AIDS due to HIV. Currently, it is estimated that one-third of the world’s population is potentially infected. However, it affects about 5–10% of infected people [5].
MicroRNA (miRNA) is attracting attention as a treatment method and biomarker for these diseases [6]. miRNA is a class of small non-coding RNA, containing around 22 nucleotides. miRNA mainly controls the turnover and translation of target mRNA by binding the 3’ untranslated region of mRNA as a target site [7]. Currently, seed theory is the mainstream in binding miRNA targets [8]. Seed theory states that the 2nd to 8th bases of miRNA bind complementarily to the target site. Various target prediction programs have been developed based on this theory, but it is known that they often have high false positive rates or false negative rates [9]. miRNAs are related to various biological processes because of their functions, and their relation to diseases is well studied [10,11]. In order to clarify the control mechanism of biological process by miRNA, a new theory for the miRNA mechanism may be necessary. For this reason, several studies have recently been conducted to clarify the relationship between diseases and miRNAs using methods that are not based on seed theory. Zou et al. combined the miRNA–miRNA network, miRNA–disease network, and disease–disease network, and performed network analysis by machine learning [12]. Luo et al., Zou et al., Lan et al., and Fu et al. applied network analysis by machine learning to miRNA–disease networks [9,13,14,15]. Ding et al. focused on miRNA clusters/family and secondary structure [16], and Zheng et al. focused on miRNA’s functional similarity and semantic similarity to diseases [17]. These methods can be expected to enable prediction of disease-related miRNA candidates that have not yet been found. For example, hsa –let –7d, hsa –mir –18a, hsa –mir –145, hsa –mir –106b, hsa –let –7e, hsa –let –7b, hsa –mir –19a, and hsa –mir –125a may be related to breast cancer [12].
Seed theory may be superior in terms of the mechanism by which miRNA suppresses translation of mRNA and the prediction of miRNA-mRNA networks. In addition, from the viewpoint of predicting disease-related miRNA candidates that have not yet been found, the above-described approach using machine learning is considered to be excellent. However, neither method may be suitable in terms of decoding the miRNA gene. This is because seed theory focuses on sequence specificity when miRNA is viewed as a translation inhibitor of mRNA, and machine learning approaches are not based on sequences. Decoding miRNA genes not only reveals the mechanism of miRNA-mRNA interactions and miRNA-non-coding RNA interactions, but may also reveal the role of short RNAs in primitive life, such as in its birth and evolution.
In order to clarify the relationship between biological phenomena and miRNA sequences, the concept of black box modeling in system identification in the engineering field was used. In black box modeling, an internal system is regarded as a black box, and system identification is performed using only Input and Output. In this study, Input was designated as miRNA sequence and Output as disease. The RNA-RNA interaction was the focus of attention when quantifying miRNA sequences. The purpose of this study was to verify whether this new approach reveals the relationship between miRNA and ST, COPD, AD, and TB.

2. Materials and Methods

2.1. Datasets

Research papers on miRNA in AD, COPD, ST, and TB using patient tissues (33, 23, 31, and 21 reports, respectively) were selected. (Table S1). When searching, “microrna Disease Name” was used as a keyword. From those papers, the name of miRNAs showing the expression level different from the control and these expression levels were extracted. If the expression level could not be obtained, the value of fold-change (FC) or log2FC was obtained. If no value was obtained in the text or in the table, the value was obtained from the graph using WebPlotDigitizer (https://apps.automeris.io/wpd/). The sequence information of miRNA was obtained from miRbase Release 22.1 (http://www.mirbase.org/).

2.2. Convert Nucleobase to Vector

Adenine (A), guanine (G), cytosine (C), and uracil (U) modelled using Avogadro (©2018 Avogadro Chemistry, v1.2.0) have undergone 25 structural optimizations. After that, it was output as a GAMESS Input file, with the calculated as single point energy, using B3LYP/6–31G (d) as the method, and water as the solvent.
The output file was calculated with Firefly QC package [18], which is partially based on the GAMESS (US) [19] source code, via MoCalc2012. Among the calculation results, MULLIKEN ATOMIC CHARGES of Atoms Used for Hydrogen Bonding (AUHB) of nucleic acids (Figure 1) was used as a vector indicating each nucleobase. When this vector is B and the element at the i−th site is e i , the vector indicating the nucleobase is represented by the following formula (Equation (1)):
B = e 1 ,   e 2 ,   e 3

2.3. Scoring of miRNAs

miRNA scoring was performed by considering miRNA as a ternary vector with quantified nucleobases (Note: not a three-dimensional vector). Here, the normal multidimensional vector DV is represented as DV , and the multi-element vector EV is represented as EV . At this time, assuming that the ternary vector in miRNA is miR , the sequence length is n, and the element at the i–th site of the j–th nucleobases is e i j , the vector indicating the nucleobase is expressed by the following formula (Equation (2)):
miR = e 1 1 ,   e 1 2 ,   e 1 3 , ,   e 1 n ,   e 2 1 ,   e 2 2 ,   e 2 3 , ,   e 2 n ,   e 3 1 ,   e 3 2 ,   e 3 3 , ,   e 3 n
In applying to disease, it was verified by four methods.
The first method was based on a mathematical method and scored with element vector size (VS). VS is expressed by the following equation from Equation (2) (Equation (3)):
VS = e 1 1 2 + e 1 2 2 + + e 1 n 2 ,   e 2 1 2 + e 2 2 2 + + e 2 n 2 ,   e 3 1 2 + e 3 2 2 + + e 3 n 2
The second method is based on the electromagnetic method and the sum of electric charges (Sum) was used (Equation (4)):
Sum = e 1 1 + e 1 2 + + e 1 n ,   e 2 1 + e 2 2 + + e 2 n ,   e 3 1 + e 3 2 + + e 3 n
The third method is based on the electromagnetic method and used the centre of the electric field vector (EV_C) of a miRNA structure. When the coordinates of the element at the i–th site of the j–th nucleobases are x i j ,   y i j ,   z i j , EV _ C is expressed by the following formula (Equation (5)). The vector element of EV_C is the x, y, and z component of the electric field vector, unlike the above two methods:
EV _ C i = 1 3 j = 1 n e i j x k x i j x k x i j 2 + y l y i j 2 + z m z i j 2   3 2 ,   i = 1 3 j = 1 n e i j y k y i j x k x i j 2 + y l y i j 2 + z m z i j 2   3 2 , i = 1 3 j = 1 n e i j z k z i j x k x i j 2 + y l y i j 2 + z m z i j 2   3 2
The fourth method is based on the electromagnetic method, and was the sum of the electric field vectors (EV_S) on the lattice points in the miRNA structure. The lattice points were included in the structure, and the maximum distance of each axis of x, y, and z was divided into 11 equal parts. When the lattice point at the a –th position among the 11 divisions is ( x a , y a , z a ), EV _ S is expressed by the following expression (Equation (6)). Note that the vector elements of EV_S are the x, y, and z components of the electric field vector as in EV _ C :
E V _ S k = 5 5 l = 5 5 m = 5 5 i = 1 3 j = 1 n e i j x k x i j x k x i j 2 + y l y i j 2 + z m z i j 2   3 2 ,   k = 5 5 l = 5 5 m = 5 5 i = 1 3 j = 1 n e i j y k y i j x k x i j 2 + y l y i j 2 + z m z i j 2   3 2 ,   k = 5 5 l = 5 5 m = 5 5 i = 1 3 j = 1 n e i j z k z i j x k x i j 2 + y l y i j 2 + z m z i j 2   3 2
In calculating EV _ C and EV _ S , the miRNA structure was regarded as a torus. This is because the existence of circular RNA is known and the topologically toric body has a high degree of freedom. Three types of miRNA structures were defined (Figure 2). In the structure X, EV _ C and EV _ S is defined as EV _ C _ X and EV _ S _ X , respectively. In each structure, the distance between bases was 12.9   , the nucleobases length was 8.6   , and the distance between sites was 2.385   . In Structure C, the angle between the nucleobases and the surface formed by the torus was 45 ° .

2.4. Application to Diseases

The miR calculated by each method was used as the miRNA score ( miR ). When FC or log2FC are x and the score of the disease obtained therefrom is Score Multiplication Expression Level (SMEL), expressed by the following formula (Equation (7):
SMEL miR = miR x

3. Results

3.1. Evaluation of Each Scoring Method

The charge at each AUHB of each nucleobase was calculated by Firefly (Table 1).
For example, since the first site of A is an amino group, there are two hydrogen atoms (H) that form hydrogen bonds. Although different values were obtained in the calculation results, the average value was defined as the H charge in this study. In the Watson–Crick pair, A and U are hydrogen bonded only at the first site and the second site. This is because the distance between H present at the third site of A and the oxygen atom (O) present at the third site of U is large. However, U is known to form a wobble base pair with G. In the wobble base pairing, the second site of U and the first site of G, and the third site of U and the second site of G, are hydrogen bonded. Therefore, the third site of U was also used for scoring miRNA. Regarding the third site of A, it was used for scoring miRNA for the following two reasons: (1) since the three scoring methods except VS are based on electromagnetic methods and (2) to unify into the same ternary vector as the other three nucleobases.
The proportion of unique values in miR of each method was verified. In addition, there are 24 duplications in the miRNA sequence and the total number of unique sequences is 2632 (Table 2).
EV_S_B and EV_S_C have a nearly 1: 1 relationship with the miRNA sequence. On the other hand, VS was the most duplicated. EV_C type and EV_S type differ greatly in the degree of overlap, and EV_C type had more overlap. In addition, the EV_C/S type had more overlap in structure A than in others.
Furthermore, single regression analysis was performed for each method (Table 3). The First element or the x component of miR was taken as the x-axis, the Second element or the y component of miR was taken as the y-axis, and the Third element or the z component of miR was taken as the z-axis. Sum had the lowest correlation compared to other methods.

3.2. Single Regression Analysis for SMEL

MiRNAs with significantly different expression levels compared to healthy controls were obtained from the paper data for four diseases (Table S1). If both healthy control and patient miRNA expression data were obtained, FC and log2FC values were calculated. When only the FC value was obtained, the log2FC value was calculated, and when only the log2FC value was obtained, the FC value was calculated. In each of the FC groups and log2FC groups for the miRNA group thus obtained, SMEL by each miR was calculated. As an example, SMEL V S of AD is shown (Figure 3).
As a result, in both FC (Figure 3a) and log2FC (Figure 3b), a very high correlation was suggested for each combination of x-axis, y-axis, and z-axis. In other words, each site that is an independent variable is considered to have multicollinearity. However, each independent variable should not be synthesized because it depends on the hydrogen bonding site of the nucleobase. Therefore, instead of multiple regression analysis, single regression analysis was performed with all independent variable combinations (Table 4). As a result, both FC and log2FC were found to be highly correlated except for Sum, EV_C_B, and EV_C_C.

3.3. Validation of the Most Suitable Method

Since SMEL is the product of miR and FC or log2FC, the intercept is zero. Therefore, the regression line for each disease is non-parallel. Thus, whether there was a significant difference in the regression line of SMEL for each disease was examined by standard deviation instead of covariance analysis. In VS, it was shown whether there was a significant difference in the inclination of the x-axis and y-axis (Table 5).
When calculated by FC, it was shown that there was a significant difference in slope in all combinations when ±2σ. On the other hand, in the case of log2FC, there was a significant difference in ST and other diseases, but there was no significant difference in other combinations. Based on Table 5, the case where there was a significant difference in all combinations of diseases was marked as ◯, and the case where there was a combination with no significant difference was marked as × (Table 6).
Then, there was an axis with a significant difference in FC, whereas there was no axis with a significant difference in log2FC. In the case of ±2σ, the number of ◯ exceeded 50% in EV_C type, EV_S_B, and EV_S_C. Since EV_C_A has only the x –axis and y –axis, the number of ◯ was interpreted as 100%. However, in the single regression analysis of each disease, EV_C_B and EV_C_C p > 0.1, and the correlation was low.

3.4. Reassessment Considering Outliers

Unlike log2FC, FC may have a very large value (Figure 3a). Therefore, when the FC values are sorted in descending order and a value that is twice as large, or larger, than the next smallest value appears, the larger value is regarded as an outlier.
Unlike the case where outliers were included, the same calculation was performed except without the outliers. However, ±3σ was all ×, so ±1σ and ±2σ were examined (Table 7). When the outliers were excluded, the only significant difference in ±2σ was in the yz axis in EV_S_B. In addition, ±1σ in EV_S_B did not change from ±2σ.

4. Discussion

Since miRNA is deeply involved in various biological reactions including diseases, it is important to know the mechanism in detail in the treatment of diseases [20,21]. The mechanism by which miRNA suppresses translation has been clarified [22], but the meaning of the sequence when miRNA is regarded as a gene has not yet been clarified. The seed region is considered important in the seed theory, but the array outside the seed region has not been fully elucidated. Even in the seed region, the relationship with the disease is unclear in terms of the base sequence. Elucidating the relationship between miRNA sequences and diseases is the first step to deciphering miRNA genes. However, it seems that conventional methods are difficult in terms of decoding miRNA genes. Therefore, we focused on RNA–RNA interactions. Nucleobase hydrogen bonds are utilized when reacting with other nucleic acids. Therefore, the base sequence was digitized by using AUHB as the miRNA score, and the relationship with the disease was verified.
In order to quantify the miRNA sequence, the charge was calculated focusing on AUHB (Table 1). The H was a positive value, and O and nitrogen atoms were negative values. In the Watson–Crick pair, the third site of A does not bond with hydrogen, and in this charge calculation, only the charge of the third site of A was smaller in absolute value than other atoms. Therefore, this calculation result was interpreted as appropriate.
In this study, miRNA was considered a torus. There are two reasons for this. First of all, it is known that there are circular forms of various RNA types. For example, circular RNAs have been reported [23,24]. In addition, tRNA is clover-shaped and can be regarded as circular in terms of topology. Pre-mRNA [25] and phytopathogenic virus viroid [26] are also circular. Although we have not found any miRNA reports that are circular, miRNA secondary structure suggests a circular shape. For example, the secondary structure of mature miRNA can be predicted with CentroidFold [27]. When predicting the secondary structure of miRNA in CentroidFold, miRNA basically has a dumbbell shape, but some miRNAs become torus (e.g., hsa-miR-1-5p; inference engine: CONTRAfold; weight of base pairs: 22). The dumbbell type is also topologically similar to a torus. Second, considering miRNA as a torus has two advantages. The first is that the clover type and dumbbell type can be considered topologically homologous, so the method of this study may be applicable to other types of RNA. Secondly, the readability of the miRNA genetic code sequence is high. In order to decipher the genetic code of miRNA sequences as the ultimate goal, this study aims to search for pre-data processing methods for machine learning. It is probably difficult to perform data preprocessing because the amount of information is too small for a one-dimensional array. On the other hand, there is too much information in an accurate secondary structure or tertiary structure. A torus can be viewed as an approximation of secondary structure.
Based on the calculation results, each miR was calculated. In order to clarify the relationship between miRNA sequences and diseases, it is desirable that miR and sequences have a 1: 1 relationship. Therefore, the ratio of unique values in each miR was verified (Table 2). The EV_C type and EV_S type differ greatly in the degree of overlap, and the EV_C type had more overlap. This is thought to be because the EV_C type is a vector with only the centre, so the influence from the charge is small, whereas the EV_S type is the sum of the electric field vectors in the torus, so it is strongly influenced by the charge. In addition, the EV_C/S type had more overlap in structure A than in others. This is presumably because the z component of the electric field vector does not exist. Surprisingly, Sum had a low overlap rate compared to the EV_C type and EV_S_A. From the viewpoint of unique values, Sum, EV_S_B, and EV_S_C were considered suitable.
A single regression analysis was performed to examine the distribution trend of each miR value (Table 3). Overall correlation was high except for Sum. B is a value uniquely determined by the nucleobase, but B itself did not show similarities in the values of each element except for the values of the second and third sites of G. Nevertheless, it was unexpected that the correlation was high. Note that the strong correlation in the order of VS, EV_S, EV_C, and Sum is thought to be due to the difference in strength affected by B . From the viewpoint of correlation, VS and EV_S types were considered suitable.
Since the correlation was high in each miR , it was suggested that SMEL could also have a high correlation. When tested in four diseases, VS and EV_S types showed a high correlation in SMEL, and EV_C type and Sum showed low correlation (Figure 3, Table 4). As described above, since the SMEL intercept is 0, the regression line of each disease is non–parallel. If there is no relationship between the miRNA sequence and the disease, it should be almost the same as the slope in miR . However, there was actually a significant difference in the slope (Table 3 and Table 4). In other words, the high correlation of SMEL is thought to be due to the correlation between the disease and miRNA, although it is influenced by the high correlation of miR .
Therefore, it was verified whether the disease can be distinguished by SMEL (Table 5 and Table 6). As a result, EV_C_A, EV_S_B, or EV_S_C was considered a suitable method. However, there may be large outliers in the case of FC (Figure 3a,b). Since these are linear regressions, outliers can have a significant impact in the case of FC. Therefore, it was verified even when outliers were not taken into account (Table 7). Based on the results with and without outliers, the most suitable miR was considered to be EV_S_B.
EV_S_B was suitable for unique values, correlations, and SMEL. The reason why EV_S_C, which showed the same tendency as EV_S_B, did not show a significant difference at 2σ of SMEL, seems to be because EV_S_B was closest to the possible structure.
It was suggested that the classification of the disease is possible by the electric field vector made by AUHB (Figure 4), but there are three main problems in this result. The first is that the structure is simplified. In the case of EV_S_B, which seemed most suitable, the First to Third sites were placed perpendicular to the plane of the torus (Figure 2b). However, the actual nucleobase of miRNA is not vertical, but seems to be angled by the π–π interaction and the phosphate–sugar chain interaction. Therefore, it is necessary to perform calculations in consideration of the angle. The second is that only four types of disease are compared. Thus, in order to prove the hypothesis that “the relationship between the miRNA sequence and disease may be revealed by focusing on hydrogen bonding sites in RNA–RNA interaction”, more diseases must be targeted. Furthermore, in this study, comparisons were made between different diseases in major classifications, but it is also necessary to verify that diseases with the same large category, such as cancer types, can be distinguished. The third is interpretation of the score. Since EV_S_B is almost sequence-specific, the results of this study are considered to be one step closer to miRNA gene decoding. However, as mentioned in the second problem, there are still only few subjects tested in the experiment. Therefore, it is not possible at this time to interpret the score. In order to interpret the score, it is necessary to collect a lot of data and perform machine learning. Therefore, it is necessary to continue research in the future.

5. Conclusions

By considering RNA as a ternary vector, it was suggested that the sequence and its expression may be related to the disease. This study is a step toward clarifying the meaning of the sequence of miRNA genes by a different approach from the conventional “method of focusing on complementary sequences of miRNA” and “analysis of gene networks”. Until now, the gene has been interpreted as a one-dimensional sequence of A/U/G/C. From now on, it can be expected that a new viewpoint for RNA-RNA interaction will be obtained by looking at a ternary vector using AUHB.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4409/8/12/1615/s1, Table S1: reference for dataset of each disease.

Author Contributions

T.O. contributed to all this work. N.Y. supervised the project.

Funding

This work was supported the grants –in –aid of MEXT, Japan [17H06105] to Naohiro Yoshida.

Acknowledgments

We are grateful to Dr. T. Ueki for comments on earlier version of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. World Health Organization. The Top 10 Causes of Death. Available online: https://www.who.int/news-room/fact-sheets/detail/the-top-10-causes-of-death (accessed on 3 September 2019).
  2. Mirzaei, H.; Momeni, F.; Saadatpour, L.; Sahebkar, A.; Goodarzi, M.; Masoudifar, A.; Kouhpayeh, S.; Salehi, H.; Mirzaei, H.R.; Jaafari, M.R. MicroRNA: Relevance to stroke diagnosis, prognosis, and therapy. J. Cell. Physiol. 2018, 233, 856–865. [Google Scholar] [CrossRef] [PubMed]
  3. Conickx, G.; Mestdagh, P.; Avila Cobos, F.; Verhamme, F.M.; Maes, T.; Vanaudenaerde, B.M.; Seys, L.J.M.; Lahousse, L.; Kim, R.Y.; Hsu, A.C.; et al. MicroRNA Profiling Reveals a Role for MicroRNA-218-5p in the Pathogenesis of Chronic Obstructive Pulmonary Disease. Am. J. Respir. Crit. Care Med. 2016, 195, 43–56. [Google Scholar] [CrossRef] [PubMed]
  4. Silvestro, S.; Bramanti, P.; Mazzon, E. Role of miRNAs in Alzheimer’s Disease and Possible Fields of Application. Int. J. Mol. Sci. 2019, 20, 3979. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Sabir, N.; Hussain, T.; Shah, S.Z.A.; Peramo, A.; Zhao, D.; Zhou, X. miRNAs in Tuberculosis: New Avenues for Diagnosis and Host-Directed Therapy. Front. Microbiol. 2018, 9, 602. [Google Scholar] [CrossRef] [Green Version]
  6. Lu, T.X.; Rothenberg, M.E. MicroRNA. J. Allergy Clin. Immunol. 2018, 141, 1202–1207. [Google Scholar] [CrossRef] [Green Version]
  7. Dong, H.; Lei, J.; Ding, L.; Wen, Y.; Ju, H.; Zhang, X. MicroRNA: Function, Detection, and Bioanalysis. Chem. Rev. 2013, 113, 6207–6233. [Google Scholar] [CrossRef]
  8. Mohr, A.M.; Mott, J.L. Overview of microRNA biology. Sem. Liver Dis. 2015, 35, 3–11. [Google Scholar] [CrossRef] [Green Version]
  9. Luo, J.; Ding, P.; Liang, C.; Cao, B.; Chen, X. Collective Prediction of Disease-Associated miRNAs Based on Transduction Learning. IEEE/ACM Trans. Comput. Biol. Bioinform. 2017, 14, 1468–1475. [Google Scholar] [CrossRef]
  10. Guo, Y. Clinical significance of serum MicroRNA-203 in patients with acute myeloid leukemia. Bioengineered 2019, 10, 345–352. [Google Scholar] [CrossRef] [Green Version]
  11. Zhou, Y.; Shen, S. MiR-520f acts as a biomarker for the diagnosis of lung cancer. Medicine 2019, 98, e16546. [Google Scholar] [CrossRef]
  12. Zou, Q.; Li, J.; Hong, Q.; Lin, Z.; Wu, Y.; Shi, H.; Ju, Y. Prediction of MicroRNA-Disease Associations Based on Social Network Analysis Methods. BioMed Res. Int. 2015, 2015, 810514. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Zou, Q.; Li, J.; Song, L.; Zeng, X.; Wang, G. Similarity computation strategies in the microRNA-disease network: A survey. Brief. Funct. Genom. 2015, 15, 55–64. [Google Scholar] [CrossRef] [PubMed]
  14. Lan, W.; Wang, J.; Li, M.; Liu, J.; Wu, F.X.; Pan, Y. Predicting MicroRNA-Disease Associations Based on Improved MicroRNA and Disease Similarities. IEEE/ACM Trans. Comput. Biol. Bioinform. 2018, 15, 1774–1782. [Google Scholar] [CrossRef] [PubMed]
  15. Fu, L.; Peng, Q. A deep ensemble model to predict miRNA-disease association. Sci. Rep. 2017, 7, 14482. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Ding, T.; Xu, J.; Sun, M.; Zhu, S.; Gao, J. Predicting microRNA biological functions based on genes discriminant analysis. Comput. Biol. Chem. 2017, 71, 230–235. [Google Scholar] [CrossRef]
  17. Zheng, K.; You, Z.H.; Wang, L.; Zhou, Y.; Li, L.P.; Li, Z.W. MLMDA: A machine learning approach to predict and validate MicroRNA-disease associations by integrating of heterogenous information sources. J. Transl. Med. 2019, 17, 260. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Granovsky, A.A. Firefly Version 8. Available online: http://classic.chem.msu.su/gran/gamess/index.html (accessed on 29 August 2019).
  19. Shmidt, M.W.; Baldrosge, K.K.; Elbert, J.A.; Gordon, M.S.; Enseh, J.H.; Koseki, S.; Matsunaga, N.; Nguyen, K.A.; Su, S.; Windus, T.L.; et al. General Atomic and Molecular Electronic Structure Systems. J. Comput. Chem. 1993, 14, 1347–1363. [Google Scholar] [CrossRef]
  20. Qin, L.; Kang, A. Epigenetic research progress in colorectal cancer. J. Cent. South Univ. Med. Sci. 2019, 44, 830–836. [Google Scholar] [CrossRef]
  21. Xu, G.; Li, N.; Zhang, Y.; Zhang, J.; Xu, R.; Wu, Y. MicroRNA-383-5p inhibits the progression of gastric carcinoma via targeting HDAC9 expression. Braz. J. Med. Biol. Res. 2019, 52, e8341. [Google Scholar] [CrossRef] [Green Version]
  22. Kura, B.; Parikh, M.; Slezak, J.; Pierce, G.N. The Influence of Diet on MicroRNAs that Impact Cardiovascular Disease. Molecules 2019, 24, 1509. [Google Scholar] [CrossRef] [Green Version]
  23. Jeck, W.R.; Sorrentino, J.A.; Wang, K.; Slevin, M.K.; Burd, C.E.; Liu, J.; Marzluff, W.F.; Sharpless, N.E. Circular RNAs are abundant, conserved, and associated with ALU repeats. RNA 2013, 19, 141–157. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Pamudurti, N.R.; Bartok, O.; Jens, M.; Ashwal-Fluss, R.; Stottmeister, C.; Ruhe, L.; Hanan, M.; Wyler, E.; Perez-Hernandez, D.; Ramberger, E.; et al. Translation of CircRNAs. Mol. Cell 2017, 66, 9–21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Welden, J.R.; Stamm, S. Pre-mRNA structures forming circular RNAs. Biochim. Biophys. Acta (BBA) Gene Regul. Mech. 2019, 1862, 194410. [Google Scholar] [CrossRef] [PubMed]
  26. Sanger, H.L.; Klotz, G.; Riesner, D.; Gross, H.J.; Kleinschmidt, A.K. Viroids are single-stranded covalently closed circular RNA molecules existing as highly base-paired rod-like structures. Proc. Natl. Acad. Sci. USA 1976, 73, 3852–3856. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Sato, K.; Hamada, M.; Asai, K.; Mituyama, T. CENTROIDFOLD: A web server for RNA secondary structure prediction. Nucleic Acids Res. 2009, 37, W277–W280. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Site where nucleobase hydrogen bonds are located. Each of the three sites was designated as the First site (red circle), the Second site (green circle), and the Third site (blue circle). In addition, third site includes Adenine C2 hydrogen that is not used for RNA–RNA interactions.
Figure 1. Site where nucleobase hydrogen bonds are located. Each of the three sites was designated as the First site (red circle), the Second site (green circle), and the Third site (blue circle). In addition, third site includes Adenine C2 hydrogen that is not used for RNA–RNA interactions.
Cells 08 01615 g001
Figure 2. Placement of nucleobases relative to the circle (a) Structure that is most susceptible to hydrogen bonding with other RNA (Structure A) (b) Structure that is most stable (Structure B) (c) Structure between A and B (Structure C) (d) Figure excerpted and enlarged from part of structure C.
Figure 2. Placement of nucleobases relative to the circle (a) Structure that is most susceptible to hydrogen bonding with other RNA (Structure A) (b) Structure that is most stable (Structure B) (c) Structure between A and B (Structure C) (d) Figure excerpted and enlarged from part of structure C.
Cells 08 01615 g002aCells 08 01615 g002b
Figure 3. 3D scatter plot of VS calculation results in AD. The X-axis was the First site score, the Y-axis was the Second site score, and the Z-axis was the Third site score. (a) Result when using fold change. (b) Figure enlarging the range of 0 to 200 in Figure 3a. (c) Result when using log2FC.
Figure 3. 3D scatter plot of VS calculation results in AD. The X-axis was the First site score, the Y-axis was the Second site score, and the Z-axis was the Third site score. (a) Result when using fold change. (b) Figure enlarging the range of 0 to 200 in Figure 3a. (c) Result when using log2FC.
Cells 08 01615 g003
Figure 4. Slope of 3D regression line for each disease.
Figure 4. Slope of 3D regression line for each disease.
Cells 08 01615 g004
Table 1. Charge at the site of hydrogen bonding in base pair.
Table 1. Charge at the site of hydrogen bonding in base pair.
NucleobaseFirstSecondThird
Adenine0.417−0.5140.228
Uracil−0.5060.408−0.471
Guanine−0.4900.4140.413
Cytosine0.379−0.558−0.476
Table 2. Duplicate in each scoring method.
Table 2. Duplicate in each scoring method.
MethodsTotal Number of Duplicate ValuesTotal Number of Unique ValuesPercentage of Unique Values
VS960169664
Sum454220283
EV_C_A832182469
EV_C_B706195073
EV_C_C715194173
EV_S_A615204177
EV_S_B78257897
EV_S_C65259198
ref*24263299
* For reference, the duplicates of miRNA sequences were presented.
Table 3. 3D regression line of miRNA in each scoring method.
Table 3. 3D regression line of miRNA in each scoring method.
Methodsxyyzzx
SlopeS.D.p-Valuer2SlopeS.D.p-Valuer2SlopeS.D.p-Valuer2
VS1.040.012.E − 1231.000.860.012.E − 1261.001.100.014.E − 1331.00
Sum−0.650.071.E − 14−0.670.450.091.E − 060.450.130.102.E − 010.11
EV_C_A3.390.014.E − 1501.00NANANANANANANANA
EV_C_B3.380.013.E − 1481.000.020.009.E − 070.462.540.502.E − 060.45
EV_C_C3.390.011.E − 1491.000.020.001.E − 120.635.360.672.E − 120.62
EV_S_A1.410.052.E − 480.93NANANANANANANANA
EV_S_B3.560.011.E − 1491.0011.080.017.E − 2161.000.030.004.E − 1531.00
EV_S_C3.560.016.E − 1501.0011.120.032.E − 1721.000.030.002.E − 1411.00
NA: The item could not be calculated because there was no value.
Table 4. 3D regression line of AD in each scoring method.
Table 4. 3D regression line of AD in each scoring method.
Methodsxyyzzx
SlopeS.D.p-Valuer2SlopeS.D.p-Valuer2SlopeS.D.p-Valuer2
log2FCVS1.040.000.E + 001.000.850.000.E + 001.001.130.000.E + 001.00
Sum−0.250.042.E − 09−0.300.120.081.E − 010.080.160.042.E − 050.21
EV_C_A3.310.010.E + 001.00NANANANANANANANA
EV_C_B3.280.020.E + 001.000.020.007.E − 150.382.240.281.E − 140.38
EV_C_C3.290.010.E + 001.000.020.008.E − 330.564.770.367.E − 330.56
EV_S_A1.250.105.E − 330.56NANANANANANANANA
EV_S_B3.520.010.E + 001.0011.100.010.E + 001.000.030.000.E + 001.00
EV_S_C3.530.010.E + 001.0011.100.010.E + 001.000.030.000.E + 001.00
FCVS1.030.000.E + 001.000.870.000.E + 001.001.110.000.E + 001.00
Sum0.060.025.E − 050.20−0.640.144.E − 06−0.240.880.042.E − 650.73
EV_C_A3.260.000.E + 001.00NANANANANANANANA
EV_C_B3.260.000.E + 001.00−0.010.001.E − 15−0.39−5.020.612.E − 15−0.39
EV_C_C3.260.000.E + 001.000.000.006.E − 010.030.611.025.E − 010.04
EV_S_A1.390.000.E + 001.00NANANANANANANANA
EV_S_B3.410.000.E + 001.0011.180.000.E + 001.000.030.000.E + 001.00
EV_S_C3.430.000.E + 001.0010.950.010.E + 001.000.030.000.E + 001.00
Table 5. Comparison of x − y slope of 3D regression line in VS.
Table 5. Comparison of x − y slope of 3D regression line in VS.
Expression ratio
COPDSTTBCOPDSTTB
log2FCAD××××
COPD × ×
ST
FCAD××
COPD ×
ST
Table 6. Comparison of slope of the 3D regression line.
Table 6. Comparison of slope of the 3D regression line.
Methods
x − yy − zz − xx − yy − zz − x
log2FCVS××××××
Sum××××××
EV_C_A×NANA×NANA
EV_C_B××××××
EV_C_C××××××
EV_S_A×NANA×NANA
EV_S_B××××××
EV_S_C××××××
FCVS×××××
Sum×××××
EV_C_ANANANANA
EV_C_B××
EV_C_C×××
EV_S_A×NANA×NANA
EV_S_B×
EV_S_C×××
Table 7. Comparison of slope of the 3D regression line without some outliers.
Table 7. Comparison of slope of the 3D regression line without some outliers.
Methods
x − yy − zz − xx − yy − zz − x
log2FCVS××××××
Sum××××××
EV_C_ANANA×NANA
EV_C_B××××××
EV_C_C×××××
EV_S_A×NANA×NANA
EV_S_B××××××
EV_S_C××××××
FCVS×××××
Sum×××××
EV_C_A×NANA×NANA
EV_C_B××××
EV_C_C××××××
EV_S_ANANA×NANA
EV_S_B××××
EV_S_C××××

Share and Cite

MDPI and ACS Style

Osone, T.; Yoshida, N. The Relationship Between the miRNA Sequence and Disease May be Revealed by Focusing on Hydrogen Bonding Sites in RNA–RNA Interactions. Cells 2019, 8, 1615. https://doi.org/10.3390/cells8121615

AMA Style

Osone T, Yoshida N. The Relationship Between the miRNA Sequence and Disease May be Revealed by Focusing on Hydrogen Bonding Sites in RNA–RNA Interactions. Cells. 2019; 8(12):1615. https://doi.org/10.3390/cells8121615

Chicago/Turabian Style

Osone, Tatsunori, and Naohiro Yoshida. 2019. "The Relationship Between the miRNA Sequence and Disease May be Revealed by Focusing on Hydrogen Bonding Sites in RNA–RNA Interactions" Cells 8, no. 12: 1615. https://doi.org/10.3390/cells8121615

APA Style

Osone, T., & Yoshida, N. (2019). The Relationship Between the miRNA Sequence and Disease May be Revealed by Focusing on Hydrogen Bonding Sites in RNA–RNA Interactions. Cells, 8(12), 1615. https://doi.org/10.3390/cells8121615

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop