title: II. Dissociation free energies in drug-receptor systems via nonequilibrium alchemical simulations: application to the FK506-related immunophilin ligands.
abstract: The recently proposed fast switching double annihilation (FS-DAM) [Cardelli et al., J. Chem. Theory Comput., 2015, 11, 423] is aimed at computing the absolute standard dissociation free energies for the chemical equilibrium RL ⇌ R + L occurring in solution through molecular dynamics (MD) simulations at the atomistic level. The technique is based on the production of fast nonequilibrium annihilation trajectories of one of the species (the ligand) in the solvated RL complex and in the bulk solvent. As detailed in the companion theoretical paper, the free energies of these two nonequilibrium annihilation processes are recovered by using an unbiased unidirectional estimate derived from the Crooks theorem exploiting the inherent Gaussian nature of the annihilation work. The FS-DAM technique was successfully applied to the evaluation of the dissociation free energy of the complexes of Zn(ii) cations with an inhibitor of the Tumor Necrosis Factor α converting enzyme. Here we apply the technique to a real drug-receptor system, by satisfactorily reproducing the experimental dissociation free energies of FK506-related bulky ligands towards the native FKBP12 enzyme and by predicting the dissociation constants for the same ligands towards the mutant I56D. The effect of such mutations on the binding affinity of FK506-related ligands is relevant for assessing the thermodynamic forces regulating molecular recognition in FKBP12 inhibition.
authors.affiliation: Department of Physics, University of Vienna, Austria.
authors.forename: Francesca
authors.initials: F
authors.lastname: Nerattini
authors.affiliation: Department of Chemistry, University of Florence, Italy. procacci@unifi.it.
authors.forename: Riccardo
authors.initials: R
authors.lastname: Chelli
authors.affiliation: Department of Chemistry, University of Florence, Italy. procacci@unifi.it.
authors.forename: Piero
authors.initials: P
authors.lastname: Procacci
mesh_terms: D065095:Calcineurin Inhibitors / Q000378:metabolism* / Q000494:pharmacology; D055808:Drug Discovery; D006801:Humans; D007166:Immunosuppressive Agents / Q000378:metabolism* / Q000494:pharmacology; D008024:Ligands; D056004:Molecular Dynamics Simulation; D017354:Point Mutation; D011485:Protein Binding; D011487:Protein Conformation; D016559:Tacrolimus / Q000378:metabolism* / Q000494:pharmacology; D022061:Tacrolimus Binding Protein 1A / Q000037:antagonists & inhibitors / Q000737:chemistry / Q000235:genetics / Q000378:metabolism*; D013816:Thermodynamics
chemical_list: D065095:Calcineurin Inhibitors; D007166:Immunosuppressive Agents; D008024:Ligands; D022061:Tacrolimus Binding Protein 1A; D016559:Tacrolimus