The effect of temperature on bacterial degradation of EDTA in pH-auxostat

The effect of temperature on the maximum specific growth rate and the cell yield was studied during cultivation of two bacterial strains (LPM-4 and Pseudomonas sp. LPM-410) on EDTA under unlimited cell growth conditions in a pH-auxostat. Both strains displayed linear dependence of reciprocal biomass yield against reciprocal specific growth rate, from which the values of rate of substrate expenditure for cell maintenance and the “maximum” yield (i.e., hypothetical yield without cell maintenance processes) were estimated. Analysis of the maximum yield values based on mass–energy balance theory suggested that oxidation of the carboxylic acid side chains of EDTA by a monooxygenase had zero or low energetic efficiency. An Arrhenius equation with different values of Arrhenius parameters within different temperature ranges gave a good fit with the temperature dependence of both growth rate and biomass yield. Specific growth rates of both strains showed a more pronounced temperature dependence than did the cell yields. A possible kinetic mechanism was suggested which might be responsible for the modes of the temperature dependences of specific growth rate and yield that were found. The mechanism is based on a hypothetical key substance governing the metabolic flows, which is formed in a zero-order reaction and destroyed in a first-order reaction, both rate constants depending on temperature according to the Arrhenius law.     

Models for contact-mediated pattern formation: cells that form parallel arrays

Kinetic continuum models are derived for cells that crawl over a 2D substrate, undergo random reorientation, and turn in response to contact with a neighbor. The integro-partial differential equations account for changes in the distribution of orientations in the population. It is found that behavior depends on parameters such as total mass, random motility, adherence, and sloughing rates, as well as on broad aspects of the contact response. Linear stability analysis, and numerical, and cellular automata simulations reveal that as parameters are varied, a bifurcation leads to loss of stability of a uniform (isotropic) steady state, in favor of an (anisotropic) patterned state in which cells are aligned in parallel arrays.     

Cannibalism as a life boat mechanism

Under certain conditions a cannibalistic population can survive when food for the adults is too scarce to support a non-cannibalistic population. Cannibalism can have this lifeboat effect if (i) the juveniles feed on a resource inaccessible to the adults; and (ii) the adults are cannibalistic and thus incorporate indirectly the inaccessible resource. Using a simple model we conclude that the mechanism works when, at low population densities, the average yield, in terms of new offspring, due to the energy provided by one cannibalized juvenile is larger than one.     

Alterations in the structure of the electrical double layer adjacent to a charged membrane arising from electromagnetically induced changes in the activity of membrane bound enzymes

Electromagnetic fields of very low amplitude have been reported to influence a number of cellular functions. Many of these effects have a high degree of frequency specificity. Herein it is suggested that some of these reported results could be explained by a fieldinduced alteration in the enzymic activity of integral membrane proteins. It is shown that such a field-induced transition from an initial nonequilibrium steady-state to a final nonequilibrium steady-state can lead to an alteration in the concentration profiles of those charged species in the cell's ambient electrolyte that comprise the so-called electrical double layer. Examples of variations in the concentration profiles of those ions that react with a membrane-bound enzyme, as well as nonreacting ionic species, are given. The modulation of such effects by systematic variations in extracellular pH and ionic strength is discussed.     

Numerical solution of a nonlinear advance-delay-differential equation from nerve conduction theory

A functional differential equation which is nonlinear and involves forward and backward deviating arguments is solved numerically. The equation models conduction in a myelinated nerve axon in which the myelin completely insulates the membrane, so that the potential change jumps from node to node. The equation is of first order with boundary values given at t=±. The problem is approximated via a difference scheme which solves the problem on a finite interval by utilizing an asymptotic representation at the endpoints, cubic interpolation and iterative techniques to approximate the delays, and a continuation method to start the procedure. The procedure is tested on a class of problems which are solvable analytically to access the scheme's accuracy and stability, then applied to the problem that models propagation in a myelinated axon. The solution's dependence on various model parameters of physical interest is studied. This is the first numerical study of myelinated nerve conduction in which the advance and delay terms are treated explicitly.Supported in part by NSF Grant MCS8301724 and by a Biomedical Research Support Grant 2SO7RR0706618 from NIH     

Resetting biological oscillators—A stochastic approach

We present a stochastic approach to phase-resetting of an ensemble of oscillators. In order to describe stimulation-induced dynamical phenomena we develop a stochastic model which consists of an ensemble of phase oscillators interacting via random forces. Every single oscillator is submitted to a phase stimulus. The ensemble's dynamics is determined by a Fokker-Planck equation. The stationary states are calculated explicitly, whereas the transients are analysed numerically. If the stimulus of a given (non-vanishing) intensity is administered at a critical initial cluster phase for a critical duration T _crit the ensemble's synchronized oscillation is annihilated. A transition from type 1 resetting to type 0 resetting occurs when the stimulation duration exceeds T _crit. Stimulation causes a shift of the mean frequency of every single oscillator. This frequency shift is explicitly calculated by deriving the mean first passage time. The model shows that there is a subcritical intensity which is connected with an enhanced vulnerability to stimulation. The desynchronized states, the so-called black holes, are typically associated with a double peak in the ensemble's phase distribution. This is important for analysing experimental data because simple peak-detection algorithms are not able to extract the underlying dynamics.Our results are discussed from the experimentator's point of view so that the insights derived from our model can improve data analysis and design of stimulation experiments.     

EFFECTS OF FLOWING WATER ON NITROGEN- AND PHOSPHORUS-LIMITED PHOTOSYNTHESIS AND OPTIMUM N:P RATIOS BY SPIROGYRA FLUVIATILIS (CHAROPHYCEAE)1,2

The effects of flowing water on net photosynthesis, dark respiration, specific growth rate, and optimum N:P ratios by Spirogyra fluviatilis Hilse were assessed. The alga was cultivated under nitrogen or phosphorus limitation in laboratory streams at three flow velocities: 3, 12, and 30 cm·s^?1. The Droop equation adequately described respiration and photosynthesis (PS_net) as a function of N or P cell quota (Q^N or Q^p). The data show that for N- or P-limited Spirogyra fluviatilis, flowing water is physiologically costly. Generally, flowing water had little effect on respiration rates; however, the proportion of gross photosynthesis devoted to dark respiration did increase with flow velocity. For photosynthesis, the minimum N and P cell quotas increased with velocity, and the theoretical PS_net maxima for N and P both appeared greatest at 12 cm·s^?1. The Droop models showed that for any given Q^N or Q^p, PS_net, was reduced by the 30-cm·s^?1 treatment. Consistent with this finding, independent estimates of specific growth rates for P-limited S. fluviatilis in the laboratory streams were inversely related to flow velocity when ambient PO₄^?3 was undetectable. However, growth was not diminished at the fastest velocity when PO₄^?3 was available for uptake. Thus, the increase in cellular phosphorus demand can be offset by flow-enhanced P uptake when conditions permit; otherwise, growth will be impaired. The optimum N:P ratios for S. fluviatilis at 3, 12, and 30 cm·s^?1 were 50, 58, and 52 by atoms, respectively, when calculated for PS_net= 0. The optimum ratios were inversely related to PS_net and decreased to approximately 20 when PS_net was near maximum. The potential for flowing water to mediate nutrient partitioning among lotic algae by altering growth rates and optimum nutrient ratios is discussed.     

Do salt bridges stabilize proteins? A continuum electrostatic analysis

The electrostatic contribution to the free energy of folding was calculated for 21 salt bridges in 9 protein X-ray crystal structures using a continuum electrostatic approach with the DELPHI computer-program package. The majority (17) were found to be electrostatically destabilizing; the average free energy change, which is analogous to mutation of salt bridging side chains to hydrophobic isosteres, was calculated to be 3.5 kcal/mol. This is fundamentally different from stability measurements using pKa shifts, which effectively measure the strength of a salt bridge relative to 1 or more charged hydrogen bonds. The calculated effect was due to a large, unfavorable desolvation contribution that was not fully compensated by favorable interactions within the salt bridge and between salt-bridge partners and other polar and charged groups in the folded protein. Some of the salt bridges were studied in further detail to determine the effect of the choice of values for atomic radii, internal protein dielectric constant, and ionic strength used in the calculations. Increased ionic strength resulted in little or no change in calculated stability for 3 of 4 salt bridges over a range of 0.1-0.9 M. The results suggest that mutation of salt bridges, particularly those that are buried, to "hydrophobic bridges" (that pack at least as well as wild type) can result in proteins with increased stability. Due to the large penalty for burying uncompensated ionizable groups, salt bridges could help to limit the number of low free energy conformations of a molecule or complex and thus play a role in determining specificity (i.e., the uniqueness of a protein fold or protein-ligand binding geometry).     

Evaluation of Fourier-transform infrared spectroscopy for data capture in predictive microbiology

Chemical changes in the medium, induced by the fermentative species Lactobacillus plantarum and Lactobacillus brevis and by the enzymatic action of a proteolytic, spoilage species, Yarrowia lipolytica, were analysed using Fourier-transform i.r. spectroscopy (FTIR). Changes in the absorbance data over time could be modelled using one of the more current predictive, mathematical models of microbial growth, such as the Gompertz equation. Moreover, a linear correlation between FTIR data (expressed as absorbance of some selected peaks) and viability data (expressed as log₁₀ c.f.u./g or ml) was observed during the fermentation process, both for L. plantarum and L. brevis.     

Coding of odor intensity in a steady-state deterministic model of an olfactory receptor neuron

The coding of odor intensity by an olfactory receptor neuron model was studied under steady-state stimulation. Our model neuron is an elongated cylinder consisting of the following three components: a sensory dendritic region bearing odorant receptors, a passive region consisting of proximal dendrite and cell body, and an axon. First, analytical solutions are given for the three main physiological responses: (1) odorant-dependent conductance change at the sensory dendrite based on the Michaelis-Menten model, (2) generation and spreading of the receptor potential based on a new solution of the cable equation, and (3) firing frequency based on a Lapicque model. Second, the magnitudes of these responses are analyzed as a function of odorant concentration. Their dependence on chemical, electrical, and geometrical parameters is examined. The only evident gain in magnitude results from the activation-to-conductance conversion. An optimal encoder neuron is presented that suggests that increasing the length of the sensory dendrite beyond about 0.3 space constant does not increase the magnitude of the receptor potential. Third, the sensivities of the responses are examined as functions of (1) the concentration at half-maximum response, (2) the lower and upper concentrations actually discriminated, and (3) the width of the dynamic range. The overall gain in sensitivity results entirely from the conductance-to-voltage conversion. The maximum conductance at the sensory dendrite appears to be the main tuning constant of the neuron because it determines the shift toward low concentrations and the increase in dynamic range. The dynamic range of the model cannot exceed 5.7 log units, for a sensitivity increase at low odor concentration is compensated by a sensitivity decrease at high odor concentration.