On the Use of Interaction Entropy and Related Methods to Estimate Binding Entropies.
Ontology highlight
ABSTRACT:
SUBMITTER: Ekberg V
PROVIDER: S-EPMC8389774 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:In the machine learning literature we can find numerous methods to solve classification problems. We propose two new performance measures to analyze such methods. These measures are defined by using the concept of proportional reduction of classification error with respect to three benchmark classifiers, the random and two intuitive classifiers which are based on how a non-expert person could realize classification simply by applying a frequentist approach. We show that these three simple methods are closely related to different aspects of the entropy of the dataset. Therefore, these measures account somewhat for entropy in the dataset when evaluating the performance of classifiers. This allows us to measure the improvement in the classification results compared to simple methods, and at the same time how entropy affects classification capacity. To illustrate how these new performance measures can be used to analyze classifiers taking into account the entropy of the dataset, we carry out an intensive experiment in which we use the well-known J48 algorithm, and a UCI repository dataset on which we have previously selected a subset of the most relevant attributes. Then we carry out an extensive experiment in which we consider four heuristic classifiers, and 11 datasets.
| S-EPMC8306704 | biostudies-literature
Project description:The hypothetical scanning molecular dynamics (HSMD) method is used here for calculating the absolute free energy of binding, ΔA(0) of the complex of the protein FKBP12 with the ligand SB2 (also denoted L8) - a system that has been studied previously for comparing the performance of different methods. Our preliminary study suggests that considering long-range electrostatics is imperative even for a hydrophobic ligand such as L8. Therefore the system is modeled by the AMBER force field using Particle Mesh Ewald (PME). HSMD consists of three stages applied to both the ligand-solvent and ligand-protein systems. (1) A small set of system configurations (frames) is extracted from an MD trajectory. (2) The entropy of the ligand in each frame is calculated by a reconstruction procedure. (3) The contribution of water and protein to ΔA(0) is calculated for each frame by gradually increasing the ligand-environment interactions from zero to their full value using thermodynamic integration (TI). Unlike the conventional methods, the structure of the ligand is kept fixed during TI, and HSMD is thus free from the end-point problem encountered with the double annihilation method (DAM); therefore, the need for applying restraints is avoided. Furthermore, unlike the conventional methods, the entropy of the ligand and water is obtained directly as a byproduct of the simulation. In this paper, in addition to the difference in the internal entropies of the ligand in the two environments, we calculate for the first time the external entropy of the ligand, which provides a measure for the size of the active site. We obtain ΔA(0) = -10.7 ±1.0 as compared to the experimental values -10.9 and -10.6 kcal/mol. However, a protein/water system treated by periodic boundary conditions grows significantly with increasing protein size and the computation of ΔA(0) would become expensive by all methods. Therefore, we also apply HSMD to FKBP12-L8 described by the GSBP/SSBP model of Roux's group (implemented in the software CHARMM) where only part of the protein and water around the active site are considered and long-range electrostatic effects are taken into account. For comparison this model was also treated by the double decoupling method (DDM). The two methods have led to comparable results for ΔA(0) which are somewhat lower than the experimental value. The ligand was found to be more confined in the active site described by GSBP/SSBP than by PME where its entropy in solvent is larger than in the active site by 1.7 and by 5.5 kcal/mol, respectively.
| S-EPMC3274823 | biostudies-literature
Project description:In our recent efforts to map protein surfaces using mixed-solvent molecular dynamics (MixMD) (Ghanakota, P.; Carlson, H. A. Moving Beyond Active-Site Detection: MixMD Applied to Allosteric Systems. J. Phys. Chem. B 2016, 120, 8685-8695), we were able to successfully capture active sites and allosteric sites within the top-four most occupied hotspots. In this study, we describe our approach for estimating the thermodynamic profile of the binding sites identified by MixMD. First, we establish a framework for calculating free energies from MixMD simulations, and we compare our approach to alternative methods. Second, we present a means to obtain a relative ranking of the binding sites by their configurational entropy. The theoretical maximum and minimum free energy and entropy values achievable under such a framework along with the limitations of the techniques are discussed. Using this approach, the free energy and relative entropy ranking of the top-four MixMD binding sites were computed and analyzed across our allosteric protein targets: Abl Kinase, Androgen Receptor, Pdk1 Kinase, Farnesyl Pyrophosphate Synthase, Chk1 Kinase, Glucokinase, and Protein Tyrosine Phosphatase 1B.
| S-EPMC6549510 | biostudies-literature
Project description:Investigators have examined whether heat mortality risk is increased in neighborhoods subject to the urban heat island (UHI) effect but have not identified degrees of difference in susceptibility to heat and cold between cool and hot areas, which we call acclimatization to the UHI.We developed methods to examine and quantify the degree of acclimatization to heat- and cold-related mortality in relation to UHI anomalies and applied these methods to London, UK.Case-crossover analyses were undertaken on 1993-2006 mortality data from London UHI decile groups defined by anomalies from the London average of modeled air temperature at a 1-km grid resolution. We estimated how UHI anomalies modified excess mortality on cold and hot days for London overall and displaced a fixed-shape temperature-mortality function ("shifted spline" model). We also compared the observed associations with those expected under no or full acclimatization to the UHI.The relative risk of death on hot versus normal days differed very little across UHI decile groups. A 1°C UHI anomaly multiplied the risk of heat death by 1.004 (95% CI: 0.950, 1.061) (interaction rate ratio) compared with the expected value of 1.070 (1.057, 1.082) if there were no acclimatization. The corresponding UHI interaction for cold was 1.020 (0.979, 1.063) versus 1.030 (1.026, 1.034) (actual versus expected under no acclimatization, respectively). Fitted splines for heat shifted little across UHI decile groups, again suggesting acclimatization. For cold, the splines shifted somewhat in the direction of no acclimatization, but did not exclude acclimatization.We have proposed two analytical methods for estimating the degree of acclimatization to the heat- and cold-related mortality burdens associated with UHIs. The results for London suggest relatively complete acclimatization to the UHI effect on summer heat-related mortality, but less clear evidence for cold-related mortality.Milojevic A, Armstrong BG, Gasparrini A, Bohnenstengel SI, Barratt B, Wilkinson P. 2016. Methods to estimate acclimatization to urban heat island effects on heat- and cold-related mortality. Environ Health Perspect 124:1016-1022;?http://dx.doi.org/10.1289/ehp.1510109.
| S-EPMC4937865 | biostudies-other
Project description:Cross-entropy was introduced in 1996 to quantify the degree of asynchronism between two time series. In 2009, a multiscale cross-entropy measure was proposed to analyze the dynamical characteristics of the coupling behavior between two sequences on multiple scales. Since their introductions, many improvements and other methods have been developed. In this review we offer a state-of-the-art on cross-entropy measures and their multiscale approaches.
| S-EPMC7516475 | biostudies-literature
Project description:Molecular dynamics (MD) simulation in the explicit water is performed to study the interaction mechanism of trypsin-ligand binding under the AMBER force field and polarized protein-specific charge (PPC) force field combined the new developed highly efficient interaction entropy (IE) method for calculation of entropy change. And the detailed analysis and comparison of the results of MD simulation for two trypsin-ligand systems show that the root-mean-square deviation (RMSD) of backbone atoms, B-factor, intra-protein and protein-ligand hydrogen bonds are more stable under PPC force field than AMBER force field. Our results demonstrate that the IE method is superior than the traditional normal mode (Nmode) method in the calculation of entropy change and the calculated binding free energy under the PPC force field combined with the IE method is more close to the experimental value than other three combinations (AMBER-Nmode, AMBER-IE and PPC-Nmode). And three critical hydrogen bonds between trypsin and ligand are broken under AMBER force field. However, they are well preserved under PPC force field. Detailed binding interactions of ligands with trypsin are further analyzed. The present work demonstrates that the polarized force field combined the highly efficient IE method is critical in MD simulation and free energy calculation.
| S-EPMC5735144 | biostudies-literature
Project description:The study of the p53-MDMX/MDM2 binding sites is a research hotspot for tumor drug design. The inhibition of p53-targeted MDMX/MDM2 has become an effective approach in anti-tumor drug development. In this paper, a theoretically rigorous and computationally accurate method, namely, the interaction entropy (IE) method, combined with the polarized protein-specific charge (PPC) force field, is used to explore the difference in the binding mechanism between p53-MDMX and p53-MDM2. The interaction of a 12mer peptide inhibitor (pDIQ), which is similar to p53 in structure, with MDMX/MDM2 is also studied. The results demonstrate that p53/pDIQ with MDM2 generates a stronger interaction than with MDMX. Compared to p53, pDIQ has larger binding free energies with MDMX and MDM2. According to the calculated binding free energies, the differences in the binding free energy among the four complexes that are obtained from the combination of PPC and IE are more consistent with the experimental values than with the results from the combination of the non-polarizable AMBER force field and IE. In addition, according to the decomposition of the binding free energy, the van der Waals (vdW) interactions are the main driving force for the binding of the four complexes. They are also the main source of the weaker binding affinity of p53/pDIQ-MDMX relative to p53/pDIQ-MDM2. Compared with p53-MDMX/MDM2, according to the analysis of the residue decomposition, the predicated total residue contributions are higher in pDIQ-MDMX/MDM2 than in p53-MDMX/MDM2, which explains why pDIQ has higher binding affinity than p53 with MDMX/MDM2. The current study provides theoretical guidance for understanding the binding mechanisms and designing a potent dual inhibitor that is targeted to MDMX/MDM2.
| S-EPMC6361799 | biostudies-literature
Project description:When constructing discrete (binned) distributions from samples of a data set, applications exist where it is desirable to assure that all bins of the sample distribution have nonzero probability. For example, if the sample distribution is part of a predictive model for which we require returning a response for the entire codomain, or if we use Kullback-Leibler divergence to measure the (dis-)agreement of the sample distribution and the original distribution of the variable, which, in the described case, is inconveniently infinite. Several sample-based distribution estimators exist which assure nonzero bin probability, such as adding one counter to each zero-probability bin of the sample histogram, adding a small probability to the sample pdf, smoothing methods such as Kernel-density smoothing, or Bayesian approaches based on the Dirichlet and Multinomial distribution. Here, we suggest and test an approach based on the Clopper-Pearson method, which makes use of the binominal distribution. Based on the sample distribution, confidence intervals for bin-occupation probability are calculated. The mean of each confidence interval is a strictly positive estimator of the true bin-occupation probability and is convergent with increasing sample size. For small samples, it converges towards a uniform distribution, i.e., the method effectively applies a maximum entropy approach. We apply this nonzero method and four alternative sample-based distribution estimators to a range of typical distributions (uniform, Dirac, normal, multimodal, and irregular) and measure the effect with Kullback-Leibler divergence. While the performance of each method strongly depends on the distribution type it is applied to, on average, and especially for small sample sizes, the nonzero, the simple "add one counter", and the Bayesian Dirichlet-multinomial model show very similar behavior and perform best. We conclude that, when estimating distributions without an a priori idea of their shape, applying one of these methods is favorable.
| S-EPMC7513126 | biostudies-literature
Project description:Multiple methods have been developed to estimate narrow-sense heritability, h2, using single nucleotide polymorphisms (SNPs) in unrelated individuals. However, a comprehensive evaluation of these methods has not yet been performed, leading to confusion and discrepancy in the literature. We present the most thorough and realistic comparison of these methods to date. We used thousands of real whole-genome sequences to simulate phenotypes under varying genetic architectures and confounding variables, and we used array, imputed, or whole genome sequence SNPs to obtain 'SNP-heritability' estimates. We show that SNP-heritability can be highly sensitive to assumptions about the frequencies, effect sizes, and levels of linkage disequilibrium of underlying causal variants, but that methods that bin SNPs according to minor allele frequency and linkage disequilibrium are less sensitive to these assumptions across a wide range of genetic architectures and possible confounding factors. These findings provide guidance for best practices and proper interpretation of published estimates.
| S-EPMC5934350 | biostudies-literature
Project description:The restriction of a small molecule's motion on binding to a protein causes a loss of configurational entropy, and thus a penalty in binding affinity. Some energy models used in computer-aided ligand design neglect this entropic penalty, whereas others account for it based on an expected drop in the number of accessible rotamers upon binding. However, the validity of the physical assumptions underlying the various approaches is largely unexamined. The present study addresses this issue by using Mining Minima calculations to analyze the association of amprenavir with HIV protease. The computed loss in ligand configurational entropy is large, contributing approximately 25 kcal/mol (4.184 kJ/kcal) to DeltaG degrees. Most of this loss results from narrower energy wells in the bound state, rather than a drop in the number of accessible rotamers. Coupling among rotation/translation and internal degrees of freedom complicates the decomposition of the entropy change into additive terms. The results highlight the potential to gain affinity by designing conformationally restricted ligands and have implications for the formulation of energy models for ligand scoring.
| S-EPMC1780070 | biostudies-literature