On the statistical equivalence of restrained-ensemble simulations with the maximum entropy method.
Ontology highlight
ABSTRACT:
SUBMITTER: Roux B
PROVIDER: S-EPMC3598863 | biostudies-other | 2013 Feb
REPOSITORIES: biostudies-other
Similar Datasets
Project description:Currently, the best existing molecular dynamics (MD) force fields cannot accurately reproduce the global free-energy minimum which realizes the experimental protein structure. As a result, long MD trajectories tend to drift away from the starting coordinates (e.g., crystallographic structures). To address this problem, we have devised a new simulation strategy aimed at protein crystals. An MD simulation of protein crystal is essentially an ensemble simulation involving multiple protein molecules in a crystal unit cell (or a block of unit cells). To ensure that average protein coordinates remain correct during the simulation, we introduced crystallography-based restraints into the MD protocol. Because these restraints are aimed at the ensemble-average structure, they have only minimal impact on conformational dynamics of the individual protein molecules. So long as the average structure remains reasonable, the proteins move in a native-like fashion as dictated by the original force field. To validate this approach, we have used the data from solid-state NMR spectroscopy, which is the orthogonal experimental technique uniquely sensitive to protein local dynamics. The new method has been tested on the well-established model protein, ubiquitin. The ensemble-restrained MD simulations produced lower crystallographic R factors than conventional simulations; they also led to more accurate predictions for crystallographic temperature factors, solid-state chemical shifts, and backbone order parameters. The predictions for (15) N R1 relaxation rates are at least as accurate as those obtained from conventional simulations. Taken together, these results suggest that the presented trajectories may be among the most realistic protein MD simulations ever reported. In this context, the ensemble restraints based on high-resolution crystallographic data can be viewed as protein-specific empirical corrections to the standard force fields.
| S-EPMC3970899 | biostudies-literature
Project description:DEER (double electron-electron resonance) spectroscopy is a powerful pulsed ESR (electron spin resonance) technique allowing the determination of spin-spin distance histograms between site-directed nitroxide label sites on a protein in their native environment. However, incorporating ESR/DEER data in structural refinement is challenging because the information from the large number of distance histograms is complex and highly coupled. Here, a novel restrained-ensemble molecular dynamics simulation method is developed to incorporate the information from multiple ESR/DEER distance histograms simultaneously. Illustrative tests on three coupled spin-labels inserted in T4 lysozyme show that the method efficiently imposes the experimental distance distribution in this system. Different rotameric states of the ?1 and ?2 dihedrals in the spin-labels are also explored by restrained ensemble simulations. Using this method, it is hoped that experimental restraints from ESR/DEER experiments can be used to refine structural properties of biological systems.
| S-EPMC3683991 | biostudies-literature
Project description:DEER (double electron-electron resonance) is a powerful pulsed ESR (electron spin resonance) technique allowing the determination of distance histograms between pairs of nitroxide spin-labels linked to a protein in a native-like solution environment. However, exploiting the huge amount of information provided by ESR/DEER histograms to refine structural models is extremely challenging. In this study, a restrained ensemble (RE) molecular dynamics (MD) simulation methodology is developed to address this issue. In RE simulation, the spin-spin distance distribution histograms calculated from a multiple-copy MD simulation are enforced, via a global ensemble-based energy restraint, to match those obtained from ESR/DEER experiments. The RE simulation is applied to 51 ESR/DEER distance histogram data from spin-labels inserted at 37 different positions in T4 lysozyme (T4L). The rotamer population distribution along the five dihedral angles connecting the nitroxide ring to the protein backbone is determined and shown to be consistent with available information from X-ray crystallography. For the purpose of structural refinement, the concept of a simplified nitroxide dummy spin-label is designed and parametrized on the basis of these all-atom RE simulations with explicit solvent. It is demonstrated that RE simulations with the dummy nitroxide spin-labels imposing the ESR/DEER experimental distance distribution data are able to systematically correct and refine a series of distorted T4L structures, while simple harmonic distance restraints are unsuccessful. This computationally efficient approach allows experimental restraints from DEER experiments to be incorporated into RE simulations for efficient structural refinement.
| S-EPMC3684008 | 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:Maximum entropy (MaxEnt) spectral reconstruction methods provide a powerful framework for spectral estimation of nonuniformly sampled datasets. Many methods exist within this framework, usually defined based on the magnitude of a Lagrange multiplier in the MaxEnt objective function. An algorithm is presented here that utilizes accelerated first-order convex optimization techniques to rapidly and reliably reconstruct nonuniformly sampled NMR datasets using the principle of maximum entropy. This algorithm - called CAMERA for Convex Accelerated Maximum Entropy Reconstruction Algorithm - is a new approach to spectral reconstruction that exhibits fast, tunable convergence in both constant-aim and constant-lambda modes. A high-performance, open source NMR data processing tool is described that implements CAMERA, and brief comparisons to existing reconstruction methods are made on several example spectra.
| S-EPMC4818692 | biostudies-literature
Project description:We describe a simulation strategy for prediction of drug targets and functional analysis of genetically defined associations. This strategy utilizes Bayesian model averaging over a sample of parameterized networks to achieve robust predictions in light of uncertainty of the network structure. We illustrate the method with an analysis of liver gene expression, serum lipid profiles and body weight measured on 120 male mice from a mouse intercross population. An ensemble of 1024 networks, gave accurate predictions of animals that were not part of the training data and explained almost twice the variance compared to quantitative trait loci alone. Additional in silico experiments identified 38 transcripts that are predicted to impact serum lipid profiles and suggest that they play a role in controlling high density lipoprotein and free triglycerides plasma concentrations.
2010-03-10 | GSE15226 | GEO
Project description:The data here reported refer to the numerical examples shown in the research article “Wave based method for flexural vibration of thin plate with general elastically restrained edges” (Liu et al., 2020 [1]). Within the examples, only the datasets regarding the plates with elastic or elastic-damping supports are provided. The datasets contain the raw data directly obtained from the forced vibration simulations. The simulations are carried out using two methods: the finite element method realized in ANSYS Mechanical APDL and the proposed wave based method (Liu et al., 2020 [1]), implemented in a MATLAB code. The data obtained from ANSYS serves as reference for the response of the plate under different boundary conditions. For each frequency, the transverse displacements of the plate at two pre-selected points are listed in the spreadsheet (e.g. MS Excel). When damping is present, they are separated into real part and imaginary part. This part of data can be used as reference when other novel methods are developed. The datasets obtained from MATLAB include the contribution factors as well as the wave functions. Based on them, one can obtain the displacement as a complex number at any point of the plate after a simple postprocessing. Postprocessing codes to obtain the frequency response function for a user-given point and the displacement field at a user-given frequency are also provided. This part of data presents much more information than the previous part as well as the corresponding results in the related research article. It makes it possible to see the responses at other points or other frequencies that are not considered in the research article, without repeating the time-consuming simulations. Moreover, if someone wants to further improve the wave based method, this part of data will be helpful, either for analysing the limitations of the proposed method or for more direct comparisons. Any research related to the flexural vibration of plate can also consider the data provided in this article.
| S-EPMC7334394 | biostudies-literature
Project description:NiP alloys are very often used in industry, due to their outstanding performance in corrosion and wear. Alloys with high phosphorus content (≥16 atom % P) are amorphous and show high corrosion resistance in both neutral and acidic solutions irrespective of the presence of chloride ions. The reason for this behavior is attributed to the formation of a "P-enriched surface layer" whose exact nature is still under debate. In this work, an iterative algorithm based on the application of maximum entropy method on nondestructive angle-resolved X-ray photoelectron spectroscopy data has been applied to the investigation of the surface layer grown on Ni18P alloys following mechanical polishing and anodic polarization in sulfate solutions. The results show that the outermost region of the examined alloy has a complex layered structure: (1) an uppermost hydrocarbon contamination layer about 1 nm thick, which includes also adsorbed water; (2) a nickel (poly)phosphate layer of about 1 nm; (3) a highly phosphorus-enriched interface being about 2 nm thick with a marked phosphorus concentration gradient, from 70 to 20 atom %; and (4) bulk alloy with the stoichiometric composition. These findings, together with the chemical state of the different phosphorus compounds, allow us to conclude that the high corrosion and wear resistance of NiP alloys might be ascribed to the presence of a thin, self-repairing nickel (poly)phosphate layer grown on a strongly P-enriched interface. Because the Auger parameter of P at the interface is similar to that of elemental P, it might be also concluded that the interface is enriched in elemental phosphorus.
| S-EPMC6645160 | biostudies-literature
Project description:Quantiles are available in various problems for developing probability distributions. In some problems quantiles are elicited from experts and used for fitting parametric models, which induce non-elicited information. In some other problems comparisons are made with a quantile of an assumed model which is noncommittal to the quantile information. The maximum entropy (ME) principle provides models that avoid these issues. However, the information theory literature has been mainly concerned about models based on moment information. This paper explores the ME models that are the minimum elaborations of the uniform and moment-based ME models by quantiles. This property provides diagnostics for the utility of elaboration in terms of the information value of each type of information over the other. The ME model with quantiles and moments is represented as the mixture of truncated distributions on consecutive intervals whose shapes and existence are determined by the moments. Elaborations of several ME distributions by quantiles are presented. The ME model based only on quantiles elicited by the fixed interval method possesses a useful property for pooling information elicited from multiple experts. The elaboration of Laplace distribution is an extension of the information theory connection with minimum risk under symmetric loss functions to the asymmetric linear loss. This extension produces a new Asymmetric Laplace distribution. Application examples compare ME priors with a parametric model fitted to elicited quantiles, illustrate measuring uncertainty and disagreement of economic forecasters based on elicited probabilities, and adjust ME models for a fundamental quantile in an inventory management problem.
| S-EPMC7414396 | biostudies-literature
Project description:Recognition of pathogens relies on families of proteins showing great diversity. Here we construct maximum entropy models of the sequence repertoire, building on recent experiments that provide a nearly exhaustive sampling of the IgM sequences in zebrafish. These models are based solely on pairwise correlations between residue positions but correctly capture the higher order statistical properties of the repertoire. By exploiting the interpretation of these models as statistical physics problems, we make several predictions for the collective properties of the sequence ensemble: The distribution of sequences obeys Zipf's law, the repertoire decomposes into several clusters, and there is a massive restriction of diversity because of the correlations. These predictions are completely inconsistent with models in which amino acid substitutions are made independently at each site and are in good agreement with the data. Our results suggest that antibody diversity is not limited by the sequences encoded in the genome and may reflect rapid adaptation to antigenic challenges. This approach should be applicable to the study of the global properties of other protein families.
| S-EPMC2851784 | biostudies-literature