Team:LZU-CHINA/Model

Model and Bioinformatic data

Model

The model can help us predict the experimental results and guide our experimental design. We designed two models, one to simulate the effect of increasing exosome secretion. It simulates the effect of different exosomes enhancing genes and exosomes booster on the experimental results, providing guidance for our experiment. Another model was used to evaluate the effectiveness of intravenous administration of tetracycline to avoid adverse effects from high concentrations of tetracycline. Therefore, we established a drug metabolism model induced by tetracycline for intravenous administration to evaluate the efficiency of microRNA expression at different concentrations.

Model A (Tetracycline iv drug model)

To activate our tetracycline-induced promoter, we administered tetracycline intravenously. However, excessive tetracycline can affect gastrointestinal tract, liver, blood system, central nervous system and kidney to different degrees. Tetracycline can also bind with phosphate in dentin and enamel, leading to teeth becoming yellow and bone dysplasia. Therefore, we need to establish a drug metabolism model to analyze and calculate the optimal drug concentration range during intravenous injection, to reduce the impact of tetracycline administration and provide certain mathematical basis for the selection of drug injection time and drug injection amount in practical application.
a. The metabolism of tetracycline in the stomach conforms to the rapid perfusion model.
 the specific calculation formula is as follows :

The parameters are described below the form:

Here comes our result:

Figure1. Over time, tetracycline concentrations gradually drop in the blood and enter the stomach tissue. It was predicted that tetracycline in the stomach reached its peak at 10hours, and then its concentration gradually decreased with the degradation of tetracycline.

Figure2. This diagram shows the functional relationship between tetracycline and absorption time and absorption coefficient in gastric tissue. The tetracycline concentration decreases gradually with the time of absorption, and the absorption efficiency is inversely proportional to the tetracycline concentration in the gastric tissue.

b. when we give the drug intravenously, the change of blood drug concentration conforms to the following formula:

Here comes our result:

Figure3. With the increase of time, tetracycline concentration in blood medicine gradually increased until reaching the plateau stage.

Model B

To increase the expression of exosomes, we added three enhanced genes to simulate the effect of increasing exosomes secretion through the model, providing guidance for our experiment.
Equations is as follow:
a.transcriptional rate simulation of three booster genes:

 b. Simulation of three booster genes transcriptional translation:

Here comes our result:

Figure4. the promotional concentrations of the three promotional genes for the secretion of exosomes. Over time, the number of exosomes secreted by cells increased.

Figure5. Comparison between the expression levels of three promoter genes and those of non-promoter exosomes. It demonstrates that exosome booster could dramatically increase exosome secrion.

Figure6. production rate of exosomes in cell. This model predicts the production rate of exosome over time in cells.

Bioinformatics data

The effects of different microRNAs on gene expression and gene pathways were analyzed and predicted in the KEGG database. We predicted miR-135b-3p, miR-769-5p and miR-942-5p, respectively, the inhibitory ability of these three microRNAs to other genes, the gene pathways that may be affected and the interactions between these gene pathways and analyzed the interaction network between the same functional proteins that are affected. We also predicted the genetic effects of all three on human cells at the same time. These bioinformatics results suggest that this is the kind of microRNA that has complex regulatory networks. Combined with our phenotypic results (cck-8 analysis, transwell test, flow cytometry analysis, scratch test), miR-135b-3p, miR-769-5p and miR-942-5p finally affected the activity, migration and DNA synthesis of gastric cancer cells through these regulatory networks. In addition, bioinformatic analysis indicates that much work remains to be done. Many regulatory mechanisms still need to be validated and interpreted experimentally. Nevertheless, as an important regulatory device in human organism, microRNA has a promising prospect in the application of cancer treatment.

1.Biological information analysis of miR-135b-3p

Figure7. According to bioinformatics prediction analysis, miR-135b-3p can inhibit 50 genes expression. The redder the color, the higher the gene expression, and the bluer the color, the lower the gene expression. As you can see from the diagram. When miR-135b-3p expression was at a high level, most genes were at a low level. As the expression level of miR-135b-3p gradually decreased, most genes showed high expression level. However, some gene expression levels are not linearly correlated, which may be due to the interaction between genes.

Figure8. The gene pathways that miR-135b-3p may affected and the interactions between these gene pathways.

Figure9.The interaction network between the same functional proteins affected by miR-135b-3p.

Figure10. GO (Gene Ontology) analysis. The molecular function affected by miR-135b-3p were analyzed by GO database, and there was a total of 20 molecular function related to miR-135b-3p. The color and size of the histogram represent the P value. The darker the color, the smaller the P value, the stronger the correlation between this function and miR-135b-3p.

2. Biological information analysis of miR-942-5p

Figure11. According to bioinformatics prediction analysis, miR-942-5p can inhibit 50 genes expression. The redder the color, the higher the gene expression, and the bluer the color, the lower the gene expression. As you can see from the diagram. When miR-942-5p expression was at a high level, most genes were at a low level. As the expression level of miR-942-5p gradually decreased, most genes showed high expression level. However, some gene expression levels are not linearly correlated, which may be due to the interaction between genes.

Figure12. The gene pathways that miR-942-5p may affected and the interactions between these genes.

Figure13. The interaction network between the same functional proteins affected by miR-942-5p.

Figure14. GO (Gene Ontology) analysis. The molecular function affected by miR-942-5p were analyzed by GO database, and there was a total of 20 molecular function related to miR-942-5p. The color and size of the histogram represent the P value. The darker the color, the smaller the P value, the stronger the correlation between this function and miR-942-5p.

3. Biological information analysis of miR-769-5p

Figure15. According to bioinformatics prediction analysis, miR-769-5p can inhibit 50 genes expression. The redder the color, the higher the gene expression, and the bluer the color, the lower the gene expression. As you can see from the diagram. When miR-769-5p expression was at a high level, most genes were at a low level. As the expression level of miR-769-5p gradually decreased, most genes showed high expression level. However, some gene expression levels are not linearly correlated, which may be due to the interaction between genes.

Figure16. The gene pathways that miR-769-5p may affected and the interactions between these genes.

Figure17. The interaction network between the same functional proteins affected by miR-769-5p.

Figure18. According to bioinformatics prediction analysis, miR-769-5p can inhibit 50 genes expression. The redder the color, the higher the gene expression, and the bluer the color, the lower the gene expression. As you can see from the diagram. When miR-769-5p expression was at a high level, most genes were at a low level. As the expression level of miR-769-5p gradually decreased, most genes showed high expression level. However, some gene expression levels are not linearly correlated, which may be due to the interaction between genes.

4. Biological information analysis of miR-135b-3p, miR-942-5p and miR-769-5p interaction

Figure19. The gene pathways that miR-135b-3p, miR-942-5p and miR-769-5p may affected and the interactions between these genes were predicted by Biological information analysis.

Figure20. The interaction network between the same functional proteins affected by miR-135b-3p, miR-942-5p and miR-769-5p.

Figure21. The functional molecules affected by miR-135b-3p, miR-942-5p and miR-769-5p were analyzed by GO database. The redder the color, the higher the gene expression, and the bluer the color, the lower the gene expression. As you can see from the diagram. When miR-135b-3p, miR-942-5p and miR-769-5p expression was at a high level, most genes were at a low level. As the expression level of miR-135b-3p, miR-942-5p and miR-769-5p gradually decreased, most genes showed high expression level. However, some gene expression levels are not linearly correlated, which may be due to the interaction between genes.

Reference

1.Sipponen, P., Kekki, M., Haapakoski, J., Ihamäki, T., & Siurala, M. (2010). Gastric cancer risk in chronic atrophic gastritis: statistical calculations of cross-sectional data. International Journal of Cancer,35(2), 173-177.
2.Zhu, Z. (2018). [2017 hotspots review and outlook on gastric cancer surgery in china]. Zhonghua Wei Chang Wai Ke Za Zhi, 21(1), 7-14.
3.Zhang, Z. Y., Zheng-Qi, W. U., & Wang, H. J. (2009). An analysis of gastric cancer detection and epidemicity in high incidental area with gastric cancer,wuwei city,gansu province. China Cancer, 19(3), 194–201.
4.Chan, S. L., & So, J. B. (2016). Gastric cancer therapy.
5.Barr, W. H., Gerbracht, L. M., Letcher, K., Plaut, M., & Strahl, N. (1972). Assessment of the biologic availability of tetracycline products in man. Clinical Pharmacology & Therapeutics, 13(1), 97–108.