Computer Simulations of Ligand-Protein Binding with Ensembles of Protein Conformations: A Monte Carlo Study of Hiv-1 Protease Binding Energy Landscapes

We present the results of molecular docking simulations with HIV-1 protease for the sb203386 and skf107457 inhibitors by Monte Carlo simulated annealing. A simplified piecewise linear energy function, the standard AMBER force field, and the AMBER force field with solvation and a soft-core smoothing component are employed in simulations with a single-protein conformation to determine the relationship between docking simulations with a simple energy function and more realistic force fields. The temperature-dependent binding free energy profiles of the inhibitors interacting with a single protein conformation provide a detailed picture of relative thermodynamic stability and a distribution of ligand binding modes in agreement with experimental crystallographic data. Using the simplified piecewise linear energy function, we also performed Monte Carlo docking simulations with an ensemble of protein conformations employing preferential biased sampling of low-energy protein conformations, and the results are analyzed in connection with the free energy profiles. (C) 1999 John Wiley & Sons, Inc.