Dataset Information

Relative Binding Enthalpies from Molecular Dynamics Simulations Using a Direct Method.

ABSTRACT: The potential for reliably predicting relative binding enthalpies, ??E, from a direct method utilizing molecular dynamics is examined for a system of three phosphotyrosyl peptides binding to a protein receptor, the Src SH2 domain. The binding enthalpies were calculated from the potential energy differences between the bound and the unbound end-states of each peptide from equilibrium simulations in explicit water. The statistical uncertainties in the ensemble-mean energy values from multiple, independent simulations were obtained using a bootstrap method. Simulations were initiated with different starting coordinates as well as different velocities. Statistical uncertainties in ??E are 2 to 3 kcal/mol based on calculations from 40, 10 ns trajectories for each system (three SH2-peptide complexes or unbound peptides). Uncertainties in relative component energies, comprising solute-solute, solute-solvent and solvent-solvent interactions, are considerably larger. Energy values were estimated from an unweighted ensemble averaging of multiple trajectories with the a priori assumption that all trajectories are equally likely. Distributions in energy-rmsd space indicate that the trajectories sample the same basin and the difference in mean energy values between trajectories is due to sampling of alternative local regions of this superbasin. The direct estimate of relative binding enthalpies is concluded to be a reasonable approach for well-ordered systems with ??E values greater than ?3 kcal/mol, although the approach would benefit from future work to determine properly distributed starting points that would enable efficient sampling of conformational space using multiple trajectories.

SUBMITTER: Roy A

PROVIDER: S-EPMC4095907 | biostudies-literature | 2014 Jul

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Relative Binding Enthalpies from Molecular Dynamics Simulations Using a Direct Method.

Roy Amitava A Hua Duy P DP Ward Joshua M JM Post Carol Beth CB

Journal of chemical theory and computation 20140612 7

The potential for reliably predicting relative binding enthalpies, ΔΔ<i>E</i>, from a direct method utilizing molecular dynamics is examined for a system of three phosphotyrosyl peptides binding to a protein receptor, the Src SH2 domain. The binding enthalpies were calculated from the potential energy differences between the bound and the unbound end-states of each peptide from equilibrium simulations in explicit water. The statistical uncertainties in the ensemble-mean energy values from multip ...[more]

PMID: 25061444

Similar Datasets

Project description:The absolute (standard) binding free energy of eight FK506-related ligands to FKBP12 is calculated using free energy perturbation molecular dynamics (FEP/MD) simulations with explicit solvent. A number of features are implemented to improve the accuracy and enhance the convergence of the calculations. First, the absolute binding free energy is decomposed into sequential steps during which the ligand-surrounding interactions as well as various biasing potentials restraining the translation, orientation, and conformation of the ligand are turned "on" and "off." Second, sampling of the ligand conformation is enforced by a restraining potential based on the root mean-square deviation relative to the bound state conformation. The effect of all the restraining potentials is rigorously unbiased, and it is shown explicitly that the final results are independent of all artificial restraints. Third, the repulsive and dispersive free energy contribution arising from the Lennard-Jones interactions of the ligand with its surrounding (protein and solvent) is calculated using the Weeks-Chandler-Andersen separation. This separation also improves convergence of the FEP/MD calculations. Fourth, to decrease the computational cost, only a small number of atoms in the vicinity of the binding site are simulated explicitly, while all the influence of the remaining atoms is incorporated implicitly using the generalized solvent boundary potential (GSBP) method. With GSBP, the size of the simulated FKBP12/ligand systems is significantly reduced, from approximately 25,000 to 2500. The computations are very efficient and the statistical error is small ( approximately 1 kcal/mol). The calculated binding free energies are generally in good agreement with available experimental data and previous calculations (within approximately 2 kcal/mol). The present results indicate that a strategy based on FEP/MD simulations of a reduced GSBP atomic model sampled with conformational, translational, and orientational restraining potentials can be computationally inexpensive and accurate.

Project description:Quinine's ability to bind DNA and potentially inhibit transcription and translation has been examined as a mode of action for its antimalarial activity. UV absorption and fluorescence-based studies have lacked the chemical specificity to develop an unambiguous molecular-level picture of the binding interaction. To address this, we use Raman spectroscopy and molecular dynamics (MD) to investigate quinine-DNA interactions. We demonstrate that quinine's strongest Raman band in the fingerprint region, which derives from a symmetric stretching mode of the quinoline ring, is highly sensitive to the local chemical environment and pH. The frequency shifts observed for this mode in solvents of varying polarity can be explained in terms of the Stark effect using a simple Onsager solvation model, indicating that the vibration reports on the local electrostatic environment. However, specific chemical interactions between the quinoline ring and its environment, such as hydrogen bonding and ?-stacking, perturb the frequency of this mode in a more complicated but predictable manner. We use this vibration as a spectroscopic probe to investigate the binding interaction between quinine and DNA. We find that, when the quinoline ring is protonated, quinine weakly intercalates into DNA by forming ?-stacking interactions with the base pairs. The Raman spectra indicate that quinine can intercalate into DNA with a ratio reaching up to roughly one molecule per 25 base pairs. Our results are confirmed by MD simulations, which also show that the quinoline ring adopts a t-shaped ?-stacking geometry with the DNA base pairs, whereas the quinuclidine head group weakly interacts with the phosphate backbone in the minor groove. We expect that the spectral correlations determined here will enable future studies to probe quinine's antimalarial activities, such as disrupting hemozoin biocrystallization, which is hypothesized to be, among other things, one of its primary modes of action against Plasmodium parasites.

Dataset Information

Relative Binding Enthalpies from Molecular Dynamics Simulations Using a Direct Method.

Publications

Relative Binding Enthalpies from Molecular Dynamics Simulations Using a Direct Method.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure