HIV-1 protease (PR) is thought to be efficient targets of anti-AIDS drug design. Molecular dynamics (MD) simulations and multiple post-processing analysis technologies were applied to decipher molecular mechanism underlying binding of three drugs Lopinavir (LPV), Nelfinavir (NFV) and Atazanavir (ATV) to the PR. Binding free energies calculated by molecular mechanics generalized Born surface area (MM-GBSA) suggest that compensation between binding enthalpy and entropy plays a vital role in binding of drugs to PR. Dynamics analyses show that binding of LPV, NFV and ATV highly affects structural flexibility, motion modes and dynamics behaviour of the PR, especially for two flaps. Computational alanine scanning and interaction network analysis verify that although three drugs have structural difference, they share similar binding modes to the PR and common interaction clusters with the PR. The current findings also confirm that residues located interaction clusters, such as Asp25/Asp25MODIFIER LETTER PRIME, Gly27/Gly27MODIFIER LETTER PRIME, Ala28/Ala28MODIFIER LETTER PRIME, Asp29, Ile47/Ile47MODIFIER LETTER PRIME, Gly49/Gly49MODIFIER LETTER PRIME, Ile50/Ile50MODIFIER LETTER PRIME, Val82/Val82MODIFIER LETTER PRIME and Ile84/Ile84, can be used as efficient targets of clinically available inhibitors towards the PR.