Partially folded states of HIV-1 protease: Molecular dynamics simulations and ligand binding

All-atom molecular dynamics simulations were performed on partially folded states (with different secondary structure content) of the dimeric enzyme HIV-1 protease in aqueous solution. The calculations were based on previous simulations of the folding process of the protein based on a Go-model. The structures turn out to be stable, and the subunit–subunit contact surface is smaller than that of the native state. Interestingly, the flexibility of the partially folded states is similar to that observed for the monomer in the native state. The intersubunit contacts are formed by conserved residues, suggesting that these residues may play a role for the folding process. Docking a large set of molecules suggests that several ligands not yet associated to HIV-1 protease may bind to these partially unfolded structures.