The Zika virus (ZIKV), a member of the Flaviviridae family, poses a major threat to human health because of the lack of effective antiviral drugs. Although the NS2B-NS3 protease of ZIKV (NS2B-NS3pro) is regarded as a major target for antiviral inhibitors, viral mutations can lead to ineffective competitive inhibitors. Allosteric inhibitors bind to highly conserved nonprotease active sites, induce conformational changes in the protease active site, and prevent substrate binding. Currently, no molecular simulation techniques are available for accurately predicting and analysing conformational changes in the protease catalytic domain. In this study, we developed a combined approach that involves blind docking, Gaussian accelerated molecular dynamics, two-dimensional potential of mean force profiling, density functional theory (DFT) calculations, and interaction region indicator (IRI) analysis and employed it to examine the allosteric inhibitor-01 molecule and its interaction with ZIKV NS2B-NS3pro. Our results indicated that the binding of inhibitor-01 to NS2B-NS3pro resulted in two major conformational states, state I and state II, which in turn changed the volume of the protease active site from 1014 Å