A multi-scaled approach for simulating chemical reaction systems

In this paper we give an overview of some very recent work, as well as presenting a new approach, on the stochastic simulation of multi-scaled systems involving chemical reactions. In many biological systems (such as genetic regulation and cellular dynamics) there is a mix between small numbers of key regulatory proteins, and medium and large numbers of molecules. In addition, it is important to be able to follow the trajectories of individual molecules by taking proper account of the randomness inherent in such a system. We describe different types of simulation techniques (including the stochastic simulation algorithm, Poisson Runge–Kutta methods and the balanced Euler method) for treating simulations in the three different reaction regimes: slow, medium and fast. We then review some recent techniques on the treatment of coupled slow and fast reactions for stochastic chemical kinetics and present a new approach which couples the three regimes mentioned above. We then apply this approach to a biologically inspired problem involving the expression and activity of LacZ and LacY proteins in E. coli, and conclude with a discussion on the significance of this work.

Gauge Poisson representations for birth/death master equations

Poisson representation techniques provide a powerful method for mapping master equations for birth/death processes -- found in many fields of physics, chemistry and biology -- into more tractable stochastic differential equations. However, the usual expansion is not exact in the presence of boundary terms, which commonly occur when the differential equations are nonlinear. In this paper, a gauge Poisson technique is introduced that eliminates boundary terms, to give an exact representation as a weighted rate equation with stochastic terms. These methods provide novel techniques for calculating and understanding the effects of number correlations in systems that have a master equation description. As examples, correlations induced by strong mutations in genetics, and the astrophysical problem of molecule formation on microscopic grain surfaces are analyzed. Exact analytic results are obtained that can be compared with numerical simulations, demonstrating that stochastic gauge techniques can give exact results where standard Poisson expansions are not able to.

Modelling the fate of chromated copper arsenate in a sandy soil

A pulse of chromated copper arsenate (CCA, a timber preservative) was applied in irrigation water to an undisturbed field soil in a laboratory column. Concentrations of various elements in the leachate from the column were measured during the experiment. Also, the remnants within the soil were measured at the end of the experiment. The geochemical modelling package, PHREEQC-2, was used to simulate the experimental data. Processes included in the CCA transport modelling were advection, dispersion, non-specific adsorption (cation exchange) and specific adsorption by clay minerals and organic matter, as well as other possible chemical reactions such as precipitation/dissolution. The modelling effort highlighted the possible complexities in CCA transport and reaction experiments. For example, the uneven dosing of CCA as well as incomplete knowledge of the soil properties resulted in simulations that gave only partial, although reasonable, agreement with the experimental data. Both the experimental data and simulations show that As and Cu are strongly adsorbed and therefore, will mostly remain at the top of the soil profile, with a small proportion appearing in leachate. On the other hand, Cr is more mobile and thus it is present in the soil column leachate. Further simulations show that both the quantity of CCA added to the soil and the pH of the irrigation water will influence CCA transport. Simulations suggest that application of larger doses of CCA to the soil will result in higher leachate concentrations, especially for Cu and As. Irrigation water with a lower pH will dramatically increase leaching of Cu. These results indicate that acidic rainfall or significant accidental spillage of CCA will increase the risk of groundwater pollution.

Exploring the potential of xenobiotic-metabolising enzymes as biocatalysts: Evolving designer catalysts from polyfunctional cytochrome P450 enzymes

1. Biological catalysts have the advantage of being able to catalyse chemical reactions with an often exquisite degree of regio- and stereospecificity in contrast with traditional methods of organic synthesis. 2. The cytochrome P450 enzymes involved in human drug metabolism are ideal starting materials for the development of designer biocatalysts by virtue of their catalytic versatility and extreme substrate diversity. Applications can be envisaged in fine chemical synthesis, such as in the pharmaceutical industry and bioremediation. 3. A variety of techniques of enzyme engineering are currently being applied to P450 enzymes to explore their catalytic potential. Although most studies to date have been performed with bacterial P450s, reports are now emerging of work with mammalian forms of the enzymes. 4. The present minireview will explore the rationale and general techniques for redesigning P450s, review the results obtained to date with xenobiotic-metabolising forms and discuss strategies to overcome some of the logistic problems limiting the full exploitation of these enzymes as industrial-scale biocatalysts.

Modern developments in hypersonic wind tunnels

The development of new methods of producing hypersonic wind-tunnel flows at increasing velocities during the last few decades is reviewed with attention to airbreathing propulsion, hypervelocity aerodynamics and superorbital aerodynamics. The role of chemical reactions in these flows leads to use of a binary scaling simulation parameter, which can be related to the Reynolds number, and which demands that smaller wind tunnels require higher reservoir pressure levels for simulation of flight phenomena. The use of combustion heated vitiated wind tunnels for propulsive research is discussed, as well as the use of reflected shock tunnels for the same purpose. A flight experiment validating shock-tunnel results is described, and relevant developments in shock tunnel instrumentation are outlined. The use of shock tunnels for hypervelocity testing is reviewed, noting the role of driver gas contamination in determining test time, and presenting examples of air dissociation effects on model flows. Extending the hypervelocity testing range into the superorbital regime with useful test times is seen to be possible by use of expansion tube/tunnels with a free piston driver.

Cobalt complexes of tripodal hexadentate ligands: Electrochemically driven rearrangements

The Co-III complexes of the hexadentate tripodal ligands HOsen (3-(2'-aminoethylamino)-2,2-bis((2 ''-aminoethylamino) methyl) propan-1-ol) and HOten (3-(2'-aminoethylthia)-2,2-bis((2 ''-aminoethylthia) methyl) propan-1-ol) have been synthesized and fully characterized. The crystal structures of [Co(HOsen)]Cl-3 center dot H2O and [Co(HOten)](ClO4)Cl-2 are reported and in both cases the ligands coordinate as tripodal hexadentate N-6 and N3S3 donors, respectively. Cyclic voltammetry of the N3S3 coordinated complex [Co(HOten)](3+) is complicated and electrode dependent. On a Pt working electrode an irreversible Co-III/II couple ( formal potential - 157 mV versus Ag-AgCl) is seen, which is indicative of dissociation of the divalent complex formed at the electrode. The free HOten released by the dissociation of [Co(HOten)](2+) can be recaptured by Hg as shown by cyclic voltammetry experiments on a static Hg drop electrode ( or in the presence of Hg2+ ions), which leads to the formation of an electroactive Hg-II complex of the N3S3 ligand (formal potential + 60 mV versus Ag-AgCl). This behaviour is in contrast to the facile and totally reversible voltammetry of the hexaamine complex [Co(HOsen)](3+) ( formal potential (Co-III/II) - 519 mV versus Ag-AgCl), which is uncomplicated by any coupled chemical reactions. Akinetic and thermodynamic analysis of the [Co(HOten)](2+)/[Hg(HOten)](2+) system is presented on the basis of digital simulation of the experimental voltammetric data.

Aspects of SAF smelting of ferrochrome

This paper provides a review of the chemical reactions occurring in the submerged arc processing of chromite ores. The excavations of industrial furnaces have shown that the charge, as it descends through the furnace, passes through a number of distinct reaction zones. Each zone is characterised by differing process conditions and reaction products. The phase equilibria, reaction steps and mechanisms occurring as the charge progresses through the furnace are examined, and the potential influences of these factors on the process outcomes are discussed.

Reactions catalyzed by bacterial cytochromes P450

The cytochromes P450 are a large family of oxidative haemoproteins that are responsible for a wide variety of oxidative transformations in a variety of organisms. This review focuses upon the reactions catalyzed specifically by bacterial enzymes, which includes aliphatic hydroxylation, alkene epoxidation, aromatic hydroxylation, oxidative phenolic coupling, heteroatom oxidation and dealkylation, and multiple oxidations including C-C bond cleavage. The potential for the practical application of the oxidizing power of these enzymes is briefly discussed.

New insights into NO-carbon and N2O-carbon reactions from quantum mechanical calculations

Density functional theory calculations were used to investigate the mechanisms of NO-carbon and N2O-carbon reactions. It was the first time that the importance of surface nitrogen groups was addressed in the kinetic behaviors of the NO-carbon reaction. It was found that the off-plane nitrogen groups that are adjacent to the zigzag edge sites and in-plane nitrogen groups that are located on the armchair sites make the bond energy of oxygen desorption even ca. 20% lower than that of the off-plane epoxy group adjacent to zigzag edge sites and in-plane o-quinone oxygen atoms on armchair sites; this may explain the reason why the experimentally obtained activation energy of the NO-carbon reaction is ca. 20% lower than that of the O-2-carbon reaction over 923 K. A higher ratio of oxygen atoms can be formed in the N2O-carbon reaction, because of the lower dissociation energy of N2O, which results in a higher ratio of off-plane epoxy oxygen atoms. The desorption energy of semiquinone with double adjacent off-plane oxygen groups is ca. 20% less than that of semiquinone with only one adjacent off-plane oxygen group. This may be the reason why the activation energy of N2O is also ca. 20% less than that of the O-2-carbon reaction. The new mechanism can also provide a good qualitative comparison for the relative reaction rates of NO-, N2O-, and O-2-carbon reactions. The anisotropic characters of these gas-carbon reactions can also be well explained.

Reactions of fac-[PtMe2(OMe)(H2O)3]+ with halide ions: effect of halide trans effect on methoxide hydrolysis

A solution of fac-[PtMe2(OMe)(H2O)(3)](+) (1) in aqueous perchloric acid underwent very slow hydrolysis of the Pt-OMe bond, over many, weeks. When chloride was added to a solution of 1, two interconverting isomers of [PtMe2(OMe)Cl(H2O)(2)] (with chloride trans to methyl) were formed, and with excess chloride, [PtMe2(OMe)Cl-2(H2O)](-) (both chloride ligands trans to methyl). This solution was stable at ambient temperature, but on heating, methanol was formed and [PtMe2Cl2(H2O)(2)] (both chloride ligands cis to methyl) was produced in the solution. It is proposed that this reaction proceeds via an intermediate complex with chloride bound trans to methoxide. Concentration gave solid [{PtMe2Cl2}n], whose identity was confirmed by conversion to [PtMe(2)Cl(2)py(2)] (pyridine, py, trans to methyl). With bromide and iodide, methoxide hydrolysis occurred at ambient temperature, more slowly with bromide than with iodide, to form solid [{PtMe2X2}(n)] without significant concentrations of [PtMe2X2(H2O)(2)] formed as an intermediate. The greater tendency for Pt-OMe bond to hydrolyse trans to halide compared with 1 was ascribed to the higher trans effect of the halide ligand compared with that of water. (C) 2003 Elsevier Science B.V. All rights reserved.

Catalytic properties of heteropolyacids supported on MCM-41 mesoporous silica for hydrocarbon cracting reactions

It is known that MCM-41 structures have very weak acid sites because of the lack of the bridging hydroxyl groups present in zeolites. Strong acidity however is required for the potential use of these materials in some specific applications such as: cracking and hydrotreating of heavy residue molecules, cracking of waste plastic, etc. The acidity enhancement of the MCM-41 materials was assessed using the n-hexane and polyethylene cracking reactions. MCM-41 samples were impregnated using heteropolyacid (HPA) such as tungestophospheric acid. The catalyst samples were characterized also by x-ray diffraction and benzene adsorption.

Nucleophilic substitution reactions of pyranose polytosylates

The 2,3,4-tri-toluenesulfonate ester derivatives of the methyl pyranosides of L-arabinose, D-ribose, D-lyxose, and D-xylose have been prepared, and their substitution reactions with various nucleophiles have been examined. For arabinose, xylose, and ribose, highly regioselective monosubstitutions were observed with benzoate, nitrite, and azide anions. These reactions have led to short and simple routes from D-xylose to L-arabinose derivatives, from L-arabinose to D-xylose derivatives, and from D-ribose to L-lyxose derivatives. The tritosylate derived from methyl alpha-D-lyxopyranoside was unreactive toward nucleophilic substitution reactions, giving instead a dihydropyran product arising from an initial E2 elimination reaction of the 2-tosylate.

