TMC-126

The binding properties of the protein-inhibitor complex of human immunodeficiency virus type 1 (HIV-1) protease with the inhibitor TMC-126 are investigated by combining computational alanine scanning (CAS) mutagenesis with binding free-energy decomposition (BFED). The calculated results demonstrate that the flap region (residues 38-58) and the active site region (residues 23-32) in HIV-1 protease contribute 63.72% of the protease to the binding of the inhibitor. In particular, the mechanisms for the interactions of key residues of these species are fully explored and analyzed. Interestingly, the regression analyses show that both CAS and BFED based on the generalized Born model yield similar results, with a correlation coefficient of 0.94. However, compared to CAS, BFED is faster and can decompose the per-residue binding free-energy contributions into backbone and side-chain contributions. The results obtained in this study are useful for studying the binding mechanism between receptor and ligand and for designing potent inhibitors that can combat diseases.

The V82F/I84V double mutation is considered as the key residue mutation of the HIV-1 protease drug resistance because it can significantly lower the binding affinity of protease inhibitors in clinical uses. In the current work, the binding of amprenavir to both of the wild-type and the drug-resistant V82F/I84V mutant of the HIV-1 protease was investigated by molecular dynamics (MD) simulations and was compared to those of two inhibitors in development, TMC126 and TMC114. Absolute binding free energies were calculated by molecular mechanics/Poisson-Boltzmann surface area (MM/ PBSA) methodology. The predicted binding affinities give a good explanation of structure-affinity relationship (SAR) of three studied inhibitors. Furthermore, in the 18 ns MD simulations on the free wild-type and the mutated proteases, we observed that the free mutated protease shows similar dynamic characteristics of the flap opening and a little higher structural stability than the free wild-type protease. This suggests that the effect of the mutations may not significantly affect the equilibrium between the semiopen and the closed conformations. Finally, decomposition analysis of binding free energies and the further structural analysis indicate that the dominating effect of the V82F/I84V double mutation is to distort the geometry of the binding site and hence weaken the interactions of inhibitors preshaped to the wild-type binding site.