skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Constructing Surrogate Models of Complex Systems with Enhanced Sparsity: Quantifying the influence of conformational uncertainty in biomolecular solvation

Journal Article · · Multiscale Modeling & Simulation
DOI:https://doi.org/10.1137/140981587· OSTI ID:1260871

Biomolecules exhibit conformational fluctuations near equilibrium states, inducing uncertainty in various biological properties in a dynamic way. We have developed a general method to quantify the uncertainty of target properties induced by conformational fluctuations. Using a generalized polynomial chaos (gPC) expansion, we construct a surrogate model of the target property with respect to varying conformational states. We also propose a method to increase the sparsity of the gPC expansion by defining a set of conformational “active space” random variables. With the increased sparsity, we employ the compressive sensing method to accurately construct the surrogate model. We demonstrate the performance of the surrogate model by evaluating fluctuation-induced uncertainty in solvent-accessible surface area for the bovine trypsin inhibitor protein system and show that the new approach offers more accurate statistical information than standard Monte Carlo approaches. Further more, the constructed surrogate model also enables us to directly evaluate the target property under various conformational states, yielding a more accurate response surface than standard sparse grid collocation methods. In particular, the new method provides higher accuracy in high-dimensional systems, such as biomolecules, where sparse grid performance is limited by the accuracy of the computed quantity of interest. Our new framework is generalizable and can be used to investigate the uncertainty of a wide variety of target properties in biomolecular systems.

Research Organization:
Pacific Northwest National Laboratory (PNNL), Richland, WA (United States)
Sponsoring Organization:
USDOE
DOE Contract Number:
AC05-76RL01830
OSTI ID:
1260871
Report Number(s):
PNNL-SA-104506; KJ0401000
Journal Information:
Multiscale Modeling & Simulation, Vol. 13, Issue 4; ISSN 1540-3459
Publisher:
SIAM
Country of Publication:
United States
Language:
English

References (48)

