Steady-state analytical approximations

To be analytically useful equation 13.16 needs to be written in terms of the concentrations of enzyme and substrate. This is accomplished by applying the steady-state approximation, in which we assume that the concentration of ES is essentially constant. After an initial period in which the enzyme-substrate complex first forms, the rate of formation of ES... 

A large number of analytical solutions of these equations appear in the literature. Mostly, however, they deal only with first order reactions. All others require solution by numerical or other approximate means. In this book, solutions of two examples are carried along analytically part way in P7.02.06 and P7.02.07. Section 7.4 considers flow through an external film, while Section 7.5 deals with diffusion and reaction in catalyst pores under steady state conditions. 

As mentioned in Section 3.1, an analytical solution can be provided for the steady-state of fully labile complexes, without needing to resort to the excess of ligand approximation ... 

Traditionally, reaction mechanisms of the kind above have been analysed based on the steady-state approximation. The differential equations for this mechanism cannot be integrated analytically. Numerical integration was not readily available and thus approximations were the only options available to the researcher. The concentrations of the catalyst and of the intermediate, activated complex B are always only very low and even more so their derivatives [Cat] and [B]. In the steady-state approach these two derivatives are set to 0. 

As flow rates decrease, the perfusion medium in the probe approaches equilibrium with the ECF (Wages et al., 1986). Therefore, the dialysate concentration of an analyte sampled at very lowflow rates more closely approximates the concentration in the extracellular environment (Menacherry et al., 1992). Like no net flux and the zero flow models, this is another steady-state analysis with limited application to transient changes based on behavior or pharmacological manipulations. However, the advent of new techniques in analytical chemistry requiring only small sample volumes from short sampling intervals may signal a potential return to the low flow method. 

Therefore, we need to find approximate methods for simultaneous reaction systems that will permit finding analytical solutions for reactants and products in simple and usable form. We use two approximations that were developed by chemists to simplify simultaneous reaction systems (1) the equilibrium step approximation and (2) the pseudo-steady-state approximation... 

Existence and uniqueness of the particular solution of (5.1) for an initial value y° can be shown under very mild assumptions. For example, it is sufficient to assume that the function f is differentiable and its derivatives are bounded. Except for a few simple equations, however, the general solution cannot be obtained by analytical methods and we must seek numerical alternatives. Starting with the known point (tD,y°), all numerical methods generate a sequence (tj y1), (t2,y2),. .., (t. y1), approximating the points of the particular solution through (tQ,y°). The choice of the method is large and we shall be content to outline a few popular types. One of them will deal with stiff differential equations that are very difficult to solve by classical methods. Related topics we discuss are sensitivity analysis and quasi steady state approximation. 

It remains now to solve Eq. (2.3). Here, there are various approaches, depending on the conditions. When a non-steady-state solution is required, one can introduce the decoupling approximation of Sumi and Marcus, if there is the difference in time scales mentioned earlier. Or one can integrate Eq. (2.3) numerically. For the steady-state approximation either Eq. (2.3) can again be solved numerically or some additional analytical approximation can be introduced. For example, one introduced elsewhere [44] is to consider the case that most of the reacting systems cross the transition state in some narrow window (X, X i jA), narrow compared with the X region of the reactant [e.g., the interval (O,Xc) in Fig. 2]. In that case the k(X) can be replaced by a delta function, fc(Xi)A5(X-Xi). Equation (2.3) is then readily integrated and the point X is obtained as the X that maximizes the rate expression. The A is obtained from the width of the distribution of rates in that system [44]. 

The rotating ring—disc electrode (RRDE) is probably the most well-known and widely used double electrode. It was invented by Frumkin and Nekrasov [26] in 1959. The ring is concentric with the disc with an insulating gap between them. An approximate solution for the steady-state collection efficiency N0 was derived by Ivanov and Levich [27]. An exact analytical solution, making the assumption that radial diffusion can be neglected with respect to radial convection, was obtained by Albery and Bruckenstein [28, 29]. We follow a similar, but simplified, argument below. 

