In-silico investigations of selective miRNA-gene targets and their validation studies in obstructive sleep apnea (OSA) patient cohorts

BackgroundObstructive sleep apnoea (OSA) is a prevalent form of sleep disordered breathing which results in sleep fragmentation and deprivation. Obesity and cardiovascular disorders are the major risk factors associated with OSA. Molecular analysis of the factors associated with OSA could demarcate the clinical analysis pattern in a population.ObjectiveThis study pertains to in-silico analyses of miRNA and their gene targets with validation for their potential role in OSA as putative biomarker candidates.MethodsmiRDB, TargetScan and miRanda databases were used to identify targets of miR-27 and let-7 that have documented role in OSA and co-related obesity and cardiovascular disorders. Quantitative PCR was used to analyze expression pattern of miR-27 and let-7 in obese and non-obese OSA patient cohorts with respective controls. In-silico analysis was done using PatchDoc to obtain atomic contact energy (ACE) scores that indicated the docked gene targets to the predicted miRNA structures. The docked structures were analysed using Maestro Suite 11 for the hydrogen and aromatic interactions.ResultsDownregulation of miR-27 and let-7 in OSA compared to controls was observed. In-silico data analysis was performed for gene targets (TGFBR1, TGFBR2, SMAD2, SMAD4, CRY2 and CNR1) of the selected miRNAs (miR-27 and let-7). Among all, CNR1 and CRY2 were found to be better targets for miR-27 and let-7 respectively as per ACE scores, ROC scores and expression fold change in OSA.ConclusionOur study gives insights to the expression profiling of miR-27 and let-7 and explore a set of potential target genes (CNR1 and CRY2) of these two miRNAs for a promising clinical relevance in OSA.     

In silico prediction of human genes as potential targets for rice miRNAs

BackgroundExogenous microRNAs (miRNAs) enter the human body through food, and their effects on metabolic processes can be considerable. It is important to determine which miRNAs from plants affect the expression of human genes and the extent of their influence.MethodThe binding sites of 738Oryza sativa miRNAs (osa-miRNAs) that interact with 17 508 mRNAs of human genes were determined using the MirTarget program.ResultThe characteristics of the binding of 46 single osa-miRNAs to 86 mRNAs of human genes with a value of free energy (ΔG) interaction equal 94%–100% from maximum ΔG were established. The findings showed that osa-miR2102-5p, osa-miR5075-3p, osa-miR2097-5p, osa-miR2919 targeted the largest number of genes at 38, 36, 23, 19 sites, respectively. mRNAs of 86 human genes were identified as targets for 93 osa-miRNAs of all family osa-miRNAs with ΔG values equal 94%–98% from maximum ΔG. Each miRNA of the osa-miR156-5p, osa-miR164-5p, osa-miR168-5p, osa-miR395-3p, osa-miR396-3p, osa-miR396-5p, osa-miR444-3p, osa-miR529-3p, osa-miR1846-3p, osa-miR2907-3p families had binding sites in mRNAs of several human target genes. The binding sites of osa-miRNAs in mRNAs of the target genes for each family of osa-miRNAs were conserved when compared to flanking nucleotide sequences.ConclusionTarget mRNA human genes of osa-miRNAs are also candidate genes of cancer, cardiovascular and neurodegenerative diseases.     

Brief survey on phytochemicals to prevent COVID-19

BackgroundThe recent pandemic by COVID-19 is a global threat to human health. The disease is caused by SARS-CoV-2 and the infection rate is increased more quickly than MERS and SARS as their rapid adaptation to varied climatic conditions through rapid mutations. It becomes more severe due to the lack of proper therapeutic drugs, insufficient diagnostic tool, scarcity of appropriate drug, life supporting medical facility and mostly lack of awareness. Therefore, preventive measure is one of the important strategies to control. In this context, herbal medicinal plants received a noticeable attention to treat COVID-19 in Indian subcontinent. Here, 44 Indian traditional plants have been discussed with their novel phytochemicals that prevent the novel corona virus. The basic of SARS-CoV-2, their common way of transmission including their effect on immune and nervous system have been discussed. We have analysed their mechanism of action against COVID-19 following in-silico analysis. Their probable mechanism and therapeutic approaches behind the activity of phytochemicals to stimulate immune response as well as inhibition of viral multiplication discussed rationally. Thus, mixtures of active secondary metabolites/phytochemicals are the only choice to prevent the disease in countries where vaccination will take long time due to overcrowded population density.     

Targeting non-structural proteins and 3CLpro in SARS-CoV-2 virus using phytochemicals from medicinal plants - In-silico approach

In the present study, the main protease 3CL^pro and non-structural protein (NSP-12 with co-factors 7 and 8) trimer complex are used to study the protein-drug interactions with the phytochemicals from Ocimum Sanctum, Tinospora Cordifolia, Glycyrrhiza Glabra, and Azadirachta Indica. Which can give insight to be used as potent antiviral drugs against SARS-CoV-2. Twenty phytochemicals, five from each plant species, known for their wide range of biological activities were chosen from the literature. The in-silico study was carried out using virtual screening tools and the top five, which showed the least binding energies, were selected. Molecular docking tools revealed that gedunin and epoxy azadiradione proved to be excellent inhibitors for 3CL^pro and so did Tinosporide for non-structural-protein complex. Further, the best-hit phytochemicals with respect to structure similarities with FDA drugs and investigatory drugs, were considered for comparative study. Molecular docking was done to check the drug-protein interactions and to check the inhibitory responses of these drugs against the viral protein. The analyses showed that the phytochemicals had similar responses on the protein complex but with exceptionally higher inhibitory responses hence which may be taken for further clinical study.     

Comparative and evolutionary studies of mammalian arylsulfatase and sterylsulfatase genes and proteins encoded on the X-chromosome

At least 19 sulfatase genes have been reported on the human genome, including four arylsulfatase (ARS) genes (ARSD; ARSE; ARSF; ARSH) and a sterylsulfatase (STS) gene located together on the X-chromosome. Bioinformatic analyses of mammalian genomes were undertaken using known human STS and ARS amino acid sequences to study the evolution of these genes and proteins encoded on eutherian and marsupial genomes. Several domain regions and key residues were conserved including signal peptides, active site residues, metal (Ca²⁺) and substrate binding sequences, transmembranes and N-glycosylation sites. Phylogenetic analyses describe the relationships and potential origins of these genes during mammalian evolution. Primate ARSH enzymes lacked signal peptide sequences which may influence their biological functions. CpG117 and CpG92 were detected within the 5′ region of the human STS and ARSD genes, respectively, and miR-205 within the 3′-UTR for the human STS gene, using bioinformatic methods A proposal is described for a primordial invertebrate STS-like gene serving as an ancestor for unequal cross over events generating the gene complex on the eutherian mammalian X-chromosome.     

Molecular characterization of pigeon torque teno virus (PTTV) in Jiangsu province

The torque teno virus (TTV) is a recently discovered DNA virus that has been detected in many different hosts, including humans, livestock and poultry. To date, there is no report of pigeon TTV (PTTV) from anywhere in the world. To investigate the distribution of PTTV in pigeons from the eastern Chinese province of Jiangsu and characterize their genomes, we employed PCR to detect PTTV in 144 samples collected from 6 pigeon plants in Jiangsu province, amplify complete genomes from representative samples and analyze genetic characteristics using bioinformatics. The results demonstrated that 71.5% (103/144) of samples were PTTV positive. The rate of sequence homology among the six PTTV complete genomes obtained from Jiangsu province ranged from 99.7% to 100%. Phylogenetic analysis suggested that PTTV genomes had a high degree of genetic similarity and were similar to chicken anemia virus that also had poultry as a host. Although with the same host, PTTV shared distant relationship with PiCV in both complete genome, Rep and Cap genes. The results of this study provided evidence that PTTV could be detected in Chinese pigeons at a high level, the evolutionary process of complete genome, Rep and Cap genes of Anelloviridae family had obvious divergence.     

Comparative study of synonymous codon usage variations between the nucleocapsid and spike genes of coronavirus, and C-type lectin domain genes of human and mouse

