Comparison of two experimental methods for the determination of Michaelis-Menten kinetics of an immobilized enzyme

For the application of immobilized enzymes, the influence of immobilization on the activity of the enzyme should be Known. This influence can be obtained by determining the intrinsic kinetic parameters of the immobilized enzyme, and by comparing them with the kinetic parameters of the suspended enzyme. This article deals with the determination of the intrinsic kinetic parameters of an agarose-gel bead immobilized oxygen-consuming enzyme: L-lactate 2-monooxygenase. The reaction rate of the enzyme can be described by Michaelis-Menten kinetics. Batch conversion experiments using a biological oxygen monitor, as well as steady-state profile measurements within the biocatalyst particles using an oxygen microsensor, were performed. Two different mathematical methods were used for the batch conversion experiments, both assuming a pseudosteady-state situation with respect to the shape of the profile inside the bead. One of the methods used an approximate relation for the effectiveness factor for Michaelis-Menten kinetics which interpolates between the analytical solutions for zero- and first-order kinetics. The other mathematical method was based on a numerical solution and combined a mass balance over the reactor with a mass balance over the bead. The main difference in the application of the two methods is the computer calculation time; the completely numerical calculation procedure was about 20 times slower than the other calculation procedure.The intrinsic kinetic parameters resulting from both experimental methods were compared to check the reliability of the methods. There was no significant difference in the intrinsic kinetic parameters obtained from the two experimental methods. By comparison of the kinetic parameters for the suspended enzyme with the intrinsic kinetic parameters for the immobilized enzyme, it appeared that immobilization caused a decrease in the value of V(m) by a factor of 2, but there was no significant difference in the values obtained for K(m).     

Theoretical evaluation of the parameters of glucose metabolism on the basis of continuous glycemia monitoring data using mathematical modeling

The aim of this paper is to construct a mathematical model that takes the main physiological parameters of blood-glucose regulation into account, in order to identify these parameters for an individual patient according to continuous glucose-monitoring data. The constructed mathematical model consists of six ordinary differential equations that describe the dynamics of changes in glucose concentrations, as well as insulin and anti-insulin factors in the blood. Estimation of the parameters of the equations was performed using an evolutionary programming method. The model predictions were fitted to the continuous glucosemonitoring data. As a result of the identification of the model parameters for two patients with type 1 diabetes mellitus, the estimated insulin secretion was close to zero and the estimated glucose utilization and insulin clearance were increased in comparison with the data for healthy donors. Here, we present a personalized model of the regulation of blood glucose, which can be used to predict the results of continuous glucose monitoring depending on modification of the prescribed glucose-lowering therapy. This approach can significantly reduce the number of iterations of the selection of medical hypoglycemic therapy and therefore increase the effectiveness of treatment according to glucose-monitoring data.

     

Measuring morphological parameters of the pelvic floor for finite element modelling purposes

The goal of this study was to obtain a complete data set needed for studying the complex biomechanical behaviour of the pelvic floor muscles using a computer model based on the finite element (FE) theory. The model should be able to predict the effect of surgical interventions and give insight into the function of pelvic floor muscles. Because there was a lack of any information concerning morphological parameters of the pelvic floor muscle structures, we performed an experimental measurement to uncover those morphological parameters. Geometric parameters as well as muscle parameters of the pelvic floor muscles were measured on an embalmed female cadaver. A three-dimensional (3D) geometric data set of the pelvic floor including muscle fibre directions was obtained using a palpator device. A 3D surface model based on the experimental data, needed for mathematical modelling of the pelvic floor, was created. For all parts of the diaphragma pelvis, the optimal muscle fibre length was determined by laser diffraction measurements of the sarcomere length. In addition, other muscle parameters such as physiological cross-sectional area and total muscle fibre length were determined. Apart from these measurements we obtained a data set of the pelvic floor structures based on nuclear magnetic resonance imaging (MRI) on the same cadaver specimen. The purpose of this experiment was to discover the relationship between the MRI morphology and geometrical parameters obtained from the previous measurements. The produced data set is not only important for biomechanical modelling of the pelvic floor muscles, but it also describes the geometry of muscle fibres and is useful for functional analysis of the pelvic floor in general. By the use of many reference landmarks all these morphologic data concerning fibre directions and optimal fibre length can be morphed to the geometrical data based on segmentation from MRI scans.These data can be directly used as an input for building a mathematical model based on FE theory.     

A stochastic model for the NMR analysis of the heteroassociation of biologically active aromatic molecules

A stochastic model for the NMR analysis of the heteroassociation of two aromatic compounds was developed, which takes into account all physically possible reactions of association of molecules in solution. Expressions for calculating the experimentally observed proton chemical shift were obtained in the general form, and an algorithm for calculating the parameters of heteroassociation using the stochastic model was proposed. The effects of limitations of the basic and general models, as compared with the stochastic model, on the model parameters of the heteroassociation of various biologically active aromatic molecules was analyzed. It was shown that the basic model can be used with a sufficient degree of accuracy for systems with a relatively small contribution of heteroassociation reactions to the total dynamic equilibrium in solution, whereas the general model describes satisfactorily the parameters of heteroassociation practically for all systems studied.     

Less invasive blood sampling in the animal laboratory: clinical chemistry and haematology of blood obtained by the Triatominae bug Dipetalogaster maximus

Dipetalogaster maximus (Dipmax), a blood-sucking bug belonging to the family Reduviidae, has been used to obtain blood samples, for example for clinical chemistry and haematology, in a variety of zoo animals and wildlife. Using this bug allows stress-free blood sampling as the bug is able to draw blood without the mammal noticing the bug. In laboratory animal science, the need for blood samples from unstressed animals may arise, especially in animal behaviour research. The use of Dipmax bugs may prove a valuable tool for this purpose. To validate the method, we compared an array of standard blood parameters sampled from New Zealand White rabbits, sampled either by the use of bugs or by the conventional method; puncture of vena auricularis caudalis. The overall hypothesis was that there was no significant difference in clinical chemistry and haematological parameters between the bug method and the conventional method. A total of 17 clinical parameters as well as 12 haematological parameters were measured and compared in New Zealand White rabbits. The results showed that for 13 of these 29 analysed parameters, the bug method and the conventional method did not give significantly different results, and the obtained results were thus directly comparable. For the remaining parameters the obtained results were significantly different. However, all parameters were measurable in the bug samples. The influences of the bug metabolism on these parameters are discussed.     

Parametric study of the swimming performance of a fish robot propelled by a flexible caudal fin

In this paper, we aim to study the swimming performance of fish robots by using a statistical approach. A fish robot employing a carangiform swimming mode had been used as an experimental platform for the performance study. The experiments conducted aim to investigate the effect of various design parameters on the thrust capability of the fish robot with a flexible caudal fin. The controllable parameters associated with the fin include frequency, amplitude of oscillation, aspect ratio and the rigidity of the caudal fin. The significance of these parameters was determined in the first set of experiments by using a statistical approach. A more detailed parametric experimental study was then conducted with only those significant parameters. As a result, the parametric study could be completed with a reduced number of experiments and time spent. With the obtained experimental result, we were able to understand the relationship between various parameters and a possible adjustment of parameters to obtain a higher thrust. The proposed statistical method for experimentation provides an objective and thorough analysis of the effects of individual or combinations of parameters on the swimming performance. Such an efficient experimental design helps to optimize the process and determine factors that influence variability.     

Influence of arterial wall-stenosis compliance on the coronary diagnostic parameters

Functional diagnostic parameters such as Fractional Flow Reserve (FFR), which is calculated from pressure measurements across stenosed arteries, are often used to determine the functional severity of coronary artery stenosis. This study evaluated the effect of arterial wall-stenosis compliance, with limiting scenarios of stenosis severity, on the diagnostic parameters. The diagnostic parameters considered in this study include an established index, FFR and two recently developed parameters: Pressure Drop Coefficient (CDP) and Lesion Flow Coefficient (LFC). The parameters were assessed for rigid artery (RR; signifying high plaque elasticity), compliant artery with calcified plaque (CC; intermediate plaque elasticity) and compliant artery with smooth muscle cell proliferation (CS; low plaque elasticity), with varying degrees of epicardial stenosis. A hyperelastic Mooney-Rivlin model was used to model the arterial wall and plaque materials. Blood was modeled as a shear thinning, non-Newtonian fluid using the Carreau model. The arterial wall compliance was evaluated using the finite element method. The present study found that, with an increase in stenosis severity, FFR decreased whereas CDP and LFC increased. The cutoff value of 0.75 for FFR was observed at 78.7% area stenosis for RR, whereas for CC and CS the cutoff values were obtained at higher stenosis severities of 81.3% and 82.7%, respectively. For a fixed stenosis, CDP value decreased and LFC value increased with a decrease in plaque elasticity (RR to CS). We conclude that the differences in diagnostic parameters with compliance at intermediate stenosis (78.7-82.7% area blockage) could lead to misinterpretation of the stenosis severity.     

On the use of one-step perturbation to investigate the dependence of NOE-derived atom–atom distance bound violations of peptides upon a variation of force-field parameters

The method of one-step perturbation can be used to predict from a single molecular dynamics simulation the values of observable quantities as functions of variations in the parameters of the Hamiltonian or biomolecular force field used in the simulation. The method is used to predict violations of nuclear overhauser effect (NOE) distance bounds measured in nuclear magnetic resonance (NMR) experiments by atom–atom distances of the NOE atom pairs when varying force-field parameters. Predictions of NOE distance bound violations between different versions of the GROMOS force field for a hexa-β-peptide in solution show that the technique works for rather large force-field parameter changes as well as for very different NOE bound violation patterns. The effect of changing individual force-field parameters on the NOE distance bound violations of the β-peptide and an α-peptide was investigated too. One-step perturbation, which in this case is equivalent to reweighting configurations, constitutes an efficient technique to predict many values of different quantities from a single conformational ensemble for a particular system, which makes it a powerful force-field development technique that easily reduces the number of required separate simulations by an order of magnitude.     

A simulation study of the effects of architectural constraints and resource translocation on population structure and competition in clonal plants

(1) Spatially explicit simulation of clonal plant growth is used to determine how ramet-level traits affect ramet density, spatial pattern of ramets and competitive ability of a clonal plant. The simulation model used combines elements of (i) an individual-based model of plant interactions, (ii) an architectural model of clonal plant growth, and (iii) a model of resource translocation within a set of physiologically integrated plant individuals. (2) The effects of two groups of parameters were studied: growth and resource acquisition parameters (resource accumulation, density-dependence of resource accumulation, resource translocation between ramets) and architectural rules (branching angle and probability of branching, internode length). The model was parameterised by values approximating those of clonally growing grasses as closely as possible. The basic parameter values were chosen from a short-turf grassland. Sensitivity analysis was carried out on relevant parameters around three basic points in the parameter space. Both single-species and two-species systems were studied. (3) It is shown that increasing resource acquisition and growth parameters increase ramet density, genet number and competitive ability. Translocation parameters and architectural parameters modify the effects of resource acquisition and growth, but their effect in single-species stands was smaller. (4) The simulations of species with fixed ramet sizes showed that ramet density in single-species stands cannot be used for predicting competitive ability. Increase in resource acquisition and growth parameters was correlated with an increase in equilibrium ramet density and competitive ability. Increasing branching angle, branching probability or internode length lead to an increased competitive ability, but did not affect equilibrium ramet density. Change of architectural parameters could therefore affect competitive ability independently of their effect on the final ramet density. (5) Spatial pattern both in single-species and two-species stands was also highly parameter-dependent. Changes in architectural parameters and in translocation usually lead to pronounced change in the spatial pattern; change in growth and resource acquisition parameters generally had little effect on spatial pattern.     

High hydrostatic pressure treatment for the inactivation of Staphylococcus aureus in human blood plasma

For the past 30years, pressure inactivation of microorganisms has been developed in biosciences, in particular for foods and more recently for biological products, including pharmaceutical ones. In many past studies, the effect of high hydrostatic pressure (HHP) processes on pathogens focused mainly on the effect of an increase of the pressure value. To assure the safety of pharmaceutical products containing fragile therapeutic components, development of new decontamination processes at the lowest pressure value is needed to maintain their therapeutic properties. The aim of this study was therefore to evaluate the impact of the process parameters characterizing high-pressure treatments [such as the pressurization rate (PR) and the application mode (AM)] on the inactivation of pathogens, in particular to determine how these parameters values could help decrease the pressure value necessary to reach the same inactivation level. The effect of these physical parameters was evaluated on the inactivation of Staphylococcus aureus ATCC 6538 which is an opportunistic pathogen of important relevance in the medical, pharmaceutical and food domains. Human blood plasma was chosen as the suspension medium because of its physiological importance in the transfusion field. It was shown that the optimization of all the selected parameters could lead to a high inactivation level (≈5log(10) decrease of the initial bacterial load) at a pressure level as low as 200MPa, underlining some synergistic effects among these parameters. Complete inactivation of the initial bacterial population was achieved for the following conditions: PR=50MPas(-1), AM=5×2min, T≈-5°C and P=300MPa.     

Development of functional variability during the motor learning process of a complex cyclic movement

Traditionally, movement variability is considered an indicator for sensorimotor malfunctioning. However, functional movement variability is also a result of compensation mechanisms e.g. to account for prior movement deviations and is, therefore, crucial for stable movements. The aim of this study was to analyze functional variability during motor learning of a complex cyclic task.Thirteen young participants practised riding a Pedalo® slalom until they were able to complete the task without errors. Since trunk movements are controlled with high priority, we analyzed trunk kinematics as a result parameter. Since lower extremities affect the result parameter, foot, thigh and pelvis kinematics are considered execution parameters. The movement variability for result and execution parameters was determined for the first (poor performance), an intermediate (medium performance) and the last (good performance) training sessions. Furthermore, the variability ratio (execution/result parameter) was calculated as a measure of functional variability.Movement variability of the result parameter decreased significantly with increasing expertise. In contrast, movement variability of all execution parameters increased significantly from measurements representing poor to medium performance. No change from medium to good performance was found. Functional variability increased over time in all execution parameters.Since the movement variability of all execution parameters did not decrease with increasing Pedalo performance, applying a traditional interpretation approach of movement variability would have led to completely wrong conclusions. Possible mechanisms explaining the increased movement variability are discussed. The variability ratio seems to be the only parameter that can reveal improved sensorimotor functioning during all analyzed stages of motor learning.     

Systematic optimization of L‐tryptophan riboswitches for efficient monitoring of the metabolite in Escherichia coli

Riboswitches form a class of genetically encoded sensor‐regulators and are considered as promising tools for monitoring various metabolites. Functional parameters of a riboswitch, like dynamic or operational range, should be optimized before the riboswitch is implemented in a specific application for monitoring the target molecule efficiently. However, optimization of a riboswitch was not straightforward and required detailed studies owing to its complex sequence‐function relationship. Here, we present three approaches for tuning and optimization of functional parameters of a riboswitch using an artificial L‐tryptophan riboswitch as an example. First, the constitutive expression level was adjusted to control the dynamic range of an L‐tryptophan riboswitch. The dynamic range increased as the constitutive expression level increased. Then, the function of a riboswitch‐encoded protein was utilized to connect the regulatory response of the riboswitch to another outcome for amplifying the dynamic range. Riboswitch‐mediated control of the host cell growth enabled the amplification of the riboswitch response. Finally, L‐tryptophan aptamers with different dissociation constants were employed to alter the operational range of the riboswitch. The dose‐response curve was shifted towards higher L‐tryptophan concentrations when an aptamer with higher dissociation constant was employed. All strategies were effective in modifying the distinct functional parameters of the L‐tryptophan riboswitch, and they could be easily applied to optimization of other riboswitches owing to their simplicity.     

Application of QbD Principles for the Evaluation of Empty Hard Capsules as an Input Parameter in Formulation Development and Manufacturing

