http://chineseinput.net/에서 pinyin(병음)방식으로 중국어를 변환할 수 있습니다.
변환된 중국어를 복사하여 사용하시면 됩니다.
Reduction of multiscale stochastic biochemical reaction networks using exact moment derivation
Kim, Jae Kyoung,Sontag, Eduardo D. Public Library of Science 2017 PLoS computational biology Vol.13 No.6
<▼1><P>Biochemical reaction networks (BRNs) in a cell frequently consist of reactions with disparate timescales. The stochastic simulations of such multiscale BRNs are prohibitively slow due to high computational cost for the simulations of fast reactions. One way to resolve this problem uses the fact that fast species regulated by fast reactions quickly equilibrate to their stationary distribution while slow species are unlikely to be changed. Thus, on a slow timescale, fast species can be replaced by their quasi-steady state (QSS): their stationary conditional expectation values for given slow species. As the QSS are determined solely by the state of slow species, such replacement leads to a reduced model, where fast species are eliminated. However, it is challenging to derive the QSS in the presence of nonlinear reactions. While various approximation schemes for the QSS have been developed, they often lead to considerable errors. Here, we propose two classes of multiscale BRNs which can be reduced by deriving an exact QSS rather than approximations. Specifically, if fast species constitute either a feedforward network or a complex balanced network, the reduced model based on the exact QSS can be derived. Such BRNs are frequently observed in a cell as the feedforward network is one of fundamental motifs of gene or protein regulatory networks. Furthermore, complex balanced networks also include various types of fast reversible bindings such as bindings between transcriptional factors and gene regulatory sites. The reduced models based on exact QSS, which can be calculated by the computational packages provided in this work, accurately approximate the slow scale dynamics of the original full model with much lower computational cost.</P></▼1><▼2><P><B>Author summary</B></P><P>Molecules inside a cell undergo various transformations via biochemical reactions with disparate rates. For instance, while transcriptional factors bind and unbind gene promoters in a time scale of seconds, mRNA transcription takes at least several minutes. For such systems regulated by both fast and slow reactions together, most of the computation time in stochastic simulations is spent on simulating the fast reactions, even if our interest is in the dynamics of slow reactions such as transcription. This problem can be resolved by deriving a reduced system, which can accurately approximate the slow dynamics of the original system without simulating fast reactions. However, when nonlinear reactions exist, the accurate reduction is often impossible. Here, we describe that two classes of nonlinear systems, both of which appear frequently as models of natural biochemical reaction networks, can be effectively reduced, and we provide a computational package to carry out the reduction. We find that the resulting reduced systems accurately capture the stochastic dynamics of the original system, while saving considerable simulation time.</P></▼2>
King, J D,Strait, E J,Boivin, R L,Taussig, D,Watkins, M G,Hanson, J M,Logan, N C,Paz-Soldan, C,Pace, D C,Shiraki, D,Lanctot, M J,La Haye, R J,Lao, L L,Battaglia, D J,Sontag, A C,Haskey, S R,Bak, J G American Institute of Physics 2014 Review of scientific instruments Vol.85 No.8
<P>The DIII-D tokamak magnetic diagnostic system [E. J. Strait, Rev. Sci. Instrum. 77, 023502 (2006)] has been upgraded to significantly expand the measurement of the plasma response to intrinsic and applied non-axisymmetric '3D' fields. The placement and design of 101 additional sensors allow resolution of toroidal mode numbers 1 n 3, and poloidal wavelengths smaller than MARS-F, IPEC, and VMEC magnetohydrodynamic model predictions. Small 3D perturbations, relative to the equilibrium field (10(-5) < δB/B0 < 10(-4)), require sub-millimeter fabrication and installation tolerances. This high precision is achieved using electrical discharge machined components, and alignment techniques employing rotary laser levels and a coordinate measurement machine. A 16-bit data acquisition system is used in conjunction with analog signal-processing to recover non-axisymmetric perturbations. Co-located radial and poloidal field measurements allow up to 14.2 cm spatial resolution of poloidal structures (plasma poloidal circumference is ~500 cm). The function of the new system is verified by comparing the rotating tearing mode structure, measured by 14 BP fluctuation sensors, with that measured by the upgraded B(R) saddle loop sensors after the mode locks to the vessel wall. The result is a nearly identical 2/1 helical eigenstructure in both cases.</P>