Abstract
This paper describes a large-scale study on the nature and the energetics of the conformational changes drug-like molecules experience upon binding. Ligand strain energies and conformational reorganization were analyzed with different computational methods on 150 crystal structures of pharmaceutically relevant protein-ligand complexes. The common knowledge that ligands rarely bind in their lowest calculated energy conformation was confirmed. Additionally, we found that over 60% of the ligands do not bind in a local minimum conformation. While approximately 60% of the ligands were calculated to bind with strain energies lower than 5 kcal/mol, strain energies over 9 kcal/mol were calculated in at least 10% of the cases regardless of the method used. A clear correlation was found between acceptable strain energy and ligand flexibility, while there was no correlation between strain energy and binding affinity, thus indicating that expensive conformational rearrangements can be tolerated in some cases without overly penalizing the tightness of binding. On the basis of the trends observed, thresholds for the acceptable strain energies of bioactive conformations were defined with consideration of the impact of ligand flexibility. An analysis of the degree of folding of the bound ligands confirmed the general tendency of small molecules to bind in an extended conformation. The results suggest that the unfolding of hydrophobic ligands during binding, which exposes hydrophobic surfaces to contact with protein residues, could be one of the factors accounting for high reorganization energies. Finally, different methods for conformational analysis were evaluated, and guidelines were defined to maximize the prevalence of bioactive conformations in computationally generated ensembles.
Keywords
Affiliated Institutions
Related Publications
Do salt bridges stabilize proteins? A continuum electrostatic analysis
Abstract The electrostatic contribution to the free energy of folding was calculated fo 21 salt bridges in 9 r X‐ray crystal structures using a continuum electrostatic approach ...
Ab initio conformational analysis of nucleic acid components: Intrinsic energetic contributions to nucleic acid structure and dynamics
In recent years, the use of high-level ab initio calculations has allowed for the intrinsic conformational properties of nucleic acid building blocks to be revisited. This has p...
Extra Precision Glide: Docking and Scoring Incorporating a Model of Hydrophobic Enclosure for Protein−Ligand Complexes
A novel scoring function to estimate protein-ligand binding affinities has been developed and implemented as the Glide 4.0 XP scoring function and docking protocol. In addition ...
State-specific protein-ligand complex structure prediction with a multi-scale deep generative model
The binding complexes formed by proteins and small molecule ligands are ubiquitous and critical to life. Despite recent advancements in protein structure prediction, existing al...
Converging free energy estimates: MM‐PB(GB)SA studies on the protein–protein complex Ras–Raf
Abstract Estimating protein–protein interaction energies is a very challenging task for current simulation protocols. Here, absolute binding free energies are reported for the c...
Publication Info
- Year
- 2004
- Type
- article
- Volume
- 47
- Issue
- 10
- Pages
- 2499-2510
- Citations
- 439
- Access
- Closed
External Links
Social Impact
Social media, news, blog, policy document mentions
Citation Metrics
Cite This
Identifiers
- DOI
- 10.1021/jm030563w