Integrative Network Analysis of Predicted miRNA-Targets Regulating Expression of Immune Response Genes in Bovine Coronavirus Infection

Dissecting the molecular mechanism of cepharanthine against COVID-19, based on a network pharmacology strategy combined with RNA-sequencing analysis, molecular docking, and molecular dynamics simulation

OBJECTIVES

Recently, it has been reported that cepharanthine (CEP) is highly likely to be an agent against Coronavirus disease 2019 (COVID-19). In the present study, a network pharmacology-based approach combined with RNA-sequencing (RNA-seq), molecular docking, and molecular dynamics (MD) simulation was performed to determine hub targets and potential pharmacological mechanism of CEP against COVID-19.

METHODS

Targets of CEP were retrieved from public databases. COVID-19-related targets were acquired from databases and RNA-seq datasets GSE157103 and GSE155249. The potential targets of CEP and COVID-19 were then validated by GSE158050. Hub targets and signaling pathways were acquired through bioinformatics analysis, including protein-protein interaction (PPI) network analysis and enrichment analysis. Subsequently, molecular docking was carried out to predict the combination of CEP with hub targets. Lastly, MD simulation was conducted to further verify the findings.

RESULTS

A total of 700 proteins were identified as CEP-COVID-19-related targets. After the validation by GSE158050, 97 validated targets were retained. Enrichment results indicated that CEP acts on COVID-19 through multiple pathways, multiple targets, and overall cooperation. Specifically, PI3K-Akt signaling pathway is the most important pathway. Based on PPI network analysis, 9 central hub genes were obtained (ACE2, STAT1, SRC, PIK3R1, HIF1A, ESR1, ERBB2, CDC42, and BCL2L1). Molecular docking suggested that the combination between CEP and 9 central hub genes is extremely strong. Noteworthy, ACE2, considered the most important gene in CEP against COVID-19, binds to CEP most stably, which was further validated by MD simulation.

CONCLUSION

Our study comprehensively illustrated the potential targets and underlying molecular mechanism of CEP against COVID-19, which further provided the theoretical basis for exploring the potential protective mechanism of CEP against COVID-19.

Deciphering inhibitory mechanism of coronavirus replication through host miRNAs-RNA-dependent RNA polymerase interactome

Despite what we know so far, Covid-19, caused by SARS-CoV-2 virus, remains a pandemic that still require urgent healthcare intervention. The frequent mutations of the SARS-CoV-2 virus has rendered disease control with vaccines and antiviral drugs quite challenging, with newer variants surfacing constantly. There is therefore the need for newer, effective and efficacious drugs against coronaviruses. Considering the central role of RNA dependent, RNA polymerase (RdRp) as an enzyme necessary for the virus life cycle and its conservation among coronaviruses, we investigated potential host miRNAs that can be employed as broad-range antiviral drugs averse to coronaviruses, with particular emphasis on BCoV, MERS-CoV, SARS-CoV and SARS-CoV-2. miRNAs are small molecules capable of binding mRNA and regulate expression at transcriptional or translational levels. Our hypothesis is that host miRNAs have the potential of blocking coronavirus replication through miRNA-RdRp mRNA interaction. To investigate this, we retrieved the open reading frame (ORF1ab) nucleotide sequences and used them to interrogate miRNA databases for miRNAs that can bind them. We employed various bioinformatics tools to predict and identify the most effective host miRNAs. In all, we found 27 miRNAs that target RdRp mRNA sequence of multiple coronaviruses, of which three - hsa-miR-1283, hsa-miR-579-3p, and hsa-miR-664b-3p target BCoV, SARS-CoV and SARS-CoV-2. Additionally, hsa-miR-374a-5p has three bovine miRNA homologs viz bta-miR-374a, bta-miR-374b, and bta-miR-374c. Inhibiting the expression of RdRp enzyme via non-coding RNA is novel and of great therapeutic importance in the control of coronavirus replication, and could serve as a broad-spectrum antiviral, with hsa-miR-1283, hsa-miR-579-3p, and hsa-miR-664b-3p as highly promising.

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.

Integrative COVID-19 biological network inference with probabilistic core decomposition

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for millions of deaths around the world. To help contribute to the understanding of crucial knowledge and to further generate new hypotheses relevant to SARS-CoV-2 and human protein interactions, we make use of the information abundant Biomine probabilistic database and extend the experimentally identified SARS-CoV-2-human protein-protein interaction (PPI) network in silico. We generate an extended network by integrating information from the Biomine database, the PPI network and other experimentally validated results. To generate novel hypotheses, we focus on the high-connectivity sub-communities that overlap most with the integrated experimentally validated results in the extended network. Therefore, we propose a new data analysis pipeline that can efficiently compute core decomposition on the extended network and identify dense subgraphs. We then evaluate the identified dense subgraph and the generated hypotheses in three contexts: literature validation for uncovered virus targeting genes and proteins, gene function enrichment analysis on subgraphs and literature support on drug repurposing for identified tissues and diseases related to COVID-19. The major types of the generated hypotheses are proteins with their encoding genes and we rank them by sorting their connections to the integrated experimentally validated nodes. In addition, we compile a comprehensive list of novel genes, and proteins potentially related to COVID-19, as well as novel diseases which might be comorbidities. Together with the generated hypotheses, our results provide novel knowledge relevant to COVID-19 for further validation.

Regulatory network of miRNA, lncRNA, transcription factor and target immune response genes in bovine mastitis

Pre- and post-transcriptional modifications of gene expression are emerging as foci of disease studies, with some studies revealing the importance of non-coding transcripts, like long non-coding RNAs (lncRNAs) and microRNAs (miRNAs). We hypothesize that transcription factors (TFs), lncRNAs and miRNAs modulate immune response in bovine mastitis and could potentially serve as disease biomarkers and/or drug targets. With computational analyses, we identified candidate genes potentially regulated by miRNAs and lncRNAs base pair complementation and thermodynamic stability of binding regions. Remarkably, we found six miRNAs, two being bta-miR-223 and bta-miR-24-3p, to bind to several targets. LncRNAs NONBTAT027932.1 and XR_003029725.1, were identified to target several genes. Functional and pathway analyses revealed lipopolysaccharide-mediated signaling pathway, regulation of chemokine (C-X-C motif) ligand 2 production and regulation of IL-23 production among others. The overarching interactome deserves further in vitro/in vivo explication for specific molecular regulatory mechanisms during bovine mastitis immune response and could lay the foundation for development of disease markers and therapeutic intervention.