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:Binding free energy calculations predict the potency of compounds to protein binding sites in a physically rigorous manner and see broad application in prioritizing the synthesis of novel drug candidates. Relative binding free energy (RBFE) calculations have emerged as an industry-standard approach to achieve highly accurate rank-order predictions of the potency of related compounds; however, this approach requires that the ligands share a common scaffold and a common binding mode, restricting the methods' domain of applicability. This is a critical limitation since complex modifications to the ligands, especially core hopping, are very common in drug design. Absolute binding free energy (ABFE) calculations are an alternate method that can be used for ligands that are not congeneric. However, ABFE suffers from a known problem of long convergence times due to the need to sample additional degrees of freedom within each system, such as sampling rearrangements necessary to open and close the binding site. Here, we report on an alternative method for RBFE, called Separated Topologies (SepTop), which overcomes the issues in both of the aforementioned methods by enabling large scaffold changes between ligands with a convergence time comparable to traditional RBFE. Instead of only mutating atoms that vary between two ligands, this approach performs two absolute free energy calculations at the same time in opposite directions, one for each ligand. Defining the two ligands independently allows the comparison of the binding of diverse ligands without the artificial constraints of identical poses or a suitable atom-atom mapping. This approach also avoids the need to sample the unbound state of the protein, making it more efficient than absolute binding free energy calculations. Here, we introduce an implementation of SepTop. We developed a general and efficient protocol for running SepTop, and we demonstrated the method on four diverse, pharmaceutically relevant systems. We report the performance of the method, as well as our practical insights into the strengths, weaknesses, and challenges of applying this method in an industrial drug design setting. We find that the accuracy of the approach is sufficiently high to rank order ligands with an accuracy comparable to traditional RBFE calculations while maintaining the additional flexibility of SepTop.

Project description:BackgroundT-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.ObjectiveTo 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.MethodsWe 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.ResultsWe 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.ConclusionThis 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

Tweets