Browse
Submit Data
Databases
API
Help

Dataset Information

14 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Ensembles Are Required to Handle Aleatoric and Parametric Uncertainty in Molecular Dynamics Simulation.

ABSTRACT: Classical molecular dynamics is a computer simulation technique that is in widespread use across many areas of science, from physics and chemistry to materials, biology, and medicine. The method continues to attract criticism due its oft-reported lack of reproducibility which is in part due to a failure to submit it to reliable uncertainty quantification (UQ). Here we show that the uncertainty arises from a combination of (i) the input parameters and (ii) the intrinsic stochasticity of the method controlled by the random seeds. To illustrate the situation, we make a systematic UQ analysis of a widely used molecular dynamics code (NAMD), applied to estimate binding free energy of a ligand-bound to a protein. In particular, we replace the usually fixed input parameters with random variables, systematically distributed about their mean values, and study the resulting distribution of the simulation output. We also perform a sensitivity analysis, which reveals that, out of a total of 175 parameters, just six dominate the variance in the code output. Furthermore, we show that binding energy calculations dampen the input uncertainty, in the sense that the variation around the mean output free energy is less than the variation around the mean of the assumed input distributions, if the output is ensemble-averaged over the random seeds. Without such ensemble averaging, the predicted free energy is five times more uncertain. The distribution of the predicted properties is thus strongly dependent upon the random seed. Owing to this substantial uncertainty, robust statistical measures of uncertainty in molecular dynamics simulation require the use of ensembles in all contexts.

SUBMITTER: Vassaux M

PROVIDER: S-EPMC8389531 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Similar Datasets

Validating solution ensembles from molecular dynamics simulation by wide-angle X-ray scattering data.

Project description:Wide-angle x-ray scattering (WAXS) experiments of biomolecules in solution have become increasingly popular because of technical advances in light sources and detectors. However, the structural interpretation of WAXS profiles is problematic, partly because accurate calculations of WAXS profiles from structural models have remained challenging. In this work, we present the calculation of WAXS profiles from explicit-solvent molecular dynamics (MD) simulations of five different proteins. Using only a single fitting parameter that accounts for experimental uncertainties because of the buffer subtraction and dark currents, we find excellent agreement to experimental profiles both at small and wide angles. Because explicit solvation eliminates free parameters associated with the solvation layer or the excluded solvent, which would require fitting to experimental data, we minimize the risk of overfitting. We further find that the influence from water models and protein force fields on calculated profiles are insignificant up to q≈15nm(-1). Using a series of simulations that allow increasing flexibility of the proteins, we show that incorporating thermal fluctuations into the calculations significantly improves agreement with experimental data, demonstrating the importance of protein dynamics in the interpretation of WAXS profiles. In addition, free MD simulations up to one microsecond suggest that the calculated profiles are highly sensitive with respect to minor conformational rearrangements of proteins, such as an increased flexibility of a loop or an increase of the radius of gyration by < 1%. The present study suggests that quantitative comparison between MD simulations and experimental WAXS profiles emerges as an accurate tool to validate solution ensembles of biomolecules.

| S-EPMC4104061 | biostudies-literature

Structural insights for designed alanine-rich helices: comparing NMR helicity measures and conformational ensembles from molecular dynamics simulation.

Project description:The temperature dependence of helical propensities for the peptides Ac-ZGG-(KAAAA)(3)X-NH(2) (Z = Y or G, X = A, K, and D-Arg) were studied both experimentally and by MD simulations. Good agreement is observed in both the absolute helical propensities as well as relative helical content along the sequence; the global minimum on the calculated free energy landscape corresponds to a single alpha-helical conformation running from K4 to A18 with some terminal fraying, particularly at the C-terminus. Energy component analysis shows that the single helix state has favorable intramolecular electrostatic energy due to hydrogen bonds, and that less-favorable two-helix globular states have favorable solvation energy. The central lysine residues do not appear to increase helicity; however, both experimental and simulation studies show increasing helicity in the series X = Ala --> Lys --> D-Arg. This C-capping preference was also experimentally confirmed in Ac-(KAAAA)(3)X-GY-NH(2) and (KAAAA)(3)X-GY-NH(2) sequences. The roles of the C-capping groups, and of lysines throughout the sequence, in the MD-derived ensembles are analyzed in detail.

| S-EPMC3816200 | biostudies-literature

Solution structure of HIV-1 protease flaps probed by comparison of molecular dynamics simulation ensembles and EPR experiments.

Project description:The introduction of multidrug treatment regimens has dramatically prolonged the progression and survival of AIDS patients. However, the success of the long-term treatment has been hindered by strains of HIV that are increasingly resistant to inhibitors of targets such as HIV protease (HIV PR). Therefore, the need for a thorough understanding of the structure and dynamics of HIV PR and how these are altered in resistant mutants is crucial for the design of more effective treatments. Crystal structures of unbound HIV PR show significant heterogeneity and often have extensive crystal packing interactions. Recent site-directed spin labeling (SDSL) and double electron-electron resonance (DEER) spectroscopy studies characterized flap conformations in HIV-1 protease in an inhibited and uninhibited form and distinguished the extent of flap opening in an unbound form. However, the correlation between EPR-measured interspin distances and structural/dynamic features of the flaps has not been established. In this report, we link EPR-based data and 900 ns of MD simulation in explicit water to gain insight into the ensemble of conformations sampled by HIV PR flaps in solution, both in the presence and in the absence of an FDA-approved HIV PR inhibitor.

| S-EPMC3390170 | biostudies-literature

