Dataset Information

The Effects of Computational Modeling Errors on the Estimation of Statistical Mechanical Variables.

ABSTRACT: Computational models used in the estimation of thermodynamic quantities of large chemical systems often require approximate energy models that rely on parameterization and cancellation of errors to yield agreement with experimental measurements. In this work, we show how energy function errors propagate when computing statistical mechanics-derived thermodynamic quantities. Assuming that each microstate included in a statistical ensemble has a measurable amount of error in its calculated energy, we derive low-order expressions for the propagation of these errors in free energy, average energy, and entropy. Through gedanken experiments we show the expected behavior of these error propagation formulas on hypothetical energy surfaces. For very large microstate energy errors, these low-order formulas disagree with estimates from Monte Carlo simulations of error propagation. Hence, such simulations of error propagation may be required when using poor potential energy functions. Propagated systematic errors predicted by these methods can be removed from computed quantities, while propagated random errors yield uncertainty estimates. Importantly, we find that end-point free energy methods maximize random errors and that local sampling of potential energy wells decreases random error significantly. Hence, end-point methods should be avoided in energy computations and should be replaced by methods that incorporate local sampling. The techniques described herein will be used in future work involving the calculation of free energies of biomolecular processes, where error corrections are expected to yield improved agreement with experiment.

SUBMITTER: Faver JC

PROVIDER: S-EPMC3568774 | biostudies-literature | 2012 Oct

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

The Effects of Computational Modeling Errors on the Estimation of Statistical Mechanical Variables.

Faver John C JC Yang Wei W Merz Kenneth M KM

Journal of chemical theory and computation 20120312 10

Computational models used in the estimation of thermodynamic quantities of large chemical systems often require approximate energy models that rely on parameterization and cancellation of errors to yield agreement with experimental measurements. In this work, we show how energy function errors propagate when computing statistical mechanics-derived thermodynamic quantities. Assuming that each microstate included in a statistical ensemble has a measurable amount of error in its calculated energy, ...[more]

PMID: 23413365

Similar Datasets

Project description:BackgroundThere is increasing interest in Positron Emission Tomography (PET) of 2-deoxy-2-[18F]flouro-D-glucose ((18)F-FDG) to evaluate pulmonary inflammation during acute lung injury (ALI). We assessed the effect of extra-vascular lung water on estimates of (18)F-FDG-kinetics parameters in experimental and simulated data using the Patlak and Sokoloff methods, and our recently proposed four-compartment model.Methodology/principal findingsEleven sheep underwent unilateral lung lavage and 4 h mechanical ventilation. Five sheep received intravenous endotoxin (10 ng/kg/min). Dynamic (18)F-FDG PET was performed at the end of the 4 h period. (18)F-FDG net uptake rate (Ki), phosphorylation rate (k(3)), and volume of distribution (F(e)) were estimated in three isogravitational regions for each method. Simulations of normal and ALI (18)F-FDG-kinetics were conducted to study the dependence of estimated parameters on the transport rate constants to (k(5)) and from (k(6)) the extra-vascular extra-cellular compartment. The four-compartment model described 85.7% of the studied (18)F-FDG-kinetics better than the Sokoloff model. Relative to the four-compartment model the Sokoloff model exhibited a consistent positive bias in Ki (3.32 [1.30-5.65] 10(-4)/min, p<0.001) and showed inaccurate estimates of the parameters composing Ki (k(3) and F(e)), even when Ki was similar for those methods. In simulations, errors in estimates of Ki due to the extra-vascular extra-cellular compartment depended on both k(5) and k(5)/k(6), with errors for the Patlak and Sokoloff methods of 0.02 [-0.01-0.18] and 0.40 [0.18-0.60] 10(-3)/min for normal lungs and of -0.47 [-0.89-0.72] and 2.35 [0.85-3.68] 10(-3)/min in ALI.Conclusions/significance(18)F-FDG accumulation in lung extra-vascular fluid, which is commonly increased during lung injury, can result in substantial estimation errors using the traditional Patlak and Sokoloff methods. These errors depend on the extra-vascular extra-cellular compartment volume and its transport rates with other compartments. The four-compartment model provides more accurate quantification of (18)F-FDG-kinetics than those methods in the presence of increased extra-vascular fluid.

Project description:Cell population heterogeneity can affect cellular response and is a major factor in drug resistance. However, there are few techniques available to represent and explore how heterogeneity is linked to population response. Recent high-throughput genomic, proteomic, and cellomic approaches offer opportunities for profiling heterogeneity on several scales. We have recently examined heterogeneity in vascular endothelial growth factor receptor (VEGFR) membrane localization in endothelial cells. We and others processed the heterogeneous data through ensemble averaging and integrated the data into computational models of anti-angiogenic drug effects in breast cancer. Here we show that additional modeling insight can be gained when cellular heterogeneity is considered. We present comprehensive statistical and computational methods for analyzing cellomic data sets and integrating them into deterministic models. We present a novel method for optimizing the fit of statistical distributions to heterogeneous data sets to preserve important data and exclude outliers. We compare methods of representing heterogeneous data and show methodology can affect model predictions up to 3.9-fold. We find that VEGF levels, a target for tuning angiogenesis, are more sensitive to VEGFR1 cell surface levels than VEGFR2; updating VEGFR1 levels in the tumor model gave a 64% change in free VEGF levels in the blood compartment, whereas updating VEGFR2 levels gave a 17% change. Furthermore, we find that subpopulations of tumor cells and tumor endothelial cells (tEC) expressing high levels of VEGFR (>35,000 VEGFR/cell) negate anti-VEGF treatments. We show that lowering the VEGFR membrane insertion rate for these subpopulations recovers the anti-angiogenic effect of anti-VEGF treatment, revealing new treatment targets for specific tumor cell subpopulations. This novel method of characterizing heterogeneous distributions shows for the first time how different representations of the same data set lead to different predictions of drug efficacy.

Dataset Information

The Effects of Computational Modeling Errors on the Estimation of Statistical Mechanical Variables.

Publications

The Effects of Computational Modeling Errors on the Estimation of Statistical Mechanical Variables.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure