Dataset Information

Application of a polarizable force field to calculations of relative protein-ligand binding affinities.

ABSTRACT: An explicitly polarizable force field based exclusively on quantum data is applied to calculations of relative binding affinities of ligands to proteins. Five ligands, differing by replacement of an atom or functional group, in complexes with three serine proteases-trypsin, thrombin, and urokinase-type plasminogen activator-with available experimental binding data are used as test systems. A special protocol of thermodynamic integration was developed and used to provide sufficiently low levels of systematic error along with high numerical efficiency and statistical stability. The calculated results are in excellent quantitative (rmsd = 1.0 kcal/mol) and qualitative (R(2) = 0.90) agreement with experimental data. The potential of the methodology to explain the observed differences in the ligand affinities is also demonstrated.

SUBMITTER: Khoruzhii O

PROVIDER: S-EPMC2492491 | biostudies-literature | 2008 Jul

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Application of a polarizable force field to calculations of relative protein-ligand binding affinities.

Khoruzhii Oleg O Donchev Alexander G AG Galkin Nikolay N Illarionov Alexei A Olevanov Mikhail M Ozrin Vladimir V Queen Cary C Tarasov Vladimir V

Proceedings of the National Academy of Sciences of the United States of America 20080724 30

An explicitly polarizable force field based exclusively on quantum data is applied to calculations of relative binding affinities of ligands to proteins. Five ligands, differing by replacement of an atom or functional group, in complexes with three serine proteases-trypsin, thrombin, and urokinase-type plasminogen activator-with available experimental binding data are used as test systems. A special protocol of thermodynamic integration was developed and used to provide sufficiently low levels o ...[more]

PMID: 18653760

Similar Datasets

Project description:Predicting relative protein-ligand binding affinities is a central pillar of lead optimization efforts in structure-based drug design. The site identification by ligand competitive saturation (SILCS) methodology is based on functional group affinity patterns in the form of free energy maps that may be used to compute protein-ligand binding poses and affinities. Presented are results obtained from the SILCS methodology for a set of eight target proteins as reported originally in Wang et al. (J. Am. Chem. Soc., 2015, 137, 2695-2703) using free energy perturbation (FEP) methods in conjunction with enhanced sampling and cycle closure corrections. These eight targets have been subsequently studied by many other authors to compare the efficacy of their method while comparing with the outcomes of Wang et al. In this work, we present results for a total of 407 ligands on the eight targets and include specific analysis on the subset of 199 ligands considered previously. Using the SILCS methodology we can achieve an average accuracy of up to 77% and 74% when considering the eight targets with their 199 and 407 ligands, respectively, for rank-ordering ligand affinities as calculated by the percent correct metric. This accuracy increases to 82% and 80%, respectively, when the SILCS atomic free energy contributions are optimized using a Bayesian Markov-chain Monte Carlo approach. We also report other metrics including Pearson's correlation coefficient, Pearlman's predictive index, mean unsigned error, and root mean square error for both sets of ligands. The results obtained for the 199 ligands are compared with the outcomes of Wang et al. and other published works. Overall, the SILCS methodology yields similar or better-quality predictions without a priori need for known ligand orientations in terms of the different metrics when compared to current FEP approaches with significant computational savings while additionally offering quantitative estimates of individual atomic contributions to binding free energies. These results further validate the SILCS methodology as an accurate, computationally efficient tool to support lead optimization and drug discovery.

Project description:Preorganizing a ligand in the conformation it adopts upon binding to a protein has long been considered to be an effective way to improve affinity by making the binding entropy more favorable. However, recent thermodynamic studies of a series of complexes of the Grb2 SH2 domain with peptide analogues having constrained and flexible replacements for a phosphotyrosine residue revealed that less favorable binding entropies may result from constraining ligands in their biologically active conformations. Toward probing the origin of this unexpected finding, we examined the complexes of four phosphotyrosine-derived analogues with the Grb2 SH2 domain using molecular dynamics simulations with a polarizable force field. Significantly, the computed values for the relative binding free energies, entropies, and enthalpies of two pairs of constrained and unconstrained ligands reproduced the trends that were determined experimentally, although the relative differences were overestimated. These calculations also revealed that a large fraction of the ligands lacking the constraining element exist in solution as compact, macrocyclic-like structures that are stabilized by interactions between the phosphate groups and the amide moieties of the C-terminal pY+2 residues. In contrast, the three-membered ring in the constrained ligands prevents the formation of such macrocyclic structures, leading instead to globally extended, less ordered conformations. Quasiharmonic analysis of these conformational ensembles suggests that the unconstrained ligands possess significantly lower entropies in solution, a finding that is consistent with the experimental observation that the binding entropies for the unconstrained ligands are more favorable than for their constrained counterparts. This study suggests that introducing local constraints in flexible molecules may have unexpected consequences, and a detailed understanding of the conformational preferences of ligands in their unbound states is a critical prerequisite to correlating changes in their chemical structure with protein binding entropies and enthalpies.

Dataset Information

Application of a polarizable force field to calculations of relative protein-ligand binding affinities.

Publications

Application of a polarizable force field to calculations of relative protein-ligand binding affinities.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets