We present a new computational strategy for calculating the absolute binding free energy for protein ligand association in the context of atomistic simulation in explicit solvent. The method is based on an appropriate definition of a solute tempering scheme enforced via Hamilton replica exchange method (HREM). The definition of "solute" includes both the ligand and the active site, with the remainder of the systems defined as "solvent". The hydrophilicity of the solute and the solute torsional plus nonbonded intrasolute interactions are increased and decreased, respectively, along the replica progression, thus favoring the extrusion of the drug form the active site in the scaled states of the generalized ensemble. The proposed technique, named "Energy Driven Undocking" (EDU-HREM), completely bypasses the need for defining and/or identifying the relevant reaction coordinates in a ligand receptor interactions and allows the calculation of the absolute binding free energy in one single generalized simulation of the drug-receptor system. The methodology is applied, with encouraging results, to the calculation of the absolute binding free energy of some FK506-related ligands of the peptidyl prolyl cis-trans isomerase protein (FKBP12) with known dissociation constants. Aspects of the binding/inhibition mechanism in FKBP12 are also analyzed and discussed.

Energy-Driven Undocking (EDU-HREM) in Solute Tempering Replica Exchange Simulations / Piero Procacci;Marco Bizzarri;Simone Marsili. - In: JOURNAL OF CHEMICAL THEORY AND COMPUTATION. - ISSN 1549-9618. - STAMPA. - 10:(2014), pp. 439-450. [10.1021/ct400809n]

Energy-Driven Undocking (EDU-HREM) in Solute Tempering Replica Exchange Simulations

PROCACCI, PIERO
;
BIZZARRI, MARCO;MARSILI, SIMONE
2014

Abstract

We present a new computational strategy for calculating the absolute binding free energy for protein ligand association in the context of atomistic simulation in explicit solvent. The method is based on an appropriate definition of a solute tempering scheme enforced via Hamilton replica exchange method (HREM). The definition of "solute" includes both the ligand and the active site, with the remainder of the systems defined as "solvent". The hydrophilicity of the solute and the solute torsional plus nonbonded intrasolute interactions are increased and decreased, respectively, along the replica progression, thus favoring the extrusion of the drug form the active site in the scaled states of the generalized ensemble. The proposed technique, named "Energy Driven Undocking" (EDU-HREM), completely bypasses the need for defining and/or identifying the relevant reaction coordinates in a ligand receptor interactions and allows the calculation of the absolute binding free energy in one single generalized simulation of the drug-receptor system. The methodology is applied, with encouraging results, to the calculation of the absolute binding free energy of some FK506-related ligands of the peptidyl prolyl cis-trans isomerase protein (FKBP12) with known dissociation constants. Aspects of the binding/inhibition mechanism in FKBP12 are also analyzed and discussed.
2014
10
439
450
Piero Procacci;Marco Bizzarri;Simone Marsili
File in questo prodotto:
Non ci sono file associati a questo prodotto.

I documenti in FLORE sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificatore per citare o creare un link a questa risorsa: https://hdl.handle.net/2158/840895
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 12
  • ???jsp.display-item.citation.isi??? 12
social impact