Browse
Submit Data
Databases
API
Help

Dataset Information

27 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Statistical Measures to Quantify Similarity between Molecular Dynamics Simulation Trajectories.

ABSTRACT: Molecular dynamics simulation is commonly employed to explore protein dynamics. Despite the disparate timescales between functional mechanisms and molecular dynamics (MD) trajectories, functional differences are often inferred from differences in conformational ensembles between two proteins in structure-function studies that investigate the effect of mutations. A common measure to quantify differences in dynamics is the root mean square fluctuation (RMSF) about the average position of residues defined by C?-atoms. Using six MD trajectories describing three native/mutant pairs of beta-lactamase, we make comparisons with additional measures that include Jensen-Shannon, modifications of Kullback-Leibler divergence, and local p-values from 1-sample Kolmogorov-Smirnov tests. These additional measures require knowing a probability density function, which we estimate by using a nonparametric maximum entropy method that quantifies rare events well. The same measures are applied to distance fluctuations between C?-atom pairs. Results from several implementations for quantitative comparison of a pair of MD trajectories are made based on fluctuations for on-residue and residue-residue local dynamics. We conclude that there is almost always a statistically significant difference between pairs of 100 ns all-atom simulations on moderate-sized proteins as evident from extraordinarily low p-values.

SUBMITTER: Farmer J

PROVIDER: S-EPMC6258182 | biostudies-literature | 2017 Dec

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Similar Datasets

Learning probabilistic models of hydrogen bond stability from molecular dynamics simulation trajectories.

Project description:BackgroundHydrogen bonds (H-bonds) play a key role in both the formation and stabilization of protein structures. They form and break while a protein deforms, for instance during the transition from a non-functional to a functional state. The intrinsic strength of an individual H-bond has been studied from an energetic viewpoint, but energy alone may not be a very good predictor.MethodsThis paper describes inductive learning methods to train protein-independent probabilistic models of H-bond stability from molecular dynamics (MD) simulation trajectories of various proteins. The training data contains 32 input attributes (predictors) that describe an H-bond and its local environment in a conformation c and the output attribute is the probability that the H-bond will be present in an arbitrary conformation of this protein achievable from c within a time duration Δ. We model dependence of the output variable on the predictors by a regression tree.ResultsSeveral models are built using 6 MD simulation trajectories containing over 4000 distinct H-bonds (millions of occurrences). Experimental results demonstrate that such models can predict H-bond stability quite well. They perform roughly 20% better than models based on H-bond energy alone. In addition, they can accurately identify a large fraction of the least stable H-bonds in a conformation. In most tests, about 80% of the 10% H-bonds predicted as the least stable are actually among the 10% truly least stable. The important attributes identified during the tree construction are consistent with previous findings.ConclusionsWe use inductive learning methods to build protein-independent probabilistic models to study H-bond stability, and demonstrate that the models perform better than H-bond energy alone.

| S-EPMC3044290 | biostudies-literature

Simulation of nitroxide electron paramagnetic resonance spectra from brownian trajectories and molecular dynamics simulations.

Project description:A simulated continuous wave electron paramagnetic resonance spectrum of a nitroxide spin label can be obtained from the Fourier transform of a free induction decay. It has been previously shown that the free induction decay can be calculated by solving the time-dependent stochastic Liouville equation for a set of Brownian trajectories defining the rotational dynamics of the label. In this work, a quaternion-based Monte Carlo algorithm has been developed to generate Brownian trajectories describing the global rotational diffusion of a spin-labeled protein. Also, molecular dynamics simulations of two spin-labeled mutants of T4 lysozyme, T4L F153R1, and T4L K65R1 have been used to generate trajectories describing the internal dynamics of the protein and the local dynamics of the spin-label side chain. Trajectories from the molecular dynamics simulations combined with trajectories describing the global rotational diffusion of the protein are used to account for all of the dynamics of a spin-labeled protein. Spectra calculated from these combined trajectories correspond well to the experimental spectra for the buried site T4L F153R1 and the helix surface site T4L K65R1. This work provides a framework to further explore the modeling of the dynamics of the spin-label side chain in the wide variety of labeling environments encountered in site-directed spin labeling studies.

| S-EPMC2367180 | biostudies-literature

Dynamics based alignment of proteins: an alternative approach to quantify dynamic similarity.

Project description:BACKGROUND:The dynamic motions of many proteins are central to their function. It therefore follows that the dynamic requirements of a protein are evolutionary constrained. In order to assess and quantify this, one needs to compare the dynamic motions of different proteins. Comparing the dynamics of distinct proteins may also provide insight into how protein motions are modified by variations in sequence and, consequently, by structure. The optimal way of comparing complex molecular motions is, however, far from trivial. The majority of comparative molecular dynamics studies performed to date relied upon prior sequence or structural alignment to define which residues were equivalent in 3-dimensional space. RESULTS:Here we discuss an alternative methodology for comparative molecular dynamics that does not require any prior alignment information. We show it is possible to align proteins based solely on their dynamics and that we can use these dynamics-based alignments to quantify the dynamic similarity of proteins. Our method was tested on 10 representative members of the PDZ domain family. CONCLUSIONS:As a result of creating pair-wise dynamics-based alignments of PDZ domains, we have found evolutionarily conserved patterns in their backbone dynamics. The dynamic similarity of PDZ domains is highly correlated with their structural similarity as calculated with Dali. However, significant differences in their dynamics can be detected indicating that sequence has a more refined role to play in protein dynamics than just dictating the overall fold. We suggest that the method should be generally applicable.

| S-EPMC2868010 | biostudies-literature

Interfacial Contact Behavior between CNTs and AgNW with Molecular Dynamics Simulation.

Project description:The behavior at an interface between carbon nanotubes (CNTs) and silver nanowire (AgNW) could hardly be observed experimentally on an atomic scale, and the interaction is difficult to accurately calculate due to nanometer size effects. In this work, the contact behavior is studied with the molecular dynamics (MD) simulation, which indicates that the CNTs and AgNW can move towards each other to form aligned structures with their interfaces in full contact. In these different composite systems, nanotubes may either keep their form of an inherent cylindrical structure or completely collapse into the nanoribbons that can tightly scroll on the AgNW periphery while wrapping it in a core-shell structure. Thus, the atomic configuration evolution that is affected by the van der Waals (vdW) interaction is closely analyzed to assist the understanding of interfacial contact behavior.

| S-EPMC7142578 | biostudies-literature

Project description:We have compared force-induced unfolding with traditional unfolding methods using apomyoglobin as a model protein. Using molecular dynamics simulation, we have investigated the structural stability as a function of the degree of mechanical perturbation. Both anisotropic perturbation by stretching two terminal atoms and isotropic perturbation by increasing the radius of gyration of the protein show the same key event of force-induced unfolding. Our primary results show that the native structure of apomyoglobin becomes destabilized against the mechanical perturbation as soon as the interhelical packing between the G and H helices is broken, suggesting that our simulation results share a common feature with the experimental observation that the interhelical contact is more important for the folding of apomyoglobin than the stability of individual helices. This finding is further confirmed by simulating both helix destabilizing and interhelical packing destabilizing mutants.

| S-EPMC1303326 | biostudies-literature

Molecular Dynamics Simulation for All.

Project description:The impact of molecular dynamics (MD) simulations in molecular biology and drug discovery has expanded dramatically in recent years. These simulations capture the behavior of proteins and other biomolecules in full atomic detail and at very fine temporal resolution. Major improvements in simulation speed, accuracy, and accessibility, together with the proliferation of experimental structural data, have increased the appeal of biomolecular simulation to experimentalists-a trend particularly noticeable in, although certainly not limited to, neuroscience. Simulations have proven valuable in deciphering functional mechanisms of proteins and other biomolecules, in uncovering the structural basis for disease, and in the design and optimization of small molecules, peptides, and proteins. Here we describe, in practical terms, the types of information MD simulations can provide and the ways in which they typically motivate further experimental work.

| S-EPMC6209097 | biostudies-literature

Direct Crystallization Resolution of Racemates Enhanced by Chiral Nanorods: Experimental, Statistical, and Quantum Mechanics/Molecular Dynamics Simulation Studies.

Project description:Three chiral nanorods of C14-l-Thea, C14-l-Phe, and C14-d-Phe were first synthesized and utilized as heterogeneous nucleants to enhance the resolution of racemic Asp via direct crystallization. Through the statistical analysis from 320 batches of nucleation experiments, we found that the apparent appearance diversity of two enantiomeric crystals of Asp existed in 80 homogeneous experiments without chiral nanorods. However, in 240 heterogeneous experiments with 4.0 wt % chiral nanorods of solute mass added, the appearance of those nuclei with the same chirality as the nanorods was apparently promoted, and that with the opposite chirality was totally inhibited. Under a supersaturation level of 1.08, the maximum ee of the initial nuclei was as high as 23.51%. When the cooling rate was 0.025 K/min, the ee of the product was up to 76.85% with a yield of 14.41%. Furthermore, the simulation results from quantum mechanics (QM) and molecular dynamics (MD) revealed that the higher chiral recognition ability of C14-l-Thea compared to C14-l-Phe that originated from the interaction difference between C14-l-Thea and Asp enantiomers was larger than that between C14-l-Phe and Asp enantiomers. Moreover, the constructed nanorods exhibited good stability and recyclability.

