Nearest-neighbor distance distributions and self-ordering in diffusion-controlled reactions. I. A+A simulations

Steady state of a two-species annihilation process with separated reactants

We describe the steady state of the annihilation process of a one-dimensional system of two initially separated reactants A and B. The parameters that define the dynamical behavior of the system are the diffusion constant, the reaction rate, and the deposition rate. Depending on the ratio between those parameters, the system exhibits a crossover between a diffusion-limited (DL) regime and a reaction-limited (RL) regime. We found that a key quantity to describe the reaction process in the system is the probability p(x_{A},x_{B}) to find the rightmost A (RMA) particle and the leftmost B (LMB) particle at the positions x_{A} and x_{B}, respectively. The statistical behavior of the system in both regimes is described using the density of particles, the gap length distribution x_{B}-x_{A}, the marginal probabilities p_{A}(x_{A}) and p_{B}(x_{B}), and the reaction kernel. For both regimes, this kernel can be approximated by using p(x_{A},x_{B}). We found an excellent agreement between the numerical and analytical results for all calculated quantities despite the reaction process being quite different in both regimes. In the DL regime, the reaction kernel can be approximated by the probability to find the RMA and LMB particles in adjacent sites. In the RL regime, the kernel depends on the marginal probabilities p_{A}(x_{A}) and p_{B}(x_{B}).

Fractal symmetry of protein interior: what have we learned?

The application of fractal dimension-based constructs to probe the protein interior dates back to the development of the concept of fractal dimension itself. Numerous approaches have been tried and tested over a course of (almost) 30 years with the aim of elucidating the various facets of symmetry of self-similarity prevalent in the protein interior. In the last 5 years especially, there has been a startling upsurge of research that innovatively stretches the limits of fractal-based studies to present an array of unexpected results on the biophysical properties of protein interior. In this article, we introduce readers to the fundamentals of fractals, reviewing the commonality (and the lack of it) between these approaches before exploring the patterns in the results that they produced. Clustering the approaches in major schools of protein self-similarity studies, we describe the evolution of fractal dimension-based methodologies. The genealogy of approaches (and results) presented here portrays a clear picture of the contemporary state of fractal-based studies in the context of the protein interior. To underline the utility of fractal dimension-based measures further, we have performed a correlation dimension analysis on all of the available non-redundant protein structures, both at the level of an individual protein and at the level of structural domains. In this investigation, we were able to separately quantify the self-similar symmetries in spatial correlation patterns amongst peptide-dipole units, charged amino acids, residues with the π-electron cloud and hydrophobic amino acids. The results revealed that electrostatic environments in the interiors of proteins belonging to 'α/α toroid' (all-α class) and 'PLP-dependent transferase-like' domains (α/β class) are highly conducive. In contrast, the interiors of 'zinc finger design' ('designed proteins') and 'knottins' ('small proteins') were identified as folds with the least conducive electrostatic environments. The fold 'conotoxins' (peptides) could be unambiguously identified as one type with the least stability. The same analyses revealed that peptide-dipoles in the α/β class of proteins, in general, are more correlated to each other than are the peptide-dipoles in proteins belonging to the all-α class. Highly favorable electrostatic milieu in the interiors of TIM-barrel, α/β-hydrolase structures could explain their remarkably conserved (evolutionary) stability from a new light. Finally, we point out certain inherent limitations of fractal constructs before attempting to identify the areas and problems where the implementation of fractal dimension-based constructs can be of paramount help to unearth latent information on protein structural properties.

Spatial fluctuations and anomalous reaction order in a reactive scheme involving a cooperative full desorption

Anomalous reaction rates have been found in the hydrogen desorption of H-terminated surfaces in semiconductor epitaxy, with a reaction order shifting from two to one, or even taking fractional values. We analyze the issue in terms of a cooperative full desorption (CFD) reaction A+A--k3-->S+S , coupled to an adsorption reaction S--k1-->A and an alternative desorption route A--k2-->S . Steady state properties of the three-step reactive scheme are analyzed in a one-dimensional lattice in the absence of diffusion. Microscopic Monte Carlo simulations show anomalous spatial distributions of reactants in the stationary state: depending on the reaction rate constants of the overall scheme, either a local "aggregation" or a local "dispersion" of A -particles is observed. The CFD reaction itself is well described by a fractional order kinetics that takes into account these anomalies and that depends on the kinetic rate constants of the overall adsorption-desorption reaction mechanism. The problem is addressed with an analytical approach for the active neighborhood of a reactant, which provides a closed expression of the reaction order as a function of the kinetic parameters. This approach is in excellent agreement with numerical simulations. Spatial correlations, as well as fluctuation correlations, are also formalized in terms of the kinetic constants. We discuss the results in the context of the hydrogen evolution reaction on silicon surfaces.

