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.

Project description:Calorimetric studies of protein-ligand binding sometimes yield thermodynamic data that are difficult to understand. Today, molecular simulations can be used to seek insight into such calorimetric puzzles, and, when simulations and experiments diverge, the results can usefully motivate further improvements in computational methods. Here, we apply near-millisecond duration simulations to estimate the relative binding enthalpies of four peptidic ligands with the Grb2 SH2 domain. The ligands fall into matched pairs, where one member of each pair has an added bond that preorganizes the ligand for binding and thus may be expected to favor binding entropically, due to a smaller loss in configurational entropy. Calorimetric studies have shown that the constrained ligands do in fact bind the SH2 domain more tightly than the flexible ones, but, paradoxically, the improvement in affinity for the constrained ligands is enthalpic, rather than entropic. The present enthalpy calculations yield the opposite trend, as they suggest that the flexible ligands bind more exothermically. Additionally, the small relative binding enthalpies are found to be balances of large differences in the energies of structural components such as ligand and the binding site residues. As a consequence, the deviations from experiment in the relative binding enthalpies represent small differences between these large numbers and hence may be particularly susceptible to error, due, for example, to approximations in the force field. We also computed first-order estimates of changes in configurational entropy on binding. These too are, arguably, paradoxical, as they tend to favor binding of the flexible ligands. The paradox is explained in part by the fact that the more rigid constrained ligands reduce the entropy of binding site residues more than their flexible analogs do, at least in the simulations. This result offers a rather general counterargument to the expectation that preorganized ligands should be associated with more favorable binding entropies, other things being equal.

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

Tweets