Understanding the product and process variable on the final product performance is an essential part of the quality-by-design (QbD) principles in pharmaceutical development. The hard capsule is an established pharmaceutical dosage form used worldwide in development and manufacturing. The empty hard capsules are supplied as an excipient that is filled by pharmaceutical manufacturers with a variety of different formulations and products. To understand the potential variations of the empty hard capsules as an input parameter and its potential impact on the finished product quality, a study was performed investigating the critical quality parameters within and in between different batches of empty hard gelatin capsules. The variability of the hard capsules showed high consistency within the specification of the critical quality parameters. This also accounts for the disintegration times, when automatic endpoint detection was used. Based on these data, hard capsules can be considered as a suitable excipient for product development using QbD principles.     

Use of an oxygen microsensor for the determination of intrinsic kinetic parameters of an immobilized oxygen reducing enzyme

An oxygen microsensor was used to measure internal oxygen profiles in biocatalyst particles of different diameter and activity. The particles were made of agarose gel and contained an oxygen reducing enzyme, L-lactate mono-oxygenase. The kinetics of the enzyme could be well described by the Michaelis-Menten equation. From the internal substrate concentration profile the intrinsic kinetic parameters were determined by means of fitting a simulated profile to the measurements, using Marquardt's algorithm. The intrinsic kinetic parameters found following this procedure appeared to be independent of particle radius or enzyme loading used, proving the method to be reliable. These parameters were also compared with the kinetic parameters of the free enzyme which were determined in a biological oxygen monitoring system. The intrinsic kinetic parameters showed a decrease with a factor 2.3 for V(m) value and with a factor 2.7 for the K(m) value compared to the parameters for the free enzyme. From this the conclusion can be drawn that the immobilization as such or the carrier material not only can have an effect on the maximum intrinsic conversion rate (V(m)) but also on the affinity of the enzyme (K(m)) for oxygen.     

Influence of pH and ionic strength on the steric mass-action model parameters around the isoelectric point of protein

The ion-exchange equilibrium and the dependence of the parameters in the steric mass-action (SMA) model on salt concentration and buffer pH around the isoelectric point of protein were studied. Bovine serum albumin (BSA, isoelectric point = 5.4) was used as a model protein and DEAE Sepharose FF as an ion exchanger. Finite batch adsorption experiments and isocratic elution chromatography were performed for the determination of the model parameters (i.e., characteristic charge, equilibrium constant, and steric factor). The results showed that pH had significant effects on the parameters. With an increase of pH from 4.5 to 6.5, the characteristic charge increased from 0.9 to 3.0 and leveled off as a plateau at pH above 5.5. The charge groups in the contact region of protein surface were considered to play a crucial role on the characteristic charge. The decrease of pH and increase of salt concentration lowered the absolute value of the zeta potential of the protein surface and led to a decrease of the equilibrium constant. The steric factor remained unchanged at about 31 at pH 5.5 and 6.0 and increased to 44.5 at pH 5.0 and 96.8 at pH 4.5, mainly as a result of the lower adsorption capacity of BSA at pH <5.5. Furthermore, the increase of the molecular volume of BSA at pH 4.5 would be an additional reason for the increase of the steric factor. Taking into account the effect of the pH and salt concentration on these parameters, the SMA model described the ion exchange equilibrium of protein more accurately.     

Influence of non-bonded parameters on the quality of NMR structures: A new force field for NMR structure calculation

The effects of different non-bonded parameters of force fields for NMR structure calculation on the quality of the resulting NMR solution structures were investigated using Interleukin 4 as a model system. NMR structure ensembles were calculated with an ab initio protocol using torsion angle dynamics. The calculations were repeated with five different non-bonded energy functions and parameters. The resulting ensembles were compared with the available X-ray structures, and their quality was assessed with common structure validation programs. In addition, the impact of torsion angle restraints and dihedral energy terms for the sidechains and the backbone was studied. The further improvement of the quality by refinement in explicit solvent was demonstrated. The optimal parameters, including those necessary for water refinement, are available in the new version of the PARALLHDG force field.     

The size of the chi-square test for the Hardy-Weinberg law

Many scientific problems can be formulated in terms of a statistical model indexed by parameters, only some of which are of scientific interest and the other parameters, called nuisance parameters, are not of interest in themselves. For testing the Hardy-Weinberg law, a relation among genotype and allele probabilities is of interest and allele probabilities are of no interest and now nuisance parameters. In this paper we investigate how the size (the maximum of the type I error rate over the nuisance parameter space) of the chi-square test for the Hardy-Weinberg law is affected by the nuisance parameters. Whether the size is well controlled or not under the nominal level has been frequently investigated as basic components of statistical tests. The size represents the type I error rate at the worst case. We prove that the size is always greater than the nominal level as the sample size increases. Extensive computations show that the size of the chi-squared test (worst type I error rate over the nuisance parameter space) deviates more upwardly from the nominal level as the sample size gets larger. The value at which the maximum of the type I error rate was found moves closer to the edges of the the nuisance parameter space with increasing sample size. An exact test is recommended as an alternative when the type I error is inflated.     

Urbanization and the Microbial Content of the North Saskatchewan River

The effect of urbanization on the microbial content of the North Saskatchewan River was determined by following the changes in the numbers of total bacteria, total eosin methylene blue (EMB) plate count, and Escherichia coli as the river flowed from its glacial source, through parklands, and out into the prairies. Changes in physical parameters such as pH, temperature, salt concentration, and the amount and nature of the suspended material were also determined to evaluate their on the microbial parameters being measured. The level of all three microbial parameters studied slowly increased as the river flowed from its glacial source out into the prairies. The major effect of small hamlets, with or without sewage treatment facilities, appears to be to supply nutrients which supports the growth of the indigenous river flora but not E. coli. In contrast, the effect of a large urban center, with a population of approximately 500,000, which utilizes primary and secondary sewage processes in disposing of sewage, is to provide the nutrients and an inoculum of E. coli which results in a marked increase in the numbers of all three microbial groups studied. The effect of this urban center was still discernible 300 miles downstream. The river was also monitored for the presence of Salmonella sp. Only one positive isolation was achieved during this study, and this isolate was characterized as being Salmonella alachua.     

Exploration of the intercellular heterogeneity of topotecan uptake into human breast cancer cells through compartmental modelling

A mathematical multi-cell model for the in vitro kinetics of the anti-cancer agent topotecan (TPT) following administration into a culture medium containing a population of human breast cancer cells (MCF-7 cell line) is described. This non-linear compartmental model is an extension of an earlier single-cell type model and has been validated using experimental data obtained using two-photon laser scanning microscopy (TPLSM). A structural identifiability analysis is performed prior to parameter estimation to test whether the unknown parameters within the model are uniquely determined by the model outputs. The full model has 43 compartments, with 107 unknown parameters, and it was found that the structural identifiability result could not be established even when using the latest version of the symbolic computation software Mathematica. However, by assuming that a priori knowledge is available for certain parameters, it was possible to reduce the number of parameters to 81, and it was found that this (Stage Two) model was globally (uniquely) structurally identifiable. The identifiability analysis demonstrated how valuable symbolic computation is in this context, as the analysis is far too lengthy and difficult to be performed by hand.     

Sources of variation for nutritional condition indices of the plasma of migratory lesser kestrels in the breeding grounds

Although published information on reference values for biochemical parameters in birds of prey has increased during the last years, little is known on their sources of variation. We used an insectivorous and small migratory raptor species, the lesser kestrel Falco naumanni, as a model. We looked for sources of variation of nutritional biochemical parameters (i.e. triglycerides, cholesterol, uric acid and urea) of both nestlings and adults. Reference values indicated that, as a rule, lesser kestrel showed more elevated triglycerides, urea and uric acid levels than other raptors. All analyzed factors except gender (i.e. year, colony, sampling time, presence/absence of a geolocator, body mass, laying date and capture date) reached significance for at least one biochemical parameter. In the morning, we found an important postprandial increase in the concentration of all biochemical parameters in nestlings, and uric acid and urea levels in adults. A positive relationship was detected between triglycerides and body mass in nestlings and adults. Although we did not find differences between blood biochemical parameters of the oldest and youngest chick of each brood, we found that cholesterol levels were lower in nestlings from larger broods. Coloration of tarsi (measured as brightness) was related to triglycerides and urea levels of nestlings and adults, respectively. The feeding habits of lesser kestrel probably explain the different levels and patterns of variation of metabolites in comparison to more carnivorous raptors eating mammals or birds.