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. (C) 2004 Elsevier Ltd. All rights reserved.

Detection of banana bunchy top virus in virus-infected plants using polymerase chain reaction

Polymerase chain reaction (PCR) was developed for the detection of Banana bunchy top virus (BBTV) at maximum after 210 min and at minimum after 90 min using Pc-1 and Pc-2, respectively. PCR detection of BBTV in crude sap indicated that the freezing of banana tissue in liquid nitrogen (LN2) before extraction was more effective than using sand as the extraction technique. BBTV was also detected using PCR assay in 69 healthy and diseased plants using Na-PO4 buffer containing 1 % SDS. PCR detection of BBTV in nucleic acid extracts using seven different extraction buffers to adapt the use of PCR in routine detection in the field was studied. Results proved that BBTV was detected with high sensitivity in nucleic acid extracts more than in infectious sap. The results also suggested the common aetiology for the BBTV by the PCR reactions of BBTV in nucleic acid extracts from Australia, Burundi, Egypt, France, Gabon, Philippines and Taiwan. Results also proved a positive relation between the Egyptian-BBTV isolate and abaca bunchy top isolate from the Philippines, but there no relation was found with the Cucumber mosaic cucumovirus (CMV) isolates from Egypt and Philippines and Banana bract mosaic virus (BBMV) were found.

A preconditioned Lanczos method for space-fractional reaction-diffusion equations on two-dimensional unstructured meshes

We consider a two-dimensional space-fractional reaction diffusion equation with a fractional Laplacian operator and homogeneous Neumann boundary conditions. The finite volume method is used with the matrix transfer technique of Ilić et al. (2006) to discretise in space, yielding a system of equations that requires the action of a matrix function to solve at each timestep. Rather than form this matrix function explicitly, we use Krylov subspace techniques to approximate the action of this matrix function. Specifically, we apply the Lanczos method, after a suitable transformation of the problem to recover symmetry. To improve the convergence of this method, we utilise a preconditioner that deflates the smallest eigenvalues from the spectrum. We demonstrate the efficiency of our approach for a fractional Fisher’s equation on the unit disk.

Addition of diazomethane to armchair single-walled carbon nanotubes and their reaction sequences : A computational study

The sidewall additions of diazomethane to (n, n), n = 3–10 armchair single-walled carbon nanotubes (SWCNTs) on two different orientations of C–C bonds have been studied using the ONIOM(B3LYP/6-31G(d):PM3) approach. The binding energies of SWCNTs complexes with CH2N2, CH2 and their transition-state structures were computed at the B3LYP/6-31G(d) level. The effects of diameters of armchair SWCNTs on their binding energies were studied. Relative reactivities of all the SWCNTs and their complexes based on their frontier orbital energies gaps are reported.

Insights into the mechanism of the reaction between ttrachloro-p-bnzoquinone and hydrogen peroxide and their implications in the catalytic role of water molecules in producing the hydroxyl radial

Detailed mechanisms for the formation of hydroxyl or alkoxyl radicals in the reactions between tetrachloro-p-benzoquinone (TCBQ) and organic hydroperoxides are crucial for better understanding the potential carcinogenicity of polyhalogenated quinones. Herein, the mechanism of the reaction between TCBQ and H2O2 has been systematically investigated at the B3LYP/6-311++G** level of theory in the presence of different numbers of water molecules. We report that the whole reaction can easily take place with the assistance of explicit water molecules. Namely, an initial intermediate is formed first. After that, a nucleophilic attack of H2O2 onto TCBQ occurs, which results in the formation of a second intermediate that contains an OOH group. Subsequently, this second intermediate decomposes homolytically through cleavage of the O-O bond to produce a hydroxyl radical. Energy analyses suggest that the nucleophilic attack is the rate-determining step in the whole reaction. The participation of explicit water molecules promotes the reaction significantly, which can be used to explain the experimental phenomena. In addition, the effects of F, Br, and CH3 substituents on this reaction have also been studied.

A combined temperature-programmed reaction spectroscopy and Fourier-transform infrared spectroscopy study of CO2/H2 and CO/CO2/H2 interactions with model ZnO/SiO2, Cu/SiO2 and Cu/ZnO/SiO2 methanol-synthesis catalysts

