Dataset Information

Self-Consistent Framework Connecting Experimental Proxies of Protein Dynamics with Configurational Entropy.

ABSTRACT: The recently developed NMR techniques enable estimation of protein configurational entropy change from the change in the average methyl order parameters. This experimental observable, however, does not directly measure the contribution of intramolecular couplings, protein main-chain motions, or angular dynamics. Here, we carry out a self-consistent computational analysis of the impact of these missing contributions on an extensive set of molecular dynamics simulations of different proteins undergoing binding. Specifically, we compare the configurational entropy change in protein complex formation as obtained by the maximum information spanning tree approximation (MIST), which treats the above entropy contributions directly, and the change in the average NMR methyl and NH order parameters. Our parallel implementation of MIST allows us to treat hard angular degrees of freedom as well as couplings up to full pairwise order explicitly, while still involving a high degree of sampling and tackling molecules of biologically relevant sizes. First, we demonstrate a remarkably strong linear relationship between the total configurational entropy change and the average change in both methyl and backbone-NH order parameters. Second, in contrast to canonical assumptions, we show that the main-chain and angular terms contribute significantly to the overall configurational entropy change and also scale linearly with it. Consequently, linear models starting from the average methyl order parameters are able to capture the contribution of main-chain and angular terms well. After applying the quantum-mechanical harmonic oscillator entropy formalism, we establish a similarly strong linear relationship for X-ray crystallographic B-factors. Finally, we demonstrate that the observed linear relationships remain robust against drastic undersampling and argue that they reflect an intrinsic property of compact proteins. Despite their remarkable strength, however, the above linear relationships yield estimates of configurational entropy change whose accuracy appears to be sufficient for qualitative applications only.

SUBMITTER: Fleck M

PROVIDER: S-EPMC9245193 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Entropy and the second law of thermodynamics are fundamental concepts that underlie all natural processes and patterns. Recent research has shown how the entropy of a landscape mosaic can be calculated using the Boltzmann equation, with the entropy of a lattice mosaic equal to the logarithm of the number of ways a lattice with a given dimensionality and number of classes can be arranged to produce the same total amount of edge between cells of different classes. However, that work seemed to also suggest that the feasibility of applying this method to real landscapes was limited due to intractably large numbers of possible arrangements of raster cells in large landscapes. Here I extend that work by showing that: (1) the proportion of arrangements rather than the number with a given amount of edge length provides a means to calculate unbiased relative configurational entropy, obviating the need to compute all possible configurations of a landscape lattice; (2) the edge lengths of randomized landscape mosaics are normally distributed, following the central limit theorem; and (3) given this normal distribution it is possible to fit parametric probability density functions to estimate the expected proportion of randomized configurations that have any given edge length, enabling the calculation of configurational entropy on any landscape regardless of size or number of classes. I evaluate the boundary limits (4) for this normal approximation for small landscapes with a small proportion of a minority class and show it holds under all realistic landscape conditions. I further (5) demonstrate that this relationship holds for a sample of real landscapes that vary in size, patch richness, and evenness of area in each cover type, and (6) I show that the mean and standard deviation of the normally distributed edge lengths can be predicted nearly perfectly as a function of the size, patch richness and diversity of a landscape. Finally, (7) I show that the configurational entropy of a landscape is highly related to the dimensionality of the landscape, the number of cover classes, the evenness of landscape composition across classes, and landscape heterogeneity. These advances provide a means for researchers to directly estimate the frequency distribution of all possible macrostates of any observed landscape, and then directly calculate the relative configurational entropy of the observed macrostate, and to understand the ecological meaning of different amounts of configurational entropy. These advances enable scientists to take configurational entropy from a concept to an applied tool to measure and compare the disorder of real landscapes with an objective and unbiased measure based on entropy and the second law.

Project description:Molecular recognition by enzymes is a complicated process involving thermodynamic energies governing protein-ligand interactions. In order to aid the estimation of inhibitory activity of compounds targeting an enzyme, several computational methods can be employed to dissect this intermolecular contact. Herein, we report a structural dynamics investigation of an epigenetic enzyme HDAC2 in differentiating its binding to various inhibitors within the sub-sites of its active site. Molecular dynamics (MD) simulation was employed to elucidate the intermolecular interactions as well as the dynamics behavior of ligand binding. MD trajectories of five distinct HDAC2-inhibitor complexes reveal that compounds lacking adequate contacts with the opening rim of the active site possess high fluctuation along the cap portion, thus weakening the overall affinity. Key intermolecular interactions determining the effective binding of inhibitors include hydrogen bonds with Gly154, Asp181, and Tyr308; hydrophobic interactions between Phe155/Phe210 and the linker region; and a pi-stacking with Arg39 at the foot pocket. Decomposition of the binding free energy calculated per-residue by MM/PBSA also indicates that the interactions within the internal foot pocket, especially with residues Met35, Leu144, Gly305, and Gly306, can contribute significantly to the ligand binding. Additionally, configurational entropy of the binding was estimated and compared to the scale of the binding free energy in order to assess its contribution to the binding and to differentiate various ligand partners. It was found that the levels of entropic contribution are comparable among a set of structurally similar carbamide ligands, while it is greatly different for the set of unrelated ligands, ranging from 2.75 to 16.38 kcal/mol for the five inhibitors examined. These findings exemplify the importance of assessing molecular dynamics as well as estimating the entropic contribution in evaluating the ligand binding mechanism.

Dataset Information

Self-Consistent Framework Connecting Experimental Proxies of Protein Dynamics with Configurational Entropy.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets