Computational analysis and molecular dynamics insights into deleterious SNPs of the gene.

Mutations in the gene, especially non-synonymous single-nucleotide polymorphisms (nsSNPs), are strongly associated with hemochromatosis, an autosomal recessive disorder characterized by intracellular iron overload, a key feature of tumor development. This study examined the structural and functional effects of deleterious nsSNPs in the gene using bioinformatics tools, gene interaction analyses, molecular docking, molecular dynamics (MD) simulations, and assessments of clinical relevance. Functional analyses identified nine deleterious nsSNPs, including C282Y, L183P, and Q283P, which disrupted disulfide bonds, hydrogen bonds, and hydrophobic interactions, destabilizing the protein. Conservation analysis revealed these mutations occur in highly conserved regions, emphasizing their structural and functional importance. Notably, five nsSNPs (R224Q, R224W, I235T, C282Y, Q283P) within the Ig-like C1-type domain were associated with cancer. Gene interaction analyses showed -related genes are linked to immunity and iron balance. Variants in interacting genes, such as and , may influence iron disorders, infection risk, and inflammation. Molecular docking showed reduced interface interactions for the C282Y mutant and altered binding to transferrin receptor 1 (TfR1), potentially destabilizing the HFE-TfR1 complex. MD simulations highlighted key differences, with the mutant showing higher RMSD, decreased compactness (Rg), increased flexibility (RMSF), and greater solvent exposure (SASA), confirming destabilization. Furthermore, expression varied across cancers, with elevated levels in twelve tumor types. Higher expression correlated with better survival in breast and gastric cancers but poorer outcomes in lung cancer. These findings highlight how deleterious nsSNPs, especially C282Y, disrupt structure and function, offering insights into disease mechanisms and guiding therapeutic strategies.