The reaction of CO2 and H2 with ZnO/SiO2 catalyst at 295 K gave predominantly hydrogencarbonate on zinc oxide and a small quantity of formate was evolved after heating at 393 K. Elevation of the reaction temperature to 503 K enhanced the rate of formation of zinc formate species. Significantly these formate species decomposed at 573 K almost entirely to CO2 and H2. Even after exposure of CO2-H2 or CO-CO2-H2 mixtures to highly defected ZnO/SiO2 catalyst, the formate species produced still decomposed to give CO2 and H2. It was concluded that carboxylate species which were formed at oxygen anion vacancies on polar Zn planes were not significantly hydrogenated to formate. Consequently it was proposed that the non-polar planes on zinc oxide contained sites which were specific for the synthesis of methanol. The interaction of CO2 and H2 with reduced Cu/ZnO/SiO2 catalyst at 393 K gave copper formate species in addition to substantial quantities of formate created at interfacial sites between copper and zinc oxide. It was deduced that interfacial formate species were produced from the hydrogenation of interfacial bidentate carbonate structures. The relevance of interfacial formate species in the methanol synthesis reaction is discussed. Experiments concerning the reaction of CO2-H2 with physical mixtures of Cu/SiO2 and ZnO/SiO2 gave results which were simply characteristic of the individual components. By careful consideration of previous data a detailed proposal regarding the role of spillover hydrogen is outlined. Admission of CO to a gaseous CO2-H2 feedstock resulted in a considerably diminished amount of formate species on copper. This was ascribed to a combination of over-reduction of the surface and site-blockage.

Ingenuity pathways analysis of urine metabonomics phenotypes toxicity of gentamicin in multiple organs

We introduce the use of Ingenuity Pathway Analysis to analyzing global metabonomics in order to characterize phenotypically biochemical perturbations and the potential mechanisms of the gentamicin-induced toxicity in multiple organs. A single dose of gentamicin was administered to Sprague Dawley rats (200 mg/kg, n = 6) and urine samples were collected at -24-0 h pre-dosage, 0-24, 24-48, 48-72 and 72-96 h post-dosage of gentamicin. The urine metabonomics analysis was performed by UPLC/MS, and the mass spectra signals of the detected metabolites were systematically deconvoluted and analyzed by pattern recognition analyses (Heatmap, PCA and PLS-DA), revealing a time-dependency of the biochemical perturbations induced by gentamicin toxicity. As result, the holistic metabolome change induced by gentamicin toxicity in the animal's organisms was characterized. Several metabolites involved in amino acid metabolism were identified in urine, and it was confirmed that gentamicin biochemical perturbations can be foreseen from these biomarkers. Notoriously, it was found that gentamicin induced toxicity in multiple organs system in the laboratory rats. The proof-of-knowledge based Ingenuity Pathway Analysis revealed gentamicin induced liver and heart toxicity, along with the previously known toxicity in kidney. The metabolites creatine, nicotinic acid, prostaglandin E2, and cholic acid were identified and validated as phenotypic biomarkers of gentamicin induced toxicity. Altogether, the significance of the use of metabonomics analyses in the assessment of drug toxicity is highlighted once more; furthermore, this work demonstrated the powerful predictive potential of the Ingenuity Pathway Analysis to study of drug toxicity and its valuable complementation for metabonomics based assessment of the drug toxicity.

Common PPARγ variants C161T and Pro12Ala are not associated with inflammatory bowel disease in an Australian cohort

Background & Aims: Peroxisome proliferator-activated receptor (PPAR) γ is a transcription factor, highly expressed in colonic epithelial cells, adipose tissue and macrophages, with an important role in the regulation of inflammatory pathways. The common PPARγ variants C161T and Pro12Ala have recently been associated with Ulcerative Colitis (UC) and an extensive UC phenotype respectively, in a Chinese population. PPARγ Pro12Ala variant homozygotes appear to be protected from the development of Crohn's disease (CD) in European Caucasians. Methods: A case-control study was performed for both variants (CD n=575, UC n=306, Controls n=360) using a polymerase chain reaction (PCR)-restriction fragment length polymorphism analysis in an Australian IBD cohort. A transmission disequilibrium test was also performed using CD trios for the PPARγ C161T variant. Genotype-phenotype analyses were also undertaken. Results: There was no significant difference in genotype distribution data or allele frequency between CD and UC patients and controls. There was no difference in allele transmission for the C161T variant. No significant relationship between the variants and disease location was observed. Conclusions: We were unable to replicate in a Caucasian cohort the recent association between PPARγ C161T and UC or between PPARγ Pro12Ala and an extensive UC phenotype in a Chinese population. There are significant ethnic differences in genetic susceptibility to IBD and its phenotypic expression.