Coronaviruses (CoVs) are single-stranded RNA viruses which contain the largest RNA genomes, and severe acute respiratory syndrome coronavirus (SARS-CoV), a newly found group 2 CoV, emerged as infectious disease with high mortality rate. In this study, we compared the synonymous codon usage patterns between the nucleocapsid and spike genes of CoVs, and C-type lectin domain (CTLD) genes of human and mouse on the codon basis. Findings indicate that the nucleocapsid genes of CoVs were affected from the synonymous codon usage bias than spike genes, and the CTLDs of human and mouse partially overlapped with the nucleocapsid genes of CoVs. In addition, we observed that CTLDs which showed the similar relative synonymous codon usage (RSCU) patterns with CoVs were commonly derived from the human chromosome 12, and mouse chromosome 6 and 12, suggesting that there might be a specific genomic region or chromosomes which show a more similar synonymous codon usage pattern with viral genes. Our findings contribute to developing the codon-optimization method in DNA vaccines, and further study is needed to determine a specific correlation between the codon usage patterns and the chromosomal locations in higher organisms.     

WNT1 Gene from WNT Signaling Pathway Is a Direct Target of miR-122 in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is an invasive form of hepatic cancer arising from the accumulation of multiple genetic alterations. In this study, the causal role of disturbed canonical Wnt/β-catenin pathway was approved, and some of HCC-driven important gene candidates were determined. MicroRNAs (miRNAs), small non-coding RNAs, are the key regulators of important cancer genes, and their participation in tumorigenesis has been shown. By reviewing literature, WNT1 gene with functional significance was selected to approve miRNAs as new subjects for targeted therapy.For proper and fast miRNA detection and also confirmation of the role of bioinformatics in obtaining practical data, we benefited from different bioinformatics tools such as TargetScan, miRanda, and DIANA. In order to use an HCC model, we used HepG2 cell line. Luciferase assay was applied to assess the ability of the selected miRNAs in targeting WNT1 3′-UTR. To overexpress the selected miRNA in HepG2 cell line, viral construct was prepared. Quantitative real-time PCR was performed to evaluate selected miRNA and target gene expression levels. miR-122 was selected according to data concerning various bioinformatics tools.miR-122 was downregulated and WNT1 gene expression was upregulated in HepG2 cell line. After viral construct transduction, miR-122 expression was elevated and WNT1 expression was notably declined. Finally, we introduced WNT1 gene as one of the important genes in HCC, and also, we showed that miR-122 can regulate WNT1 gene expression.Moreover, our study determines the potential of bioinformatics analyses in providing accurate and reliable data for miRNA: messenger RNA (mRNA) prediction.     

Identification of novel human USP2 inhibitor and its putative role in treatment of COVID-19 by inhibiting SARS-CoV-2 papain-like (PLpro) protease

Human ubiquitin carboxyl-terminal hydrolase-2 (USP2) inhibitors, such as thiopurine analogs, have been reported to inhibit SARS-CoV papain-like proteases (PLpro). The PLpro have significant functional implications in the innate immune response during SARS-CoV-2 infection and considered an important antiviral target. Both proteases share strikingly similar USP fold with right-handed thumb–palm–fingers structural scaffold and conserved catalytic triad Cys-His-Asp/Asn. In this urgency situation of COVID-19 outbreak, there is a lack of in-vitro facilities readily available to test SARS-CoV-2 inhibitors in whole-cell assays. Therefore, we adopted an alternate route to identify potential USP2 inhibitor through integrated in-silico efforts. After an extensive virtual screening protocol, the best compounds were selected and tested. The compound Z93 showed significant IC₅₀ value against Jurkat (9.67 μM) and MOTL-4 cells (11.8 μM). The binding mode of Z93 was extensively analyzed through molecular docking, followed by MD simulations, and molecular interactions were compared with SARS-CoV-2. The relative binding poses of Z93 fitted well in the binding site of both proteases and showed consensus π-π stacking and H-bond interactions with histidine and aspartate/asparagine residues of the catalytic triad. These results led us to speculate that compound Z93 might be the first potential chemical lead against SARS-CoV-2 PLpro, which warrants in-vitro evaluations.     

Epstein–Barr virus latent genes

