{"title":"自动前向模型参数化与构象群体的贝叶斯推断","authors":"Robert M. Raddi, Tim Marshall, Vincent A. Voelz","doi":"arxiv-2405.18532","DOIUrl":null,"url":null,"abstract":"To quantify how well theoretical predictions of structural ensembles agree\nwith experimental measurements, we depend on the accuracy of forward models.\nThese models are computational frameworks that generate observable quantities\nfrom molecular configurations based on empirical relationships linking specific\nmolecular properties to experimental measurements. Bayesian Inference of\nConformational Populations (BICePs) is a reweighting algorithm that reconciles\nsimulated ensembles with ensemble-averaged experimental observations, even when\nsuch observations are sparse and/or noisy. This is achieved by sampling the\nposterior distribution of conformational populations under experimental\nrestraints as well as sampling the posterior distribution of uncertainties due\nto random and systematic error. In this study, we enhance the algorithm for the\nrefinement of empirical forward model (FM) parameters. We introduce and\nevaluate two novel methods for optimizing FM parameters. The first method\ntreats FM parameters as nuisance parameters, integrating over them in the full\nposterior distribution. The second method employs variational minimization of a\nquantity called the BICePs score that reports the free energy of `turning on`\nthe experimental restraints. This technique, coupled with improved likelihood\nfunctions for handling experimental outliers, facilitates force field\nvalidation and optimization, as illustrated in recent studies (Raddi et al.\n2023, 2024). Using this approach, we refine parameters that modulate the\nKarplus relation, crucial for accurate predictions of J-coupling constants\nbased on dihedral angles between interacting nuclei. We validate this approach\nfirst with a toy model system, and then for human ubiquitin, predicting six\nsets of Karplus parameters. This approach, which does not rely on predetermined\nparameters, enhances predictive accuracy and can be used for many applications.","PeriodicalId":501065,"journal":{"name":"arXiv - PHYS - Data Analysis, Statistics and Probability","volume":"86 1","pages":""},"PeriodicalIF":0.0000,"publicationDate":"2024-05-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"Automatic Forward Model Parameterization with Bayesian Inference of Conformational Populations\",\"authors\":\"Robert M. Raddi, Tim Marshall, Vincent A. Voelz\",\"doi\":\"arxiv-2405.18532\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"To quantify how well theoretical predictions of structural ensembles agree\\nwith experimental measurements, we depend on the accuracy of forward models.\\nThese models are computational frameworks that generate observable quantities\\nfrom molecular configurations based on empirical relationships linking specific\\nmolecular properties to experimental measurements. Bayesian Inference of\\nConformational Populations (BICePs) is a reweighting algorithm that reconciles\\nsimulated ensembles with ensemble-averaged experimental observations, even when\\nsuch observations are sparse and/or noisy. This is achieved by sampling the\\nposterior distribution of conformational populations under experimental\\nrestraints as well as sampling the posterior distribution of uncertainties due\\nto random and systematic error. In this study, we enhance the algorithm for the\\nrefinement of empirical forward model (FM) parameters. We introduce and\\nevaluate two novel methods for optimizing FM parameters. The first method\\ntreats FM parameters as nuisance parameters, integrating over them in the full\\nposterior distribution. The second method employs variational minimization of a\\nquantity called the BICePs score that reports the free energy of `turning on`\\nthe experimental restraints. This technique, coupled with improved likelihood\\nfunctions for handling experimental outliers, facilitates force field\\nvalidation and optimization, as illustrated in recent studies (Raddi et al.\\n2023, 2024). Using this approach, we refine parameters that modulate the\\nKarplus relation, crucial for accurate predictions of J-coupling constants\\nbased on dihedral angles between interacting nuclei. We validate this approach\\nfirst with a toy model system, and then for human ubiquitin, predicting six\\nsets of Karplus parameters. This approach, which does not rely on predetermined\\nparameters, enhances predictive accuracy and can be used for many applications.\",\"PeriodicalId\":501065,\"journal\":{\"name\":\"arXiv - PHYS - Data Analysis, Statistics and Probability\",\"volume\":\"86 1\",\"pages\":\"\"},\"PeriodicalIF\":0.0000,\"publicationDate\":\"2024-05-28\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"arXiv - PHYS - Data Analysis, Statistics and Probability\",\"FirstCategoryId\":\"1085\",\"ListUrlMain\":\"https://doi.org/arxiv-2405.18532\",\"RegionNum\":0,\"RegionCategory\":null,\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"\",\"JCRName\":\"\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"arXiv - PHYS - Data Analysis, Statistics and Probability","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/arxiv-2405.18532","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
Automatic Forward Model Parameterization with Bayesian Inference of Conformational Populations
To quantify how well theoretical predictions of structural ensembles agree
with experimental measurements, we depend on the accuracy of forward models.
These models are computational frameworks that generate observable quantities
from molecular configurations based on empirical relationships linking specific
molecular properties to experimental measurements. Bayesian Inference of
Conformational Populations (BICePs) is a reweighting algorithm that reconciles
simulated ensembles with ensemble-averaged experimental observations, even when
such observations are sparse and/or noisy. This is achieved by sampling the
posterior distribution of conformational populations under experimental
restraints as well as sampling the posterior distribution of uncertainties due
to random and systematic error. In this study, we enhance the algorithm for the
refinement of empirical forward model (FM) parameters. We introduce and
evaluate two novel methods for optimizing FM parameters. The first method
treats FM parameters as nuisance parameters, integrating over them in the full
posterior distribution. The second method employs variational minimization of a
quantity called the BICePs score that reports the free energy of `turning on`
the experimental restraints. This technique, coupled with improved likelihood
functions for handling experimental outliers, facilitates force field
validation and optimization, as illustrated in recent studies (Raddi et al.
2023, 2024). Using this approach, we refine parameters that modulate the
Karplus relation, crucial for accurate predictions of J-coupling constants
based on dihedral angles between interacting nuclei. We validate this approach
first with a toy model system, and then for human ubiquitin, predicting six
sets of Karplus parameters. This approach, which does not rely on predetermined
parameters, enhances predictive accuracy and can be used for many applications.