DOE PAGES 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

Abstract

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. Finally, our new framework is generalizablemore » and can be used to investigate the uncertainty of a wide variety of target properties in biomolecular systems.« less

Authors:
 [1];  [1];  [1];  [2];  [1]
  1. Pacific Northwest National Lab. (PNNL), Richland, WA (United States)
  2. Purdue Univ., West Lafayette, IN (United States). Dept. of Mathematics
Publication Date:
Research Org.:
Pacific Northwest National Laboratory (PNNL), Richland, WA (United States)
Sponsoring Org.:
USDOE
OSTI Identifier:
1255409
Report Number(s):
PNNL-SA-110085
Journal ID: ISSN 1540-3459; KJ0401000
Grant/Contract Number:  
AC05-76RL01830
Resource Type:
Accepted Manuscript
Journal Name:
Multiscale Modeling & Simulation
Additional Journal Information:
Journal Volume: 13; Journal Issue: 4; Journal ID: ISSN 1540-3459
Publisher:
SIAM
Country of Publication:
United States
Language:
English
Subject:
59 BASIC BIOLOGICAL SCIENCES; uncertainty quantification; coarse graining; stochastic; biomolecular conformation fluctuation; polynomial chaos; compressive sensing method; model reduction

Citation Formats

Lei, Huan, Yang, Xiu, Zheng, Bin, Lin, Guang, and Baker, Nathan A. Constructing Surrogate Models of Complex Systems with Enhanced Sparsity: Quantifying the Influence of Conformational Uncertainty in Biomolecular Solvation. United States: N. p., 2015. Web. doi:10.1137/140981587.
Lei, Huan, Yang, Xiu, Zheng, Bin, Lin, Guang, & Baker, Nathan A. Constructing Surrogate Models of Complex Systems with Enhanced Sparsity: Quantifying the Influence of Conformational Uncertainty in Biomolecular Solvation. United States. https://doi.org/10.1137/140981587
Lei, Huan, Yang, Xiu, Zheng, Bin, Lin, Guang, and Baker, Nathan A. Thu . "Constructing Surrogate Models of Complex Systems with Enhanced Sparsity: Quantifying the Influence of Conformational Uncertainty in Biomolecular Solvation". United States. https://doi.org/10.1137/140981587. https://www.osti.gov/servlets/purl/1255409.
@article{osti_1255409,
title = {Constructing Surrogate Models of Complex Systems with Enhanced Sparsity: Quantifying the Influence of Conformational Uncertainty in Biomolecular Solvation},
author = {Lei, Huan and Yang, Xiu and Zheng, Bin and Lin, Guang and Baker, Nathan A.},
abstractNote = {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. Finally, our new framework is generalizable and can be used to investigate the uncertainty of a wide variety of target properties in biomolecular systems.},
doi = {10.1137/140981587},
journal = {Multiscale Modeling & Simulation},
number = 4,
volume = 13,
place = {United States},
year = {Thu Nov 05 00:00:00 EST 2015},
month = {Thu Nov 05 00:00:00 EST 2015}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 23 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

Molecular Dynamics:  Survey of Methods for Simulating the Activity of Proteins
journal, May 2006

  • Adcock, Stewart A.; McCammon, J. Andrew
  • Chemical Reviews, Vol. 106, Issue 5
  • DOI: 10.1021/cr040426m

Progress in the prediction of pKa values in proteins
journal, October 2011

  • Alexov, Emil; Mehler, Ernest L.; Baker, Nathan
  • Proteins: Structure, Function, and Bioinformatics, Vol. 79, Issue 12
  • DOI: 10.1002/prot.23189

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

  • Brooks, B.; Karplus, M.
  • Proceedings of the National Academy of Sciences, Vol. 80, Issue 21
  • DOI: 10.1073/pnas.80.21.6571

From Sparse Solutions of Systems of Equations to Sparse Modeling of Signals and Images
journal, February 2009

  • Bruckstein, Alfred M.; Donoho, David L.; Elad, Michael
  • SIAM Review, Vol. 51, Issue 1, p. 34-81
  • DOI: 10.1137/060657704

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

  • Constantine, Paul G.; Dow, Eric; Wang, Qiqi
  • SIAM Journal on Scientific Computing, Vol. 36, Issue 4
  • DOI: 10.1137/130916138

Stable recovery of sparse overcomplete representations in the presence of noise
journal, January 2006

  • Donoho, D. L.; Elad, M.; Temlyakov, V. N.
  • IEEE Transactions on Information Theory, Vol. 52, Issue 1
  • DOI: 10.1109/TIT.2005.860430

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

  • Foo, Jasmine; Wan, Xiaoliang; Karniadakis, George Em
  • Journal of Computational Physics, Vol. 227, Issue 22
  • DOI: 10.1016/j.jcp.2008.07.009

Sparse grid collocation schemes for stochastic natural convection problems
journal, July 2007

  • Ganapathysubramanian, Baskar; Zabaras, Nicholas
  • Journal of Computational Physics, Vol. 225, Issue 1
  • DOI: 10.1016/j.jcp.2006.12.014

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

  • Go, N.; Noguti, T.; Nishikawa, T.
  • Proceedings of the National Academy of Sciences, Vol. 80, Issue 12
  • DOI: 10.1073/pnas.80.12.3696

The Split Bregman Method for L1-Regularized Problems
journal, January 2009

  • Goldstein, Tom; Osher, Stanley
  • SIAM Journal on Imaging Sciences, Vol. 2, Issue 2
  • DOI: 10.1137/080725891

Gaussian Dynamics of Folded Proteins
journal, October 1997


Influence of Grid Spacing in Poisson–Boltzmann Equation Binding Energy Estimation
journal, July 2013

  • Harris, Robert C.; Boschitsch, Alexander H.; Fenley, Marcia O.
  • Journal of Chemical Theory and Computation, Vol. 9, Issue 8
  • DOI: 10.1021/ct300765w

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

  • Li, Xin
  • IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems, Vol. 29, Issue 11
  • DOI: 10.1109/TCAD.2010.2061292

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

  • Narayan, Akil; Xiu, Dongbin
  • SIAM Journal on Scientific Computing, Vol. 34, Issue 3
  • DOI: 10.1137/110854059

Constructing Nested Nodal Sets for Multivariate Polynomial Interpolation
journal, January 2013

  • Narayan, Akil; Xiu, Dongbin
  • SIAM Journal on Scientific Computing, Vol. 35, Issue 5
  • DOI: 10.1137/12089613X

An Anisotropic Sparse Grid Stochastic Collocation Method for Partial Differential Equations with Random Input Data
journal, January 2008

  • Nobile, F.; Tempone, R.; Webster, C. G.
  • SIAM Journal on Numerical Analysis, Vol. 46, Issue 5
  • DOI: 10.1137/070680540

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

  • Ren, Pengyu; Chun, Jaehun; Thomas, Dennis G.
  • Quarterly Reviews of Biophysics, Vol. 45, Issue 4
  • DOI: 10.1017/S003358351200011X

Solvent accessible surface area and excluded volume in proteins
journal, September 1984


Uncertainty Quantification in MD Simulations. Part I: Forward Propagation
journal, January 2012

  • Rizzi, F.; Najm, H. N.; Debusschere, B. J.
  • Multiscale Modeling & Simulation, Vol. 10, Issue 4
  • DOI: 10.1137/110853169

Uncertainty Quantification in MD Simulations. Part II: Bayesian Inference of Force-Field Parameters
journal, January 2012

  • Rizzi, F.; Najm, H. N.; Debusschere, B. J.
  • Multiscale Modeling & Simulation, Vol. 10, Issue 4
  • DOI: 10.1137/110853170

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

  • Tama, F.; Sanejouand, Y. -H.
  • Protein Engineering, Design and Selection, Vol. 14, Issue 1
  • DOI: 10.1093/protein/14.1.1

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

  • van den Berg, Ewout; Friedlander, Michael P.
  • SIAM Journal on Scientific Computing, Vol. 31, Issue 2
  • DOI: 10.1137/080714488

Structure of bovine pancreatic trypsin inhibitor
journal, December 1984


High-Order Collocation Methods for Differential Equations with Random Inputs
journal, January 2005

  • Xiu, Dongbin; Hesthaven, Jan S.
  • SIAM Journal on Scientific Computing, Vol. 27, Issue 3
  • DOI: 10.1137/040615201

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

  • Yin, Wotao; Osher, Stanley; Goldfarb, Donald
  • SIAM Journal on Imaging Sciences, Vol. 1, Issue 1
  • DOI: 10.1137/070703983

Error Estimates for the ANOVA Method with Polynomial Chaos Interpolation: Tensor Product Functions
journal, January 2012

  • Zhang, Zhongqiang; Choi, Minseok; Karniadakis, George Em
  • SIAM Journal on Scientific Computing, Vol. 34, Issue 2
  • DOI: 10.1137/100788859

A Recursive Sparse Grid Collocation Method for Differential Equations with White Noise
journal, January 2014

  • Zhang, Z.; Tretyakov, M. V.; Rozovskii, B.
  • SIAM Journal on Scientific Computing, Vol. 36, Issue 4
  • DOI: 10.1137/130938906

Works referencing / citing this record:

Calibration of reduced-order model for a coupled Burgers equations based on PC-EnKF
journal, January 2018

  • Wang, Yuepeng; Cheng, Yue; Zhang, Zongyuan
  • Mathematical Modelling of Natural Phenomena, Vol. 13, Issue 2
  • DOI: 10.1051/mmnp/2018023

Probabilistic surrogate models for uncertainty analysis: Dimension reduction‐based polynomial chaos expansion
journal, November 2019

  • Son, Jeongeun; Du, Yuncheng
  • International Journal for Numerical Methods in Engineering, Vol. 121, Issue 6
  • DOI: 10.1002/nme.6262