Dataset Information

Ranking Reversible Covalent Drugs: From Free Energy Perturbation to Fragment Docking.

ABSTRACT: Reversible covalent inhibitors have drawn increasing attention in drug design, as they are likely more potent than noncovalent inhibitors and less toxic than covalent inhibitors. Despite those advantages, the computational prediction of reversible covalent binding presents a formidable challenge because the binding process consists of multiple steps and quantum mechanics (QM) level calculation is needed to estimate the covalent binding free energy. It has been shown that the dissociation rates and the equilibrium dissociation constants vary significantly even with similar warheads, due to noncovalent interactions. We have previously used a simplistic two-state model for predicting the relative binding selectivity of reversible covalent inhibitors ( J. Am. Chem. Soc. 2017, 139 , 17945 ). Here we go beyond binding selectivity and demonstrate that it is possible to use free energy perturbation (FEP) molecular dynamics (MD) to calculate the overall reversible covalent binding using a specially designed thermodynamic cycle. We show that FEP can predict the varying binding free energies of the analogs sharing a common warhead. More importantly, our results revealed that the chemical modification away from warhead alters the binding affinity at both noncovalent and covalent binding states, and the computational prediction can be improved by considering the binding free energy of both states. Furthermore, we explored the possibility of using a more rapid computational method, site-identification by ligand competitive saturation (SILCS), to rank the same set of reversible covalent inhibitors. We found that the fragment docking to a set of precomputed fragment maps produces a reasonable ranking. In conclusion, two independent approaches provided consistent results that the covalent binding state is suitable for the initial ranking of the reversible covalent drug candidates. For lead-optimization, the FEP approach designed here can provide more rigorous and detailed information regarding how much the covalent and noncovalent binding states are contributing to the overall binding affinity, thus offering a new avenue for fine-tuning the noncovalent interactions for optimizing reversible covalent drugs.

SUBMITTER: Zhang H

PROVIDER: S-EPMC6610880 | biostudies-literature | 2019 May

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Ranking Reversible Covalent Drugs: From Free Energy Perturbation to Fragment Docking.

Zhang Han H Jiang Wenjuan W Chatterjee Payal P Luo Yun Y

Journal of chemical information and modeling 20190227 5

Reversible covalent inhibitors have drawn increasing attention in drug design, as they are likely more potent than noncovalent inhibitors and less toxic than covalent inhibitors. Despite those advantages, the computational prediction of reversible covalent binding presents a formidable challenge because the binding process consists of multiple steps and quantum mechanics (QM) level calculation is needed to estimate the covalent binding free energy. It has been shown that the dissociation rates a ...[more]

PMID: 30763080

Similar Datasets

Project description:Free-energy calculations play an important role in the application of computational chemistry to a range of fields, including protein biochemistry, rational drug design, or materials science. Importantly, the free-energy difference is directly related to experimentally measurable quantities such as partition and adsorption coefficients, water activity, and binding affinities. Among several techniques aimed at predicting free-energy differences, perturbation approaches, involving the alchemical transformation of one molecule into another through intermediate states, stand out as rigorous methods based on statistical mechanics. However, despite the importance of free-energy calculations, the applicability of the perturbation approaches is still largely impeded by a number of challenges, including the definition of the perturbation path, i.e., alchemical changes leading to the transformation of one molecule to the other. To address this, an automatic perturbation topology builder based on a graph-matching algorithm is developed, which can identify the maximum common substructure (MCS) of two or multiple molecules and provide the perturbation topologies suitable for free-energy calculations using the GROMOS and the GROMACS simulation packages. Various MCS search options are presented leading to alternative definitions of the perturbation pathway. Moreover, perturbation topologies generated using the default multistate MCS search are used to calculate the changes in free energy between lysine and its two post-translational modifications, 3-methyllysine and acetyllysine. The pairwise free-energy calculations performed on this test system led to a cycle closure of 0.5 ± 0.3 and 0.2 ± 0.2 kJ mol-1, with GROMOS and GROMACS simulation packages, respectively. The same relative free energies between the three states are obtained by employing the enveloping distribution sampling (EDS) approach when compared to the pairwise perturbations. Importantly, this toolkit is made available online as an open-source Python package (https://github.com/drazen-petrov/SMArt).

Dataset Information

Ranking Reversible Covalent Drugs: From Free Energy Perturbation to Fragment Docking.

Publications

Ranking Reversible Covalent Drugs: From Free Energy Perturbation to Fragment Docking.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets