Screening of miRNAs as prognostic biomarkers and their associated hub targets across Hepatocellular carcinoma using survival-based bioinformatics approach.

BACKGROUND: The hepatocellular carcinoma (HCC) incident rate is gradually increasing yearly despite all the research and efforts taken by scientific communities and governing bodies. Approximately 90% of all liver cancer cases belong to HCC. Usually, HCC patients approach the treatment in the late stages of this malignancy which becomes the primary cause of high mortality rate. The knowledge about molecular pathogenesis of HCC is limited and needs more attention from researchers to identify the driver genes and miRNAs, which causes to translate this information into clinical practice. Therefore, the key regulators identification of miRNA-mRNA regulatory network is essential to identify HCC-associated genes. METHODOLOGY: We extracted microRNA (miRNA) and messenger RNA (mRNA) expression datasets of normal and tumor HCC patient samples from UCSC Xena followed by identifying differentially expressed genes (DEGs) and differentially expressed miRNAs (DEMs). Univariate and multivariate cox-proportional hazard models were utilized to identify DEMs having significant association with overall survival (OS). Kaplan-Meier (KM) plotter was used to validate the presence of prognostic DEMs. A risk-score model was used to evaluate the effectiveness of KM-plotter validated DEMs combination on risk of samples. Target DEGs of prognostic miRNAs were identified via sources such as miRTargetLink and miRWalk followed by their validation in an external microarray cohort and enrichment analysis. RESULTS: 562 DEGs and 388 DEMs were identified followed by seven prognostic miRNAs (i.e., miR-19a, miR-19b, miR-30d-5p, miR-424-5p, miR-3677-5p, miR-3913-5p, miR-7705) post univariate, multivariate, risk-score model evaluation and KM-plotter analyses. ANLN, MRO, CPEB3 were their targets and were also validated in GSE84005 dataset. CONCLUSIONS: The findings of this study decipher that most significant miRNAs and their identified target genes have association with apoptosis, inflammation, cell cycle regulation and cancer-related pathways, which appear to contribute to HCC pathogenesis and therefore, the discovery of new targets.

Identification of potential therapeutic targets for COVID-19 through a structural-based similarity approach between SARS-CoV-2 and its human host proteins.

Background: The COVID-19 pandemic caused by SARS-CoV-2 has led to millions of deaths worldwide, and vaccination efficacy has been decreasing with each lineage, necessitating the need for alternative antiviral therapies. Predicting host-virus protein-protein interactions (HV-PPIs) is essential for identifying potential host-targeting drug targets against SARS-CoV-2 infection. Objective: This study aims to identify therapeutic target proteins in humans that could act as virus-host-targeting drug targets against SARS-CoV-2 and study their interaction against antiviral inhibitors. Methods: A structure-based similarity approach was used to predict human proteins similar to SARS-CoV-2 ("hCoV-2"), followed by identifying PPIs between hCoV-2 and its target human proteins. Overlapping genes were identified between the protein-coding genes of the target and COVID-19-infected patient's mRNA expression data. Pathway and Gene Ontology (GO) term analyses, the construction of PPI networks, and the detection of hub gene modules were performed. Structure-based virtual screening with antiviral compounds was performed to identify potential hits against target gene-encoded protein. Results: This study predicted 19,051 unique target human proteins that interact with hCoV-2, and compared to the microarray dataset, 1,120 target and infected group differentially expressed genes (TIG-DEGs) were identified. The significant pathway and GO enrichment analyses revealed the involvement of these genes in several biological processes and molecular functions. PPI network analysis identified a significant hub gene with maximum neighboring partners. Virtual screening analysis identified three potential antiviral compounds against the target gene-encoded protein. Conclusion: This study provides potential targets for host-targeting drug development against SARS-CoV-2 infection, and further experimental validation of the target protein is required for pharmaceutical intervention.

Structural and thermodynamic properties of conserved water molecules in Mpro native: A combined approach by MD simulation and Grid Inhomogeneous Solvation Theory.

