Dataset Information

Interaction of I50V mutant and I50L/A71V double mutant HIV-protease with inhibitor TMC114 (darunavir): molecular dynamics simulation and binding free energy studies.

ABSTRACT: In the present work, the binding of inhibitor TMC114 (darunavir) to wild-type (WT), single (I50V) as well as double (I50L/A71V) mutant HIV-proteases (HIV-pr) was investigated with all-atom molecular dynamics (MD) simulations as well as molecular mechanic-Poisson-Boltzmann surface area (MM-PBSA) calculation. For both the apo and complexed HIV-pr, many intriguing effects due to double mutant, I50L/A71V, are observed. For example, the flap-flap distance and the distance from the active site to the flap residues in the apo I50L/A71V-HIV-pr are smaller than those of WT- and I50V-HIV-pr, probably making the active site smaller in volume and closer movement of flaps. For the complexed HIV-pr with TMC114, the double mutant I50L/A71V shows a less curling of the flap tips and less flexibility than WT and the single mutant I50V. As for the other previous studies, the present results also show that the single mutant I50V decreases the binding affinity of I50V-HIV-pr to TMC, resulting in a drug resistance; whereas the double mutant I50L/A71V increases the binding affinity, and as a result of the stronger binding, the I50L/A71V may be well adapted by the TMC114. The energy decomposition analysis suggests that the increase of the binding for the double mutant I50L/A71V-HIV-pr can be mainly attributed to the increase in electrostatic energy by -5.52 kacl/mol and van der Waals by -0.42 kcal/mol, which are canceled out in part by the increase of polar solvation energy of 1.99 kcal/mol. The I50L/A71V mutant directly increases the binding affinity by approximately -0.88 (Ile50 to Leu50) and -0.90 (Ile50' to Leu50') kcal/mol, accounting 45% for the total gain of the binding affinity. Besides the direct effects from the residues Leu50 and Leu50', the residue Gly49' increases the binding affinity of I50L/A71V-HIV-pr to the inhibitor by -0.74 kcal/mol, to which the electrostatic interaction of Leu50's backbone contributes by -1.23 kcal/mol. Another two residues Ile84 and Ile47' also increase the binding affinity by -0.22 and -0.29 kcal/mol, respectively, which can be mainly attributed to van der Waals terms (?T(vdw) = -0.21 and -0.39 kcal/mol).

SUBMITTER: Meher BR

PROVIDER: S-EPMC3288396 | biostudies-literature | 2012 Feb

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Interaction of I50V mutant and I50L/A71V double mutant HIV-protease with inhibitor TMC114 (darunavir): molecular dynamics simulation and binding free energy studies.

Meher Biswa Ranjan BR Wang Yixuan Y

The journal of physical chemistry. B 20120203 6

In the present work, the binding of inhibitor TMC114 (darunavir) to wild-type (WT), single (I50V) as well as double (I50L/A71V) mutant HIV-proteases (HIV-pr) was investigated with all-atom molecular dynamics (MD) simulations as well as molecular mechanic-Poisson-Boltzmann surface area (MM-PBSA) calculation. For both the apo and complexed HIV-pr, many intriguing effects due to double mutant, I50L/A71V, are observed. For example, the flap-flap distance and the distance from the active site to the ...[more]

PMID: 22239286

Similar Datasets

Project description:We report the binding free energy calculation and its decomposition for the complexes of alpha-lytic protease and its protein inhibitors using molecular dynamics simulation. Standard mechanism serine protease inhibitors eglin C and OMTKY3 are known to have strong binding affinity for many serine proteases. Their binding loops have significant similarities, including a common P1 Leu as the main anchor in the binding interface. However, recent experiments demonstrate that the two inhibitors have vastly different affinity towards alpha-lytic protease (ALP), a bacterial serine protease. OMTKY3 inhibits the enzyme much more weakly (by approximately 10(6) times) than eglin C. Moreover, a variant of OMTKY3 with five mutations, OMTKY3M, has been shown to inhibit 10(4) times more strongly than the wild-type inhibitor. The underlying mechanisms for the unusually large difference in binding affinities and the effect of mutation are not well understood. Here we use molecular dynamics simulation with molecular mechanics-Poisson Boltzmann/surface area method (MM-PB/SA) to investigate quantitatively the binding specificity. The calculated absolute binding free energies correctly differentiate the thermodynamic stabilities of these protein complexes, but the magnitudes of the binding affinities are systematically overestimated. Analysis of the binding free energy components provides insights into the molecular mechanism of binding specificity. The large DeltaDeltaG(bind) between eglin C and wild type OMTKY3 towards ALP is mainly attributable to the stronger nonpolar interactions in the ALP-eglin C complex, arising from a higher degree of structural complementarity. Here the electrostatic interaction contributes to a lesser extent. The enhanced inhibition in the penta-mutant OMTKY3M over its wild type is entirely due to an overall improvement in the solvent-mediated electrostatic interactions in the ALP-OMTKY3M complex. The results suggest that for these protein-complexes and similar enzyme-inhibitor systems (1) the binding is driven by nonpolar interactions, opposed by overall electrostatic and solute entropy contributions; (2) binding specificity can be tuned by improving the complementarity in electrostatics between two associating proteins. Binding free energy decomposition into contributions from individual protein residues provides additional detailed information on the structural determinants and subtle conformational changes responsible for the binding specificity.

Dataset Information

Interaction of I50V mutant and I50L/A71V double mutant HIV-protease with inhibitor TMC114 (darunavir): molecular dynamics simulation and binding free energy studies.

Publications

Interaction of I50V mutant and I50L/A71V double mutant HIV-protease with inhibitor TMC114 (darunavir): molecular dynamics simulation and binding free energy studies.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets