Constructing Surrogate Models of Complex Systems with Enhanced Sparsity: Quantifying the Influence of Conformational Uncertainty in Biomolecular Solvation
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 fluctuationinduced uncertainty in solventaccessible 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 highdimensional 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 »
 Authors:

^{[1]};
^{[1]};
^{[1]};
^{[2]};
^{[1]}
 Pacific Northwest National Lab. (PNNL), Richland, WA (United States)
 Purdue Univ., West Lafayette, IN (United States). Dept. of Mathematics
 Publication Date:
 Report Number(s):
 PNNLSA110085
Journal ID: ISSN 15403459; KJ0401000
 Grant/Contract Number:
 AC0576RL01830
 Type:
 Accepted Manuscript
 Journal Name:
 Multiscale Modeling & Simulation
 Additional Journal Information:
 Journal Volume: 13; Journal Issue: 4; Journal ID: ISSN 15403459
 Publisher:
 SIAM
 Research Org:
 Pacific Northwest National Lab. (PNNL), Richland, WA (United States)
 Sponsoring Org:
 USDOE
 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
 OSTI Identifier:
 1255409
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.,
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. doi:10.1137/140981587.
Lei, Huan, Yang, Xiu, Zheng, Bin, Lin, Guang, and Baker, Nathan A. 2015.
"Constructing Surrogate Models of Complex Systems with Enhanced Sparsity: Quantifying the Influence of Conformational Uncertainty in Biomolecular Solvation". United States.
doi: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 fluctuationinduced uncertainty in solventaccessible 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 highdimensional 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 = {2015},
month = {11}
}