The new viral strains of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) are continuously rising, becoming more virulent, and transmissible. Therefore, the development of new antiviral drugs is essential. Due to its significant role in the viral life cycle of SARS-CoV-2, the main protease (Mpro) enzyme is a leading target for antiviral drug design. The Mpro monomer consists of domain DI, DII, and DI-DII interface. Twenty-one conserved water molecules (W4-W24) are occupied at these domains according to multiple crystal structure analyses. The crystal and MD structures reveal the presence of eight conserved water sites in domain DI, DII and remaining in the DI-DII interface. Grid-based inhomogeneous fluid solvation theory (GIST) was employed on MD structures of Mpro native to predict structural and thermodynamic properties of each conserved water site for focusing to identify the specific conserved water molecules that can easily be displaced by proposed ligands. Finally, MD water W13 is emerged as a promising candidate for water mimic drug design due to its low mean interaction energy, loose binding character with the protein, and its involvement in a water-mediated H-bond with catalytic His41 via the interaction Thr25(OG)---W13---W---His41(NE2). In this context, water occupancy, relative interaction energy, entropy, and topologies of W13 are thermodynamically acceptable for the water displacement method. Therefore, the strategic use of W13's geometrical position in the DI domain may be implemented for drug discovery against COVID disease by designing new ligands with appropriately oriented chemical groups to mimic its structural, electronic, and thermodynamic properties.

New biochemical insights into the dynamics of water molecules at the GMP or IMP binding site of human GMPR enzyme: A molecular dynamics study.

Human GMP reductase (hGMPR) enzyme is involved in a cellular metabolic pathway, converting GMP into IMP, and also it is an important target for anti-leukemic agents. Present computational investigations explain dynamical behavior of water molecules during the conformational transition process from GMP to IMP using molecular dynamics simulations. Residues at substrate-binding site of cancerous protein (PDB Id. 2C6Q) are mostly more dynamic in nature than the normal protein (PDB Id. 2BLE). Nineteen conserved water molecules are identified at the GMP/IMP binding site and are classified as (i) conserved stable dynamic and (ii) infrequent dynamic. Water molecules W11, W14, and W16 are classified as conserved stable dynamic due to their immobile character, whereas remaining water molecules (W1, W2, W3, W4, W5, W7, W8, W9, W10, W12, W13, W15, W17, W18, and W19) are infrequent with dynamic nature. Entrance or displacement of these infrequent water molecules at GMP/IMP sites may occur due to forward and backward movement of reference residues involving ligands. Four water molecules of hGMPR-I and nine water molecules of hGMPR-II are observed in repetitive transitions from GMP to IMP pathway, which indicates discrimination between two isoforms of hGMPRs. Water molecules in cancerous protein are more dynamic and unstable compared to normal protein. These water molecules execute rare dynamical events at GMP binding site and could assist in detailed understanding of conformational transitions that influence the hGMPR's biological functionality. The present study should be of interest to the experimental community engaged in leukemia research and drug discovery for CML cancer.

Integrative Multiomics and Regulatory Network Analyses Uncovers the Role of OAS3, TRAFD1, miR-222-3p, and miR-125b-5p in Hepatitis E Virus Infection.

The hepatitis E virus (HEV) is a long-ignored virus that has spread globally with time. It ranked 6th among the top risk-ranking viruses with high zoonotic spillover potential; thus, considering its viral threats is a pressing priority. The molecular pathophysiology of HEV infection or the underlying cause is limited. Therefore, we incorporated an unbiased, systematic methodology to get insights into the biological heterogeneity associated with the HEV. Our study fetched 93 and 2016 differentially expressed genes (DEGs) from chronic HEV (CHEV) infection in kidney-transplant patients, followed by hub module selection from a weighted gene co-expression network (WGCN). Most of the hub genes identified in this study were associated with interferon (IFN) signaling pathways. Amongst the genes induced by IFNs, the 2'-5'-oligoadenylate synthase 3 (OAS3) protein was upregulated. Protein-protein interaction (PPI) modular, functional enrichment, and feed-forward loop (FFL) analyses led to the identification of two key miRNAs, i.e., miR-222-3p and miR-125b-5p, which showed a strong association with the OAS3 gene and TRAF-type zinc finger domain containing 1 (TRAFD1) transcription factor (TF) based on essential centrality measures. Further experimental studies are required to substantiate the significance of these FFL-associated genes and miRNAs with their respective functions in CHEV. To our knowledge, it is the first time that miR-222-3p has been described as a reference miRNA for use in CHEV sample analyses. In conclusion, our study has enlightened a few budding targets of HEV, which might help us understand the cellular and molecular pathways dysregulated in HEV through various factors. Thus, providing a novel insight into its pathophysiology and progression dynamics.

