The aspartyl dyad of free HIV-1 protease has apparent pK(a)s of approximately 3 and approximately 6, but recent NMR studies indicate that the aspartyl dyad is fixed in the doubly protonated form over a wide pH range when cyclic urea inhibitors are bound, and in the monoprotonated form when the inhibitor KNI-272 is bound. We present computations and measurements related to these changes in protonation and to the thermodynamic linkage between protonation and inhibition. The Poisson-Boltzmann model of electrostatics is used to compute the apparent pK(a)s of the aspartyl dyad in the free enzyme and in complexes with four different inhibitors. The calculations are done with two parameter sets. One assigns epsilon = 4 to the solute interior and uses a detailed model of ionization; the other uses epsilon = 20 for the solute interior and a simplified representation of ionization. For the free enzyme, both parameter sets agree well with previously measured apparent pK(a)s of approximately 3 and approximately 6. However, the calculations with an internal dielectric constant of 4 reproduce the large pKa shifts upon binding of inhibitors, but the calculations with an internal dielectric constant of 20 do not. This observation has implications for the accurate calculation of pK(a)s in complex protein environments. Because binding of a cyclic urea inhibitor shifts the pK(a)s of the aspartyl dyad, changing the pH is expected to change its apparent binding affinity. However, we find experimentally that the affinity is independent of pH from 5.5 to 7.0. Possible explanations for this discrepancy are discussed.

The structure and dynamics of water in the vicinity of the hemispherical hydrophobic pocket of 8 A radius were examined via molecular dynamics simulations in NVT ensemble. Density, hydrogen bonding properties, and residence times of water molecules were projected on two-dimensional planes providing a spatial description of water behavior. We found that the average water density is significantly depleted relative to bulk value. A detailed analysis of pocket occupancy revealed fluctuations between states of completely empty pocket and a pocket filled with a bulklike fluid, which seem to result from collective behavior of water molecules. Free energy differences accompanying these fluctuations are rather small, suggesting that the given pocket radius is close to the critical one for transition between gas and liquid phases in the considered system. We show that the situation is different in the case of a simple Lennard-Jones fluid. These results indicate that changing the surface curvature from flat to concave may lead to qualitative difference in water behavior in its vicinity. We think that our studies may also put some light on binding site desolvation process which is necessary to understand to make correct predictions of binding energies.

No abstract

The mechanism of the first steps of the reaction catalyzed by HIV-1 protease was studied through molecular dynamics simulations. The potential energy surface in the active site was generated using the approximate valence bond method. The approximate valence bond (AVB) method was parameterized based on density functional calculations. The surrounding protein and explicit water environment was modeled with conventional, classical force field. The calculations were performed based on HIV-1 protease complexed with the MVT-101 inhibitor that was modified to a model substrate. The protonation state of the catalytic aspartates was determined theoretically. Possible reaction mechanisms involving the lytic water molecule are accounted for in this study. The modeled steps include the dissociation of the lytic water molecule and proton transfer onto Asp-125, the nucleophilic attack followed by a proton transfer onto peptide nitrogen. The simulations show that in the active site most preferable energetically are structures consisting of ionized or polarized molecular fragments that are not accounted for in conventional molecular dynamics. The mobility of the lytic water molecule, the dynamics of the hydrogen bond network, and the conformation of the aspartates in the active center were analyzed.

No abstract