
Citation: David K. Y. Lim, Peer M. Schenk. Microalgae selection and improvement as oil crops: GM vs non-GM strain engineering[J]. AIMS Bioengineering, 2017, 4(1): 151-161. doi: 10.3934/bioeng.2017.1.151
[1] | Aslak Tveito, Glenn T. Lines, Pan Li, Andrew McCulloch . Defining candidate drug characteristics for Long-QT (LQT3) syndrome. Mathematical Biosciences and Engineering, 2011, 8(3): 861-873. doi: 10.3934/mbe.2011.8.861 |
[2] | Virginia González-Vélez, Amparo Gil, Iván Quesada . Minimal state models for ionic channels involved in glucagon secretion. Mathematical Biosciences and Engineering, 2010, 7(4): 793-807. doi: 10.3934/mbe.2010.7.793 |
[3] | Hongtao Liu, Shuqin Liu, Xiaoxu Ma, Yunpeng Zhang . A numerical model applied to the simulation of cardiovascular hemodynamics and operating condition of continuous-flow left ventricular assist device. Mathematical Biosciences and Engineering, 2020, 17(6): 7519-7543. doi: 10.3934/mbe.2020384 |
[4] | Siqi Chen, Ran Su . An autonomous agent for negotiation with multiple communication channels using parametrized deep Q-network. Mathematical Biosciences and Engineering, 2022, 19(8): 7933-7951. doi: 10.3934/mbe.2022371 |
[5] | Jiuju Shen . EM-TSA: An ensemble machine learning-based transient stability assessment approach for operation of power systems. Mathematical Biosciences and Engineering, 2023, 20(5): 8226-8240. doi: 10.3934/mbe.2023358 |
[6] | Yue Li, Wusheng Xu, Wei Li, Ang Li, Zengjin Liu . Research on hybrid intrusion detection method based on the ADASYN and ID3 algorithms. Mathematical Biosciences and Engineering, 2022, 19(2): 2030-2042. doi: 10.3934/mbe.2022095 |
[7] | Cong Zhao, Na Deng . An actor-critic framework based on deep reinforcement learning for addressing flexible job shop scheduling problems. Mathematical Biosciences and Engineering, 2024, 21(1): 1445-1471. doi: 10.3934/mbe.2024062 |
[8] | T. J. Newman . Modeling Multicellular Systems Using Subcellular Elements. Mathematical Biosciences and Engineering, 2005, 2(3): 613-624. doi: 10.3934/mbe.2005.2.613 |
[9] | Jimin Yu, Jiajun Yin, Shangbo Zhou, Saiao Huang, Xianzhong Xie . An image super-resolution reconstruction model based on fractional-order anisotropic diffusion equation. Mathematical Biosciences and Engineering, 2021, 18(5): 6581-6607. doi: 10.3934/mbe.2021326 |
[10] | Na Pang . Nonlinear neural networks adaptive control for a class of fractional-order tuberculosis model. Mathematical Biosciences and Engineering, 2023, 20(6): 10464-10478. doi: 10.3934/mbe.2023461 |
Since advanced computing technology makes the computational models of cardiac electrophysiology getting mature, the computational models become an important basis for clinical applications [1,2]. Moreover, cardiac cells up to an order of one billion and the internal processes of each cell are also quite complex in the human heart [3]. Therefore, the computational speed, accuracy, and stability are key hotspots in whole heart simulation. Firstly, advanced computer parallelization greatly reduces computing cost of tissue-scale cardiac simulations [4,5,6,7,8]. Secondly, the adaptive algorithms of time-steps and space-steps are a powerful tool for improving the accuracy of computation [6,9,10]. Thirdly, stiffness is ubiquitous in cardiac cell models, which easily leads to unstable computation [11]. Aiming to solve the stiff model of a single cell, the advanced numerical methods of ODEs system, Exponential Adams-Bashforth integrators [12] and Gating-enhanced implicit-explicit linear multistep methods [13] were proposed. In addition, Sundials CVODE [14] is a useful solver for ODEs system. The backward-differentiation-formula (BDF) method in CVODE with low tolerances can obtain a highly converged solution for stiff problems. Due to linear multistep methods consuming high computation cost for high precision solutions, using one-step methods to study the tradeoff between speedup and stability has been a new hot-spot in recent years.
The stability of stiff Markov model for ion channels is our main concern. In recent years, Markovian model for ion channels has been widely used in cardiac arrhythmia [15,16,17] and pharmacologic screening for anti-arrhythmic drugs [18,19]. The introduction of Markov model makes the cell model become quite stiff, resulting in unstable simulation as large time-step is used to accelerate the simulation. Therefore, the confrontation between the computational speed and stability of a single cell is a noticeable problem. To ensure high simulation stability, Johnny Moreira et al. [20] compared UNI with Forward Euler (FE), Rush-Larsen, First-Order Sundnes et al. (SAST1), Second-Order Sundnes et al. (SAST2), and Second-Order Runge-Kutta (RK2) methods and found that UNI could substantially increase the stability of simple first-order (fixed time-step) methods but with very low speedup. Stary et al. [21] proposed the HOS method in comparison to FE and Matrix Rush-Larsen (MRL) methods to show that HOS method could also increase stability with moderate speedup. While, CCL algorithm [22] has been verified to be an efficient method for the stability of the action potential simulation. And our previous work [23] confirmed that CCL algorithm combined with HOS or UNI method could save markedly computation cost and increased the stability of adaptive methods for time-step from 0.001 to 1 millisecond.
For the simulation of cardiac action potential, Marsh et al. [24] proposed that Mixed Root Mean Square(MRMS) errors should be less than 5% to meet the requirement of clinical applications. The simulation results by Chen et al. [23] even with high speed could not meet the accuracy requirement of clinical applications (less than 5% errors). Due to CCL algorithm only adjusts time-steps according to the first and second derivatives of membrane potential, during the period of slow voltage change (from phase 3 to phase 4), the time-step increases rapidly, causing simulation instability of the state variables for the WT-9 model as the state variables still rapidly change their values, resulting in over 5% MRMS errors at simulations of membrane potential and sodium current. For instance, using fixed time-step HOS method as a reference, CCL+HOS solver can get 295 times speedup, but with 22% and 6.9% MRMS errors for potential and sodium current, respectively. The big errors come from large time-step causing the state occupancies of WT-9 model to be simulated unstably at fast inactivation, closed-inactivation and the close states. To eliminate this problem, we propose a multi-variable CCL method that takes the first and second derivatives of membrane potential and state occupancies as input parameters, and then generates the appropriate time-step based on the quadratic equation solution of time-step.
The Markov model of the Wild-Type (WT) sodium channel [18,25] is integrated into the dynamic Luo-Rudy (LRd) model [26]. To simulate the ventricular action potential, the formula is given by,
dV/dt=−(1/Cm)(Iion+Ist) | (2.1) |
Iion=INa+INa,b+ICa(L),Na+Ins,Na+3INaK+3INaCa+IK+IKi+IKp+ICa(L),K−2INaK+Ins,K+ICa(L)+ICa,b+Ip(Ca)−2INaCa | (2.2) |
where Ist is an external stimulus current, with time-step 0.01 or 0.001 ms with duration 0.5 ms for WT-9 or WT-8 model, respectively. INa,INa,b,ICa(L),Na,Ins,Na,3INaK,3INaCa are sodium channel currents. IK,IKi,IKp,ICa(L),K,2INaK,Ins,K are potassium channel currents. ICa(L),ICa,b,Ip(Ca),2INaCa are calcium channel currents. For the detailed descriptions, please see LRd model [26].
The WT-9 Markovian model of the sodium channel includes nine states [25]. Each state is represented by the proportion of channels with respect to the total number of channels called state occupancy. A Markovian model of the WT sodium channel (Figure 1) is composed of one conducting open state (O), three closed states (C1,C2,C3), two closed-inactivation states (IC3,IC2), one fast inactivation state(IF) and two intermediate inactivation states (IM1,IM2). According to the principle of microscopic reversibility[27], the change in one state occupancy with continuous-time is equal to the net flow for the transition state. This change can be represented by an ordinary differential equation (ODE). The WT-9 Markovian sodium channel has nine equations that form a linear system known as Kolmogorov equations.
d→u/dt=Q(V)→u | (2.3) |
where Q(V) is the 9×9 matrix of transition rates (af11,af12,af13,af2,af3,af4,af5,bt11,bt12,bt13,bt2,bt3,bt4,bt5). Q(V) is a real matrix of nine rows by nine columns, and it contains elements qij. →uis a vector with nine states [O,C1,C2,C3,IC3,IC2,IF,IM1,IM2]T.
The formulation of the sodium channel is described by,
INa=¯gNa×O×(V−ENa) | (2.4) |
where ¯gNa is the conductance constant of the sodium channel and ENa is the equilibrium potential.
A WT-8 Markovian model with eight states of the sodium channel [18] is composed of one conducting open state (O), three closed states (C1,C2,C3), two closed-inactivation states (IC3,IC2), one fast inactivation state(IF) and slow inactivation state (IS) shown in Figure 2.
Q(V) in Eq (2.3) is the 8×8 matrix of transition rates (af11,af12,af13,af2,af3,afx,bt11,bt12,bt13,bt2,bt3,btx). Q(V) is a real matrix of eight rows by eight columns, and it contains elements qij. The vector →u includes eight states [O,C1,C2,C3,IC3,IC2,IF,IS]T.
CCL method is applied to change time-step with membrane potential in simulations. The maximum time-step Δtmax=1 ms and the minimum time-step Δtmin=0.001 ms. Time-step is changed in the scope of Δtmin∼Δtmax. To create a fine adaptive time-step, the new approximated membrane potential with second-order Taylor expansion is given by, V(tn+1)≅V(tn)+dV/dt×Δt+1/2×d2V/dt2×Δt2. Which can be a quadratic Δt equation used to update time-steps adaptively. For a detailed description, please see Chen et al. [22]. The equation is given by,
1/2×a×Δt2+b×Δt−c≅0,where a=d2V/dt2,b=dV/dt,c=V(tn+1)−V(tn)=±0.1mV | (2.5) |
Let D=b2+4×a/2×c, The solutions of Eq (2.5) are shown,
Δt=(−b+√D)/a, if b≥0,D≥0, where Δt∈(0.001∼1)ms,c=0.1mV | (2.6) |
Δt=(−b−√D)/a, if b<0,D≥0, where Δt∈(0.001∼1)ms,c=−0.1mV | (2.7) |
Δt=−b/a, if D<0, where Δt∈(0.001∼1)ms | (2.8) |
Δtnew=Minimum(Δt,Δtmax), where Δtmax=Minimum(2×Δtold,Δtmax) | (2.9) |
In addition to membrane potential, the variables control in MCCL method is expanded to multiple variables given by,
1/2×a×Δt2i+b×Δti−c≅0, where a=d2Xi/dt2,b=dXi/dt | (2.10) |
c is the offset of variable Xi, Xi is one of the variables Vn,O,C1,C2,C3,IC3,IC2,IF,IS,IM1,IM2 at ith time. Vn is the normalized membrane potential given by, Vn=(V−Vrest)/NV,where NV=|Vrest−Vreversal|. Vrest is the membrane potential of rest state, Vreversal is the reversal potential of sodium channel. Vrest=−88.65 mV and Vreversal=70.50 mV are set respectively. Let c=Xi(tn+1)−Xi(tn)=±0.1/NV and D=b2+4×a/2×c, then The solutions of Eq (2.10) are shown,
Δti=(−b+√D)/a, if b≥0,D≥0,where Δti∈(Δtmin∼Δtmax)ms,c=0.1/NV | (2.11) |
Δti=(−b−√D)/a, if b<0,D≥0,where Δti∈(Δtmin∼Δtmax)ms,c=−0.1/NV | (2.12) |
Δti=−b/a, if D<0,where Δti∈(Δtmin∼Δtmax)ms | (2.13) |
Δti=Minimum(Δti,Δtmax),where Δtmax=Minimum(2×Δtold,Δtmax) | (2.14) |
For WT−9 model,Δtnew=Minimum(ΔtVn,ΔtO,ΔtC1,ΔtC2,ΔtC3,ΔtIC3,ΔtIC2,ΔtIF,ΔtIM1,ΔtIM2) | (2.15) |
Or, for WT−8 model,Δtnew=Minimum(ΔtVn,ΔtO,ΔtC1,ΔtC2,ΔtC3,ΔtIC3,ΔtIC2,ΔtIF,ΔtIS) | (2.16) |
Four methods [23] are mainly adopted to solve transient solutions of Eq (2.3) as shown below,
1) Forward Euler method ("FE") and FE method is given by,
→u(tn+Δt)=→u(tn)+ΔtQ(V(tn))→u(tn) | (2.17) |
2) Taylor series expansion ("TAL") and second-order Taylor series method is given by,
→u(tn+Δt)=→u(tn)+ΔtQ(V(tn))1!→u(tn)+(ΔtQ(V(tn)))22!→u(tn) | (2.18) |
3) Uniformization ("UNI") and the uniformization method by Sidje et al.[28] is given by,
→u(tn+Δt)=exp(ΔtQ(V(tn)))→u(tn)≈→u(tn)+ΔtQ(V(tn))1!→u(tn)+⋅⋅⋅+(ΔtQ(V(tn)))ll!→u(tn) | (2.19) |
→u(tn+Δt)=exp(qΔtQ(V(tn))∗−I)→u(tn)=exp(−qΔt)exp(qΔtQ(V(tn))∗)→u(tn) | (2.20) |
ˉ→u(tn+Δt)=l∑k=0exp(−qΔt)(qΔt)kk!(Q(V(tn))∗)k→u(tn) | (2.21) |
l∑k=0exp(−qΔt)(qΔt)kk!≥1−εtol | (2.22) |
ˉ→u(tn+Δt)=(exp(−qΔt)exp(qΔtQ(V(tn))∗))m→u(tn) | (2.23) |
where m=⌈qΔtθ⌉,h=Δtm, Q(V(tn))∗=I+Q(V(tn))q,q=max|qii|. q is the maximum diagonal element, in absolute value, of matrix Q. εtol denotes the prescribed error tolerance, θ and m are parameters to partition the integral domain for overflow prevention. Two parameters are set in advance, the default εtol is set 10−11 and θ=1 unless otherwise specified. The truncation point l is determined by Eq (2.22).
4) Hybrid Operator Splitting ("HOS") and HOS method is given by,
Q(V(tn))=Q0(V(tn))+Q1(V(tn))+Q2(V(tn)) | (2.24) |
→un+1/3=exp(ΔtQ0(V(tn)))→un | (2.25) |
→un+2/3=exp(ΔtQ1(V(tn)))→un+1/3 | (2.26) |
→un+1=exp(ΔtQ2(V(tn)))→un+2/3 | (2.27) |
For a detailed description, please see Stary et al. [21].
To compare the performance and accuracy between CCL and MCCL methods combined with Markov solver, we have tested eight hybrid solvers, in which each of four methods (FE, TAL, UNI, HOS) is combined with CCL or MCCL method as listed in Table 1. The reference solver is fourth-order four-stage Runge-Kutta (RK4) method used to solve Markov model with fixed time-step (Δt=0.001 ms). For the very sharply stiff model, BDF (backward-differentiation formula) with tight tolerance and fixed time-step 0.001 ms in Sundials CVODE as a benchmark is an alternative good choice due to high stability even with high computation cost. RK4 with fixed step 0.001 ms has less than 0.6% MRMSE deviation and 10 times smaller computation cost in comparison to BDF with tight tolerance and fixed time-step 0.001 ms in Sundials CVODE.
Solver | Time-steps | Solver | Time-steps |
CCL+FE | (0.001~1 ms) | MCCL+FE | (0.001~1 ms) |
CCL+TAL | (0.001~1 ms) | MCCL+TAL | (0.001~1 ms) |
CCL+UNI | (0.001~1 ms) | MCCL+UNI | (0.001~1 ms) |
CCL+HOS | (0.001~1 ms) | MCCL+HOS | (0.001~1 ms) |
We implement the numerical experiments with Visual Studio 2013 C/C++ compiler. All the numerical tests are performed on a desktop computer equipped with an Intel(R) Core (TM) i7-7700 CPU running 64-bit Windows10 × 64. The fixed time-step scheme with RK4 is used as reference solution to evaluate the performance of adaptive time-step schemes with Markov solvers. The Mixed Root Mean Square Error (MRMSE) [24] is used as a metric of quantitative error. It is given by,
MRMSE≐√N∑i=1(Xi−ˆXi1+|ˆXi|)2N | (2.28) |
Where Xi is the numerical solution, ˆXi is the reference solution, N is the number of temporal points. All of the numerical schemes run 3000 basic cycle length (BCL) or beats with 1000 ms per BCL or beat. The simulation data are taken from the 1st, 500th, 1000th, 1500th, 2500th, and 3000th beats.
Table 2 shows the MRMSE and computational performance of the eight schemes for the WT-9 model.
Solvers | MRMSE-V | MRMSE-I | Speedup | Cost(ms) |
DT+RK4 | 0.0 | 0.0 | 1.0X | 1.0×107±9.8×104 |
CCL+FE | 0.016 | 0.588 | 109.9X | 9.5×104±4.0×102 |
CCL+TAL | 0.745 | 1.158 | 33.0X | 3.2×105±1.2×103 |
CCL+UNI | 0.056 | 0.012 | 16.4X | 6.4×105±8.4×102 |
CCL+HOS | 0.236 | 0.064 | 412.3X | 2.6×104±3.2×102 |
MCCL+FE | 0.022 | 0.065 | 96.9X | 1.1×105±2.7×102 |
MCCL+TAL | 0.015 | 0.124 | 29.8X | 3.5×105±1.4×103 |
MCCL+UNI | 0.021 | 0.009 | 14.8X | 7.0×105±3.4×103 |
MCCL+HOS | 0.006 | 0.002 | 17.2X | 6.1×105±3.7×103 |
Note: MRMSE-V: MRMSE of voltage in comparison to the reference. MRMSE-I: MRMSE of sodium current in comparison to the reference. Speedup: The computation cost ratio of one of eight schemes to the reference solution DT+RK4. X, times speedup. Ten simulations for each solver. Cost in Mean ± Standard deviation. |
To meet the requirement of clinical applications, Marsh et al. [24] suggested that MRMSE of simulations should be less than 5%. Table 2 shows that all of CCL combined any Markov solver (FE, TAL, UNI, HOS) do not satisfy the requirement, nor MCCL+FE or MCCL+TAL solver. Only MCCL+UNI and MCCL+HOS solvers can meet the 5% requirement. Particularly, MCCL+HOS solver has less than 1% MRMSEs in voltage and sodium current with 17.2 times (17.2X) speedup which is greater than 14.8X for MCCL+UNI solver. It is suggested that the optimal solver is MCCL+HOS solver among eight solvers in Table 2 for the WT-9 model.
Figure 3 shows the long-term stability (3000 beats) of CCL method or MCCL method combined Markov solvers for Vmax and APD. Figure 3a, e show that both Vmax and APD of MCCL+FE solver have some unstable fine fluctuation, but they are closer to the reference solution than those of CCL+FE solver. Figure 3b, f show that Vmax and APD of MCCL+TAL solver are close to the reference solution with tiny unstable fluctuations, while Vmax and APD of CCL+TAL solver not only have some tiny unstable oscillation but also deviates far from the reference solution. Importantly, Figure 3c, g show that Vmax of CCL+UNI and that of MCCL+UNI solver have a slight error, but remain stable. Figure 3d, h show that the Vmax and APD of MCCL+HOS solver are closer to the reference solution than that of CCL+HOS solver. It is indicated that MCCL method can guarantee long-term stability and obtain more accurate solutions about Vmax and APD than CCL method.
Figure 4 shows the occupancy of nine states of Markov model for CCL or MCCL combined HOS solver at first beat. Figure 4a, c, f–i show that all these state occupancies (O,C1,C3,IC3,IM1,IM2) are close to the reference solution very well for CCL+HOS and MCCL+HOS solvers. While, in Figure 4b shows that the state occupancy IF of MCCL+HOS solver has tiny unstable fluctuation (inset in Figure 4b), but IF of CCL+HOS solver shows a small deviation near 200 ms. Similarly, in Figure 4d, e state C2 and IC2 occupancies of MCCL+HOS solver only have very minor fluctuation (insets in Figure 4d, e), while that of CCL+HOS solver presents a deviation error with rapidly descending towards abscissa near 200 ms. It is suggested that MCCL+HOS solver enhances the instability of CCL+HOS solver.
Table 3 shows the MRMSE and computational performance of the eight simulate schemes for the WT-8 model.
Solvers | MRMSE-V | MRMSE-I | Speedup | Cost(ms) |
DT+RK4 | 0.0 | 0.0 | 1.0X | 1.1×107±6.0×104 |
CCL+FE | - - | - - | - - | failure |
CCL+TAL | - - | - - | - - | failure |
CCL+UNI | 0.019 | 0.008 | 1.2X | 8.9×106±1.8×104 |
CCL+HOS | 0.250 | 0.114 | 266.6X | 3.3×104±1.3×102 |
MCCL+FE | 0.103 | 3.345 | 19.8X | 5.3×105±2.5×103 |
MCCL+TAL | - - | - - | - - | failure |
MCCL+UNI | 0.076 | 0.042 | 1.2X | 8.9×106±1.1×104 |
MCCL+HOS | 0.032 | 0.013 | 21.1X | 5.0×105±1.8×103 |
CCL+FE, CCL+TAL and MCCL+TAL solvers are the failure, while MCCL+FE solver is success but with MRMSE higher than 5%. Only CCL+UNI and MCCL+HOS solvers have MRMSE-V and MRMSE-I less than 5%. CCL+UNI solver almost takes no speedup, while MCCL+HOS solver still obtains 21.1 times speedup. For WT-8 model, MCCL+HOS solver is the best solver among eight solvers in Table 3.
Due to an enormous maximum time-step size, Δtmax, used by CCL and MCCL algorithms, the simulation is prone to instability and even divergence to failure for CCL and MCCL combined with FE or TAL, especially for WT-8 model. It is improper to use FE or TAL combined with adaptive time-step methods unless a smaller maximum time-step size, Δtmax, is used but then the adaptivity of the method is limited.
In Tables 2 and 3, the speedup performance of MCCL+UNI and CCL+UNI solvers is similar since the UNI algorithm consumes most of the computation time. The formula (2.23) is iterated 260 times with q=260 for the WT-8 model to reach the preset tolerant error even time-step is set as 1ms. As a result, there is almost no speedup of MCCL or CCL combined with UNI for WT-8 model.
While, the speedup of MCCL+HOS and CCL+HOS solvers is significant due to HOS being lightweight. CCL+HOS solver gets high over 250 times speedup for WT-8 and WT-9 models while MCCL+HOS only gets about 20 times. However, CCL+HOS solver has large MRMSE for WT-9 model and even 25% MRMSE for WT-8 failing the clinical application requirement. While MCCL+HOS solver has MRMSEs less than 5% clinical application requirement even with lower speedup.
In our previous work [23], we find that CCL+UNI solver is the most stable but very time-consuming while CCL+HOS solver is less stable but more efficient. It is suggested the optimal solver should take tradeoff between stability and acceleration. In this work, the MCCL+HOS solver is more stable, more accurate and faster than MCCL+UNI solver. For the WT-8 model, CCL+UNI solver is the most accurate but has poor speedup performance (only 1.2X). It is concluded that MCCL+HOS solver is a good solver, especially for high stiff Markov model such as WT-8 sodium channel and MCCL instead of CCL can be a good speedup algorithm tradeoff for stability against acceleration.
Both CCL and MCCL have their limitations in which the speedup of MCCL+HOS solver is far less than CCL+HOS while the accuracy of CCL+HOS is far worse than MCCL+HOS. To achieve better accuracy, MCCL+HOS is highly recommended but its speedup performance is only moderate and, importantly, the speedup performance is about the same for moderate stiff WT-9 and high stiff WT-8 models. It is highly recommended to use MCCL+HOS for a high stiff model, but for moderate stiff model MCCL+HOS should be revised to gain higher speedup for the logic of lower stiff higher speedup performance in the near future.
This work was supported by Sun Yat-sen University, China, under Scientific Initiation Project [No.67000-18821109] for High-level Experts and by the Ministry of Science and Technology, Taiwan, under grant 107-2115-M-006-013 separately.
The authors declare that they have no conflict of interest.
[1] |
Schenk PM, Thomas-Hall SR, Stephens E, et al. (2008) Second generation biofuels: high-efficiency microalgae for biodiesel production. Bioenerg Res 1: 20-43. doi: 10.1007/s12155-008-9008-8
![]() |
[2] |
Singh B, Guildhe A, Rawat I, et al. (2014) Towards a sustainable approach for development of biodiesel from plant and microalgae. Renew Sust Energ Rev 29: 216-245. doi: 10.1016/j.rser.2013.08.067
![]() |
[3] |
Razeghifard R (2013) Algal biofuels. Photosynth Res 117: 207-219. doi: 10.1007/s11120-013-9828-z
![]() |
[4] |
Chisti Y (2008) Biodiesel from microalgae beats bioethanol. Trends Biotechnol 26: 126-131. doi: 10.1016/j.tibtech.2007.12.002
![]() |
[5] |
Dermirbas A, Dermirbas MF (2011) Importance of algae oil as a source of biodiesel. Energ Convers Manag e 52: 163-170. doi: 10.1016/j.enconman.2010.06.055
![]() |
[6] |
Duong VT, Li Y, Nowak E, et al. (2012) Microalgae isolation and selection for prospective biodiesel production. Energies 5: 1835-1849. doi: 10.3390/en5061835
![]() |
[7] | Davey HM, Kell DB (1996) Flow cytometry and cell sorting of heterogenous microbial populations: the importance of single-cel analyses. Microbiol Res 60: 641-696. |
[8] | Mutanda T, Ramesh D, Karthukeyan S, et al. (2011) Bioprospecting for hyper-lipid producing microalgal strains for sustainable biofuel production. Bioresour Technol 102: 57-70. |
[9] | Sieracki M, Poulton N, Crosbie N (2005) Automated isolation techniques for microalgae. Algal Cult Tech, 101-106. |
[10] | Mendoza H, De La JA, Presmanes KF, et al. (2008) Characterization of dunaliella salina strains by flow cytometry: a new approach to select carotenoid hyperproducing strains. Electron J Biotech nol 11: 5-6. |
[11] |
Larkum AWD, Ross IL, Kruse O, et al. (2012) Selection, breeding and engineering of microalgae for bioenergy and biofuel production. Trends Biotechnol 30: 198-205. doi: 10.1016/j.tibtech.2011.11.003
![]() |
[12] |
Sharma KS, Schuhmann H, Schenk PM (2012) High lipid induction in microalage for biodiesel production. Energies 5: 1532-1553. doi: 10.3390/en5051532
![]() |
[13] |
Mata TM, Martins AA, Caetano NS (2010) Microalgae for biodiesel production and other applications: a review. Renew Sust Energ Rev 14: 217-232. doi: 10.1016/j.rser.2009.07.020
![]() |
[14] |
Chisti Y (2007) Biodiesel from microalgae. Biotechnol Adv 25: 294-306. doi: 10.1016/j.biotechadv.2007.02.001
![]() |
[15] |
Huerlimann R, De Nys R, Heimann K (2010) Growth, lipid content, productivity, and fatty acid composition of tropical microalgae for scale-up production. Biotechnol Bioeng 107: 245-257. doi: 10.1002/bit.22809
![]() |
[16] |
Lim DKY, Garg S, Timmins M, et al. (2012) Isolation and evaluation of oil-producing microalgae from subtropical coastal and brackish waters. Plos One 7: e40751-e40763. doi: 10.1371/journal.pone.0040751
![]() |
[17] |
Rodolfi L, Zittelli GC, Bassi N, et al. (2009) Microalgae for oil: strain selection, induction of lipid synthesis and outdoor mass cultivation in a low-cost photobioreactor. Biotechnol Bioeng 102: 100-112. doi: 10.1002/bit.22033
![]() |
[18] |
Mendoza H, De La JA, Presmanes KF, et al. (2012) Quick estimation of intraspeciific variation of fatty acid composition in Dunaliella salina using flow cytometry and nile red. J Adminppl Phycol 24: 1237-1243. doi: 10.1007/s10811-011-9768-y
![]() |
[19] | Nascimento IA, Marques SSI, Cabanelas ITD, et al. (2012) Screening microalgae strains for biodiesel production: lipid productivity and estimation of fuel quality based on fatty acids profiles as selective criteria. Bioenerg Res 6: 1-13. |
[20] | Christenson L, Sims R (2011) Production and harvesting of microalgae for wastewater treatment, biofuels, and bioproducts. Biotechnol Adv 29: 686-702. |
[21] |
Sharma KS, Garg S, Li Y, et al. (2013) Critical analysis of current microalgae dewatering techniques. Biofuels 4: 397-404. doi: 10.4155/bfs.13.25
![]() |
[22] |
Chen W, Zhang C, Song L, et al. (2009) A high throughput nile red method for quantitative measurement of neutral lipids in microalgae. J Microbiol Meth 77: 41-71. doi: 10.1016/j.mimet.2009.01.001
![]() |
[23] | Montero MF, Aristizabal M, Reina GG (2011) Isolation of high-lipid content strains of the marine microala tetraselmis suecica for biodiesel prodction by flow cytometry and single-cell sorting. J App l Phycol 23: 1053-1057. |
[24] | Anandarajah K, Mahendraperumal G, Sommerfeld M, et al. (2012) Characterization of microalga Nannochloropsis sp. mutants for improved production of biofuels. Appl Energ 96: 371-377. |
[25] |
Vigeolas H, Duby F, Kaymak E, et al. (2012) Isolation and partial characterization of mutants with elecated lipid content in Chlorella sorokiana and Scenedesmus obliquus. J Biotechnol 162: 3-12. doi: 10.1016/j.jbiotec.2012.03.017
![]() |
[26] | Alonso DL, Castillo CIS, Grima EM, et al. (1996) First insights into improvement of eicosapentaenoic acid content in Phaeodactylum tricornutum (bacillariophyceae) by induced mutagenesis. J Phycol 32: 553-558. |
[27] |
Chaturvedi R, Fujita Y (2006) Isolation of enhanced eicosapentanoic acid producing mutant of Nannochloropsis oculata ST-6 using ethyl methane sulfonate induced mutagenesis techniques and their characterization at mRNA transcript level. Phycol Res 54: 208-219. doi: 10.1111/j.1440-1835.2006.00428.x
![]() |
[28] |
Chaturvedi R, Uppalapati SR, Alamsjah MA, et al. (2004) Isolation of quizalofop-resistant mutants of Nannochloropsis oculata (eustigmatophyceae) with high eicosapentanoic acid following N-methyl-N-nitrosourea-induced mutagenesis. J Appl Phycol 16: 135-144. doi: 10.1023/B:JAPH.0000044826.70360.8e
![]() |
[29] |
Meireles LA, Guesdes AC, Malcata FX (2003) Increase of the yields of eicosapentanoic and docosahexanoic acids by the microalga Pavlova lutheri following random mutagenesis. Biotechnol Bioeng 81: 50-55. doi: 10.1002/bit.10451
![]() |
[30] |
Nymark M, Sharma AK, Sparstad T, et al. (2016) A CRISPR/Cas9 system adapted for gene editing in marine algae. Sci Rep 6: 24951-24956. doi: 10.1038/srep24951
![]() |
[31] |
Schuhmann H, Lim DKY, Schenk PM (2012) Perspectives on metabolic engineering for increased lipid contents in microalgae. Biofuels 3: 71-86. doi: 10.4155/bfs.11.147
![]() |
[32] |
Rismani-Yazdi H, Haznedaroglu BZ, Bibby K, et al. (2011) Transcriptome sequencing and annotation of the microalgae Dunaliella tertiolecta: pathway description and gene discovery for production of next-generation biofuels. BMC Genomics 12: 148-164. doi: 10.1186/1471-2164-12-148
![]() |
[33] |
Lei A, Chen H, Shenhu Z, et al. (2012) Expression of fatty acid synthesis genes and fatty acid accumulation in Haematococcus pluvialis under different stressors. Biotechnol Biofuels 5: 18-28. doi: 10.1186/1754-6834-5-18
![]() |
[34] |
Valenzuela J, Mazurie A, Carlson RP, et al. (2012) Potential role of multiple carbon fixation pathways during lipid accumulation in Phaeodactylum tricornutum. Biotechnol Biofuels 5: 40-56. doi: 10.1186/1754-6834-5-40
![]() |
[35] |
Rismani-Yazdi H, Haznedaroglu BZ, Hsin C, et al. (2012) Transcriptomic analysis of the oleaginous microalga Neochloris oleoabundans reveals metabolic insights into triacylglyceride accumulation. Biotechnol Biofuels 5: 1-16. doi: 10.1186/1754-6834-5-1
![]() |
[36] |
Guarnieri MT, Nag A, Smolinski SL, et al. (2011) Examination of triacylglycerol biosynthetic pathways via de novo transcriptomic and proteomic analyses in an unsequenced microalga. Plos One 6: e25851-e25863. doi: 10.1371/journal.pone.0025851
![]() |
[37] |
Lv H, Qu G, Qi X, et al. (2013) Transcriptome analysis of Chlamydomonas reinhardtii during the process of lipid accumulation. Genomics 101: 229-237. doi: 10.1016/j.ygeno.2013.01.004
![]() |
[38] |
Miller R, Wu G, Deshpande RR, et al. (2010) Changes in transcript abundance in Chlamydomonas reinhardtii following nitrogen deprivation predict diversion of metabolism. Plant Physiol 154: 1737-1752. doi: 10.1104/pp.110.165159
![]() |
[39] | Dunhay TG, Jarvis EE, Dais SS, et al. (1996) Manipulation of microalgal lipid production using genetic engineering. Appl Biochem Biotechnol 58: 223-231. |
[40] |
Deng X, Cai J, Fei X (2013) Effect of the expression and knockdown of citrate synthase gene on carbon flux during triacylglycerol biosynthesis by green algae Chlamydomonas reinhardtii. BMC Genomics 14: 38-48. doi: 10.1186/1471-2164-14-38
![]() |
[41] |
Niu Y, Zhang M, Li D, et al. (2013) Effect of the expression and knockdown of citrate synthase gene on carbon flux during triacylglycerol biosynthesis by green algae Chlamydomonas reinhardtii. Mar Drugs 11: 4558-4569. doi: 10.3390/md11114558
![]() |
[42] |
Mussgnug JH, Thomas-Hall SR, Rupprecht J, et al. (2007) Engineering photosynthetic light capture: impacts on improved solar energy to biomass conversion. Plant Biotechnol J 5: 802-814. doi: 10.1111/j.1467-7652.2007.00285.x
![]() |
[43] |
Beckmann J, Lehr F, Finazzi G, et al. (2009) Improvement of light to biomass conversion by de-regulation of light-harvesting protein translation in Chlamydomonas reinhardtii. J Biotechnol 142: 70-77. doi: 10.1016/j.jbiotec.2009.02.015
![]() |
[44] | Bougaran G, Rouxel C, Dubois N, et al. (2012) Enhancement of neutral lipid productivity in the microalga Isochrysis affinis galbana (T-Iso) by a mutation-selection procedure. Biotechnol Bioeng 11: 2737-2745. |
[45] | Doan YTT, Obbard JP (2011) Enhanced intracellular lipid in Nannochloropsis sp. via random mutagenesis and flow cytometric cell sorting. Glob Change Biol Bioener g 3: 264-270. |
[46] | Kamath BE, Vidhyavathi R, Sarada R, et al. (2008) Enhancement of carotenoids by mutation and stress induced carotenogenic genes in Haematococcus pluvialis mutants. Bioresour Technol 99: 8867-8673. |
[47] | Lian M, Huang H, Ren L, et al. (2010) Increase of docosahexaenoic acid production by Schizochytrium sp. through mutagenesis and enzyme assay. Appl Microbiol Biot 162: 935-941. |
[48] | Lim D, Schuhmann H, Sharma K, et al. (2015) Isolation of high-lipid Tetraselmis suecica strains following repeated UV-C mutagenesis, FACS, and high-throughput growth selection. Bio energ Res 8: 750-759. |
1. | Ning He, Changyou Tang, Xianchao Zhou, Shuolin Feng, M Pallikonda Rajasekaran, Diagnostic Value of Magnetic Resonance Image Feature Analysis under Reconstruction Algorithm for Knee Epiphyseal Injury, 2022, 2022, 1875-919X, 1, 10.1155/2022/1783975 |
Solver | Time-steps | Solver | Time-steps |
CCL+FE | (0.001~1 ms) | MCCL+FE | (0.001~1 ms) |
CCL+TAL | (0.001~1 ms) | MCCL+TAL | (0.001~1 ms) |
CCL+UNI | (0.001~1 ms) | MCCL+UNI | (0.001~1 ms) |
CCL+HOS | (0.001~1 ms) | MCCL+HOS | (0.001~1 ms) |
Solvers | MRMSE-V | MRMSE-I | Speedup | Cost(ms) |
DT+RK4 | 0.0 | 0.0 | 1.0X | 1.0×107±9.8×104 |
CCL+FE | 0.016 | 0.588 | 109.9X | 9.5×104±4.0×102 |
CCL+TAL | 0.745 | 1.158 | 33.0X | 3.2×105±1.2×103 |
CCL+UNI | 0.056 | 0.012 | 16.4X | 6.4×105±8.4×102 |
CCL+HOS | 0.236 | 0.064 | 412.3X | 2.6×104±3.2×102 |
MCCL+FE | 0.022 | 0.065 | 96.9X | 1.1×105±2.7×102 |
MCCL+TAL | 0.015 | 0.124 | 29.8X | 3.5×105±1.4×103 |
MCCL+UNI | 0.021 | 0.009 | 14.8X | 7.0×105±3.4×103 |
MCCL+HOS | 0.006 | 0.002 | 17.2X | 6.1×105±3.7×103 |
Note: MRMSE-V: MRMSE of voltage in comparison to the reference. MRMSE-I: MRMSE of sodium current in comparison to the reference. Speedup: The computation cost ratio of one of eight schemes to the reference solution DT+RK4. X, times speedup. Ten simulations for each solver. Cost in Mean ± Standard deviation. |
Solvers | MRMSE-V | MRMSE-I | Speedup | Cost(ms) |
DT+RK4 | 0.0 | 0.0 | 1.0X | 1.1×107±6.0×104 |
CCL+FE | - - | - - | - - | failure |
CCL+TAL | - - | - - | - - | failure |
CCL+UNI | 0.019 | 0.008 | 1.2X | 8.9×106±1.8×104 |
CCL+HOS | 0.250 | 0.114 | 266.6X | 3.3×104±1.3×102 |
MCCL+FE | 0.103 | 3.345 | 19.8X | 5.3×105±2.5×103 |
MCCL+TAL | - - | - - | - - | failure |
MCCL+UNI | 0.076 | 0.042 | 1.2X | 8.9×106±1.1×104 |
MCCL+HOS | 0.032 | 0.013 | 21.1X | 5.0×105±1.8×103 |
Solver | Time-steps | Solver | Time-steps |
CCL+FE | (0.001~1 ms) | MCCL+FE | (0.001~1 ms) |
CCL+TAL | (0.001~1 ms) | MCCL+TAL | (0.001~1 ms) |
CCL+UNI | (0.001~1 ms) | MCCL+UNI | (0.001~1 ms) |
CCL+HOS | (0.001~1 ms) | MCCL+HOS | (0.001~1 ms) |
Solvers | MRMSE-V | MRMSE-I | Speedup | Cost(ms) |
DT+RK4 | 0.0 | 0.0 | 1.0X | 1.0×107±9.8×104 |
CCL+FE | 0.016 | 0.588 | 109.9X | 9.5×104±4.0×102 |
CCL+TAL | 0.745 | 1.158 | 33.0X | 3.2×105±1.2×103 |
CCL+UNI | 0.056 | 0.012 | 16.4X | 6.4×105±8.4×102 |
CCL+HOS | 0.236 | 0.064 | 412.3X | 2.6×104±3.2×102 |
MCCL+FE | 0.022 | 0.065 | 96.9X | 1.1×105±2.7×102 |
MCCL+TAL | 0.015 | 0.124 | 29.8X | 3.5×105±1.4×103 |
MCCL+UNI | 0.021 | 0.009 | 14.8X | 7.0×105±3.4×103 |
MCCL+HOS | 0.006 | 0.002 | 17.2X | 6.1×105±3.7×103 |
Note: MRMSE-V: MRMSE of voltage in comparison to the reference. MRMSE-I: MRMSE of sodium current in comparison to the reference. Speedup: The computation cost ratio of one of eight schemes to the reference solution DT+RK4. X, times speedup. Ten simulations for each solver. Cost in Mean ± Standard deviation. |
Solvers | MRMSE-V | MRMSE-I | Speedup | Cost(ms) |
DT+RK4 | 0.0 | 0.0 | 1.0X | 1.1×107±6.0×104 |
CCL+FE | - - | - - | - - | failure |
CCL+TAL | - - | - - | - - | failure |
CCL+UNI | 0.019 | 0.008 | 1.2X | 8.9×106±1.8×104 |
CCL+HOS | 0.250 | 0.114 | 266.6X | 3.3×104±1.3×102 |
MCCL+FE | 0.103 | 3.345 | 19.8X | 5.3×105±2.5×103 |
MCCL+TAL | - - | - - | - - | failure |
MCCL+UNI | 0.076 | 0.042 | 1.2X | 8.9×106±1.1×104 |
MCCL+HOS | 0.032 | 0.013 | 21.1X | 5.0×105±1.8×103 |