A potential role for isothermal calorimetry in studies of the effects of thermodynamic non-ideality in enzyme-catalyzed reactions

Attention is drawn to the feasibility of using isothermal calorimetry for the characterization of enzyme reactions under conditions bearing greater relevance to the crowded biological environment, where kinetic parameters are likely to differ significantly from those obtained by classical enzyme kinetic studies in dilute solution. An outline of the application of isothermal calorimetry to the determination of enzyme kinetic parameters is followed by considerations of the nature and consequences of crowding effects in enzyme catalysis. Some of those effects of thermodynamic non-ideality are then illustrated by means of experimental results from calorimetric studies of the effect of molecular crowding on the kinetics of catalysis by rabbit muscle pyruvate kinase. This review concludes with a discussion of the potential of isothermal calorimetry for the experimental determination of kinetic parameters for enzymes either in biological environments or at least in media that should provide reasonable approximations of the crowded conditions encountered in vivo. Copyright (C) 2004 John Wiley Sons, Ltd.

Transition-state-theory calculations for reactions of O(3P) with halogenated olefins

The potential energy surfaces for the reactions of atomic oxygen in its ground electronic state, O(P-3), with the olefins: CF2=CCl2 and CF2=CF - CF3, have been characterized using ab initio molecular orbital calculations. Geometry optimization and vibrational frequency calculations were performed for reactants, transition states and products at the MP2 and QCISD levels of theory using the 6-31G(d) basis set. This database was then used to calculate the rate constants by means of Transition-State-Theory. To obtain a better reference and to test the reliability of the activation barriers we have also carried out computations using the CCSD(T)(fc)/6-311Gdagger, MP4(SDQ)(fc)/CBSB4 and MP2(fc)/CBSB3 single point energy calculations at both of the above levels of theory, as well as with the composite CBS-RAD procedure ( P. M. Mayer, C. J. Parkinson, D. M. Smith and L. Radom, J. Chem. Phys., 1998, 108, 604) and a modi. cation of this approach, called: CBS-RAD( MP2, MP2). It was found that the kinetic parameters obtained in this work particularly with the CBS-RAD ( MP2, MP2) procedure are in reasonable agreement with the experimental values. For both reactions it is found that the channels leading to the olefin double-bond addition predominates with respect to any other reaction pathway. However, on account of the different substituents in the alkenes we have located, at all levels of theory, two transition states for each reaction. Moreover, we have found that, for the reactions studied, a correlation exists between the activation energies and the electronic structure of the transition states which can explain the influence of the substituent effect on the reactivity of the halo-olefins.

Intra- vs intermolecular photoinduced electron transfer reactions of a macrocyclic donor-acceptor dyad

The synthesis, structural characterization, and photophysical behavior of a 14-membered tetraazamacrocycle with pendant 4-dimethylaminobenzyl (DMAB) and 9-anthracenylmethyl groups is reported (L-3, 6-((9-anthracenylmethyl)amino)-trans-6,13-dimethyl-13-((4-dimethylaminobenzyl)amino)-1,4,8,11-tetraaza-cyclotetradecane). In its free base form, this compound displays rapid intramolecular photoinduced electron transfer (PET) quenching of the anthracene emission, with both the secondary amines and the DMAB group capable of acting as electron donors. When complexed with Zn(II), the characteristic fluorescence of the anthracene chromophore is restored as the former of these pathways is deactivated by coordination. Importantly, it is shown that the DMAB group, which remains uncoordinated and PET active, acts only very weakly to quench emission, by comparison to the behavior of a model Zn complex lacking the pendant DMAB group, [ZnL2](2+) (Chart 1). By contrast, Stern-Volmer analysis of intermolecular quenching of [ZnL2](2+) by N,N-dimethylaniline (DMA) has shown that this reaction is diffusion limited. Hence, the pivotal role of the bridge in influencing intramolecular PET is highlighted.