Investigation of structural analogs of hydroxychloroquine for SARS-CoV-2 main protease (Mpro): A computational drug discovery study.

The main protease (Mpro) is the key enzyme of nCOVID-19 and plays a decisive role that makes it an attractive drug target. Multiple analysis of crystal structures reveals the presence of W1, W2, and W3 water locations in the active site pocket of Mpro; W1 and W2 are unstable and are weakly bonded with protein in comparison to W3 of Mpro-native. So, we adopt the water displacement method to occupy W1 or W2 sites by triggering HCQ or its analogs to inactivate the enzyme. Virtual screening is employed to find out best analogs of HCQ, molecular docking is used for water displacement from catalytic region of Mpro, and finally, MD simulations are conducted for validation of these findings. The docking study reveals that W1 and W2 are occupied by respective atoms of ZINC28706440 whereas W2 by HCQ and indacaterol. Finally, MD results demonstrate (i) HCQ occupies W1 and W2 positions, but its analogs (indacaterol and ZINC28706440) are inadequate to retain either W1 or W2 (ii) His41 and Asp187 are stabilized by W3 in Mpro-native and His41, Cys145 and HCQ by W7 in ZINC28706440, and W4, W5, and W6 make water mediated bridge between indacaterol with His41. The structural, dynamical, and thermodynamic (WFP and J value) profiling parameters suggest that W3, W4, and W7 are prominent in their corresponding positions in comparison with W5 and W6. The final results conclude that ZINC28706440 may act as a best analog of HCQ with acceptable physico-chemical and toxicological scores and may further be synthesized for experimental validation.

Structural and Dynamical Impact of Water Molecules at Substrate- or Product-Binding Sites in Human GMPR Enzyme: A Study by Molecular Dynamics Simulations.

Human guanosine monophosphate reductase (hGMPR) enzyme maintains the intracellular balance between adenine and guanine nucleotide pools, and it is an excellent target for the design of isoform-specific antileukemic agents. In the present study, we have investigated solvation properties of substrate GMP or product inosine-5'-monophosphate (IMP)-binding pocket of hGMPR by employing molecular dynamics simulations on conformations A (substrate GMP), B [substrate GMP with cofactor nicotinamide adenine dinucleotide phosphate (NDP)], C (product IMP with cofactor NDP), and D (product IMP). Nineteen water sites are identified precisely; they are responsible for the catalytic activity of this site, control structural and dynamical integrity, and electronic consequences of GMP or IMP in the binding site of hGMPR. The water sites of category-1 (W1, W4, W5, W6, W13, and W15) in normal protein and category-2 (W2, W3, W7, W8, W10, W17, and W18) in cancerous protein are unique and stabilize the guanosine or inosine group of GMP or IMP for participation in the enzymatic reaction, whereas the remaining water centers either stabilize pentose sugar ribose or the phosphate group of GMP or IMP. Furthermore, water sites of category-4 (W11, W14, and W16) appear to be conserved in all conformations during the entire simulation. The GMP-binding site in cancerous protein 2C6Q is significantly expanded, and its dynamics are very different from normal protein 2BLE. Furthermore, unique interactions of GMP(N1)···W2···Asp129/Asn158, IMP(N1)···W3···Glu289, and IMP(O6)···W10···Ser270 might be used in a water mimic drug design for hGMPR-II. In this context, water finding probability, relative interaction energy (J) associated with water site W, entropy, and topologies of these three water sites are thermodynamically acceptable for the water displacement method by the modified ligand. Hence, their positions in the catalytic pocket may also facilitate future drug discovery for chronic myelogenous leukemia by the design of appropriately oriented chemical groups that may displace these water molecules to mimic their structural, electronic, and thermodynamic properties.