The University of Southampton
University of Southampton Institutional Repository

The development of free energy methods for protein-ligand complexes

The development of free energy methods for protein-ligand complexes
The development of free energy methods for protein-ligand complexes

Methods for the calculation of the relative binding free energies of ligands to a protein are investigated and developed.  The aim of these investigations was to improve the reliability and speed of free energy methods, such that they become practical tools for commercial rational drug design. To this end, the relative hydration free energy of water methane, and the relative binding free energies of halides to a calix[4]pyrrole derivative were investigated by three established free energy methods (Free Energy Perturbation (FEP). Finite Difference Thermodynamic Integration (FDTI) and Adaptive Umbrella WHAM (Ad-UmWHAM). The results of these applications showed that inconsistencies in sampling led to unreliable free energy predictions. To overcome these problems, a series of four new free energy methods were developed (Bivariate Multicanonical WHAM (BMW), Parallel Tempering Thermodynamic Integration (PTTI), Replica Exchange Free Energy Perturbation (REFEP) and Replica Exchange Thermodynamic Integration (RETI). These approaches all combined traditional free energy methods with generalised ensembles. Testing of these methods revealed that Replica Exchange Thermodynamic Integration was the superior of all seven methods. FDTI and RETI were then tested by calculating the relative binding free energies of a group of SB1-like ligands to p38 MAP kinase.  The results of this test showed that RETI was still the superior method. This test also revealed that there were still sampling issues that needed to be resolved. A new Monte Carlo code was developed to run the tests on p38. The optimised data structure of the code led to a ten to twelve fold speed up compared to an established MC code.  This, combined with the use of a large Linux Beowulf cluster, enabled each protein-ligand free energy calculation to be run within 1.5 days. We predict that, using the latest computers, these calculations could take less than 12 hours.

University of Southampton
Woods, Christopher J
54c21a03-755f-4977-aa2c-72c1b073e734
Woods, Christopher J
54c21a03-755f-4977-aa2c-72c1b073e734

Woods, Christopher J (2003) The development of free energy methods for protein-ligand complexes. University of Southampton, Doctoral Thesis.

Record type: Thesis (Doctoral)

Abstract

Methods for the calculation of the relative binding free energies of ligands to a protein are investigated and developed.  The aim of these investigations was to improve the reliability and speed of free energy methods, such that they become practical tools for commercial rational drug design. To this end, the relative hydration free energy of water methane, and the relative binding free energies of halides to a calix[4]pyrrole derivative were investigated by three established free energy methods (Free Energy Perturbation (FEP). Finite Difference Thermodynamic Integration (FDTI) and Adaptive Umbrella WHAM (Ad-UmWHAM). The results of these applications showed that inconsistencies in sampling led to unreliable free energy predictions. To overcome these problems, a series of four new free energy methods were developed (Bivariate Multicanonical WHAM (BMW), Parallel Tempering Thermodynamic Integration (PTTI), Replica Exchange Free Energy Perturbation (REFEP) and Replica Exchange Thermodynamic Integration (RETI). These approaches all combined traditional free energy methods with generalised ensembles. Testing of these methods revealed that Replica Exchange Thermodynamic Integration was the superior of all seven methods. FDTI and RETI were then tested by calculating the relative binding free energies of a group of SB1-like ligands to p38 MAP kinase.  The results of this test showed that RETI was still the superior method. This test also revealed that there were still sampling issues that needed to be resolved. A new Monte Carlo code was developed to run the tests on p38. The optimised data structure of the code led to a ten to twelve fold speed up compared to an established MC code.  This, combined with the use of a large Linux Beowulf cluster, enabled each protein-ligand free energy calculation to be run within 1.5 days. We predict that, using the latest computers, these calculations could take less than 12 hours.

Text
914698.pdf - Version of Record
Available under License University of Southampton Thesis Licence.
Download (18MB)

More information

Published date: 2003

Identifiers

Local EPrints ID: 465115
URI: http://eprints.soton.ac.uk/id/eprint/465115
PURE UUID: 5fde7bc4-9d49-488e-ac3f-4fec5a597c95

Catalogue record

Date deposited: 05 Jul 2022 00:24
Last modified: 16 Mar 2024 19:57

Export record

Contributors

Author: Christopher J Woods

Download statistics

Downloads from ePrints over the past year. Other digital versions may also be available to download e.g. from the publisher's website.

View more statistics

Atom RSS 1.0 RSS 2.0

Contact ePrints Soton: eprints@soton.ac.uk

ePrints Soton supports OAI 2.0 with a base URL of http://eprints.soton.ac.uk/cgi/oai2

This repository has been built using EPrints software, developed at the University of Southampton, but available to everyone to use.

We use cookies to ensure that we give you the best experience on our website. If you continue without changing your settings, we will assume that you are happy to receive cookies on the University of Southampton website.

×