Stochastic simulation of chemical reactions with spatial resolution and single molecule detail

Methods are presented for simulating chemical reaction networks with a spatial resolution that is accurate to nearly the size scale of individual molecules. Using an intuitive picture of chemical reaction systems, each molecule is treated as a point-like particle that diffuses freely in three-dimensional space. When a pair of reactive molecules collide, such as an enzyme and its substrate, a reaction occurs and the simulated reactants are replaced by products. Achieving accurate bimolecular reaction kinetics is surprisingly difficult, requiring a careful consideration of reaction processes that are often overlooked. This includes whether the rate of a reaction is at steady-state and the probability that multiple reaction products collide with each other to yield a back reaction. Inputs to the simulation are experimental reaction rates, diffusion coefficients and the simulation time step. From these are calculated the simulation parameters, including the 'binding radius' and the 'unbinding radius', where the former defines the separation for a molecular collision and the latter is the initial separation between a pair of reaction products. Analytic solutions are presented for some simulation parameters while others are calculated using look-up tables. Capabilities of these methods are demonstrated with simulations of a simple bimolecular reaction and the Lotka-Volterra system.     

Reaction rate calculations for cosubstrates diffusing into catalyst layer from opposite sides

The effects of diffusion on a reaction taking place in a permeable catalyst are examined theoretically for the case where the reaction has two substrates supplied from opposite sides of a catalytic slab. The solutions of the reaction-diffusion equation for combinations of zeroth-and first-order kinetics are given in terms of an effectiveness factor and a parameter describing the position in the layer where the reaction occurs. In these terms, the results vary only weakly with reaction order. The use of the exact solutions for a reaction that is zeroth order in both sub strates is proposed as a general rule for estimating the reaction rate and the reaction position.     

Inter-chromophore interactions in pigment-modified and dimer-less bacterial photosynthetic reaction centers

The magnitudes of inter-chromophore interactions in bacterial photosynthetic reaction centers are investigated by measuring absorption and Stark spectra of reaction centers in which monomeric chromophores are modified and in a novel triplet mutant which lacks the special pair. The circular dichroism spectrum of the triple mutant reaction center was also measured. Only small changes in the spectroscopic properties are observed, as has also been found for several types of reaction centers in which the absorption or chemical properties of a chromophore are altered by site-specific mutations. We conclude that the electronic absorption, circular dichroism and Stark features of the special pair and the monomeric chromophores in the reaction center are relatively insensitive to inter-chromophore interactions.     

An autocatalytic step in the reaction cycle of hydrogenase from Thiocapsa roseopersicina can explain the special characteristics of the enzyme reaction

A moving front has been observed as a special pattern during the hydrogenase-catalyzed reaction of hydrogen uptake with benzyl viologen as electron acceptor in a thin-layer reaction chamber. Such fronts start spontaneously and at random times at different points of the reaction chamber; blue spheres are seen expanding at constant speed and amplitude. The number of observable starting points depends on the hydrogenase concentration. Fronts can be initiated by injecting either a small amount of completed reaction mixture or activated hydrogenase, but not by injecting a low concentration of reduced benzyl viologen. These characteristics are consistent with an autocatalytic reaction step in the enzyme reaction. The special characteristics of the hydrogen-uptake reaction in the bulk reaction (a long lag phase, and the enzyme concentration dependence of the lag phase) support the autocatalytic nature. We conclude that there is at least one autocatalytic reaction step in the hydrogenase-catalyzed reaction. The two possible autocatalytic schemes for hydrogenase are prion-type autocatalysis, in which two enzyme forms interact, and product-activation autocatalysis, where a reduced electron acceptor and an inactive enzyme form interact. The experimental results strongly support the occurrence of prion-type autocatalysis.     

Influence of the conserved active site residues of histidyl tRNA synthetase on the mechanism of aminoacylation reaction

The relation between the conservation of active site residues and the molecular mechanism of aminoacylation reaction is an unexplored problem. In the present paper, the influences of the conserved active site residues on the reaction mechanism as well as the electrostatic potential near the reaction center are analyzed for Histidyl tRNA synthetase from Escherichia coli, Thermus thermophilus and Staphylococcus aureus. While the primary structures show both convergence as well as divergence, the secondary level structures of the active sites of the three species show considerable conservation in the respective structural organizations. The conserved active site residues near the reaction center, which have a major role in the reaction mechanism and catalysis, retain their specific position and orientation relative to the substrate in the three species. In order to understand the influence of different conserved and nonconserved residues near the reaction center, two different models are considered. First, a large model of active site with the substrates, Mg²⁺ ions and water is constructed in which the first shell residues (including both conserved as well as nonconserved) near the reaction center are studied. From the large model, a smaller model is constructed for reaction path modeling individually for three species. Validation of the smaller model is carried out by comparing the energy surfaces of large and small models as a function of reaction coordinates. Further, the electrostatic potential near the reaction center for the large and small model are compared. The transition state structures of the activation step of aminoacylation reaction for E. coli, T. thermophilus and S. aureus are calculated using the combined ab-initio/semi-empirical calculation. The similarity of the energy profiles as a function of the relevant reaction coordinate and the orientation of the catalytic residue, Arg259, indicate that the reaction mechanisms are identical which are guided by the strikingly similar structural pattern formed by conserved residues for three species. The energy surfaces have close resemblance in three species and present a clear perspective that how the reaction proceeds with the aid of different conserved residues. The study of electrostatic potential confirms this view. The present study provides an understanding of the relationship between the conservation of residues and the efficient reaction mechanism of aminoacylation reaction.     

Inversion of the anomeric configuration of the transferred sugar during inactivation of the macrolide antibiotic oleandomycin catalyzed by a macrolide glycosyltransferase

Macrolides are a group of antibiotics structurally characterized by a macrocyclic lactone to which one or several deoxy-sugar moieties are attached. The sugar moieties are transferred to the different aglycones by glycosyltransferases (GTF). The OleI GTF of an oleandomycin producer, Streptomyces antibioticus, catalyzes the inactivation of this macrolide by glycosylation. The product of this reaction was isolated and its structure elucidated. The donor substrate of the reaction was UDP-alpha-D-glucose, but the reaction product showed a beta-glycosidic linkage. The inversion of the anomeric configuration of the transferred sugar and other data about the kinetics of the reaction and primary structure analysis of several GTFs are compatible with a reaction mechanism involving a single nucleophilic substitution at the sugar anomeric carbon in the catalytic center of the enzyme.     

Investigation of Embedded Si/C System Exposed to a Hybrid Reaction of Centrifugal-Assisted Thermite Method

A novel method is proposed to study the behavior and phase formation of a Si+C compacted pellet under centrifugal acceleration in a hybrid reaction. Si+C as elemental mixture in the form of a pellet is embedded in a centrifugal tube. The pellet assembly and tube are exposed to the sudden thermal energy of a thermite reaction resulted in a hybrid reaction. The hybrid reaction of thermite and Si+C produced unique phases. X-ray diffraction pattern (XRD) as well as microstructural and elemental analyses are then investigated. XRD pattern showed formation of materials with possible electronic and magnetic properties. The cooling rate and the molten particle viscosity mathematical model of the process are meant to assist in understanding the physical and chemical phenomena took place during and after reaction. The results analysis revealed that up to 85% of materials converted into secondary products as ceramics-matrix composite.     

Effect of dynamic channelling on the transient-state kinetics of coupled enzyme reactions.

Evidence is presented to show that the presence of a channelled reaction branch in a coupled two-enzyme reaction leads to a slower overall relaxation of the system towards steady state following perturbations of the reaction conditions. A predominant part of the total reaction flux change resulting from perturbations may be governed by transients that are more rapid in the presence than in the absence of a channel, but only under specific conditions that are not likely to be of biological relevance.     

Oscillatory kinetics of the peroxidase-oxidase reaction in an open system. Experimental and theoretical studies.

1. The oscillations in the peroxidase (donor: hydrogen-peroxide oxidoreductase, EC 1.11.1.7)-catalyzed reaction between NADH and O2 are undamped when the reaction is carried out in a system open to both substrates and when 2,4-dichlorophenol and methylene blue are present in the solution. 2. The waveform of the oscillations changes when the concentration of peroxidase is varied. 3. The waveforms obtained experimentally can be simulated by a branched chain reaction model in which the branching is quadratic. 4. A correlation between the present knowledge of the reaction and the model can be made by combining well established and hypothetical reaction steps into a few reaction schemes. A selection among schemes however, is not possible at the present time. 5. Compound III participates in the reaction as an active intermediate. This is possible because dichlorophenol stimulates the break down of compound III.     

Oscillations in peroxidase-catalyzed reactions and their potential function in vivo

The peroxidase-oxidase reaction has become a model system for the study of oscillations and complex dynamics in biochemical systems. In the present paper we give an overview of previous experimental and theoretical studies of the peroxidase-oxidase reaction. Recent in vitro experiments have raised the question whether the reaction also exhibits oscillations and complex dynamics in vivo. To investigate this possibility further we have undertaken new experimental studies of the reaction, using horseradish extracts and phenols which are widely distributed in plants. The results are discussed in light of the occurrence and a possible functional role of oscillations and complex dynamics of the peroxidase-oxidase reaction in vivo.     

Freeze fracture studies of reaction centers from Rhodospirillum rubrum in chromatophores and liposomes

In freeze-fractures of chromatophores of Rhodospirillum rubrum the reaction centers are seen as hexagonal arranged particles of 13 nm diameter with a density of around 5,500 particles per m². Similar regions on the cytoplasmic membrane suggest that these parts are the prospective invagination sites.Isolated reaction centers are easily incorporated into liposomes. In freeze fractures of liposomes particles similar in shape and size, although less dense as in chromatophores are observed. In negative staining much smaller units of only 5 nm in diameter are found indicating that reaction centers occur in the membrane as tri- or tetramers. There is a strong correlation between particle density in chromatophores and titratable reaction centers remaining in these membranes after extraction of reaction centers by detergents; both values are in good agreement with the yield of reaction centers at a given detergent concentration.Abbreviations LDAO Lauryldimethylamine oxide - PF protoplasmic fracture face - EF exoplasmic fracture face     

Real time measurement and control of thermodynamic water activities for enzymatic catalysis in hexane

The esterification reaction of geraniol with acetic acid catalyzed by immobilized Candida antarctica lipase B was studied in hexane using a pervaporation-assisted batch reactor. The effect of thermodynamic water activity (a(w)) on the initial reaction rate was investigated at a(w) ranging from 0.02 to 1.0. The a(w) was monitored on-line in real time. a(w) was actively controlled throughout the reaction by using highly water-selective membrane pervaporation. This novel combination of a(w) sensing and control eliminates changes in a(w) during the reaction even in the initial phase of relatively rapid water release during an esterification. No chemicals are introduced for a(w) control, and no purge gases or liquids are needed. A maximum in the initial reaction rate was found approximately at a(w)=0.1. The initial reaction rate declined quickly at higher a(w), and dropped precipitously at lower a(w).     

Utilization of peroxide and its relevance in oxygen insertion reactions catalyzed by chloroperoxidase.

Chloroperoxidase (CPO) catalyzed oxygen insertions are highly enantioselective and hence of immense biotechnological potential. A peroxide activation step is required to give rise to the compound I species that catalyzes this chiral reaction. A side reaction, a catalase type peroxide dismutation, is another feature of CPO's versatility. This work systematically investigates the utilization of different peroxides for the two reactions, i.e. the catalase type reaction and the oxygen insertion reaction. For the oxygen insertion reaction, indene and phenylethyl sulfide were chosen as substrate models for epoxidation and sulfoxidation respectively. The results clearly show that CPO is stable towards hydrogen peroxide and has a total number of turnovers near one million prior to deactivation. The epoxidation reactions terminate before completion because the enzyme functioning in its catalatic mode quickly removes all of the hydrogen peroxide from the reaction mixture. Sulfoxidation reactions are much faster than epoxidation reactions and thus are better able to compete with the catalase reaction for hydrogen peroxide utilization. A preliminary study towards optimizing the reaction system components for a laboratory scale synthetic epoxidation is reported.     

Multicomponent reactions in crop protection chemistry

An overview is given of the significance of multicomponent reactions in the synthesis of agrochemicals. The most important applications of multicomponent condensations, such as the Biginelli reaction, Bucherer-Bergs reaction, Hantzsch dihydropyridine synthesis, Kabachnik-Fields reaction, Mannich reaction, Passerini reaction, Petasis reaction, Strecker reaction, Ugi reaction and Willgerodt-Kindler reaction, to the synthesis of herbicidally, fungicidally and insecticidally active compounds are presented. Also the mode of action and biological activity of these multicomponent reaction products are reported.     

