Extending the kinetic solution of the classic Michaelis-Menten model of enzyme action

The principal aim of studies of enzyme-mediated reactions has been to provide comparative and quantitative information on enzyme-catalyzed reactions under distinct conditions. The classic Michaelis-Menten model (Biochem Zeit 49:333, 1913) for enzyme kinetic has been widely used to determine important parameters involved in enzyme catalysis, particularly the Michaelis-Menten constant (K (M) ) and the maximum velocity of reaction (V (max) ). Subsequently, a detailed treatment of the mechanisms of enzyme catalysis was undertaken by Briggs-Haldane (Biochem J 19:338, 1925). These authors proposed the steady-state treatment, since its applicability was constrained to this condition. The present work describes an extending solution of the Michaelis-Menten model without the need for such a steady-state restriction. We provide the first analysis of all of the individual reaction constants calculated analytically. Using this approach, it is possible to accurately predict the results under new experimental conditions and to characterize and optimize industrial processes in the fields of chemical and food engineering, pharmaceuticals and biotechnology.

Numerical solution of the two-phase expansion of a metastable flashing liquid jet using the dispersion-controlled dissipative scheme

This paper presents a study of the stationary phenomenon of superheated or metastable liquid jets, flashing into a two-dimensional axisymmetric domain, while in the two-phase region. In general, the phenomenon starts off when a high-pressure, high-temperature liquid jet emerges from a small nozzle or orifice expanding into a low-pressure chamber, below its saturation pressure taken at the injection temperature. As the process evolves, crossing the saturation curve, one observes that the fluid remains in the liquid phase reaching a superheated condition. Then, the liquid undergoes an abrupt phase change by means of an oblique evaporation wave. Across this phase change the superheated liquid becomes a two-phase high-speed mixture in various directions, expanding to supersonic velocities. In order to reach the downstream pressure, the supersonic fluid continues to expand, crossing a complex bow shock wave. The balance equations that govern the phenomenon are mass conservation, momentum conservation, and energy conservation, plus an equation-of-state for the substance. A false-transient model is implemented using the shock capturing scheme: dispersion-controlled dissipative (DCD), which was used to calculate the flow conditions as the steady-state condition is reached. Numerical results with computational code DCD-2D vI have been analyzed. Copyright (C) 2009 John Wiley & Sons, Ltd.

Numerical solution of the eXtended Pom-Pom model for viscoelastic free surface flows

In this paper we present a finite difference method for solving two-dimensional viscoelastic unsteady free surface flows governed by the single equation version of the eXtended Pom-Pom (XPP) model. The momentum equations are solved by a projection method which uncouples the velocity and pressure fields. We are interested in low Reynolds number flows and, to enhance the stability of the numerical method, an implicit technique for computing the pressure condition on the free surface is employed. This strategy is invoked to solve the governing equations within a Marker-and-Cell type approach while simultaneously calculating the correct normal stress condition on the free surface. The numerical code is validated by performing mesh refinement on a two-dimensional channel flow. Numerical results include an investigation of the influence of the parameters of the XPP equation on the extrudate swelling ratio and the simulation of the Barus effect for XPP fluids. (C) 2010 Elsevier B.V. All rights reserved.

Numerical solution of the PTT constitutive equation for unsteady three-dimensional free surface flows

This work deals with the development of a numerical technique for simulating three-dimensional viscoelastic free surface flows using the PTT (Phan-Thien-Tanner) nonlinear constitutive equation. In particular, we are interested in flows possessing moving free surfaces. The equations describing the numerical technique are solved by the finite difference method on a staggered grid. The fluid is modelled by a Marker-and-Cell type method and an accurate representation of the fluid surface is employed. The full free surface stress conditions are considered. The PTT equation is solved by a high order method, which requires the calculation of the extra-stress tensor on the mesh contours. To validate the numerical technique developed in this work flow predictions for fully developed pipe flow are compared with an analytic solution from the literature. Then, results of complex free surface flows using the FIT equation such as the transient extrudate swell problem and a jet flowing onto a rigid plate are presented. An investigation of the effects of the parameters epsilon and xi on the extrudate swell and jet buckling problems is reported. (C) 2010 Elsevier B.V. All rights reserved.

