Abstract:Water molecules play a pivotal role in modulating protein-ligand interactions, wielding a substantial influence on ligand binding affinity and orientation in many systems. Water molecules can facilitate the interaction between a protein and ligand through hydrogen-bond networks. Conversely, water molecules occupying a ligand binding site may need to be displaced at an energetic cost prior to ligand binding. As a result, appropriate treatment of water contributions to protein-ligand interactions is important for ligand design. However, this is a very challenging problem in the context of adequately determining the number of waters to consider and undertaking the conformational sampling of the ligands, the waters, and the surrounding protein. In the present study, an extension of the Site Identification by Ligand Competitive Saturation-Monte Carlo (SILCS-MC) docking approach is presented that enables determination of the location of water molecules in the binding pocket and their impact on the predicted ligand binding orientation and affinities. The approach, termed SILCS-WATER, involves initial exhaustive MC docking of the ligand and 15 water molecules in a binding site followed by selection of a subset of waters for additional MC sampling. In the approach, the Metropolis criteria for the MC is based on the overlap of the ligand and water with the SILCS FragMaps and the interaction energy between the waters and ligand. Results show that the SILCS-WATER methodology is able to capture important waters and improve ligand positions and orientation. For 6 of 10 multiple-ligand protein systems the method improved the relative binding affinity prediction against experimental results to varying degrees when compared to standard SILCS-MC. The method identifies waters interacting with ligands that occupy unfavorable locations with respect to the protein whose displacement through the appropriate ligands modifications should improve ligands binding affinity. Results are consistent with the binding affinity being modeled as a ligand-water complex interacting with the protein. The presented approach offers new possibilities in revealing water networks and their contributions to the binding affinity of a ligand to a protein and is anticipated to be of utility for computer-aided drug design.

Identification of Protein-Ligand Binding Sites by the Level-Set Variational Implicit-Solvent Approach.

Coupling Monte Carlo, Variational Implicit Solvation, and Binary Level-Set for Simulations of Biomolecular Binding.

Level-Set Variational Implicit-Solvent Modeling of Biomolecules with the Coulomb-Field Approximation

Variational Implicit-Solvent Predictions Of The Dry-Wet Transition Pathways For Ligand-Receptor Binding And Unbinding Kinetics

Identification of Cavities on Protein Surface Using Multiple Computational Approaches for Drug Binding Site Prediction

Combined physics- and machine-learning-based method to identify druggable binding sites using SILCS-Hotspots

Identification of Protein-Ligand Binding Sites by Sequence Information and Ensemble Classifier.

Identifying ligand binding sites and poses using GPU-accelerated Hamiltonian replica exchange molecular dynamics

Optimizing the Affinity and Specificity of Ligand Binding with the Inclusion of Solvation Effect.

Detection of Putative Ligand Dissociation Pathways in Proteins using Site-Identification by Ligand Competitive Saturation

PointSite: A Point Cloud Segmentation Tool for Identification of Protein Ligand Binding Atoms

Modeling Ligand Binding Site Water Networks with Site-Identification by Ligand Competitive Saturation: Impact on Ligand Binding Orientations and Relative Binding Affinities

Application of Site-Identification by Ligand Competitive Saturation in Computer-Aided Drug Design

Variational Implicit-Solvent Modeling of Host-Guest Binding: A Case Study on Cucurbit[7]uril|

Binding site detection and druggability prediction of protein targets for structure-based drug design.

Protein-protein interface analysis and hot spots identification for chemical ligand design.

Protein-Ligand Binding Site Recognition Using Complementary Binding-Specific Substructure Comparison And Sequence Profile Alignment

Virtual screening using molecular simulations

Predicting Protein-Ligand Binding Sites Based on an Improved Geometric Algorithm

Protein-protein binding site identification by enumerating the configurations

“Martinizing” the Variational Implicit Solvent Method (VISM): Solvation Free Energy for Coarse-Grained Proteins