Stochastic self-tuning hybrid algorithm for reaction-diffusion systems.

Many biochemical phenomena involve reactants with vastly different concentrations, some of which are amenable to continuum-level descriptions, while the others are not. We present a hybrid self-tuning algorithm to model such systems. The method combines microscopic (Brownian) dynamics for diffusion with mesoscopic (Gillespie-type) methods for reactions and remains efficient in a wide range of regimes and scenarios with large variations of concentrations. Its accuracy, robustness, and versatility are balanced by redefining propensities and optimizing the mesh size and time step. We use a bimolecular reaction to demonstrate the potential of our method in a broad spectrum of scenarios: from almost completely reaction-dominated systems to cases where reactions rarely occur or take place very slowly. The simulation results show that the number of particles present in the system does not degrade the performance of our method. This makes it an accurate and computationally efficient tool to model complex multireaction systems.

Assessment of Insecticidal Efficacy of Diatomaceous Earth and Powders of Common Lavender and Field Horsetail against Bean Weevil Adults.

In the search for an effective and sustainable control method against the bean weevil Acanthoscelides obtectus (Say), an important insect pest affecting stored common beans and other legumes, three different powders were tested against adult been weevils under laboratory conditions. The three powders were diatomaceous earth (DE) (commercial product SilicoSec®), common lavender (Lavandula angustifolia) powder and field horsetail (Equisetum arvense) powder. The substances were tested at five temperatures (15, 20, 25, 30, and 35°C), two relative humidity levels (RH) (55 and 75%), and four concentrations (100, 300, 500, and 900 ppm). The mortality of adults was measured after the 1st, 2nd, 4th, and 7th days of exposure. The efficacy of the powders increased with the temperature, whereas in general, RH did not have a significant effect on the adults' survival. According to common practice of storing common beans, we recommend the use of DE against the pest in question, as this inert powder showed the highest efficacy at lower temperatures and concentrations. Concerning the wider use of common lavender and field horsetail powders, we suggest studying their combined use with other environmentally friendly methods with the aim of achieving the highest synergistic effect possible.

Miniature endplate current rise times less than 100 microseconds from improved dual recordings can be modeled with passive acetylcholine diffusion from a synaptic vesicle.

We recorded miniature endplate currents (mEPCs) using simultaneous voltage clamp and extracellular methods, allowing correction for time course measurement errors. We obtained a 20-80% rise time (tr) of approximately 80 micros at 22 degrees C, shorter than any previously reported values, and tr variability (SD) with an upper limit of 25-30 micros. Extracellular electrode pressure can increase tr and its variability by 2- to 3-fold. Using Monte Carlo simulations, we modeled passive acetylcholine diffusion through a vesicle fusion pore expanding radially at 25 nm x ms(-1) (rapid, from endplate omega figure appearance) or 0.275 nm x ms(-1) (slow, from mast cell exocytosis). Simulated mEPCs obtained with rapid expansion reproduced tr and the overall shape of our experimental mEPCs, and were similar to simulated mEPCs obtained with instant acetylcholine release. We conclude that passive transmitter diffusion, coupled with rapid expansion of the fusion pore, is sufficient to explain the time course of experimentally measured synaptic currents with trs of less than 100 micros.

Monte Carlo simulation of miniature endplate current generation in the vertebrate neuromuscular junction.

A Monte Carlo method for modeling the neuromuscular junction is described in which the three-dimensional structure of the synapse can be specified. Complexities can be introduced into the acetylcholine kinetic model used with only a small increase in computing time. The Monte Carlo technique is shown to be superior to differential equation modeling methods (although less accurate) if a three-dimensional representation of synaptic geometry is desired. The conceptual development of the model is presented and the accuracy estimated. The consequences of manipulations such as varying the spacing of secondary synaptic folds or that between the release of multiple quantal packets of acetylcholine, are also presented. Increasing the spacing between folds increases peak current. Decreased spacing of adjacent quantal release sites increases the potentiation of peak current.