Deconvoluting Protein (Un)folding Structural Ensembles Using X-Ray Scattering, Nuclear Magnetic Resonance Spectroscopy and Molecular Dynamics Simulation.

Project description:The folding and unfolding of protein domains is an apparently cooperative process, but transient intermediates have been detected in some cases. Such (un)folding intermediates are challenging to investigate structurally as they are typically not long-lived and their role in the (un)folding reaction has often been questioned. One of the most well studied (un)folding pathways is that of Drosophila melanogaster Engrailed homeodomain (EnHD): this 61-residue protein forms a three helix bundle in the native state and folds via a helical intermediate. Here we used molecular dynamics simulations to derive sample conformations of EnHD in the native, intermediate, and unfolded states and selected the relevant structural clusters by comparing to small/wide angle X-ray scattering data at four different temperatures. The results are corroborated using residual dipolar couplings determined by NMR spectroscopy. Our results agree well with the previously proposed (un)folding pathway. However, they also suggest that the fully unfolded state is present at a low fraction throughout the investigated temperature interval, and that the (un)folding intermediate is highly populated at the thermal midpoint in line with the view that this intermediate can be regarded to be the denatured state under physiological conditions. Further, the combination of ensemble structural techniques with MD allows for determination of structures and populations of multiple interconverting structures in solution.

| S-EPMC4422743 | biostudies-literature

Bayesian regression explains how human participants handle parameter uncertainty.

Project description:Accumulating evidence indicates that the human brain copes with sensory uncertainty in accordance with Bayes' rule. However, it is unknown how humans make predictions when the generative model of the task at hand is described by uncertain parameters. Here, we tested whether and how humans take parameter uncertainty into account in a regression task. Participants extrapolated a parabola from a limited number of noisy points, shown on a computer screen. The quadratic parameter was drawn from a bimodal prior distribution. We tested whether human observers take full advantage of the given information, including the likelihood of the quadratic parameter value given the observed points and the quadratic parameter's prior distribution. We compared human performance with Bayesian regression, which is the (Bayes) optimal solution to this problem, and three sub-optimal models, which are simpler to compute. Our results show that, under our specific experimental conditions, humans behave in a way that is consistent with Bayesian regression. Moreover, our results support the hypothesis that humans generate responses in a manner consistent with probability matching rather than Bayesian decision theory.

| S-EPMC7259793 | biostudies-literature

Uncertainty in ensembles of global biodiversity scenarios.

Project description:While there is a clear demand for scenarios that provide alternative states in biodiversity with respect to future emissions, a thorough analysis and communication of the associated uncertainties is still missing. Here, we modelled the global distribution of ~11,500 amphibian, bird and mammal species and project their climatic suitability into the time horizon 2050 and 2070, while varying the input data used. By this, we explore the uncertainties originating from selecting species distribution models (SDMs), dispersal strategies, global circulation models (GCMs), and representative concentration pathways (RCPs). We demonstrate the overwhelming influence of SDMs and RCPs on future biodiversity projections, followed by dispersal strategies and GCMs. The relative importance of each component varies in space but also with the selected sensitivity metrics and with species' range size. Overall, this means using multiple SDMs, RCPs, dispersal assumptions and GCMs is a necessity in any biodiversity scenario assessment, to explicitly report associated uncertainties.

| S-EPMC6441032 | 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

Molecular nonlinear dynamics and protein thermal uncertainty quantification.

Project description:This work introduces molecular nonlinear dynamics (MND) as a new approach for describing protein folding and aggregation. By using a mode system, we show that the MND of disordered proteins is chaotic while that of folded proteins exhibits intrinsically low dimensional manifolds (ILDMs). The stability of ILDMs is found to strongly correlate with protein energies. We propose a novel method for protein thermal uncertainty quantification based on persistently invariant ILDMs. Extensive comparison with experimental data and the state-of-the-art methods in the field validate the proposed new method for protein B-factor prediction.

| S-EPMC3899061 | biostudies-literature

Mapping conformational ensembles of a? oligomers in molecular dynamics simulations.

Project description:Although the oligomers formed by A? peptides appear to be the primary cytotoxic species in Alzheimer's disease, detailed information about their structures appears to be lacking. In this article, we use exhaustive replica exchange molecular dynamics and an implicit solvent united-atom model to study the structural properties of A? monomers, dimers, and tetramers. Our analysis suggests that the conformational ensembles of A? dimers and tetramers are very similar, but sharply distinct from those sampled by the monomers. The key conformational difference between monomers and oligomers is the formation of ?-structure in the oligomers occurring together with the loss of intrapeptide interactions and helix structure. Our simulations indicate that, independent of oligomer order, the A? aggregation interface is largely confined to the sequence region 10-23, which forms the bulk of interpeptide interactions. We show that the fractions of ? structure computed in our simulations and measured experimentally are in good agreement.

| S-EPMC2941018 | biostudies-literature

Semi-parametric Regression under Model Uncertainty: Economic Applications.

Project description:Economic theory does not always specify the functional relationship between dependent and explanatory variables, or even isolate a particular set of covariates. This means that model uncertainty is pervasive in empirical economics. In this paper, we indicate how Bayesian semi-parametric regression methods in combination with stochastic search variable selection can be used to address two model uncertainties simultaneously: (i) the uncertainty with respect to the variables which should be included in the model and (ii) the uncertainty with respect to the functional form of their effects. The presented approach enables the simultaneous identification of robust linear and nonlinear effects. The additional insights gained are illustrated on applications in empirical economics, namely willingness to pay for housing, and cross-country growth regression.

| S-EPMC6850473 | 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