Fourier transform infrared spectroscopy of special pair bacteriochlorophylls in homodimeric reaction centers of heliobacteria and green sulfur bacteria

Heliobacteria and green sulfur bacteria have type I homodimeric reaction centers analogous to photosystem I. One remaining question regarding these homodimeric reaction centers is whether the structures and electron transfer reactions are truly symmetric or not. This question is relevant to the origin of the heterodimeric reaction centers, such as photosystem I and type II reaction centers. In this mini-review, Fourier transform infrared studies on the special pair bacteriochlorophylls, P798 in heliobacteria and P840 in green sulfur bacteria, are summarized. The data are reinterpreted in the light of the X-ray crystallographic structure of photosystem I and the sequence alignments of type I reaction center proteins, and discussed in terms of hydrogen bonding interactions and the symmetry of charge distribution over the dimer.     

Effect of calcium ions on creatine kinase systems of myocardial cells

The kinetics of the creatine phosphokinase reaction catalyzed by different creatine phosphokinase (CPK) isoenzymes from the heart and the mechanism of the regulatory action of calcium ions on this reaction have been studied. It has been shown that the kinetic parameters of the reaction in the presence of calcium are similar to those in the presence of magnesium with the exception of the maximal rate of the reverse reaction, which is about three times lower in the presence of calcium. Calcium ions are able to exert a significant regulatory action on the CPK reaction. Simultaneous regulation of the CPK reaction by calcium and magnesium ions can be quantitatively described by a kinetic model, which takes into account the formation of complexes of adenine nucleotides with calcium and magnesium. The mechanism of the CPK reaction regulation by metal ions involves changes in concentrations of the metal -- adenine nucleotide compexes. The regulatory action of calcium on the creatine kinase reaction is the same for all CPK isoenzymes studied.     

A novel method to highly versatile monomeric PNA building blocks by multi component reactions.

A novel approach to monomeric PNA building blocks by a solution phase Ugi multi component reaction (MCR) is described. The reaction is easily performed in 96 well plates. The products precipitate from the reaction solution and are thus obtained in high yields and purity. Those products are not amenable by another multi step synthesis in such a diversity.     

One-electron reduction of flavodoxin. A fast kinetic study.

The reaction kinetics of fully oxidized flavodoxin from Clostridium MP with the hydrated electron have been investigated by the pulse radiolysis method. Four spectrally distinct processes have been observed with the ultimate formation of the singly reduced flavin form of the protein. The last two species obtained in the reaction sequence are spectrally similar, and are connected through a reaction which is first order. It is proposed that this reaction involves a protein conformational alteration.     

Effects of pressure and temperature on enzymatic reactions in supercritical fluids

Supercritical fluids (SCFs) are receiving increasing attention as reaction media because they permit higher reaction rates compared with the conventional solvents. The ease of manipulating the physical properties of the SCFs enables easier control of the reaction conditions and easier solvent removal after the reaction. This review focuses on effects of pressure, temperature and the properties of the SCFs, on enzymatic reactions. Phase behavior, reaction rate and activation volume in SCFs are discussed. Within the ranges of pressure (10-40 MPa) and temperature (35-60 degrees C) that typically characterize the supercritical region, an increase in pressure and/or a decrease in temperature lead to a decrease in the enzyme turnover because the diffusion coefficients of the substrates migrating to the active sites of enzymes are affected.     

Flavin-photosensitized reactions of retinal and stilbene

1. Retinol and stilbene are both isomerized when they are illuminated anaerobically in the presence of flavins. 2. Triplet quenchers (e.g. oxygen, potassium iodide and paramagnetic ions) inhibit the reaction more efficiently than they quench flavin fluorescence. At 77 degrees C in a diethyl ether-isopentane-ethanol (5:5:2) glass retinol quenches flavin phosphorescence, but not its fluorescence. 3. For the stilbene reaction cis/trans photostationary-state mixtures are obtained with different flavins and these are linearly related to the phosphorescence transition energies of the flavins used. 4. The reaction involves the triplet state of flavin and a scheme for the reaction is suggested. 5. The dependence of the rate of reaction on substrate concentration is explicable in terms of this scheme. 6. The photobleaching of rhodopsin sensitized by flavin is also demonstrated.