Dataset Information

On the use of orientational restraints and symmetry corrections in alchemical free energy calculations.

ABSTRACT: Alchemical free energy calculations are becoming a useful tool for calculating absolute binding free energies of small molecule ligands to proteins. Here, we find that the presence of multiple metastable ligand orientations can cause convergence problems when distance restraints alone are used. We demonstrate that the use of orientational restraints can greatly accelerate the convergence of these calculations. However, even with this acceleration, we find that sufficient sampling requires substantially longer simulations than are used in many published protocols. To further accelerate convergence, we introduce a new method of configuration space decomposition by orientation which reduces required simulation lengths by at least a factor of 5 in the cases examined. Our method is easily parallelizable, well suited for cases where a ligand cocrystal structure is not available, and can utilize initial orientations generated by docking packages.

SUBMITTER: Mobley DL

PROVIDER: S-EPMC3583553 | biostudies-literature | 2006 Aug

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

On the use of orientational restraints and symmetry corrections in alchemical free energy calculations.

Mobley David L DL Chodera John D JD Dill Ken A KA

The Journal of chemical physics 20060801 8

Alchemical free energy calculations are becoming a useful tool for calculating absolute binding free energies of small molecule ligands to proteins. Here, we find that the presence of multiple metastable ligand orientations can cause convergence problems when distance restraints alone are used. We demonstrate that the use of orientational restraints can greatly accelerate the convergence of these calculations. However, even with this acceleration, we find that sufficient sampling requires substa ...[more]

PMID: 16965052

Similar Datasets

Project description:Progress in the development of GPU-accelerated free energy simulation software has enabled practical applications on complex biological systems and fueled efforts to develop more accurate and robust predictive methods. In particular, this work re-examines concerted (a.k.a., one-step or unified) alchemical transformations commonly used in the prediction of hydration and relative binding free energies (RBFEs). We first classify several known challenges in these calculations into three categories: endpoint catastrophes, particle collapse, and large gradient-jumps. While endpoint catastrophes have long been addressed using softcore potentials, the remaining two problems occur much more sporadically and can result in either numerical instability (i.e., complete failure of a simulation) or inconsistent estimation (i.e., stochastic convergence to an incorrect result). The particle collapse problem stems from an imbalance in short-range electrostatic and repulsive interactions and can, in principle, be solved by appropriately balancing the respective softcore parameters. However, the large gradient-jump problem itself arises from the sensitivity of the free energy to large values of the softcore parameters, as might be used in trying to solve the particle collapse issue. Often, no satisfactory compromise exists with the existing softcore potential form. As a framework for solving these problems, we developed a new family of smoothstep softcore (SSC) potentials motivated by an analysis of the derivatives along the alchemical path. The smoothstep polynomials generalize the monomial functions that are used in most implementations and provide an additional path-dependent smoothing parameter. The effectiveness of this approach is demonstrated on simple yet pathological cases that illustrate the three problems outlined. With appropriate parameter selection, we find that a second-order SSC(2) potential does at least as well as the conventional approach and provides vast improvement in terms of consistency across all cases. Last, we compare the concerted SSC(2) approach against the gold-standard stepwise (a.k.a., decoupled or multistep) scheme over a large set of RBFE calculations as might be encountered in drug discovery.

Dataset Information

On the use of orientational restraints and symmetry corrections in alchemical free energy calculations.

Publications

On the use of orientational restraints and symmetry corrections in alchemical free energy calculations.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets