Browse
Submit Data
Databases
API
Help

Dataset Information

21 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Use of Bayesian Inference in Crystallographic Structure Refinement via Full Diffraction Profile Analysis.

ABSTRACT: A Bayesian inference method for refining crystallographic structures is presented. The distribution of model parameters is stochastically sampled using Markov chain Monte Carlo. Posterior probability distributions are constructed for all model parameters to properly quantify uncertainty by appropriately modeling the heteroskedasticity and correlation of the error structure. The proposed method is demonstrated by analyzing a National Institute of Standards and Technology silicon standard reference material. The results obtained by Bayesian inference are compared with those determined by Rietveld refinement. Posterior probability distributions of model parameters provide both estimates and uncertainties. The new method better estimates the true uncertainties in the model as compared to the Rietveld method.

SUBMITTER: Fancher CM

PROVIDER: S-EPMC4994022 | biostudies-literature | 2016 Aug

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Publications

Use of Bayesian Inference in Crystallographic Structure Refinement via Full Diffraction Profile Analysis.

Fancher Chris M CM Han Zhen Z Levin Igor I Page Katharine K Reich Brian J BJ Smith Ralph C RC Wilson Alyson G AG Jones Jacob L JL

Scientific reports 20160823

A Bayesian inference method for refining crystallographic structures is presented. The distribution of model parameters is stochastically sampled using Markov chain Monte Carlo. Posterior probability distributions are constructed for all model parameters to properly quantify uncertainty by appropriately modeling the heteroskedasticity and correlation of the error structure. The proposed method is demonstrated by analyzing a National Institute of Standards and Technology silicon standard referenc ...[more]

PMID: 27550221

Similar Datasets

Towards automated crystallographic structure refinement with phenix.refine.

Project description:phenix.refine is a program within the PHENIX package that supports crystallographic structure refinement against experimental data with a wide range of upper resolution limits using a large repertoire of model parameterizations. It has several automation features and is also highly flexible. Several hundred parameters enable extensive customizations for complex use cases. Multiple user-defined refinement strategies can be applied to specific parts of the model in a single refinement run. An intuitive graphical user interface is available to guide novice users and to assist advanced users in managing refinement projects. X-ray or neutron diffraction data can be used separately or jointly in refinement. phenix.refine is tightly integrated into the PHENIX suite, where it serves as a critical component in automated model building, final structure refinement, structure validation and deposition to the wwPDB. This paper presents an overview of the major phenix.refine features, with extensive literature references for readers interested in more detailed discussions of the methods.

| S-EPMC3322595 | biostudies-literature

Structure of the full-length Shaker potassium channel Kv1.2 by normal-mode-based X-ray crystallographic refinement.

Project description:Voltage-dependent potassium channels (Kv) are homotetramers composed of four voltage sensors and one pore domain. Because of high-level structural flexibility, the first mammalian Kv structure, Kv1.2 at 2.9 A, has about 37% molecular mass of the transmembrane portion not resolved. In this study, by applying a novel normal-mode-based X-ray crystallographic refinement method to the original diffraction data and structural model, we established the structure of full-length Kv1.2 in its native form. This structure offers mechanistic insights into voltage sensing. Particularly, it shows a hydrophobic layer of about 10 A at the midpoint of the membrane bilayer, which is likely the molecular basis for the observed "focused electric field" of Kv1.2 between the internal and external solutions. This work also demonstrated the potential of the refinement method in bringing up large chunks of missing densities, thus beneficial to structural refinement of many difficult systems.

| S-EPMC2895106 | biostudies-literature

Robust Bayesian inference via coarsening.

Project description:The standard approach to Bayesian inference is based on the assumption that the distribution of the data belongs to the chosen model class. However, even a small violation of this assumption can have a large impact on the outcome of a Bayesian procedure. We introduce a novel approach to Bayesian inference that improves robustness to small departures from the model: rather than conditioning on the event that the observed data are generated by the model, one conditions on the event that the model generates data close to the observed data, in a distributional sense. When closeness is defined in terms of relative entropy, the resulting "coarsened" posterior can be approximated by simply tempering the likelihood-that is, by raising the likelihood to a fractional power-thus, inference can usually be implemented via standard algorithms, and one can even obtain analytical solutions when using conjugate priors. Some theoretical properties are derived, and we illustrate the approach with real and simulated data using mixture models and autoregressive models of unknown order.

| S-EPMC6961963 | biostudies-literature

Hirshfeld atom refinement and dynamical refinement of hexagonal ice structure from electron diffraction data.

Project description:Reaching beyond the commonly used spherical atomic electron density model allows one to greatly improve the accuracy of hydrogen atom structural parameters derived from X-ray data. However, the effects of atomic asphericity are less explored for electron diffraction data. In this work, Hirshfeld atom refinement (HAR), a method that uses an accurate description of electron density by quantum mechanical calculation for a system of interest, was applied for the first time to the kinematical refinement of electron diffraction data. This approach was applied here to derive the structure of ordinary hexagonal ice (Ih). The effect of introducing HAR is much less noticeable than in the case of X-ray refinement and it is largely overshadowed by dynamical scattering effects. It led to only a slight change in the O-H bond lengths (shortening by 0.01 Å) compared with the independent atom model (IAM). The average absolute differences in O-H bond lengths between the kinematical refinements and the reference neutron structure were much larger: 0.044 for IAM and 0.046 Å for HAR. The refinement results changed considerably when dynamical scattering effects were modelled - with extinction correction or with dynamical refinement. The latter led to an improvement of the O-H bond length accuracy to 0.021 Å on average (with IAM refinement). Though there is a potential for deriving more accurate structures using HAR for electron diffraction, modelling of dynamical scattering effects seems to be a necessary step to achieve this. However, at present there is no software to support both HAR and dynamical refinement.

| S-EPMC11364029 | biostudies-literature

Automating crystallographic structure solution and refinement of protein-ligand complexes.

Project description:High-throughput drug-discovery and mechanistic studies often require the determination of multiple related crystal structures that only differ in the bound ligands, point mutations in the protein sequence and minor conformational changes. If performed manually, solution and refinement requires extensive repetition of the same tasks for each structure. To accelerate this process and minimize manual effort, a pipeline encompassing all stages of ligand building and refinement, starting from integrated and scaled diffraction intensities, has been implemented in Phenix. The resulting system is able to successfully solve and refine large collections of structures in parallel without extensive user intervention prior to the final stages of model completion and validation.

| S-EPMC3919266 | biostudies-literature

Crystallographic refinement of ligand complexes.

Project description:Model building and refinement of complexes between biomacromolecules and small molecules requires sensible starting coordinates as well as the specification of restraint sets for all but the most common non-macromolecular entities. Here, it is described why this is necessary, how it can be accomplished and what pitfalls need to be avoided in order to produce chemically plausible models of the low-molecular-weight entities. A number of programs, servers, databases and other resources that can be of assistance in the process are also discussed.

| S-EPMC2483469 | biostudies-literature

Bayesian inference of protein structure from chemical shift data.

Project description:Protein chemical shifts are routinely used to augment molecular mechanics force fields in protein structure simulations, with weights of the chemical shift restraints determined empirically. These weights, however, might not be an optimal descriptor of a given protein structure and predictive model, and a bias is introduced which might result in incorrect structures. In the inferential structure determination framework, both the unknown structure and the disagreement between experimental and back-calculated data are formulated as a joint probability distribution, thus utilizing the full information content of the data. Here, we present the formulation of such a probability distribution where the error in chemical shift prediction is described by either a Gaussian or Cauchy distribution. The methodology is demonstrated and compared to a set of empirically weighted potentials through Markov chain Monte Carlo simulations of three small proteins (ENHD, Protein G and the SMN Tudor Domain) using the PROFASI force field and the chemical shift predictor CamShift. Using a clustering-criterion for identifying the best structure, together with the addition of a solvent exposure scoring term, the simulations suggests that sampling both the structure and the uncertainties in chemical shift prediction leads more accurate structures compared to conventional methods using empirical determined weights. The Cauchy distribution, using either sampled uncertainties or predetermined weights, did, however, result in overall better convergence to the native fold, suggesting that both types of distribution might be useful in different aspects of the protein structure prediction.

| S-EPMC4375973 | biostudies-literature

Phasing electron diffraction data by molecular replacement: strategy for structure determination and refinement.

Project description:Electron crystallography is arguably the only electron cryomicroscopy (cryo EM) technique able to deliver atomic resolution data (better then 3 Å) for membrane proteins embedded in a membrane. The progress in hardware improvements and sample preparation for diffraction analysis resulted in a number of recent examples where increasingly higher resolutions were achieved. Other chapters in this book detail the improvements in hardware and delve into the intricate art of sample preparation for microscopy and electron diffraction data collection and processing. In this chapter, we describe in detail the protocols for molecular replacement for electron diffraction studies. The use of a search model for phasing electron diffraction data essentially eliminates the need of acquiring image data rendering it immune to aberrations from drift and charging effects that effectively lower the attainable resolution.

| S-EPMC3851287 | biostudies-literature

Analyzing single-molecule time series via nonparametric Bayesian inference.

Project description:The ability to measure the properties of proteins at the single-molecule level offers an unparalleled glimpse into biological systems at the molecular scale. The interpretation of single-molecule time series has often been rooted in statistical mechanics and the theory of Markov processes. While existing analysis methods have been useful, they are not without significant limitations including problems of model selection and parameter nonidentifiability. To address these challenges, we introduce the use of nonparametric Bayesian inference for the analysis of single-molecule time series. These methods provide a flexible way to extract structure from data instead of assuming models beforehand. We demonstrate these methods with applications to several diverse settings in single-molecule biophysics. This approach provides a well-constrained and rigorously grounded method for determining the number of biophysical states underlying single-molecule data.

| S-EPMC4317543 | biostudies-literature

Unbiased Bayesian inference for population Markov jump processes via random truncations.

Project description:We consider continuous time Markovian processes where populations of individual agents interact stochastically according to kinetic rules. Despite the increasing prominence of such models in fields ranging from biology to smart cities, Bayesian inference for such systems remains challenging, as these are continuous time, discrete state systems with potentially infinite state-space. Here we propose a novel efficient algorithm for joint state/parameter posterior sampling in population Markov Jump processes. We introduce a class of pseudo-marginal sampling algorithms based on a random truncation method which enables a principled treatment of infinite state spaces. Extensive evaluation on a number of benchmark models shows that this approach achieves considerable savings compared to state of the art methods, retaining accuracy and fast convergence. We also present results on a synthetic biology data set showing the potential for practical usefulness of our work.

| S-EPMC5477715 | biostudies-other

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