Efficient and accurate experimental design for enzyme kinetics: Bayesian studies reveal a systematic approach

In areas such as drug development, clinical diagnosis and biotechnology research, acquiring details about the kinetic parameters of enzymes is crucial. The correct design of an experiment is critical to collecting data suitable for analysis, modelling and deriving the correct information. As classical design methods are not targeted to the more complex kinetics being frequently studied, attention is needed to estimate parameters of such models with low variance. We demonstrate that a Bayesian approach (the use of prior knowledge) can produce major gains quantifiable in terms of information, productivity and accuracy of each experiment. Developing the use of Bayesian Utility functions, we have used a systematic method to identify the optimum experimental designs for a number of kinetic model data sets. This has enabled the identification of trends between kinetic model types, sets of design rules and the key conclusion that such designs should be based on some prior knowledge of K-M and/or the kinetic model. We suggest an optimal and iterative method for selecting features of the design such as the substrate range, number of measurements and choice of intermediate points. The final design collects data suitable for accurate modelling and analysis and minimises the error in the parameters estimated. (C) 2003 Elsevier Science B.V. All rights reserved.

Gas-phase kinetics

Commentary on: 'Rates of pyrolysis and bond energies of substituted organic iodides' (Part 1), E.T. Butler and M. Polanyi, Trans Faraday Society, 1943, 39, 19.

A temperature-dependent relative-rate study of the OH initiated oxidation of n-butane: the kinetics of the reactions of the 1-and 2-butoxy radicals

The kinetics of the reactions of 1-and 2-butoxy radicals have been studied using a slow-flow photochemical reactor with GC-FID detection of reactants and products. Branching ratios between decomposition, CH3CH(O-.)CH2CH3 CH3CHO + C2H5, reaction (7), and reaction with oxygen, CH3CH(O-.)CH2CH3 + O-2 -> CH3C(O)C2H5 + HO2, reaction (6), for the 2-butoxy radical and between isomerization, CH3CH2CH2CH2O. -> CH2CH2CH2CH2OH, reaction (9), and reaction with oxygen, CH3CH2CH2CH2O. + O-2 -> C3H7CHO + HO2, reaction (8), for the 1-butoxy radical were measured as a function of oxygen concentration at atmospheric pressure over the temperature range 250-318 K. Evidence for the formation of a small fraction of chemically activated alkoxy radicals generated from the photolysis of alkyl nitrite precursors and from the exothermic reaction of 2-butyl peroxy radicals with NO was observed. The temperature dependence of the rate constant ratios for a thermalized system is given by k(7)/k(6) = 5.4 x 1026 exp[(-47.4 +/- 2.8 kJ mol(-1))/RT] molecule cm(-3) and k(9)/k(8) = 1.98 x 10(23) exp[(-22.6 +/- 3.9 kJ mol(-1))/RT] molecule cm(-3). The results agree well with the available experimental literature data at ambient temperature but the temperature dependence of the rate constant ratios is weaker than in current recommendations.

A pressure-jump time-resolved x-ray diffraction study of cubic-cubic transition kinetics in monoolein

In the past two decades, the geometric pathways involved in the transformations between inverse bicontinuous cubic phases in amphiphilic systems have been extensively theoretically modeled. However, little experimental data exists on the cubic-cubic transformation in pure lipid systems. We have used pressure-jump time-resolved X-ray diffraction to investigate the transition between the gyroid Q(II)(G) and double-diamond Q(II)(D) phases in mixtures of 1-monoolein in 30 wt% water. We find for this system that the cubic-cubic transition occurs without any detectable intermediate structures. In addition, we have determined the kinetics of the transition, in both the forward and reverse directions, as a function of pressure-jump amplitude, temperature, and water content. A recently developed model allows (at least in principle) the calculation of the activation energy for lipid phase transitions from such data. The analysis is applicable only if kinetic reproducibility is achieved, at least within one sample, and achievement of such kinetic reproducibility is shown here, by carrying out prolonged pressure-cycling. The rate of transformation shows clear and consistent trends with pressure-jump amplitude, temperature, and water content, all of which are shown to be in agreement with the effect of the shift in the position of the cubic-cubic phase boundary following a change in the thermodynamic parameters.

The kinetics of the 2π+2π photodimerisation reactions of single-crystalline derivatives of trans-cinnamic acid: a study by infrared microspectroscopy

The kinetics of the photodimerisation reactions of the 2- and 4-β-halogeno-derivatives of trans-cinnamic acid (where the halogen is fluorine, chlorine or bromine) have been investigated by infrared microspectroscopy. It is found that none of the reactions proceed to 100% yield. This is in line with a reaction mechanism developed by Wernick and his co-workers that postulates the formation of isolated monomers within the solid, which cannot react. β-4-Bromo and β-4-chloro-trans-cinnamic acids show approximately first order kinetics, although in both cases the reaction accelerates somewhat as it proceeds. First order kinetics is explained in terms of a reaction between one excited- and one ground-state monomer molecule, while the acceleration of the reaction implies that it is promoted as defects are formed within the crystal. By contrast β-2-chloro-trans-cinnamic acid shows a strongly accelerating reaction which models closely to the contracting cube equation. β-2-Fluoro- and β-4-fluoro-trans-cinnamic acids show a close match to first order kinetics. The 4-fluoro-derivative, however, shows a reaction that proceeds via a structural intermediate. The difference in behaviour between the 2-fluoro- and 4-fluoro-derivative may be due to different C–HF hydrogen bonds observed within these single-crystalline starting materials.

The kinetics of the gas-phase decomposition of the 2-methyl-2-butoxyl and 2-methyl-2-pentoxyl radicals

The kinetics of the title reactions have been studied by relative-rate methods as a function of temperature. Relative-rate coefficients for the two decomposition channels of 2-methyl-2-butoxyl have been measured at five different temperatures between 283 and 345 K and the observed temperature dependence is consistent with the results of some previous experimental studies. The kinetics of the two decomposition channels of 2-methyl-2-pentoxyl have also been investigated, as a function of temperature, relative to the estimated rate of isomerisation of this radical. Room-temperature rate coefficient data for the two decomposition channels of both 2-methyl-2-pentoxyl and 2-methyl-2-butxoyl (after combining the relative rate coefficient for this latter with a value for the rate coefficient of the major channel, extrapolated from the data presented by Batt et al., Int. J. Chem. Kinet., 1978, 10, 931) are shown to be consistent with a non-linear kinetic correlation, for alkoxyl radical decomposition rate data, previously presented by this laboratory (Johnson et al., Atmos. Environ., 2004, 38, 1755-1765).

Correlation-type structure activity relationships for the kinetics of gas-phase RO2 self-reactions and reaction with HO2

Structure activity relationships (SARs) are presented for the gas-phase reactions of RO2 with HO2, and the self- and cross-reactions of RO2. For RO2+HO2 the SAR is based upon a correlation between the logarithm of the measured rate coefficient and a calculated ionisation potential for the molecule R-CH=CH2, R being the same group in both the radical and molecular analogue. The correlation observed is strong and only for one RO2 species does the measured rate coefficient deviate by more than a factor of two from the linear least-squares regression line. For the self- and cross-reactions of RO2 radicals, the SAR is based upon a correlation between the logarithm of the measured rate coefficient and the calculated electrostatic potential (ESP) at the equivalent carbon atom in the RH molecule to which oxygen is attached in RO2, again R being the same group in the molecule and the radical. For cases where R is a simple alkyl-group, a strong linear correlation observed. For RO2 radicals which contain lone pair-bearing substituents and for which the calculated ESP<-0.05 self-reaction rate coefficients appear to be insensitive to the value of the ESP. For RO2 of this type with ESP>-0.05 a linear relationship between log k and the ESP is again observed. Using the relationships, 84 out of the 85 rate coefficients used to develop the SARs are predicted to within a factor of three of their measured values. A relationship is also presented that allows the prediction of the Arrhenius parameters for the self-reactions of simple alkyl RO2 radicals. On the basis of the correlations, predictions of room-temperature rate coefficients are made for a number of atmospherically important peroxyl-peroxyl radical reactions. (C) 2003 Elsevier Ltd. All rights reserved.

Oxovanadium(V) complexes of bis(phenolate) ligands with acetylacetone as co-ligand: synthesis, crystal structure, electrochemical and kinetics studies on the oxidation of ascorbic acid

Two vanadium(V) complexes, [VO(L-1)]acac)] (1) and [VO(L-2)(acac)] (2), where H2L1 = N,N-bis(2-hydroxy-3-5-di-tert-butyl-benzyl)propylamine and H2L2 = 2,2'-selenobis(4,6-di-tert-butylphenol), have been synthesized and characterized by elemental analyses, IR, V-51 NMR, both in the solid and in solution, and cyclic voltammetric studies. Single crystal X-ray studies reveal that in complex 1 the vanadium atom is octahedrally coordinated with an O5N donor environment, where the oxygen atom of the V-V=O moiety and the N atom of the ONO ligand occupy the axial sites while two oxygen atoms (O1 and O2) from the bisphenolate ligand and two oxygen atoms (O3 and O4) from the acac ligand occupy the equatorial plane. A similar bonding pattern has also been encountered for 2 with the exception that a Se atom instead of N is involved in weak bonding to the metal center. Both complexes showed reversible cyclic voltammeric responses and E-1/2 appears at -0.18 and 0.10 V versus NHE for complexes 1 and 2, respectively. The kinetics of oxidation of ascorbic acid by complex 1 were carried out in 50% MeCN-50% HO (v/v) at 25 degrees C. The high formation constant value, Q = 63 +/- 7 M-1, reveals that the reaction proceeds through the rapid formation of a H-bonded intermediate. The low k(2)Q(2)/k(1)Q(1) ratio (13.4) for 1 points out that there is extensive H-bonding between the oxygen atom of the V-V=O group and the OH group of ascorbic acid. (c) 2007 Published by Elsevier Ltd.

Self-nucleation and crystallization kinetics of double crystalline poly(p-dioxanone)-b-poly(epsilon-caprolactone) diblock copolymers

The crystallization kinetics of each constituent of poly(p-dioxanone)-b-poly(epsilon-caprolactone) diblock copolymers (PPDX-b-PCL) has been determined in a wide composition range by differential scanning calorimetry and compared to that of the equivalent homopolymers. Spherulitic growth rates were also measured by polarized optical microscopy while atomic force microscopy was employed to reveal the morphology of one selected diblock copolymer. It was found that crystallization drives structure formation and both components form lamellae within mixed spherulitic superstructures. The overall isothermal crystallization kinetics of the PPDX block at high temperatures, where the PCL is molten, was determined by accelerating the kinetics through a previous self-nucleation procedure. The application of the Lauritzen and Ho. man theory to overall growth rate data yielded successful results for PPDX and the diblock copolymers. The theory was applied to isothermal overall crystallization of previously self-nucleated PPDX ( where growth should be the dominant factor if self-nucleation was effective) and the energetic parameters obtained were perfectly matched with those obtained from spherulitic growth rate data of neat PPDX. A quantitative estimate of the increase in the energy barrier for crystallization of the PPDX block, caused by the covalently bonded molten PCL as compared to homo-PPDX, was thus determined. This energy increase can dramatically reduce the crystallization rate of the PPDX block as compared to homo-PPDX. In the case of the PCL block, both the crystallization kinetics and the self-nucleation results indicate that the PPDX is able to nucleate the PCL within the copolymers and heterogeneous nucleation is always present regardless of composition. Finally, preliminary results on hydrolytic degradation showed that the presence of relatively small amounts of PCL within PPDX-bPCL copolymers substantially retards hydrolytic degradation of the material in comparison to homo-PPDX. This increased resistance to hydrolysis is a complex function of composition and its knowledge may allow future prediction of the lifetime of the material for biomedical applications.

Quantitative model for the kinetics of lyotropic phase transitions involving changes in monolayer curvature

We present a kinetic model for transformations between different self-assembled lipid structures. The model shows how data on the rates of phase transitions between mesophases of different geometries can be used to provide information on the mechanisms of the transformations and the transition states involved. This can be used, for example, to gain an insight into intermediate structures in cell membrane fission or fusion. In cases where the monolayer curvature changes on going from the initial to the final mesophase, we consider the phase transition to be driven primarily by the change in the relaxed curvature with pressure or temperature, which alters the relative curvature elastic energies of the two mesophase structures. Using this model, we have analyzed previously published kinetic data on the inter-conversion of inverse bicontinuous cubic phases in the 1-monoolein-30 wt% water system. The data are for a transition between QII(G) and QII(D) phases, and our analysis indicates that the transition state more closely resembles the QII(D) than the QII(G) phase. Using estimated values for the monolayer mean curvatures of the QII(G) and QII(D) phases of -0.123 nm(-1) and -0.133 nm(-1), respectively, gives values for the monolayer mean curvature of the transition state of between -0.131 nm(-1) and -0.132 nm(-1). Furthermore, we estimate that several thousand molecules undergo the phase transition cooperatively within one "cooperative unit", equivalent to 1-2 unit cells of QII(G) or 4-10 unit cells of QII(D).

Thermochemical kinetics: does it still give insights?

This tutorial review revisits the subject of the seminal book written by Sidney Benson in 1968. A short summary of the nature of the subject is presented, including its place in the wider world of quantitative chemistry. A number of themes are selected to illustrate its previous and continuing usefulness in evaluating numerical values of important quantities, and probing ideas of reaction mechanism. These include strain enthalpies for biradical combination, chain reactions, why some reactions don't occur and the involvement of carbenes in hydrocarbon rearrangements.

Formation of crystalline macrocyclic phases during electrophilic precipitation-polycondensation syntheses of poly(arylene ether ketone)s

Elongated crystalline particles formed as by-products during poly(arylene ether ketone) synthesis by electrophilic precipitation-polycondensation of 4,4'-diphenoxybenzophenone with terephthaloyl chloride or isophthaloyl chloride, thought previously to be polymer-whiskers, have now been identified as macrocyclic phases. Single crystal X-ray analysis of the needle-like particles formed in the reaction with terephthaloyl chloride, using the microdiffraction technique with synchrotron radiation, revealed that they consist of a macrocylic compound containing ten phenylene units, i.e. the [2 + 2] cyclic dimer. An analogous structure has also been demonstrated for the corresponding macrocycle derived from the reaction of 4,4-diphenoxybenzophenone with isophthaloyl chloride. Chloroform extraction of the products of the two polycondensations dissolved the macrocyclic material (but not the linear polymer), and analysis of the extracts by MALDI-TOF mass spectrometry demonstrated the presence in both cases of homologous families of macrocyclic products. Higher yields of macrocycles were obtained under pseudo-high dilution conditions, enabling the [2 + 2] cyclodimers from reactions of 4,4'-diphenoxybenzophenone with both terephthaloyl and isophthaloyl chloride to be isolated as pure compounds and fully characterised. (C) 2003 Published by Elsevier Ltd.

A comparison of two time-domain analysis procedures in the determination of VO2 kinetics by pseudorandom binary sequence exercise testing

The purpose of this study was to apply and compare two time-domain analysis procedures in the determination of oxygen uptake (VO2) kinetics in response to a pseudorandom binary sequence (PRBS) exercise test. PRBS exercise tests have typically been analysed in the frequency domain. However, the complex interpretation of frequency responses may have limited the application of this procedure in both sporting and clinical contexts, where a single time measurement would facilitate subject comparison. The relative potential of both a mean response time (MRT) and a peak cross-correlation time (PCCT) was investigated. This study was divided into two parts: a test-retest reliability study (part A), in which 10 healthy male subjects completed two identical PRBS exercise tests, and a comparison of the VO2 kinetics of 12 elite endurance runners (ER) and 12 elite sprinters (SR; part B). In part A, 95% limits of agreement were calculated for comparison between MRT and PCCT. The results of part A showed no significant difference between test and retest as assessed by MRT [mean (SD) 42.2 (4.2) s and 43.8 (6.9) s] or by PCCT [21.8 (3.7) s and 22.7 (4.5) s]. Measurement error (%) was lower for MRT in comparison with PCCT (16% and 25%, respectively). In part B of the study, the VO2 kinetics of ER were significantly faster than those of SR, as assessed by MRT [33.4 (3.4) s and 39.9 (7.1) s, respectively; P<0.01] and PCCT [20.9 (3.8) s and 24.8 (4.5) s; P < 0.05]. It is possible that either analysis procedure could provide a single test measurement Of VO2 kinetics; however, the greater reliability of the MRT data suggests that this method has more potential for development in the assessment Of VO2 kinetics by PRBS exercise testing.

Modelling the uptake and growth kinetics of Penicillium glabrum in a tannic acid-containing solid-state fermentation for tannase production

A mathematical growth model for the batch solid-state fermentation process for fungal tannase production was developed and tested experimentally. The unstructured model describes the uptake and growth kinetics of Penicillium glabrum in an impregnated polyurethane foam substrate system. In general, good agreement between the experimental data and model simulations was obtained. Biomass, tannase and spore production are described by logistic kinetics with a time delay between biomass production and tannase and spore formation. Possible induction mechanisms for the latter are proposed. Hydrolysis of tannic acid, the main carbon source in the substrate system, is reasonably well described with Michaelis-Menten kinetics with time-varying enzyme concentration but a more complex reaction mechanism is suspected. The metabolism of gallic acid, a tannase-hydrolysis product of tannic acid, was shown to be growth limiting during the main growth phase. (c) 2004 Elsevier Ltd. All rights reserved.

Kinetics of high pressure facilitated starch gelatinisation in Thai glutinous rice

The combined effect of pressure and temperature on the rate of gelatinisation of starch present in Thai glutinous rice was investigated. Pressure was found to initiate gelatinisation when its value exceeded 200 MPa at ambient temperature. On the other hand, complete gelatinisation was observed at 500 and 600 MPa at 70 degrees C, when the rice was soaked in water under these conditions for 120 min. A first-order kinetic model describing the rate of gelatinisation was developed to estimate the values of the rate constants as a function of pressure and temperature in the range: 0.1-600 MPa and 20-70 degrees C. The model, based on the well-known Arrhenius and Eyring equations, assumed the form [GRAPHICS] The constants k(0), E-a, and Delta V were found to take values: 31.19 s(-1), 37.89 kJ mol(-1) and -9.98 cm(3) mol(-1), respectively. It was further noted that the extent of gelatinisation occurring at any time, temperature and pressure, could be exclusively correlated with the grain moisture content. (c) 2006 Elsevier Ltd. All rights reserved.