Latent Epstein–Barr virus (EBV) infection has a substantial role in causing many human disorders. The persistence of these viral genomes in all malignant cells, yet with the expression of limited latent genes, is consistent with the notion that EBV latent genes are important for malignant cell growth. While the EBV-encoded nuclear antigen-1 (EBNA-1) and latent membrane protein-2A (LMP-2A) are critical, the EBNA-leader proteins, EBNA-2, EBNA-3A, EBNA-3C and LMP-1, are individually essential for in vitro transformation of primary B cells to lymphoblastoid cell lines. EBV-encoded RNAs and EBNA-3Bs are dispensable. In this review, the roles of EBV latent genes are summarized.     

The regulation of microRNA in each of cancer stage from two different ethnicities as potential biomarker for breast cancer

miRNA has recently emerged as a potential biomarker for breast cancer. Even though many studies have identified ethnic variation affecting miRNA regulation, the effect of cancer stage within specific ethnicities on miRNA epigenetic remains unclear. The present study is designed to investigate miRNA regulation from two distinct ethnicities in specific cancer stages (non-Hispanic white and non-Hispanic black) using the TCGA dataset. Differentially expressed miRNAs were calculated by using the edgeR package. miRNAs with the highest or lowest log fold Change from each cancer stage were selected as a potential biomarker. miRNA-gene interaction was analyzed by using spearman correlation analysis, CLUEGO, and DIANA-mirpath. The association of biomarker candidates with diagnostic and prognostic performance was assessed using ROC and Kaplan-Meier survival analysis. miRNA-gene interaction analysis revealed the involvement of selected miRNAs in cancer progression. From eleven selected aberrant miRNAs, four of the miRNAs (hsa-mir-495, hsa-mir-592, hsa-mir-6501, and hsa-mir-937) are significantly detrimental to breast cancer diagnosis and prognosis. Hence, our result provides valuable information to explore miRNA’s role in each cancer stage between non-Hispanic white and non-Hispanic black.     

Multiplexed and amplified chemiluminescence resonance energy transfer (CRET) detection of genes and microRNAs using dye-loaded hemin/G-quadruplex-modified UiO-66 metal–organic framework nanoparticles

Dye-loaded UiO-66 metal–organic framework nanoparticles (NMOFs) modified with catalytic hemin/G-quadruplex DNAzyme labels act as functional hybrid modules for the chemiluminescence resonance energy transfer (CRET) analysis of miRNAs (miRNA-155 or miRNA-21) or genes (p53 or BRCA1). The dye-loaded NMOFs (dye = fluorescein (Fl) or rhodamine 6G (Rh 6G)) are modified with hairpin probes that are engineered to include in their loop domains recognition sequences for the miRNAs or genes, and in their stem regions caged G-quadruplex domains. In the presence of the analytes miRNAs or genes, the hairpin structures are opened, leading, in the presence of hemin, to the self-assembly of hemin/G-quadruplex DNAzyme labels linked to the dye-loaded NMOFs. In the presence of luminol and H₂O₂, the hemin/G-quadruplex DNAzyme labels catalyze the generation of chemiluminescence that provides radiative energy to stimulate the process of CRET to the dye loaded in the NMOFs, resulting in the luminescence of the loaded dye without external excitation. The resulting CRET signals relate to the concentrations of the miRNAs or the genes and allow the sensitive analysis of miRNAs and genes. In addition, the DNA hairpin-functionalized dye-loaded NMOF sensing modules were further applied to develop amplified miRNA or gene CRET-based sensing platforms. The dye-loaded NMOFs were modified with hairpin probes that include in their loop domain the recognition sequences for miRNA-155 or miRNA-21 or the recognition sequences for the p53 or BRCA1 genes. Subjecting the hairpin-modified NMOFs to the respective miRNAs or genes, in the presence of two hairpins H_i and H_j that include in their stem regions caged G-quadruplex subunit domains, results in the analyte-triggered opening of the probe hairpin linked to the NMOFs, and the opened hairpin tethers induce the cross-opening of the hairpins H_i and H_j by the hybridization chain reaction, HCR, resulting in the assembly of G-quadruplex wires tethered to the NMOFs. The binding of hemin to the HCR-generated chains yields hemin/G-quadruplex DNAzyme wires that enhance, in the presence of luminol/H₂O₂, the CRET processes in the hybrid nanostructures. These amplification platforms lead to the amplified sensing of miRNAs and genes. By mixing the Fl- and Rh 6G-loaded hairpin-functionalized UiO NMOFs, the multiplexed CRET detection of miRNA-155, miRNA-21 and the p53 and BRCA1 genes is demonstrated.

Hemin/G-quadruplex DNAzyme-modified metal–organic framework nanoparticles act as functional hybrids for the catalyzed oxidation of luminol by H₂O₂, causing chemiluminescence and activation of chemiluminescence resonance energy transfer to the dye loads.     

Sequence-based analysis of 5′UTR and coding regions of CASP3 in terms of miRSNPs and SNPs in targetting miRNAs

Apoptosis is described as a mechanism of cell death occurring after adequate cellular harm. Deregulation of apoptosis occurs in many human conditions such as autoimmune disorders, ischemic damage, neurodegenerative diseases and different cancer types. Information relating miRNAs to cancer is increasing. miRNAs can affect development of cancer via many different pathways, including apoptosis. Polymorphisms in miRNA genes or miRNA target sites (miRSNPs) can change miRNA activity. Although polymorphisms in miRNA genes are very uncommon, SNPs in miRNA-binding sites of target genes are quite common. Many researches have revealed that SNPs in miRNA target sites improve or decrease the efficacy of the interaction between miRNAs and their target genes. Our aim was to specify miRSNPs on CASP3 gene (caspase-3) and SNPs in miRNA genes targeting 5′UTR and coding exons of CASP3, and evaluate the effect of these miRSNPs and SNPs of miRNA genes with respect to apoptosis. We detected 141 different miRNA binding sites (126 different miRNAs) and 7 different SNPs in binding sites of miRNA in 5′UTR and CDS of CASP3 gene. Intriguingly, miR-339-3p’s binding site on CASP3 has a SNP (rs35372903, G/A) on CASP3 5′UTR and its genomic sequence has a SNP (rs565188493, G/A) at the same nucleotide with rs35372903. Also, miR-339-3p has two other SNPs (rs373011663, C/T rs72631820, A/G) of which the first is positioned at the binding site. Here, miRSNP (rs35372903) at CASP3 5′UTR and SNP (rs565188493) at miR-339-3p genomic sequence cross-matches at the same site of binding region. Besides, miR-339-3p targets many apoptosis related genes (ZNF346, TAOK2, PIM2, HIP1, BBC3, TNFRSF25, CLCF1, IHPK2, NOL3) although it had no apoptosis related interaction proven before. This means that miR-339-3p may also have a critical effect on apoptosis via different pathways other than caspase-3. Hence, we can deduce that this is the first study demonstrating a powerful association between miR-339-3p and apoptosis upon computational analysis.     

Application of an integrated cheminformatics-molecular docking approach for discovery for physicochemically similar analogs of fluoroquinolones as putative HCV inhibitors

BackgroundHepatitis C Virus (HCV) infection is a major public health concern across the globe. At present, direct-acting antivirals are the treatment of choice. However, the long-term effect of this therapy has yet to be ascertained. Previously, fluoroquinolones have been reported to inhibit HCV replication by targeting NS3 protein. Therefore, it is logical to hypothesize that the natural analogs of fluoroquinolones will exhibit NS3 inhibitory activity with substantially lesser side effects.MethodIn this study, we tested the application of a recently devised integrated in-silico Cheminformatics-Molecular Docking approach to identify physicochemically similar natural analogs of fluoroquinolones from the available databases (Ambinter, Analyticon, Indofines, Specs, and TimTec). Molecular docking and ROC curve analyses were performed, using PatchDock and Graphpad software, respectively, to compare and analyze drug-protein interactions between active natural analogs, Fluoroquinolones, and HCV NS3 protein.ResultIn our analysis, we were able to shortlist 18 active natural analogs, out of 10,399, that shared physicochemical properties with the template drugs (fluoroquinolones). These analogs showed comparable binding efficacy with fluoroquinolones in targeting 32 amino acids in the HCV NS3 active site that are crucial for NS3 activity. Our approach had around 80 % sensitivity and 70 % specificity in identifying physicochemically similar analogs of fluoroquinolones.ConclusionOur current data suggest that our approach can be efficiently applied to identify putative HCV drug inhibitors that can be taken for in vitro testing. This approach can be applied to discover physicochemically similar analogs of virtually any drug, thus providing a speedy and inexpensive approach to complement drug discovery and design, which can tremendously economize on time and money spent on the screening of putative drugs.     