Analysis of Sdhd and Mmp12 in an affected solar keratosis and control cohort

The incidence of Squamous Cell Carcinoma (SCG) is growing in certain populations to the extent that it is now the most common skin lesion in young men and women in high ultraviolet exposure regions such as Queensland. In terms of incidence up to 40% of the Australian population over 40 years of age is thought to possess the precancerous Solar Keratosis (SK) lesion and with a small, but significant, chance of progression into SCC, understanding the genetic events that play a role in this process is essential. The major aims of this study were to analyse whole blood derived samples for DNA aberrations in genes associated with tumour development and cellular maintenance, with the ultimate aim of identifying genes associated with non-melanoma skin cancer development. More specifically the first aim of this project was to analyse the SDHD and MMP12 genes via Dual-Labelled Probe Real-Time PCR for copy number aberrations in an affected Solar Keratosis and control cohort. It was found that 12 samples had identifiable copy-number aberrations in either the SDHD or MMP12 gene (this means that a genetic section of either of these two genes is aberrantly amplified or deleted), with five of the samples exhibiting aberrations in both genes. The significance of this study is the contribution to the knowledge of the genetic pathways that are malformed in the progression and development of the pre-cancerous skin lesion Solar Keratosis. © 2008 Springer Science+Business Media, LLC.

Cyclin A/CDK2 regulates multiple G2 phase pathways