Numerical Solution of the Upper-Convected Maxwell Model for Three-Dimensional Free Surface Flows

This work presents a finite difference technique for simulating three-dimensional free surface flows governed by the Upper-Convected Maxwell (UCM) constitutive equation. A Marker-and-Cell approach is employed to represent the fluid free surface and formulations for calculating the non-Newtonian stress tensor on solid boundaries are developed. The complete free surface stress conditions are employed. The momentum equation is solved by an implicit technique while the UCM constitutive equation is integrated by the explicit Euler method. The resulting equations are solved by the finite difference method on a 3D-staggered grid. By using an exact solution for fully developed flow inside a pipe, validation and convergence results are provided. Numerical results include the simulation of the transient extrudate swell and the comparison between jet buckling of UCM and Newtonian fluids.

Approximate solution of the autocatalytic hydrolysis of cellulose

Depolymerization of cellulose in homogeneous acidic medium is analyzed on the basis of autocatalytic model of hydrolysis with a positive feedback of acid production from the degraded biopolymer. The normalized number of scissions per cellulose chain, S(t)/nA degrees A = 1 - C(t)/C(0), follows a sigmoid behavior with reaction time t, and the cellulose concentration C(t) decreases exponentially with a linear and cubic time dependence, C(t) = C(0)exp[-at - bt (3)], where a and b are model parameters easier determined from data analysis.

Characterization of the interdependency between residues that bind the substrate in a β-glycosidase

The manner by which effects of simultaneous mutations combine to change enzymatic activity is not easily predictable because these effects are not always additive in a linear manner. Hence, the characterization of the effects of simultaneous mutations of amino acid residues that bind the substrate can make a significant contribution to the understanding of the substrate specificity of enzymes. In the β-glycosidase from Spodoptera frugiperda (Sfβgly), both residues Q39 and E451 interact with the substrate and this is essential for defining substrate specificity. Double mutants of Sfβgly (A451E39, S451E39 and S451N39) were prepared by site-directed mutagenesis, expressed in bacteria and purified using affinity chromatography. These enzymes were characterized using p-nitrophenyl β-galactoside and p-nitrophenyl β-fucoside as substrates. The k cat/Km ratio for single and double mutants of Sfβgly containing site-directed mutations at positions Q39 and E451 was used to demonstrate that the effect on the free energy of ES‡ (enzyme-transition state complex) of the double mutations (∆∆G‡xy) is not the sum of the effects resulting from the single mutations (∆∆G‡x and ∆∆G‡y). This difference in ∆∆G‡ indicates that the effects of the single mutations partially overlap. Hence, this common effect counts only once in ∆∆G‡xy. Crystallographic data on β-glycosidases reveal the presence of a bidentate hydrogen bond involving residues Q39 and E451 and the same hydroxyl group of the substrate. Therefore, both thermodynamic and crystallographic data suggest that residues Q39 and E451 exert a mutual influence on their respective interactions with the substrate.

Cyclic variability of the circumstellar disk of the Be star zeta Tauri II. Testing the 2D global disk oscillation model

Context. About 2/3 of the Be stars present the so-called V/R variations, a phenomenon characterized by the quasi-cyclic variation in the ratio between the violet and red emission peaks of the HI emission lines. These variations are generally explained by global oscillations in the circumstellar disk forming a one-armed spiral density pattern that precesses around the star with a period of a few years. Aims. This paper presents self-consistent models of polarimetric, photometric, spectrophotometric, and interferometric observations of the classical Be star zeta Tauri. The primary goal is to conduct a critical quantitative test of the global oscillation scenario. Methods. Detailed three-dimensional, NLTE radiative transfer calculations were carried out using the radiative transfer code HDUST. The most up-to-date research on Be stars was used as input for the code in order to include a physically realistic description for the central star and the circumstellar disk. The model adopts a rotationally deformed, gravity darkened central star, surrounded by a disk whose unperturbed state is given by a steady-state viscous decretion disk model. It is further assumed that this disk is in vertical hydrostatic equilibrium. Results. By adopting a viscous decretion disk model for zeta Tauri and a rigorous solution of the radiative transfer, a very good fit of the time-average properties of the disk was obtained. This provides strong theoretical evidence that the viscous decretion disk model is the mechanism responsible for disk formation. The global oscillation model successfully fitted spatially resolved VLTI/AMBER observations and the temporal V/R variations in the H alpha and Br gamma lines. This result convincingly demonstrates that the oscillation pattern in the disk is a one-armed spiral. Possible model shortcomings, as well as suggestions for future improvements, are also discussed.