Cross-Kingdom Gene regulation via miRNAs of Hypericum perforatum (St. John’s wort) flower dietetically absorbed: An in silico approach to define potential biomarkers for prostate cancer

Prostate cancer (PCa) is the most frequent type of cancer in men. Hypericum perforatum (H. Perforatum) extract (HPE) administration provides remarkable decrease of PCa development. H. perforatum contains 7 conserved miRNAs (Hyp-miR-156a, Hyp-miR-156b, Hyp-miR-166, Hyp-miR-390, Hyp-miR-394, Hyp-miR-396 and Hyp-miR-414) with different targets. In this study, we aimed to investigate cross-kingdom gene regulation via miRNAs of H. perforatum flower dietetically absorbed in manner of an in silico approach to define potential biomarkers for PCa. psRNATarget database was used to find human genes targeted by 7 pre-defined H. perforatum miRNAs. We defined the mostly affected gene families from these miRNAs as ZNF, TMEM, SLC and FAM gene families. GeneMANIA database was used to define the most affected genes (TMEM41B and SLC4A7) from these 7 miRNAs. cBioPortal database was used to define alteration frequencies of TMEM41B and SLC4A7 on different types of PCa and to measure the mutual interaction potency and significance of co-occurence in PCa. This analysis showed that neuroendocrine prostate cancer (NEPC) had the highest total mutation frequency (22%) of TMEM41B and SLC4A7 genes. Also, TMEM41B and SLC4A7 genes had an average 2.1% pathway change potential among all different types of PCa. Moreover, TMEM41B and SLC4A7 gene pair was found significantly co-occurrent in PCa (p < 0.001). Finally, via GEPIA database, we used Spearman correlation analysis to measure the correlation degree of TMEM41B and SLC4A7 genes in PCa and found their significant correlation with PCa (p = 1.2 × 10⁻¹², R = 0.28). All in all, it was proved in silico and supported with previously known clinical data that SLC4A7 and TMEM41B potentially have a significant and critical tumor suppressive role for PCa, and show this effect combinatorily working together. This is the first study correlating SLC4A7 and TMEM41B with PCa significantly.     

In Silico Mining of Conserved miRNAs of Indian Catfish <Emphasis Type="Italic">Clarias batrachus</Emphasis> (Linnaeus, 1758) from the Contigs,ESTs, and BAC End Sequences

MicroRNAs are small non-coding RNAs which play significant role in RNA interference. The present work deals with mining of the conserved miRNA and their target genes from the contigs, ESTs, and BAC end sequences of commercially important catfish, Clarias batrachus, from India. A total of 138, 1 and 1 conserved pre-miRNA sequences, were mined from the contigs, ESTs, and BAC end sequences, respectively. The analysis of families of the conserved pre-miRNA revealed conservation of the fish-specific family mir-430 and other important families, such as mir-455, let-7, mir-133, and mir-137. The mir-455 is involved in hypoxia signaling, let-7 family represents potential anti-tumor molecules involved in human cancer therapy, whereas mir-133 and mir-137 have high therapeutic potentials. Using an alternate computational in silico approach, mining of mature miRNAs resulted in identification of 210 mature miRNAs from contigs, 1 from EST, and 2 each from forward as well as reverse BAC end sequences. Target prediction of these putative miRNAs resulted in the identification of 66,758 and 18,747 target genes in C. batrachus and Danio rerio, respectively. Functional annotation of these miRNAs indicated their involvement in diverse biological functions. The findings of the present study can serve as a valuable resource for further functional genomics studies in C. batrachus.