The cell cycle is a carefully choreographed series of phases that when executed successfully will allow the complete replication of the genome and the equal division of the genome and other cellular content into two independent daughter cells. The inability of the cell to execute cell division successfully can result in either checkpoint activation to allow repair and/or apoptosis and/or mutations/errors that may or may not lead to tumourgenesis. Cyclin A/CDK2 is the primary cyclin/CDK regulating G2 phase progression of the cell cycle. Cyclin A/CDK2 activity peaks in G2 phase and its inhibition causes a G2 phase delay that we have termed 'the cyclin A/CDK2 dependent G2 delay'. Understanding the key pathways that are involved in the cyclin A/CDK2 dependent G2 delay has been the primary focus of this study. Characterising the cyclin A/CDK2 dependent G2 delay revealed accumulated levels of the inactive form of the mitotic regulator, cyclin B/CDK1. Surprisingly, there was also increased microtubule nucleation at the centrosomes, and the centrosomes stained for markers of cyclin B/CDK1 activity. Both microtubule nucleation at the centrosomes and phosphoprotein markers were lost with short-term treatment of CDK1/2 inhibition. Cyclin A/CDK2 localised at the centrosomes in late G2 phase after separation of the centrosomes but before the start of prophase. Thus G2 phase cyclin A/CDK2 controls the timing of entry into mitosis by controlling the subsequent activation of cyclin B/CDK1, but also has an unexpected role in coordinating the activation of cyclin B/CDK1 at the centrosome and in the nucleus. In addition to regulating the timing of cyclin B/CDK1 activation and entry into mitosis in the unperturbed cell cycle, cyclin A/CDK2 also was shown to have a role in G2 phase checkpoint recovery. Known G2 phase regulators were investigated to determine whether they had a role in imposing the cyclin A/ CDK2 dependent G2 delay. Examination of the critical G2 checkpoint arrest protein, Chk1, which also has a role during unperturbed G2/M phases revealed the presence of activated Chk1 in G2 phase, in a range of cell lines. Activated Chk1 levels were shown to accumulate in cyclin A/CDK2 depleted/inhibited cells. Further investigations revealed that Chk1, but not Chk2, depletion could reverse the cyclin A/CDK2 dependent G2 delay. It was confirmed that the accumulative activation of Chk1 was not a consequence of DNA damage induced by cyclin A depletion. The potential of cyclin A/CDK2 to regulate Chk1 revealed that the inhibitory phosphorylations, Ser286 and Ser301, were not directly catalysed by cyclin A/CDK2 in G2 phase to regulate mitotic entry. It appeared that the ability of cyclin A/CDK2 to regulate cyclin B/CDK1 activation impacted cyclin B/CDK1s phosphorylation of Chk1 on Ser286 and Ser301, thereby contributing to the delay in G2/M phase progression. Chk1 inhibition/depletion partially abrogated the cyclin A/CDK2 dependent G2 delay, and was less effective in abrogating G2 phase checkpoint suggesting that other cyclin A/CDK2 dependent mechanisms contributed to these roles of cyclin A/CDK2. In an attempt to identify these other contributing factors another G2/M phase regulator known to be regulated by cyclin A/CDK2, Cdh1 and its substrates Plk1 and Claspin were examined. Cdh1 levels were reduced in cyclin A/CDK2 depleted/inhibited cells although this had little effect on Plk1, a known Cdh1 substrate. However, the level of another substrate, Claspin, was increased. Cdh1 depletion mimicked the effect of cyclin A depletion but to a weaker extent and was sufficient at increasing Claspin levels similar to the increase caused by cyclin A depletion. Co-depletion of cyclin A and Claspin blocked the accumulation of activated Chk1 normally seen with cyclin A depletion alone. However Claspin depletion alone did not reduce the cyclin A/CDK2 dependent G2 delay but this is likely to be a result of inhibition of S phase roles of Claspin. Together, these data suggest that cyclin A/CDK2 regulates a number of different mechanisms that contribute to G2/M phase progression. Here it has been demonstrated that in normal G2/M progression and possibly to a lesser extent in G2 phase checkpoint recovery, cyclin A/CDK2 regulates the level of Cdh1 which in turn affects at least one of its substrates, Claspin, and consequently results in the increased level of activated Chk1 observed. However, the involvement of Cdh1 and Claspin alone does not explain the G2 phase delay observed with cyclin A/CDK2 depletion/inhibition. It is likely that other mechanisms, possibly including cyclin A/CDK2 regulation of Wee1 and FoxM1, as reported by others, combine with the mechanism described here to regulate normal G2/M phase progression and G2 phase checkpoint recovery. These findings support the critical role for cyclin A/CDK2 in regulating progression into mitosis and suggest that upstream regulators of cyclin A/CDK2 activation will also be critical controllers of this cell cycle transition. The pathways that work to co-ordinate cell cycle progression are very intricate and deciphering these pathways, required for normal cell cycle progression, is key to understanding tumour development. By understanding cell cycle regulatory pathways it will allow the identification of the pathway/s and their mechanism/s that become affected in tumourgenesis. This will lead to the development of better targeted therapies, inferring better efficacy with fewer side effects than commonly seen with the use of traditional therapies, such as chemotherapy. Furthermore, this has the potential to positively impact the development of personalised medicines and the customisation of healthcare.

A simplified approach for calculating the moments of action for linear reaction-diffusion equations

The mean action time is the mean of a probability density function that can be interpreted as a critical time, which is a finite estimate of the time taken for the transient solution of a reaction-diffusion equation to effectively reach steady state. For high-variance distributions, the mean action time under-approximates the critical time since it neglects to account for the spread about the mean. We can improve our estimate of the critical time by calculating the higher moments of the probability density function, called the moments of action, which provide additional information regarding the spread about the mean. Existing methods for calculating the nth moment of action require the solution of n nonhomogeneous boundary value problems which can be difficult and tedious to solve exactly. Here we present a simplified approach using Laplace transforms which allows us to calculate the nth moment of action without solving this family of boundary value problems and also without solving for the transient solution of the underlying reaction-diffusion problem. We demonstrate the generality of our method by calculating exact expressions for the moments of action for three problems from the biophysics literature. While the first problem we consider can be solved using existing methods, the second problem, which is readily solved using our approach, is intractable using previous techniques. The third problem illustrates how the Laplace transform approach can be used to study coupled linear reaction-diffusion equations.

Effects of mobile phone distraction on drivers’ reaction times

