Sorafenib is the only approved targeted drug for hepatocellular carcinoma (HCC), but its effect on patients survival gain is limited and varies over a wide range depending on patho-genetic conditions. Thus, enhancing the efficacy of sorafenib and finding a reliable predictive biomarker are crucuial to achieve efficient control of HCCs. In this study, we employed a systems approach by combining transcriptome analysis of the mRNA changes in HCC cell lines in response to sorafenib with network analysis to investigate the action and resistance mechanism of sorafenib. Gene ontology and gene set analysis revealed that proteotoxic stress and apoptosis modules are activated in the presence of sorafenib. Further analysis of the endoplasmic reticulum (ER) stress network model combined with in vitro experiments showed that introducing an additional stress by treating the orally active protein disulfide isomerase (PDI) inhibitor (PACMA 31) can synergistically increase the efficacy of sorafenib in vitro and in vivo, which was confirmed using a mouse xenograft model. We also found that HCC patients with high PDI expression show resistance to sorafenib and poor clinical outcomes, compared to the low PDI expression group. These results suggest that PDI is a promising therapeutic target for enhancing the efficacy of sorafenib and can also be a biomarker for predicting sorafenib responsiveness.