Molecular Dynamics:  Survey of Methods for Simulating the Activity of Proteins journal May 2006
Progress in the prediction of pKa values in proteins journal October 2011
Anisotropy of Fluctuation Dynamics of Proteins with an Elastic Network Model journal January 2001
Biomolecular Applications of Poisson-Boltzmann Methods book January 2005
Harmonic dynamics of proteins: normal modes and fluctuations in bovine pancreatic trypsin inhibitor. journal November 1983
From Sparse Solutions of Systems of Equations to Sparse Modeling of Signals and Images journal February 2009
The restricted isometry property and its implications for compressed sensing journal May 2008
Decoding by Linear Programming journal December 2005
Solvent-accessible surfaces of proteins and nucleic acids journal August 1983
Active Subspace Methods in Theory and Practice: Applications to Kriging Surfaces journal January 2014
Stable recovery of sparse overcomplete representations in the presence of noise journal January 2006
A non-adapted sparse approximation of PDEs with stochastic inputs journal April 2011
Biomolecular Simulation: A Computational Microscope for Molecular Biology journal June 2012
Multi-element probabilistic collocation method in high dimensions journal March 2010
The multi-element probabilistic collocation method (ME-PCM): Error analysis and applications journal November 2008
Sparse grid collocation schemes for stochastic natural convection problems journal July 2007
Fully symmetric interpolatory rules for multiple integrals over infinite regions with Gaussian weight journal July 1996
Dynamics of a small globular protein in terms of low-frequency vibrational modes. journal June 1983
The Split Bregman Method for L1-Regularized Problems journal January 2009
Gaussian Dynamics of Folded Proteins journal October 1997
Influence of Grid Spacing in Poisson–Boltzmann Equation Binding Energy Estimation journal July 2013
The interpretation of protein structures: Estimation of static accessibility journal February 1971
Finding Deterministic Solution From Underdetermined Equation: Large-Scale Performance Variability Modeling of Analog/RF Circuits journal November 2010
An adaptive hierarchical sparse grid collocation algorithm for the solution of stochastic differential equations journal May 2009
An adaptive high-dimensional stochastic model representation technique for the solution of stochastic partial differential equations journal May 2010
Stochastic Collocation Methods on Unstructured Grids in High Dimensions via Interpolation journal January 2012
Constructing Nested Nodal Sets for Multivariate Polynomial Interpolation journal January 2013
An Anisotropic Sparse Grid Stochastic Collocation Method for Partial Differential Equations with Random Input Data journal January 2008
Sparse Legendre expansions via <mml:math altimg="si12.gif" display="inline" overflow="scroll" xmlns:xocs="http://www.elsevier.com/xml/xocs/dtd" xmlns:xs="http://www.w3.org/2001/XMLSchema" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns="http://www.elsevier.com/xml/ja/dtd" xmlns:ja="http://www.elsevier.com/xml/ja/dtd" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:tb="http://www.elsevier.com/xml/common/table/dtd" xmlns:sb="http://www.elsevier.com/xml/common/struct-bib/dtd" xmlns:ce="http://www.elsevier.com/xml/common/dtd" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:cals="http://www.elsevier.com/xml/common/cals/dtd"><mml:msub><mml:mrow><mml:mi>ℓ</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:math>-minimization journal May 2012
Biomolecular electrostatics and solvation: a computational perspective journal November 2012
Solvent accessible surface area and excluded volume in proteins journal September 1984
Uncertainty Quantification in MD Simulations. Part I: Forward Propagation journal January 2012
Uncertainty Quantification in MD Simulations. Part II: Bayesian Inference of Force-Field Parameters journal January 2012
Implicit solvent models journal April 1999
Environment and exposure to solvent of protein atoms. Lysozyme and insulin journal September 1973
Conformational change of proteins arising from normal mode calculations journal January 2001
Basis adaptation in homogeneous chaos spaces journal February 2014
Large Amplitude Elastic Motions in Proteins from a Single-Parameter, Atomic Analysis journal August 1996
Probing the Pareto Frontier for Basis Pursuit Solutions journal January 2009
Structure of bovine pancreatic trypsin inhibitor journal December 1984
High-Order Collocation Methods for Differential Equations with Random Inputs journal January 2005
The Wiener--Askey Polynomial Chaos for Stochastic Differential Equations journal January 2002
STOCHASTIC COLLOCATION ALGORITHMS USING l1-MINIMIZATION journal January 2012
Adaptive ANOVA decomposition of stochastic incompressible and compressible flows journal February 2012
Reweighted minimization method for stochastic elliptic differential equations journal September 2013
Bregman Iterative Algorithms for $\ell_1$-Minimization with Applications to Compressed Sensing journal January 2008
Error Estimates for the ANOVA Method with Polynomial Chaos Interpolation: Tensor Product Functions journal January 2012
A Recursive Sparse Grid Collocation Method for Differential Equations with White Noise journal January 2014

Cited By (3)

A comparative study of coarse-graining methods for polymeric fluids: Mori-Zwanzig vs. iterative Boltzmann inversion vs. stochastic parametric optimization journal July 2016
Calibration of reduced-order model for a coupled Burgers equations based on PC-EnKF journal January 2018
Probabilistic surrogate models for uncertainty analysis: Dimension reduction‐based polynomial chaos expansion journal November 2019

Similar Records

Constructing Surrogate Models of Complex Systems with Enhanced Sparsity: Quantifying the Influence of Conformational Uncertainty in Biomolecular Solvation
Journal Article · Thu Nov 05 00:00:00 EST 2015 · Multiscale Modeling & Simulation · OSTI ID:1260871

A data-driven framework for sparsity-enhanced surrogates with arbitrary mutually dependent randomness
Journal Article · Sat Jun 15 00:00:00 EDT 2019 · Computer Methods in Applied Mechanics and Engineering · OSTI ID:1260871

Atomic Radius and Charge Parameter Uncertainty in Biomolecular Solvation Energy Calculations
Journal Article · Tue Jan 16 00:00:00 EST 2018 · Journal of Chemical Theory and Computation · OSTI ID:1260871