Distraction resulting from mobile phone use whilst driving has been shown to increase the reaction times of drivers, thereby increasing the likelihood of a crash. This study compares the effects of mobile phone conversations on reaction times of drivers responding to traffic events that occur at different points in a driver’s field of view. The CARRS-Q Advanced Driving Simulator was used to test a group of young drivers on various simulated driving tasks including a traffic event that occurred within the driver’s central vision—a lead vehicle braking suddenly—and an event that occurred within the driver’s peripheral—a pedestrian entering a zebra crossing from a footpath. Thirty-two licensed drivers drove the simulator in three phone conditions: baseline (no phone conversation), and while engaged in hands-free and handheld phone conversations. The drivers were aged between 21 to 26 years and split evenly by gender. Differences in reaction times for an event in a driver’s central vision were not statistically significant across phone conditions, probably due to a lower speed selection by the distracted drivers. In contrast, the reaction times to detect an event that originated in a distracted driver’s peripheral vision were more than 50% longer compared to the baseline condition. A further statistical analysis revealed that deterioration of reaction times to an event in the peripheral vision was greatest for distracted drivers holding a provisional licence. Many critical events originate in a driver’s periphery, including vehicles, bicyclists, and pedestrians emerging from side streets. A reduction in the ability to detect these events while distracted presents a significant safety concern that must be addressed.

A parametric duration model of the reaction times of drivers distracted by mobile phone conversations

The use of mobile phones while driving is more prevalent among young drivers—a less experienced cohort with elevated crash risk. The objective of this study was to examine and better understand the reaction times of young drivers to a traffic event originating in their peripheral vision whilst engaged in a mobile phone conversation. The CARRS-Q Advanced Driving Simulator was used to test a sample of young drivers on various simulated driving tasks, including an event that originated within the driver’s peripheral vision, whereby a pedestrian enters a zebra crossing from a sidewalk. Thirty-two licensed drivers drove the simulator in three phone conditions: baseline (no phone conversation), hands-free and handheld. In addition to driving the simulator each participant completed questionnaires related to driver demographics, driving history, usage of mobile phones while driving, and general mobile phone usage history. The participants were 21 to 26 years old and split evenly by gender. Drivers’ reaction times to a pedestrian in the zebra crossing were modelled using a parametric accelerated failure time (AFT) duration model with a Weibull distribution. Also tested where two different model specifications to account for the structured heterogeneity arising from the repeated measures experimental design. The Weibull AFT model with gamma heterogeneity was found to be the best fitting model and identified four significant variables influencing the reaction times, including phone condition, driver’s age, license type (Provisional license holder or not), and self-reported frequency of usage of handheld phones while driving. The reaction times of drivers were more than 40% longer in the distracted condition compared to baseline (not distracted). Moreover, the impairment of reaction times due to mobile phone conversations was almost double for provisional compared to open license holders. A reduction in the ability to detect traffic events in the periphery whilst distracted presents a significant and measurable safety concern that will undoubtedly persist unless mitigated.

Galvanic replacement mediated synthesis of hollow Pt nanocatalysts : significance of residual Ag for the H2 evolution reaction

With the increasing popularity of the galvanic replacement approach towards the development of bimetallic nanocatalysts, special emphasis has been focused on minimizing the use of expensive metal (e.g. Pt), in the finally formed nanomaterials (e.g. Ag/Pt system as a possible catalyst for fuel cells). However, the complete removal of the less active sacrificial template is generally not achieved during galvanic replacement, and its residual presence may significantly impact on the electrocatalytic properties of the final material. Here, we investigate the hydrogen evolution reaction (HER) activity of Ag nanocubes replaced with different amounts of Pt, and demonstrate how the bimetallic composition significantly affects the activity of the alloyed nanomaterial.

A study of the galvanic replacement reaction at surfaces and the role of lateral charge propagation

The galvanic replacement of isolated nanostructures of copper and silver on conducting supports as well as continuous films of copper with gold is reported. The surface morphology was characterized by scanning electron microscopy and the replacement with gold was confirmed by EDX analysis. It was found that lateral charge propagation during the replacement reaction had a significant effect in all cases. For the isolated nanostructures the deposition of gold was observed not only at the sacrificial template but also at the surrounding unmodified areas of the conducting substrate. In the case of copper films the role of lateral charge propagation was also confirmed by connecting it to an ITO electrode through an external circuit upon which gold deposition was also observed to occur. Interestingly, by inhibiting the rate of charge propagation, through the introduction of a series resistor, the morphology of gold on the copper substrate could be changed from discrete surface decoration with cube like nanoparticles to a more porous rough surface.