What values of relative numerical tolerance should be chosen in simulation of a deterministic model of a biochemical reaction is unclear, which impairs the modeling effort since the simulation outcomes of a model may depend on the relative numerical tolerance values. In an attempt to provide a guideline to selecting appropriate numerical tolerance values in simulation of in vivo biochemical reactions, reasonable numerical tolerance values were estimated based on the uncertainty principle and assumptions of related cellular parameters. The calculations indicate that relative numerical tolerance values can be reasonably set at or around 10^4 for the concentrations expressed in ng/L. This work also suggests that further reducing relative numerical values may result in erroneous simulation results.Plant Biology, Ecology, and Evolutio

Yang, Louis Z.

Yang, Ming

English

SHAREOK repository

1  Determining a reasonable range of relative numerical tolerance values for simulating deterministic models of biochemical reactions  Ming Yanga,*, Louis Z. Yangb a Department of Plant Biology, Ecology, and Evolution, Oklahoma State University, 301 Physical Sciences, Stillwater, OK 74078, USA b Department of Finance and Business Economics, University of Southern California, HOH Hall - 231, MC-1422, 701 Exposition Boulevard, Ste. 231, Los Angeles, CA 90089, USA  * Corresponding author Email address: ming.yang@okstate.edu (M. Yang).    .CC-BY-NC-ND 4.0 International licenseacertified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under The copyright holder for this preprint (which was notthis version posted January 29, 2018. ; https://doi.org/10.1101/256214doi: bioRxiv preprint 2  ABSTRACT  What values of relative numerical tolerance should be chosen in simulation of a deterministic model of a biochemical reaction is unclear, which impairs the modeling effort since the simulation outcomes of a model may depend on the relative numerical tolerance values. In an attempt to provide a guideline to selecting appropriate numerical tolerance values in simulation of in vivo biochemical reactions, reasonable numerical tolerance values were estimated based on the uncertainty principle and assumptions of related cellular parameters. The calculations indicate that relative numerical tolerance values can be reasonably set at or around 410  for the concentrations expressed in ng/L. This work also suggests that further reducing relative numerical values may result in erroneous simulation results.      Keywords: Chemical concentration  Uncertainty principle Standard deviation Subcellular compartment System behavior    .CC-BY-NC-ND 4.0 International licenseacertified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under The copyright holder for this preprint (which was notthis version posted January 29, 2018. ; https://doi.org/10.1101/256214doi: bioRxiv preprint 3  A deterministic model of a biological process in the form of differential equations presents mathematical relationships among variables of precise values, but in simulation of such a model, only approximate values of the variables can be employed. The approximation is governed by either an absolute or a relative numerical tolerance value. An absolute numerical tolerance value is the absolute difference between the actual and expected values. Dividing such an absolute difference by the expected value gives rise to a relative numerical tolerance value. Because a model of a biological system may produce different simulation outcomes with different ranges of numerical tolerance values, it is important to determine a reasonable range of numerical tolerance values. Currently a reasonable range of numerical tolerance values for modeling a biological process such as that involving biochemical reactions has not been explored.   Numerical tolerance issues have been studied before, which led to the general conclusion that decreasing numerical tolerance values to a sufficiently small and yet practical level may effectively correct erroneous outcomes in computer simulation (Byrne and Thompson, 2013). The erroneous outcomes arise from numerical stiffness, a phenomenon of unacceptably large errors in dependent variables as a result of small errors in independent variables (DiStefano III, 2015). Here we present a case that a stable numerical outcome may be achieved only when the numerical tolerance values approach zero. This case was previously described in an attempt to model oscillations of biochemical reactions in a negative feedback system (Yang and Yang, 2015). In this case, the magnitudes of the oscillations of chemicals x and y decrease along with the reduction in numerical tolerance values, but chemical concentrations remain oscillatory even when the relative and absolute numerical tolerance values are reduced to 10-10 and 10-12, respectively (Fig. 1). It seems that the chemical concentrations will remain oscillatory unless the numerical tolerance values approach zero or are much smaller than those in Fig. 1C. A question then arises: which simulation outcome should be accepted in this case, an outcome with a certain range of numerical tolerance values far from zero or one with the numerical values approaching or close to zero?    We would argue that the answer to the above question is an outcome with the numerical tolerance values far from zero. We think that at the fundamental level, chemical concentrations cannot be precisely determined and the inherent variance of a chemical concentration is governed by the uncertainly principle. The uncertainty principle was typically applied to a subatomic particle. Because each molecule can be considered as a population of subatomic particles, the uncertainty level of such a population in position and energy state should be equal to or more than that of a subatomic particle, depending on the complexity of the molecule. The uncertainty level of a subatomic particle, therefore, can serve as a guidance value for the estimation of the lower limit of the uncertainty level of molecules. Precisely calculating the uncertainty level of a molecule, or a population of molecules in the case of determining the concentration of a molecular species, is challenging due to the complexity arising from the large number of subatomic particles involved and the dynamics of the molecule dictated by its internal structure and external environment. Towards finding the lower limit of the uncertainty level of a population of molecules, we start from the uncertainty level of a subatomic particle, which is typically presented as the following:  4hpx            (1)  .CC-BY-NC-ND 4.0 International licenseacertified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under The copyright holder for this preprint (which was notthis version posted January 29, 2018. ; https://doi.org/10.1101/256214doi: bioRxiv preprint 4  where x  and p are the standard deviations of the position and momentum of a particle, respectively, and h is the Planck constant. Now let’s consider a biochemical reaction as a representative case. The kinetics of a chemical reaction is typically described by an equation with variables of the concentrations of the reactant(s) and product(s), and the standard deviation of the concentration c  is, in effect, the relative numerical tolerance value, which exerts its influence on the simulation when the concentration values are reasonably large (e.g. >0.1). The standard deviation of the concentration of a chemical, c , can be calculated as the following:   nvwni vwiic2001)(           (2)  where 0w , 0v , and n  are the mean of the chemical weight measurement, mean of the volume measurement, and the number of the measurements, respectively.  To simply the analysis, it is assumed that the weight measurements are quite accurate, i.e., iw ≈ 0w , therefore,      nnivvvvwiic 1)( 2000          (3)  For a given shape of the volume,   3ii axv  .           (4)  Therefore,  nxxxxawni iic 123033030)( .         (5)  Because )(3)(3 02000303 xxxxxxxxx iiii   when 0xxi  ,   xniicaxwnxxaxw400120400 3)(3 .       (6)  For a subcellular compartment, assuming its shape is a sphere, the above equation is then simplified as .CC-BY-NC-ND 4.0 International licenseacertified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under The copyright holder for this preprint (which was notthis version posted January 29, 2018. ; https://doi.org/10.1101/256214doi: bioRxiv preprint 5   40049xw xc  .          (7)  It is proposed that in Eq. 1, when px   , it strikes a balance between the measurement accuracy of the position and that of the momentum of a particle. Under this condition,   4hx  .           (8)  Assuming that x in Eq. 7 is equivalent to x in Eq. 8 since they both represent a standard deviation of distance measurements, and based on the early justification that the uncertainty level of a population of molecules is equal to or more than that of a subatomic particle,   340089hxwc  .         (9)  To calculate values of c in Eq. 9, biologically relevant values of 0w and 0x should be determined. Based on the study by Milo and Phillips (2015), it is assumed that, in a subcellular compartment,180 109817.4 w g (the approximate weight of 100 molecules of a protein of 30000 Daltons), and 10 x  µm 610 m, then   13101928.8 c g m-3 4101928.8  ng/L.  The above results indicate that when the in vivo concentration values are expressed in ng/L (likely being reasonably large absolute numerical values) in simulation of biochemical reactions, the relative numerical tolerance can be reasonably set at or around 410  . For chemical concentrations expressed in a unit other than ng/L, the corresponding relative numerical tolerance values in simulation may be adjusted based on the above c value. The proposed range of relative numerical tolerance values suggests that the simulation results in Fig. 1A and B are likely closer to reality than that in Fig. 1C is.  Based on the above analysis, reducing relative numerical tolerance values for chemical concentration values in simulation will inevitably be met with increased standard deviation values in determining the energy (momentum) state of the molecules. This inverse relationship between the chemical concentrations and the momentum of the molecules likely means that the mathematical relationship among variables in a deterministic model will increasingly suffer inaccuracy with the increase in the accuracy of the concentration values. Furthermore, what is presented here may also be relevant in general to mathematical models that need to account for positions of molecules.     .CC-BY-NC-ND 4.0 International licenseacertified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under The copyright holder for this preprint (which was notthis version posted January 29, 2018. ; https://doi.org/10.1101/256214doi: bioRxiv preprint 6  References  Byrne, G. D., Thompson, S., 2013. Error control matters. http://www.radford.edu/~thompson/vodef90web/ErrorControlMatters.pdf.  DiStefano III, J., 2015. Dynamic Systems Biology Modeling and Simulation.  Elsevier, Waltham, MA, USA.  Milo, R., Phillips, B., 2016. Cell Biology by the Numbers. Garland Science, Taylor & Francis Group, New York, New York, USA.  Yang, L. Z., Yang, M., 2015. Modeling biological oscillations: integration of short reaction pauses into a stationary model of a negative feedback loop generates sustained long oscillations. PeerJ Preprint, doi:10.7287/peerj.preprints.1272v1.  Figure legend Fig. 1. Simulation outcomes of a negative feedback system of biochemical reactions at three numerical tolerance levels. This is a two-ODE system numerically integrated as described for Fig. 2B in Yang and Yang (2015). (a) Relative numerical tolerance = 10-3, and absolute numerical tolerance = 10-6. (b) Relative numerical tolerance = 10-5, and absolute numerical tolerance = 10-8. (c) Relative numerical tolerance = 10-10, and absolute numerical tolerance = 10-12.  .CC-BY-NC-ND 4.0 International licenseacertified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under The copyright holder for this preprint (which was notthis version posted January 29, 2018. ; https://doi.org/10.1101/256214doi: bioRxiv preprint .CC-BY-NC-ND 4.0 International licenseacertified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under The copyright holder for this preprint (which was notthis version posted January 29, 2018. ; https://doi.org/10.1101/256214doi: bioRxiv preprint 

Determining a reasonable range of relative numerical tolerance values for simulating deterministic models of biochemical reactions

https://shareok.org/bitstream/11244/336862/3/oksd_yang_determining_a_reasonable_range_2018.pdf

Determining a reasonable range of relative numerical tolerance values for simulating deterministic models of biochemical reactions

Abstract

Similar works

Full text

Available Versions

SHAREOK repository