HIV-1 has a large genetic diversity. Subtype B HIV-1 is commonly found in patients in developed countries. In contrast, an increasing number of patients are infected with the non-B subtype viruses, especially with subtype C HIV-1, in developing countries. It remains to be clarified how mutations or polymorphisms in non-B subtype HIV-1 influence the efficacy of the approved inhibitors. In this study, we have performed molecular dynamics simulations on clinically isolated subtype C HIV-1 proteases in complex with three kinds of approved inhibitors. From the structural and energetic viewpoints, we identified the polymorphisms influencing on the binding of the inhibitors. The effect of the V82I mutation on the association with chemicals and the reason for rare appearance of the D30N mutation in subtype C HIV-1 were discussed in terms of the change of geometry of the residues in HIV-1 protease.