Can an optimization/scoring procedure in ligand-protein docking be employed to probe drug-resistant mutations in...

A simple ligand-protein structural optimization and binding evaluation procedure has been routinely used in high-speed ligand-protein docking studies. In this work, we examine whether such an optimization/scoring procedure is useful in indicating possible drug-resistant mutations in proteins. Crystal structures of three wild-type enzymes (HIV-1 protease, HIV-1 reverse transcriptase, and Mycobacterium tuberculosis H37Rv enoyl-ACP reductase) complexed to a variety of inhibitors are studied. Mutations are introduced into these structures by using the molecular modeling software, SYBYL. Structural optimization and scoring of a mutant complex is conducted by a procedure similar to that used in a recent docking study (Wang et al., 1999). The computed results are compared with observed drug resistance data and the profile of nonresistant mutations. Most mutations studied show an energy change in the same direction as those indicated by observed resistance data. 50% of the polar to polar or nonpolar to nonpolar mutations are found to correlate qualitatively with observed drug resistance data. Van der Waals interactions account for most of these changes, which is in agreement with conclusions from structural studies. Substantially larger deviations are found between computed results and observed data for most polar to nonpolar or nonpolar to polar mutations, which result from deficiency in modelling and scoring ligand-protein interactions in our procedure. Our results suggest that an optimization/docking scoring procedure is useful for qualitatively probing polar to polar or nonpolar to nonpolar resistant mutations in addition to its application in screening active compounds. More accurate description of ligand-protein interactions and the use of methods such as free energy perturbation and Poisson-Boltzmann may be needed to further improve the quality of prediction.