Binding free energies and free energy components from molecular dynamics and Poisson-Boltzmann calculations. Application to amino acid recognition by aspartyl-tRNA synthetase11Edited by A. R. Fersht

2001 ◽  
Vol 306 (2) ◽  
pp. 307-327 ◽  
Author(s):  
Georgios Archontis ◽  
Thomas Simonson ◽  
Martin Karplus
2021 ◽  
Author(s):  
Alexander Wade ◽  
Agastya Bhati ◽  
Shunzhou Wan ◽  
Peter Coveney

The binding free energy between a ligand and its target protein is an essential quantity to know at all stages of the drug discovery pipeline. Assessing this value computationally can offer insight into where efforts should be focused in the pursuit of effective therapeutics to treat myriad diseases. In this work we examine the computation of alchemical relative binding free energies with an eye to assessing reproducibility across popular molecular dynamics packages and free energy estimators. The focus of this work is on 54 ligand transformations from a diverse set of protein targets: MCL1, PTP1B, TYK2, CDK2 and thrombin. These targets are studied with three popular molecular dynamics packages: OpenMM, NAMD2 and NAMD3. Trajectories collected with these packages are used to compare relative binding free energies calculated with thermodynamic integration and free energy perturbation methods. The resulting binding free energies show good agreement between molecular dynamics packages with an average mean unsigned error between packages of 0.5 $kcal/mol$ The correlation between packages is very good with the lowest Spearman's, Pearson's and Kendall's tau correlation coefficient between two packages being 0.91, 0.89 and 0.74 respectively. Agreement between thermodynamic integration and free energy perturbation is shown to be very good when using ensemble averaging.


Author(s):  
Shunzhou Wan ◽  
Peter V Coveney ◽  
Darren R Flower

The binding to the T cell receptor of wild-type and variant HTLV-1 Tax peptide complexed to the major histocompatibility complex has been investigated by means of molecular dynamics simulations. The binding free energy difference is calculated using the molecular mechanics Poisson–Boltzmann surface area and linear interaction energy methods. These methods extract useful information on the binding energetics from simulations of the physical states of the ligands, which are more computationally expedient than the commonly used thermodynamic integration method. The successful reproduction of the relative binding free energies shows that these methods can be useful for free energy calculations and the rational design of drugs and vaccines.


2008 ◽  
Vol 7 (3) ◽  
pp. 103-116 ◽  
Author(s):  
Atsushi SUENAGA ◽  
Osamu UMEZU ◽  
Tadashi ANDO ◽  
Ichiro YAMATO ◽  
Takeshi MURATA ◽  
...  

2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Germano Heinzelmann ◽  
Michael K. Gilson

AbstractAbsolute binding free energy calculations with explicit solvent molecular simulations can provide estimates of protein-ligand affinities, and thus reduce the time and costs needed to find new drug candidates. However, these calculations can be complex to implement and perform. Here, we introduce the software BAT.py, a Python tool that invokes the AMBER simulation package to automate the calculation of binding free energies for a protein with a series of ligands. The software supports the attach-pull-release (APR) and double decoupling (DD) binding free energy methods, as well as the simultaneous decoupling-recoupling (SDR) method, a variant of double decoupling that avoids numerical artifacts associated with charged ligands. We report encouraging initial test applications of this software both to re-rank docked poses and to estimate overall binding free energies. We also show that it is practical to carry out these calculations cheaply by using graphical processing units in common machines that can be built for this purpose. The combination of automation and low cost positions this procedure to be applied in a relatively high-throughput mode and thus stands to enable new applications in early-stage drug discovery.


Molecules ◽  
2021 ◽  
Vol 26 (6) ◽  
pp. 1586
Author(s):  
Leonor Contreras ◽  
Ignacio Villarroel ◽  
Camila Torres ◽  
Roberto Rozas

Doxorubicin (DOX), a recognized anticancer drug, forms stable associations with carbon nanotubes (CNTs). CNTs when properly functionalized have the ability to anchor directly in cancerous tumors where the release of the drug occurs thanks to the tumor slightly acidic pH. Herein, we study the armchair and zigzag CNTs with Stone–Wales (SW) defects to rank their ability to encapsulate DOX by determining the DOX-CNT binding free energies using the MM/PBSA and MM/GBSA methods implemented in AMBER16. We investigate also the chiral CNTs with haeckelite defects. Each haeckelite defect consists of a pair of square and octagonal rings. The armchair and zigzag CNT with SW defects and chiral nanotubes with haeckelite defects predict DOX-CNT interactions that depend on the length of the nanotube, the number of present defects and nitrogen doping. Chiral nanotubes having two haeckelite defects reveal a clear dependence on the nitrogen content with DOX-CNT interaction forces decreasing in the order 0N > 4N > 8N. These results contribute to a further understanding of drug-nanotube interactions and to the design of new drug delivery systems based on CNTs.


1998 ◽  
Vol 275 (5) ◽  
pp. 823-846 ◽  
Author(s):  
Georgios Archontis ◽  
Thomas Simonson ◽  
Dino Moras ◽  
Martin Karplus

2021 ◽  
Author(s):  
Safak OZHAN KOCAKAYA

Abstract Recently, protein tyrosine phosphatase 1B (PTP1B) inhibitors have become the frontier as possible targeting for anti-cancer and antidiabetic drugs. The contemporary observe represents a pc assisted version to investigate the importance of precise residues within the binding web site of PTP1B with numerous Sanggenon derivatives remoted from nature. Molecular dynamics (MD) simulations were performed to estimate the dynamics of the complexes, and absolute binding unfastened energies have been calculated with exclusive additives, and carried out through the usage of the Molecular Mechanics-Poisson-Boltzmann floor region (MM-PB/SA) and Generalized Born surface vicinity (MM-GB/SA) strategies. The effects show that the expected free energies of the complexes are normally constant with the available experimental statistics. MM/GBSA free energy decomposition analysis shows that the residues Asp29, Arg24, Met258, and , Arg254 in the second active site in PTP1B are crucial for the excessive selectivity of the inhibitors.


2021 ◽  
Author(s):  
Yuriy Khalak ◽  
Gary Tresdern ◽  
Matteo Aldeghi ◽  
Hannah Magdalena Baumann ◽  
David L. Mobley ◽  
...  

The recent advances in relative protein-ligand binding free energy calculations have shown the value of alchemical methods in drug discovery. Accurately assessing absolute binding free energies, although highly desired, remains...


Sign in / Sign up

Export Citation Format

Share Document