Reaction front in an A+B-->C reaction-subdiffusion process

We study the reaction front for the process A+B-->C in which the reagents move subdiffusively. Our theoretical description is based on a fractional reaction-subdiffusion equation in which both the motion and the reaction terms are affected by the subdiffusive character of the process. We design numerical simulations to check our theoretical results, describing the simulations in some detail because the rules necessarily differ in important respects from those used in diffusive processes. Comparisons between theory and simulations are on the whole favorable, with the most difficult quantities to capture being those that involve very small numbers of particles. In particular, we analyze the total number of product particles, the width of the depletion zone, the production profile of product and its width, as well as the reactant concentrations at the center of the reaction zone, all as a function of time. We also analyze the shape of the product profile as a function of time, in particular, its unusual behavior at the center of the reaction zone.

Monte carlo simulations of enzyme reactions in two dimensions: fractal kinetics and spatial segregation

Conventional equations for enzyme kinetics are based on mass-action laws, that may fail in low-dimensional and disordered media such as biological membranes. We present Monte Carlo simulations of an isolated Michaelis-Menten enzyme reaction on two-dimensional lattices with varying obstacle densities, as models of biological membranes. The model predicts that, as a result of anomalous diffusion on these low-dimensional media, the kinetics are of the fractal type. Consequently, the conventional equations for enzyme kinetics fail to describe the reaction. In particular, we show that the quasi-stationary-state assumption can hardly be retained in these conditions. Moreover, the fractal characteristics of the kinetics are increasingly pronounced as obstacle density and initial substrate concentration increase. The simulations indicate that these two influences are mainly additive. Finally, the simulations show pronounced S-P segregation over the lattice at obstacle densities compatible with in vivo conditions. This phenomenon could be a source of spatial self organization in biological membranes.

Reaction front structure in the diffusion-limited A+B model with initially randomized reactants

Subtle features of the reaction front formation in the A+B-->0 reaction are reported for the initially random and equal A+B reactant distribution. Three nonclassical parameters (initial linewidth, minimum, and maximum), for each interparticle gap and nearest neighbor distance distributions, are derived, as a function of time, using Monte Carlo simulations. These empirical front measures and their temporal scaling exponents are compared with the previously studied ones for the reactant interparticle distributions.

Exciton microscopy and reaction kinetics in restricted spaces

We describe the development of a new biologically non-invasive ultraresolution light microscopy, based on combining the energy transfer "spectral ruler" method with the micro-movement technology employed in scanning tunneling microscopy (STM). We use near-field scanning optical microscopy, with micropipettes containing crystals of energy packaging donor molecules in the tips that can have apertures below 5 nm. The excitation of these tips extends near field microscopy well beyond the 50 nm limit. The theoretical resolution limit for this spectrally sensitive light microscopy is well below 1 nm. Exciton microscopy is ideally suited for kinetic studies that are spatially resolved on the molecular scale, i.e., at a single molecule site. Moreover, the successful operation of the scanning exciton tip depends on an understanding of reaction kinetics in restricted spaces. In contrast to the many recent reviews on scanning tip microscopies, there is no adequate review of the recent revolutionary developments in the area of reaction kinetics in confined geometries. We thus attempt such a review in this paper. Reactions in restricted spaces rarely get stirred vigorously by convection and are thus often controlled by diffusion. Furthermore, the compactness of the Brownian motion leads to both anomalous diffusion and anomalous reaction kinetics. Elementary binary reactions of the type A + A----Products, A + B----Products and A + C----C + Products are discussed theoretically for both batch and steady-state conditions. The anomalous reaction orders and time exponents (for the rate coefficients) are discussed for various situations. Global and local rate laws are related to particle distribution functions. Only Poissonian distributions guarantee the classical rate laws. Reactant self-organization leads to interesting new phenomena. These are demonstrated by theory, simulations, and experiments. The correlation length of reactant production affects the self-ordering length-scale. These effects are demonstrated experimentally, including the stability of reactant segregation observed in chemical reactions in one-dimensional spaces, e.g., capillaries and microcapillaries. The gap between the reactant A (cation) and B (anion) actually increases in time, and extends over millimeters. Excellent agreement is found among theory, simulation, and experiment for the various scaling exponents.