| S-EPMC9202296 | biostudies-literature

Interaction between β-lactoglobulin and EGCG under high-pressure by molecular dynamics simulation.

Project description:The binding between β-lactoglobulin and epigallocatechin gallate (EGCG) under the pressure of 600 MPa was explored using molecular docking and molecular dynamics (MD) simulation. EGCG bound mainly in two regions with site 1 in internal cavity of the β-barrel and site 2 on the surface of protein. 150 ns MD was performed starting from the structure with the optimal binding energy at the two sites in molecular docking, respectively. It was found that the protein fluctuated greatly when small molecule bound to site 2 at 0.1 MPa, and the protein fluctuation and solvent accessible surface area became smaller under high-pressure. The binding of small molecules made the protein structure more stable with increasing of α-helix and β-sheet, while high-pressure destroyed α-helix of protein. The binding energy of small molecules at site 1was stronger than that at site 2 under 0.1 MPa, with stronger van der Waals and hydrophobic interaction at site 1 while more hydrogen bonds were present at site 2. The binding energy of both sites weakened under high-pressure, especially at site 1, causing the binding force to be weaker at site 1 than that at site 2 under high-pressure.

| S-EPMC8691620 | biostudies-literature

Interaction between Curcumin and β-Casein: Multi-Spectroscopic and Molecular Dynamics Simulation Methods.

Project description:Effect of temperature and pH on the interaction of curcumin with β-casein was explored by fluorescence spectroscopy, ultraviolet-visible spectroscopy and molecular dynamics simulation. The spectroscopic results showed that curcumin could bind to β-casein to form a complex which was driven mainly by electrostatic interaction. The intrinsic fluorescence of β-casein was quenched by curcumin through static quenching mechanism. The binding constants of curcumin to β-casein were 6.48 × 104 L/mol (298 K), 6.17 × 104 L/mol (305 K) and 5.73 × 104 L/mol (312 K) at pH 2.0, which was greater than that (3.98 × 104 L/mol at 298 K, 3.90 × 104 L/mol at 305 K and 3.41 × 104 L/mol at 312 K) at pH 7.4. Molecular docking study showed that binding energy of β-casein-curcumin complex at pH 2.0 (-7.53 kcal/mol) was lower than that at pH 7.4 (-7.01 kcal/mol). The molecular dynamics simulation study showed that the binding energy (-131.07 kJ/mol) of β-casein-curcumin complex was relatively low at pH 2.0 and 298 K. α-Helix content in β-casein was decreased and random coil content was increased in the presence of curcumin. These results can promote a deep understanding of interaction between curcumin and β-casein and provide a reference for improving the bioavailability of curcumin.

| S-EPMC8398086 | biostudies-literature

Molecular dynamics simulation of the interactions between EHD1 EH domain and multiple peptides.

Project description:ObjectiveTo provide essential information for peptide inhibitor design, the interactions of Eps15 homology domain of Eps15 homology domain-containing protein 1 (EHD1 EH domain) with three peptides containing NPF (asparagine-proline-phenylalanine), DPF (aspartic acid-proline-phenylalanine), and GPF (glycine-proline-phenylalanine) motifs were deciphered at the atomic level. The binding affinities and the underlying structure basis were investigated.MethodsMolecular dynamics (MD) simulations were performed on EHD1 EH domain/peptide complexes for 60 ns using the GROMACS package. The binding free energies were calculated and decomposed by molecular mechanics/generalized Born surface area (MM/GBSA) method using the AMBER package. The alanine scanning was performed to evaluate the binding hot spot residues using FoldX software.ResultsThe different binding affinities for the three peptides were affected dominantly by van der Waals interactions. Intermolecular hydrogen bonds provide the structural basis of contributions of van der Waals interactions of the flanking residues to the binding.Conclusionsvan der Waals interactions should be the main consideration when we design peptide inhibitors of EHD1 EH domain with high affinities. The ability to form intermolecular hydrogen bonds with protein residues can be used as the factor for choosing the flanking residues.

| S-EPMC4609540 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data