A note on the likelihood and moments of the skew-normal distribution

In this paper an alternative approach to the one in Henze (1986) is proposed for deriving the odd moments of the skew-normal distribution considered in Azzalini (1985). The approach is based on a Pascal type triangle, which seems to greatly simplify moments computation. Moreover, it is shown that the likelihood equation for estimating the asymmetry parameter in such model is generated as orthogonal functions to the sample vector. As a consequence, conditions for a unique solution of the likelihood equation are established, which seem to hold in more general setting.

Performance of the Equivalent Reservoir Modelling Technique for Multi-Reservoir Hydropower Systems

This paper investigates the validity of a simplified equivalent reservoir representation of a multi-reservoir hydroelectric system for modelling its optimal operation for power maximization. This simplification, proposed by Arvanitidis and Rosing (IEEE Trans Power Appar Syst 89(2):319-325, 1970), imputes a potential energy equivalent reservoir with energy inflows and outflows. The hydroelectric system is also modelled for power maximization considering individual reservoir characteristics without simplifications. Both optimization models employed MINOS package for solution of the non-linear programming problems. A comparison between total optimized power generation over the planning horizon by the two methods shows that the equivalent reservoir is capable of producing satisfactory power estimates with less than 6% underestimation. The generation and total reservoir storage trajectories along the planning horizon obtained by equivalent reservoir method, however, presented significant discrepancies as compared to those found in the detailed modelling. This study is motivated by the fact that Brazilian generation system operations are based on the equivalent reservoir method as part of the power dispatch procedures. The potential energy equivalent reservoir is an alternative which eliminates problems with the dimensionality of state variables in a dynamic programming model.

New Analytic Solution of Boltzmann Transform for Horizontal Water Infiltration into Sand

The use of the Boltzmann transform function, lambda(theta), to solve the Richards equation when the diffusivity, D, is a function of only soil water content,., is now commonplace in the literature. Nevertheless, a new analytic solution of the Boltzmann transform lambda(h) as a function of matric potential for horizontal water infiltration into a sand was derived without invoking the concept or use of D(theta). The derivation assumes that a similarity exists between the soil water retention function and the Boltzmann transform lambda(theta). The solution successfully described soil water content profiles experimentally measured for different infiltration times into a homogeneous sand and agrees with those presented by Philip in 1955 and 1957. The applicability of this solution for all soils remains open, but it is anticipated to hold for soils whose air-filled pore-size distribution before wetting is sufficiently narrow to yield a sharp increase of water content at the wetting front during infiltration. It also improves and provides a versatile alternative to the well-known analysis pioneered by Green and Ampt in 1911.

QUADRANT-SEARCHING: A NEW TECHNIQUE TO REDUCE THE SEARCH SPACE OF THE INVERSE PROBLEM OF EIT USING BEM TO THE DIRECT PROBLEM

The image reconstruction using the EIT (Electrical Impedance Tomography) technique is a nonlinear and ill-posed inverse problem which demands a powerful direct or iterative method. A typical approach for solving the problem is to minimize an error functional using an iterative method. In this case, an initial solution close enough to the global minimum is mandatory to ensure the convergence to the correct minimum in an appropriate time interval. The aim of this paper is to present a new, simple and low cost technique (quadrant-searching) to reduce the search space and consequently to obtain an initial solution of the inverse problem of EIT. This technique calculates the error functional for four different contrast distributions placing a large prospective inclusion in the four quadrants of the domain. Comparing the four values of the error functional it is possible to get conclusions about the internal electric contrast. For this purpose, initially we performed tests to assess the accuracy of the BEM (Boundary Element Method) when applied to the direct problem of the EIT and to verify the behavior of error functional surface in the search space. Finally, numerical tests have been performed to verify the new technique.