Wiedemeier and coworkers (1996) have suggested two methods to approximate biodegradation rates in groundwater field studies (a) use a biologically recalcitrant tracer (e.g., three isomers of trimethylbenzene) in the groundwater to correct for dilution, sorption, and/or volatilization and calculate the rate constant by using the downgradient travel time or (b) assume that the plume has evolved to a dynamic steady-state equilibrium and develop a one-dimensional analytical solution to the advection-dispersion equation. 

The solution of the above problem applies to both transient and steady-state feedback experiments. Since transient SECM measurements are somewhat less accurate and harder to perform, most quantitative studies were carried out under steady-state conditions. The non-steady-state SECM response depends on too many parameters to allow presentation of a complete set of working curves, which would cover all experimental possibilities. The steady-state theory is simpler and often can be expressed in the form of dimensionless working curves or analytical approximations. 

Liljeroth et al. [80] used SECM in the feedback mode to study the electronic conductivity of a film of gold nanoparticles deposited at various pressures on a nonconductive substrate. They were able to observe an insulator-to-metal transition associated with a change in surface pressure. Unwin Whitworth et al. [83] have also developed a method to determine the electronic conductivity of ultrathin films using SECM under steady-state conditions. They obtained analytical approximations for the fitting of approach curves. The usefulness of their approach was demonstrated by investigating the effect of surface pressure on conductivity of a polyaniline monolayer at the water-air interface. 

For most macroscopic dynamic systems, the neglect of correlations and fluctuations is a legitimate approximation [383]. For these cases the deterministic and stochastic approaches are essentially equivalent, and one is free to use whichever approach turns out to be more convenient or efficient. If an analytical solution is required, then the deterministic approach will always be much easier than the stochastic approach. For systems that are driven to conditions of instability, correlations and fluctuations will give rise to transitions between nonequihbrium steady states and the usual deterministic approach is incapable of accurately describing the time behavior. On the other hand, the stochastic simulation algorithm is directly applicable to these studies. 

The shape of the growth spirals may be quite complex and only approximated by analytical functions such as the Archimedean spiral employed by Burton et al. [24] and shown in Fig. 12. The relationship between the velocity of a straight step and curved step may be derived using the approach of Nielsen [27]. This enables the steady-state shape of a spiral in contact with a fixed supersaturation solution to be calculated. 

Often the key entity one is interested in obtaining in modeling enzyme kinetics is the analytical expression for the turnover flux in quasi-steady state. Equations (4.12) and (4.38) are examples. These expressions are sometimes called Michaelis-Menten rate laws. Such expressions can be used in simulation of cellular biochemical systems, as is the subject of Chapters 5, 6, and 7 of this book. However, one must keep in mind that, as we have seen, these rates represent approximations that result from simplifications of the kinetic mechanisms. We typically use the approximate Michaelis-Menten-type flux expressions rather than the full system of equations in simulations for several reasons. First, often the quasi-steady rate constants (such as Ks and K in Equation (4.38)) are available from experimental data while the mass-action rate constants (k+i, k-i, etc.) are not. In fact, it is possible for different enzymes with different detailed mechanisms to yield the same Michaelis-Menten rate expression, as we shall see below. Second, in metabolic reaction networks (for example), reactions operate near steady state in vivo. Kinetic transitions from one in vivo steady state to another may not involve the sort of extreme shifts in enzyme binding that have been illustrated in Figure 4.7. Therefore the quasi-steady approximation (or equivalently the approximation of rapid enzyme turnover) tends to be reasonable for the simulation of in vivo systems. 

The same considerations made before are valid for this case and it is very important to have an available validated reaction mechanism. It can be obtained from three main sources (Blelski et al., 1985 Buxton et al., 1988 Stefan and Bolton, 1998) and it is shown in Table 5. With the available information about the constant k2, k, k, fcg, and k-j, it could be possible to solve a system of four differential equations and extract from the experimental data, the missing constants 4> and k (that in real terms is k /Co2)-This method would provide good information about the kinetic constants, but it is not the best result for studying temperature effects if the same information is not available for the pre-exponential factors and the activation energies. Then, it is better to look for an analytical expression even if it is necessary to make some approximations. This is particularly true in this case, where the direct application of the micro steady-state approximation (MSSA) is more difficult due to the existence of a recombination step that includes the two free radicals formed in the reaction. From the available information, it is possible to know that to calculate the pseudo-steady-state... 

---