Dataset Information

Using the fast fourier transform in binding free energy calculations.

ABSTRACT: According to implicit ligand theory, the standard binding free energy is an exponential average of the binding potential of mean force (BPMF), an exponential average of the interaction energy between the unbound ligand ensemble and a rigid receptor. Here, we use the fast Fourier transform (FFT) to efficiently evaluate BPMFs by calculating interaction energies when rigid ligand configurations from the unbound ensemble are discretely translated across rigid receptor conformations. Results for standard binding free energies between T4 lysozyme and 141 small organic molecules are in good agreement with previous alchemical calculations based on (1) a flexible complex ( R?0.9 for 24 systems) and (2) flexible ligand with multiple rigid receptor configurations ( R?0.8 for 141 systems). While the FFT is routinely used for molecular docking, to our knowledge this is the first time that the algorithm has been used for rigorous binding free energy calculations. © 2017 Wiley Periodicals, Inc.

SUBMITTER: Nguyen TH

PROVIDER: S-EPMC5834390 | biostudies-literature | 2018 Apr

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Using the fast fourier transform in binding free energy calculations.

Nguyen Trung Hai TH Zhou Huan-Xiang HX Minh David D L DDL

Journal of computational chemistry 20171222 11

According to implicit ligand theory, the standard binding free energy is an exponential average of the binding potential of mean force (BPMF), an exponential average of the interaction energy between the unbound ligand ensemble and a rigid receptor. Here, we use the fast Fourier transform (FFT) to efficiently evaluate BPMFs by calculating interaction energies when rigid ligand configurations from the unbound ensemble are discretely translated across rigid receptor conformations. Results for stan ...[more]

PMID: 29270990

Similar Datasets

Project description:The absolute (standard) binding free energy of eight FK506-related ligands to FKBP12 is calculated using free energy perturbation molecular dynamics (FEP/MD) simulations with explicit solvent. A number of features are implemented to improve the accuracy and enhance the convergence of the calculations. First, the absolute binding free energy is decomposed into sequential steps during which the ligand-surrounding interactions as well as various biasing potentials restraining the translation, orientation, and conformation of the ligand are turned "on" and "off." Second, sampling of the ligand conformation is enforced by a restraining potential based on the root mean-square deviation relative to the bound state conformation. The effect of all the restraining potentials is rigorously unbiased, and it is shown explicitly that the final results are independent of all artificial restraints. Third, the repulsive and dispersive free energy contribution arising from the Lennard-Jones interactions of the ligand with its surrounding (protein and solvent) is calculated using the Weeks-Chandler-Andersen separation. This separation also improves convergence of the FEP/MD calculations. Fourth, to decrease the computational cost, only a small number of atoms in the vicinity of the binding site are simulated explicitly, while all the influence of the remaining atoms is incorporated implicitly using the generalized solvent boundary potential (GSBP) method. With GSBP, the size of the simulated FKBP12/ligand systems is significantly reduced, from approximately 25,000 to 2500. The computations are very efficient and the statistical error is small ( approximately 1 kcal/mol). The calculated binding free energies are generally in good agreement with available experimental data and previous calculations (within approximately 2 kcal/mol). The present results indicate that a strategy based on FEP/MD simulations of a reduced GSBP atomic model sampled with conformational, translational, and orientational restraining potentials can be computationally inexpensive and accurate.

Project description:T-wave alternans (TWA) has been associated with increased vulnerability to ventricular tachyarrhythmias and sudden cardiac death. However, both random (white) noise and (patho)physiologic processes (ie, premature ventricular contractions and heart and respiration rates) may hamper TWA estimation and therefore lessen its clinical utility for risk stratification.To investigate the effect of random noise and certain (patho)physiologic processes on the estimation of TWA by using the fast Fourier transform method and to develop methods to overcome these potential sources of error.We used a combination of human electrocardiogram data and computer simulations to assess the effects of a premature ventricular contraction and random and colored noise on the accuracy of TWA estimation.We quantitatively demonstrate that replacing a "bad" beat with an odd/even median beat is a more accurate approach than replacing it with the overall average or the overall median beat. We also show that phase resetting may have a significant effect on alternans estimation and that estimation of alternans by using frequencies >0.4922 cycles/beat in a 128-point fast Fourier transform provides the most accurate approach for estimating the alternans when phase resetting is likely to occur. In addition, our data demonstrate that the number of indeterminate TWA tests due to high levels of noise can be reduced when the alternans voltage exceeds a new higher threshold. Furthermore, the amplitude of random noise has a significant effect on alternans estimation and the alternans voltage threshold should be adjusted for noise levels >1.8 ?V. Finally, we quantitatively demonstrate that colored noise may lead to a false-positive or a false-negative result. We propose methods to estimate the effect of these (patho)physiologic processes on the alternans estimation in order to determine whether a TWA test is likely to be a true positive or a true negative.This study introduces novel methods to overcome potential sources of error in the estimation of TWA. These methods may improve the utility of TWA either for ambulatory monitoring or for clinical risk stratification for ventricular arrhythmias and sudden cardiac death.

Dataset Information

Using the fast fourier transform in binding free energy calculations.

Publications

Using the fast fourier transform in binding free energy calculations.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure