Modeling studies were performed on HCV NS5B Polymerase in an effort to design new inhibitors. The binding models of five different scaffold inhibitors were investigated and compared by using molecular dynamics simulations, free energy calculation and decomposition. Our results show Tyr448 plays the most critical role in the binding of most inhibitors. In addition, favorable contributions of residues Pro197, Arg200, Cys366, Met414 and Tyr448 in a deep hydrophobic pocket prove to be important for the selectivity of inhibitors. Furthermore, an optimized docking protocol was presented based on cross-docking the five inhibitors in the palm binding site of this enzyme using the Autodock program. This protocol was used later to virtually screen NCI and Maybridge diversity set libraries. The binding site was profiled via the statistics and analysis of the hydrogen bond networks formed between the receptor and the top-ranked diversity set compounds. Based on our detailed binding site analysis two useful rules were proposed to guide the selection of promising hits.Response to Reviewers: Question1: It is a pity that nothing in detail is discussed about found structures resulting from the virtual screening. Did the authors find new promising ligands with high (higher than the template ligands?) affinity and are they structurally diverse or similar to the used ones?Answer: We have found some promising compounds after the virtual screening. They are structurally diverse. Most of them form hydrogen bonds with the critical residues such as Tyr448 and have strong hydrophobic interaction with the residues in the deep hydrophobic pocket mentioned in our article. These compounds are under biological test.Question2: Additionally to figure 5 at least one figure should be added, showing the amino acid residues of the active site with the discussed most important interactions with one (the best?) ligand. This would considerably help to understand the discussion without looking in the original pdb-files.
AbstractModeling studies were performed on HCV NS5B Polymerase in an effort to design new inhibitors. The binding models of five different scaffold inhibitors were investigated and compared by using molecular dynamics simulations, free energy calculation and decomposition. Our results show Tyr448 plays the most critical role in the binding of most inhibitors. In addition, favorable contributions of residues Pro197, Arg200, Cys366, Met414 and Tyr448 in a deep hydrophobic pocket prove to be important for the selectivity of inhibitors.Furthermore, an optimized docking protocol was presented based on cross-docking the five inhibitors in the palm binding site of this enzyme using the Autodock program. This protocol was used later to virtually screen NCI and Maybridge diversity set libraries. The binding site was profiled via the statistics and analysis of the hydrogen bond networks formed between the receptor and the top-ranked diversity set compounds. Based on our detailed binding site analysis two useful rules were proposed to ...