Network-base approaches to identify therapeutic biomarkers in hepatocellular carcinoma and search for drug hunting utilizing molecular dynamics simulations

Abstract

The presence of conditions like Alpha-1 antitrypsin deficiency, hemochromatosis, non-alcoholic fatty liver diseases and metabolic syndrome can elevate the susceptibility to hepatic cellular carcinoma (HCC). Utilizing network-based gene expression profiling via network analyst tools, presents a novel approach for drug target discovery. The significance level (p-score) obtained through Cytoscape in the intended center gene survival assessment confirms the identification of all target center genes, which play a fundamental role in disease formation and progression in HCC. A total of 1064 deferential expression genes were found. These include MCM2 with the highest degree, followed by 4917 MCM6 and MCM4 with a 3944-degree score. We investigated the regulatory kinases involved in establishing the protein-protein interactions network using X2K web tool. The docking approach yields a favorable binding affinity of −8.7 kcal/mol against the target MCM2 using Auto-Dock Vina. Interestingly after simulating the complex system via AMBER16 package, results showed that the root mean square deviation values remained within 4.74 Å for a protein and remains stable throughout the time intervals. Additionally, the ligand’s fit to the protein exhibited fluctuations at some intervals but remains stable. Finally, Gibbs free energy was found to be at its lowest at 1 kcal/mol which presents the real time interactive binding of the atomic residues among inhibitor and protein. The displacement of the ligand was measured showing stable movement and displacement along the active site. These findings increased our understanding for potential biomarkers in hepatocellular carcinoma and an experimental approach will further enhance our outcomes in future.