Evaluation of the aquatic macroinvertebrate community as a tool for monitoring a reservoir in the Pitangui river basin, Parana, Brazil

Evaluation of the aquatic macroinvertebrate community as a tool for monitoring a reservoir in the Pitangui river basin, Parana, Brazil. Benthic and nektonic macroinvertebrates play an important role in the structure and function of aquatic ecosystems and their distribution is influenced by chemical features of the substrate, vegetation composition, and water depth. Knowledge on the fauna contributes to the assessment of water quality and development of biodiversity conservation activities. Different biotic factors affecting the invertebrate community were evaluated in the Alagados reservoir, the main water source of the city of Ponta Grossa, Parana. In five different sampling points, 18,473 specimens of aquatic or semiaquatic macroinvertebrates were collected, belonging to 46 taxa of the phylla Annelida (Hirudinea and Oligochaeta), Mollusca (Gastropoda), Platyhelminthes (Turbellaria), Nematoda and Arthropoda (Arachnida, Crustacea and Insecta). This community was composed mainly of predators (45.7% of the taxa sampled), collectors and/or filterers (23.9%), scrapers (15.2%), shredders (13.0%) and detritivores (2.2%). Diversity (H`) and evenness (J) indices were significantly low for the sites examined, and H` ranged between 0.3301 and 1.0396. Regarding tolerance of organisms to organic pollution, more sensitive taxa were very rare (Plecoptera) or unusual (Trichoptera and Ephemeroptera). Among the more resistant groups are Chironomidae and Hirudinea, both fairly common in the samples. This study corroborates the importance of bioindicators as a tool to assess water quality for human consumption and for the conservation of aquatic environments, integrating physical, chemical and biological factors in monitoring programs.

The role of an insulating surface layer on the relaxation time of the ionic redistribution in an electrolytic cell

We analyze the influence of a surface dielectric layer on the transient phenomena related to the ionic redistribution in an electrolytic cell submitted to a step-like external voltage. The adsorption-desorption phenomenon is taken into account in the famework of the Gouy-Chapman approximation, where the ions are assumed dimensionless. In the limit of small amplitude of the applied voltage, where the equations of the problem can be linearized, we obtain an analytical solution for the surface densities of ions, for the electrical potential and for the relaxation time for the transient phenomena. In the general case, when the linearized analysis is no longer valid, the solution of the problem is obtained numerically. The role of the thickness of the dielectric layer on the relaxation time is also discussed.

Kinetic and mechanistic investigation of the ozonolysis of 2,4-xylidine (2,4-dimethyl-aniline) in acidic aqueous solution

The ozonolysis of 2,4-xylidine (2,4-dimethyl-aniline) in acidic aqueous solution was investigated by determining the major reaction products and their evolution as a function of the reaction time and their dependence on the pH of the reaction system. 2,4-Dimethyl-nitrobenzene and 2,4-dimethyl-phenol were found to be primary reaction products; their formation might be explained by electron transfer and substitution reactions. 2,4-Dimethyl-phenol was further oxidized yielding 2,4-dimethyl- and/or 4,6-dimethyl-resorcinol by electrophilic addition of HO(center dot) radicals. The best fitting phenomenological kinetic model and the good convergence of calculated and experimentally determined rate constants imply two additional competitive pathways of substrate oxidation: (i) electrophilic addition of HO(center dot) radicals and fast subsequent substitution would also yield the resorcinol derivatives. (ii) Substrate and isolated products are thought to be oxidized by hydrogen abstraction at the benzylic sites, but the corresponding products (alcohols, aldehydes, and carboxylic acids) could not be identified. Fe(II) was added to probe for the presence of H(2)O(2), but had no or only a minor effect on the kinetics of the ozonolysis. (c) 2009 Elsevier B.V. All rights reserved.