Theoretical studies of relaxation of a monomeric subunit of HIV-1 protease in water using molecular dynamics

The dynamic behavior of one 99-residue subunit of the dimeric aspartyl protease of HIV-1 was studied in a 160 psec molecular dynamics simulation at 300 K in water. The crystal structure of one of the identical subunits of the dimer was the starting point, with the aqueous phase modeled by 4,331 explicit waters in a restrained spherical droplet. Analysis of the simulations showed that the monomer displayed considerable flexibility in the interfacial portions of the flap (the region which folds over the substrate), the N- and C-termini, and, to a lesser extent, the active site. The flap undergoes significant motion as an independent rigid finger, but without the cantilever previously reported in a simulation of the dimer. The N-terminus displayed the greatest fluctuational disorder whereas the C-terminus exhibited the greatest root mean square movement from the crystal structure. The central core of the monomer had a heavy-atom root mean square deviation from the initial structure of about 3.0 A during the latter half of the simulation. Although this is larger than the 1.6 A found for comparable simulations of typical globular proteins, the general features of the tertiary structure were preserved over the course of the simulation. Overall, these results indicate that the relaxed structure obtained in these simulations may provide a better model for the tertiary structure of the solvated HIV-1 protease monomer than the subunit conformation seen in the X-ray crystallographic structure of the dimer. Except in the flap region, the design of compounds intended to interfere with dimerization should take this relaxation and the flexibility of the solvated monomer, especially at the termini, into account.