Unifying Time to Contact Estimation and Collision Avoidance across Species

The [Formula: see text]-function and the [Formula: see text]-function are phenomenological models that are widely used in the context of timing interceptive actions and collision avoidance, respectively. Both models were previously considered to be unrelated to each other: [Formula: see text] is a decreasing function that provides an estimation of time-to-contact (ttc) in the early phase of an object approach; in contrast, [Formula: see text] has a maximum before ttc. Furthermore, it is not clear how both functions could be implemented at the neuronal level in a biophysically plausible fashion. Here we propose a new framework - the corrected modified Tau function - capable of predicting both [Formula: see text]-type ("[Formula: see text]") and [Formula: see text]-type ("[Formula: see text]") responses. The outstanding property of our new framework is its resilience to noise. We show that [Formula: see text] can be derived from a firing rate equation, and, as [Formula: see text], serves to describe the response curves of collision sensitive neurons. Furthermore, we show that [Formula: see text] predicts the psychophysical performance of subjects determining ttc. Our new framework is thus validated successfully against published and novel experimental data. Within the framework, links between [Formula: see text]-type and [Formula: see text]-type neurons are established. Therefore, it could possibly serve as a model for explaining the co-occurrence of such neurons in the brain.     

Trans-Theta Logistics: A New Family of Population Growth Sigmoid Functions

Sigmoid functions have been applied in many areas to model self limited population growth. The most popular functions; General Logistic (GL), General von Bertalanffy (GV), and Gompertz (G), comprise a family of functions called Theta Logistic ([Formula: see text] L). Previously, we introduced a simple model of tumor cell population dynamics which provided a unifying foundation for these functions. In the model the total population (N) is divided into reproducing (P) and non-reproducing/quiescent (Q) sub-populations. The modes of the rate of change of ratio P/N was shown to produce GL, GV or G growth. We now generalize the population dynamics model and extend the possible modes of the P/N rate of change. We produce a new family of sigmoid growth functions, Trans-General Logistic (TGL), Trans-General von Bertalanffy (TGV) and Trans-Gompertz (TG)), which as a group we have named Trans-Theta Logistic (T [Formula: see text] L) since they exist when the [Formula: see text] L are translated from a two parameter into a three parameter phase space. Additionally, the model produces a new trigonometric based sigmoid (TS). The [Formula: see text] L sigmoids have an inflection point size fixed by a single parameter and an inflection age fixed by both of the defining parameters. T [Formula: see text] L and TS sigmoids have an inflection point size defined by two parameters in bounding relationships and inflection point age defined by three parameters (two bounded). While the Theta Logistic sigmoids provided flexibility in defining the inflection point size, the Trans-Theta Logistic sigmoids provide flexibility in defining the inflection point size and age. By matching the slopes at the inflection points we compare the range of values of inflection point age for T [Formula: see text] L versus [Formula: see text] L for model growth curves.     

Limits to the rate of adaptive substitution in sexual populations

In large populations, many beneficial mutations may be simultaneously available and may compete with one another, slowing adaptation. By finding the probability of fixation of a favorable allele in a simple model of a haploid sexual population, we find limits to the rate of adaptive substitution, [Formula: see text], that depend on simple parameter combinations. When variance in fitness is low and linkage is loose, the baseline rate of substitution is [Formula: see text], where [Formula: see text] is the population size, [Formula: see text] is the rate of beneficial mutations per genome, and [Formula: see text] is their mean selective advantage. Heritable variance [Formula: see text] in log fitness due to unlinked loci reduces [Formula: see text] by [Formula: see text] under polygamy and [Formula: see text] under monogamy. With a linear genetic map of length [Formula: see text] Morgans, interference is yet stronger. We use a scaling argument to show that the density of adaptive substitutions depends on [Formula: see text], [Formula: see text], [Formula: see text], and [Formula: see text] only through the baseline density: [Formula: see text]. Under the approximation that the interference due to different sweeps adds up, we show that [Formula: see text], implying that interference prevents the rate of adaptive substitution from exceeding one per centimorgan per 200 generations. Simulations and numerical calculations confirm the scaling argument and confirm the additive approximation for [Formula: see text]; for higher [Formula: see text], the rate of adaptation grows above [Formula: see text], but only very slowly. We also consider the effect of sweeps on neutral diversity and show that, while even occasional sweeps can greatly reduce neutral diversity, this effect saturates as sweeps become more common-diversity can be maintained even in populations experiencing very strong interference. Our results indicate that for some organisms the rate of adaptive substitution may be primarily recombination-limited, depending only weakly on the mutation supply and the strength of selection.     

Suicide inactivation of tyrosinase in its action on tetrahydropterines

Tetrahydrobiopterin (BH(4)), methyl-tetrahydropterin (MBH(4)) and dimethyl-tetrahydropterin (DMBH(4)) are oxidized by tyrosinase in a process during which the suicide inactivation of tyrosinase may occur. From the kinetic study of this process, [Formula: see text] (apparent maximum constant for the suicide inactivation), [Formula: see text] (Michaelis constant for the substrate) and r (number of turnovers that the enzyme makes before the inactivation) can be obtained. From the results obtained, it can be deduced that the velocity of the inactivation governed by ([Formula: see text]) and the potency of the same ([Formula: see text]) follow the order: BH(4) > MBH(4) > DMBH(4).     

Fracture Faces in the Cell Envelope of Escherichia coli

Freeze-fracturing of Escherichia coli cells in the presence of 30% (v/v) glycerol resulted in a double cleavage of the cell envelope exposing two convex and two concave fracture faces ([Formula: see text], [Formula: see text] and [Formula: see text], [Formula: see text]) with characteristic patterns. Complementary replicas revealed the relationship of the fracture faces to their corresponding fracture planes. The inner fracture plane splits the plasma membrane at one particular level. Apparently the outer fracture plane was located in the outer part of the wall, as it was separated by a layer ([Formula: see text]) from the fractured profile (CW1) presumably corresponding to the murein layer. The outer fracture plane did alternate toward the cell periphery, exposing complementary smooth areas ([Formula: see text] and [Formula: see text]). When cells were freeze-fractured in the absence of glycerol, the outer cell surface appeared as an etching face rather than a fracture face. A schematic representation of the relative location of the different fracture faces in the E. coli cell envelope is given.     

Architecture of the Outer Membrane of Escherichia coli III. Protein-Lipopolysaccharide Complexes in Intramembraneous Particles

In a previous paper (A. Verkleij, L. van Alphen, J. Bijvelt, and B. Lugtenberg, Biochim. Biophys. Acta 466:269-282, 1977) we have hypothesized that particles on the outer fracture face of the outer membrane ([Formula: see text]), with corresponding pits on the inner fracture face of the outer membrane ([Formula: see text]), consist of lipopolysaccharide (LPS) aggregates stabilized by divalent cations and that they might contain protein and/or phospholipid. In the present paper the roles of LPS, cations, and proteins in these [Formula: see text] particles are described more extensively, using a strain that lacks the major outer membrane proteins, b, c, and d (b(-) c(-) d(-)), and has a reduction in the number of [Formula: see text] particles of 75%. To study the role of divalent cations in the formation of [Formula: see text] particles, these b(-) c(-) d(-) cells were grown or incubated with Ca(2+), Mg(2+), or putrescine. The presence of Ca(2+) resulted in the appearance of many [Formula: see text] particles and [Formula: see text] pits. Mg(2+) and putrescine were less effective than Ca(2+). Introduction of these particles was not accompanied by alterations in the relative amounts of LPS and cell envelope proteins. Ca(2+) treatment of a heptoseless derivative of a b(-) c(-) d(-) strain did not result in morphological changes. Incubation of Ca(2+)-treated cells with ethylenediaminetetraacetate caused the disappearance of the introduced particles as well as the release of more than 60% of the cellular LPS. These results strongly support the hypothesis that LPS is involved in the formation of [Formula: see text] particles and [Formula: see text] pits. The roles of various outer membrane proteins in the formation of [Formula: see text] particles were studied by comparing the freeze-fracture morphology of b(-) c(-) d(-) cells with that of cells which contain one of the outer membrane proteins b, c, d, and e or the receptor protein for bacteriophage lambda. The results showed that the presence of any of these five proteins in a b(-) c(-) d(-) background resulted in a large increase in the number of [Formula: see text] particles and [Formula: see text] pits, indicating that these proteins are, independent of each other, involved in the formation of [Formula: see text] particles and [Formula: see text] pits. The simplest explanation for the results is that in wild-type cells each particle consists of LPS complexed with some molecules of a single protein species, stabilized by either divalent cations or polyamines. It is hypothesized that the outer membrane of the wild-type cell contains a heterogeneous population of particles, of which 75% consists of protein b-LPS, protein c-LPS, and protein d-LPS particles. A function of these particles as aqueous pores is proposed.     

Mapping of redox state of mitochondrial cytochromes in live cardiomyocytes using Raman microspectroscopy

This paper presents a nonivasive approach to study redox state of reduced cytochromes [Formula: see text], [Formula: see text] and [Formula: see text] of complexes II and III in mitochondria of live cardiomyocytes by means of Raman microspectroscopy. For the first time with the proposed approach we perform studies of rod- and round-shaped cardiomyocytes, representing different morphological and functional states. Raman mapping and cluster analysis reveal that these cardiomyocytes differ in the amounts of reduced cytochromes [Formula: see text], [Formula: see text] and [Formula: see text]. The rod-shaped cardiomyocytes possess uneven distribution of reduced cytochromes [Formula: see text], [Formula: see text] and [Formula: see text] in cell center and periphery. Moreover, by means of Raman spectroscopy we demonstrated the decrease in the relative amounts of reduced cytochromes [Formula: see text], [Formula: see text] and [Formula: see text] in the rod-shaped cardiomyocytes caused by H(2)O(2)-induced oxidative stress before any visible changes. Results of Raman mapping and time-dependent study of reduced cytochromes of complexes II and III and cytochrome [Formula: see text] in cardiomyocytes are in a good agreement with our fluorescence indicator studies and other published data.     

A Mathematical Model of Rift Valley Fever with Human Host

Rift Valley Fever is a vector-borne disease mainly transmitted by mosquito. To gain some quantitative insights into its dynamics, a deterministic model with mosquito, livestock, and human host is formulated as a system of nonlinear ordinary differential equations and analyzed. The disease threshold [Formula: see text] is computed and used to investigate the local stability of the equilibria. A sensitivity analysis is performed and the most sensitive model parameters to the measure of initial disease transmission [Formula: see text] and the endemic equilibrium are determined. Both [Formula: see text] and the disease prevalence in mosquitoes are more sensitive to the natural mosquito death rate, d(m). The disease prevalence in livestock and humans are more sensitive to livestock and human recruitment rates, [Formula: see text] and [Formula: see text], respectively, suggesting isolation of livestock from humans is a viable preventive strategy during an outbreak. Numerical simulations support the analytical results in further exploring theoretically the long-term dynamics of the disease at the population level.     

Theoretical analysis of simultaneous mixed irradiation with multiple types of radiation.

A survival model Eq.1 was presented for cells irradiated simultaneously with multiple types of radiation using the extended Zaider-Rossi model, which is model for mixed irradiation with two types of radiation. [equation : see text] Eq.1. Where q(t)=2t0/t-2(t0/t)2 ?1-exp(-t0/t)? Eq.2. Eq.1 was proved by mathematical induction using the concept that mixed irradiation with n types of radiation is considered as mixed irradiation with two types of radiation regarding n-1 types as one type of radiation. The model is not limited by the dose rate of radiation, because its effect is corrected by reduction factor Eq.2. The problem of the model is that Eq.2 was led assuming repair function to be exponential given by Eq.3. tau(t)=exp(-t/t0) Eq.3. However, the repair function is usually expressed by biphasic Eq.4 rather than monophasic Eq.3. tau(t)=Aexp(-t/01)+(1-A)exp(-t/t02) Eq.4. It is, therefore, important to keep in mind that Eq. 4 should be used instead of Eq.3.     

Uptake of Different Species of Iodine by Water Spinach and Its Effect to Growth

A hydroponic experiment has been carried out to study the influence of iodine species [iodide (I(-)), iodate ([Formula: see text]), and iodoacetic acid (CH(2)ICOO(-))] and concentrations on iodine uptake by water spinach. Results show that low levels of iodine in the nutrient solution can effectively stimulate the growth of biomass of water spinach. When iodine levels in the nutrient solution are from 0 to 1.0 mg/l, increases in iodine levels can linearly augment iodine uptake rate by the leafy vegetables from all three species of iodine, and the uptake effects are in the following order: [Formula: see text]. In addition, linear correlation was observed between iodine content in the roots and shoots of water spinach, and their proportion is 1:1. By uptake of I(-), vitamin C (Vit C) content in water spinach increased, whereas uptake of [Formula: see text] and CH(2)ICOO(-) decreased water spinach Vit C content. Furthermore, through uptake of I(-) and [Formula: see text], the nitrate content in water spinach was increased by different degrees.     

Complete determination of the Pin1 catalytic domain thermodynamic cycle by NMR lineshape analysis

The phosphorylation-specific peptidyl-prolyl isomerase Pin1 catalyzes the isomerization of the peptide bond preceding a proline residue between cis and trans isomers. To best understand the mechanisms of Pin1 regulation, rigorous enzymatic assays of isomerization are required. However, most measures of isomerase activity require significant constraints on substrate sequence and only yield rate constants for the cis isomer, [Formula: see text] and apparent Michaelis constants, [Formula: see text]. By contrast, NMR lineshape analysis is a powerful tool for determining microscopic rates and populations of each state in a complex binding scheme. The isolated catalytic domain of Pin1 was employed as a first step towards elucidating the reaction scheme of the full-length enzyme. A 24-residue phosphopeptide derived from the amyloid precurser protein intracellular domain (AICD) phosphorylated at Thr668 served as a biologically-relevant Pin1 substrate. Specific (13)C labeling at the Pin1-targeted proline residue provided multiple reporters sensitive to individual isomer binding and on-enzyme catalysis. We have performed titration experiments and employed lineshape analysis of phosphopeptide (13)C-(1)H constant time HSQC spectra to determine [Formula: see text], [Formula: see text], [Formula: see text], and [Formula: see text] for the catalytic domain of Pin1 acting on this AICD substrate. The on-enzyme equilibrium value of [E·trans]/[E·cis]?=?3.9 suggests that the catalytic domain of Pin1 is optimized to operate on this substrate near equilibrium in the cellular context. This highlights the power of lineshape analysis for determining the microscopic parameters of enzyme catalysis, and demonstrates the feasibility of future studies of Pin1-PPIase mutants to gain insights on the catalytic mechanism of this important enzyme.     

Robustness and information propagation in attractors of random boolean networks

Attractors represent the long-term behaviors of Random Boolean Networks. We study how the amount of information propagated between the nodes when on an attractor, as quantified by the average pairwise mutual information ([Formula: see text]), relates to the robustness of the attractor to perturbations ([Formula: see text]). We find that the dynamical regime of the network affects the relationship between [Formula: see text] and [Formula: see text]. In the ordered and chaotic regimes, [Formula: see text] is anti-correlated with [Formula: see text], implying that attractors that are highly robust to perturbations have necessarily limited information propagation. Between order and chaos (for so-called "critical" networks) these quantities are uncorrelated. Finite size effects cause this behavior to be visible for a range of networks, from having a sensitivity of 1 to the point where [Formula: see text] is maximized. In this region, the two quantities are weakly correlated and attractors can be almost arbitrarily robust to perturbations without restricting the propagation of information in the network.     

Critical motor number for fractional steps of cytoskeletal filaments in gliding assays

In gliding assays, filaments are pulled by molecular motors that are immobilized on a solid surface. By varying the motor density on the surface, one can control the number [Formula: see text] of motors that pull simultaneously on a single filament. Here, such gliding assays are studied theoretically using Brownian (or Langevin) dynamics simulations and taking the local force balance between motors and filaments as well as the force-dependent velocity of the motors into account. We focus on the filament stepping dynamics and investigate how single motor properties such as stalk elasticity and step size determine the presence or absence of fractional steps of the filaments. We show that each gliding assay can be characterized by a critical motor number, [Formula: see text]. Because of thermal fluctuations, fractional filament steps are only detectable as long as [Formula: see text]. The corresponding fractional filament step size is [Formula: see text] where [Formula: see text] is the step size of a single motor. We first apply our computational approach to microtubules pulled by kinesin-1 motors. For elastic motor stalks that behave as linear springs with a zero rest length, the critical motor number is found to be [Formula: see text], and the corresponding distributions of the filament step sizes are in good agreement with the available experimental data. In general, the critical motor number [Formula: see text] depends on the elastic stalk properties and is reduced to [Formula: see text] for linear springs with a nonzero rest length. Furthermore, [Formula: see text] is shown to depend quadratically on the motor step size [Formula: see text]. Therefore, gliding assays consisting of actin filaments and myosin-V are predicted to exhibit fractional filament steps up to motor number [Formula: see text]. Finally, we show that fractional filament steps are also detectable for a fixed average motor number [Formula: see text] as determined by the surface density (or coverage) of the motors on the substrate surface.     

Hydroxylation of p-substituted phenols by tyrosinase: further insight into the mechanism of tyrosinase activity

A study of the monophenolase activity of tyrosinase by measuring the steady state rate with a group of p-substituted monophenols provides the following kinetic information: k(cat)(m) and the Michaelis constant, K(M)(m). Analysis of these data taking into account chemical shifts of the carbon atom supporting the hydroxyl group (δ) and σ(p)(+), enables a mechanism to be proposed for the transformation of monophenols into o-diphenols, in which the first step is a nucleophilic attack on the copper atom on the form E(ox) (attack of the oxygen of the hydroxyl group of C-1 on the copper atom) followed by an electrophilic attack (attack of the hydroperoxide group on the ortho position with respect to the hydroxyl group of the benzene ring, electrophilic aromatic substitution with a reaction constant ρ of -1.75). These steps show the same dependency on the electronic effect of the substituent groups in C-4. Furthermore, a study of a solvent deuterium isotope effect on the oxidation of monophenols by tyrosinase points to an appreciable isotopic effect. In a proton inventory study with a series of p-substituted phenols, the representation of [Formula: see text] / [Formula: see text] against n (atom fractions of deuterium), where [Formula: see text] is the catalytic constant for a molar fraction of deuterium (n) and [Formula: see text] is the corresponding kinetic parameter in a water solution, was linear for all substrates. These results indicate that only one of the proton transfer processes from the hydroxyl groups involved the catalytic cycle is responsible for the isotope effects. We suggest that this step is the proton transfer from the hydroxyl group of C-1 to the peroxide of the oxytyrosinase form (E(ox)). After the nucleophilic attack, the incorporation of the oxygen in the benzene ring occurs by means of an electrophilic aromatic substitution mechanism in which there is no isotopic effect.     

Fully relativistic rovibrational energies and spectroscopic constants of the lowest {\text{X}}:(1)0_g^{ + } , A′:(1)2 u , A:(1)1 u , {\text{B}}\prime :(1)0_u^{ - } and {\text{B}}:(1)0_u^{ + } states of molecular chlorine

The main goal of this paper is to present the rovibrational energies and spectroscopic constants of the Cl(2) molecular system in the relativistic states [Formula: see text], A':(1)2( u ), A:(1)1( u ), [Formula: see text] and [Formula: see text]. More precisely, we have evaluated the Cl(2) ω ( e ), ω ( e ) x ( e ), ω ( e ) y ( e ), α ( e ), γ ( e ) and B ( e ) rovibrational spectroscopic constants using two different procedures. The first was obtained by combining the rovibrational energies, calculated through solving Schr?dinger's nuclear equation and the diatomic rovibrational energy equation. The second was obtained by using the Dunham method. The calculated properties are in good agreement with available experimental data.     

Specificity of Hpa II and Hae III DNA methylases

The methylases M.HaeIII and M.HpaII recognize the tetranucleotide sequences [Formula: see text] and [Formula: see text] respectively, in DNA, and transfer a methyl group from S-adenosylmethionine to the 5-position of cytosine on each strand as indicated by the asterisks. Restriction endonuclease R.HaeIII does not cleave the methylated sequence [Formula: see text] but can cleave [Formula: see text] in which methylation is introduced on the unnatural external cytosine positions. Similarly, R.HpaII does not cleave [Formula: see text] but can cleave [Formula: see text].Images     

Is thiocyanate peroxidation at equilibrium in vivo?

The peroxidase-catalyzed oxidation of SCN- by H2O2 is an important in vivo reaction because it limits the accumulation of toxic H2O2 and provides significant concentrations of the antimicrobial agents, HOSCN and OSCN-. Data presented in this report suggest that the reaction: (Formula: see text) is in a state of dynamic equilibrium in vivo. Since OSCN- can form the weak acid HOSCN (pKa = 5.3), the equilibrium constant expression (Kox) for thiocyanate peroxidation is dependent on the concentration of hydrogen ions as well as the concentrations of H2O2, SCN-, HOSCN, OSCN- and water, and on the HOSCN ionization constant, Ka: (Formula: see text). The concentration of water is assumed to be constant and unaffected by the other components and is omitted from the Kox equation. The value of Kox was estimated from in vitro data to be 3.7 X 10(3) M-1 (S.D. = 0.8 X 10(3) M-1, n = 8). Using this value for Kox and observations of salivary concentrations of SCN- and HOSCN + OSCN- from several previous reports, the equilibrium concentrations of H2O2 in whole saliva were calculated to range from 8 to 13 microM. This range is consistent with reported estimates of 10 microM as the hydrogen peroxide tolerance limit for human cells.     

Complexity-entropy causality plane as a complexity measure for two-dimensional patterns

Complexity measures are essential to understand complex systems and there are numerous definitions to analyze one-dimensional data. However, extensions of these approaches to two or higher-dimensional data, such as images, are much less common. Here, we reduce this gap by applying the ideas of the permutation entropy combined with a relative entropic index. We build up a numerical procedure that can be easily implemented to evaluate the complexity of two or higher-dimensional patterns. We work out this method in different scenarios where numerical experiments and empirical data were taken into account. Specifically, we have applied the method to [Formula: see text] fractal landscapes generated numerically where we compare our measures with the Hurst exponent; [Formula: see text] liquid crystal textures where nematic-isotropic-nematic phase transitions were properly identified; [Formula: see text] 12 characteristic textures of liquid crystals where the different values show that the method can distinguish different phases; [Formula: see text] and Ising surfaces where our method identified the critical temperature and also proved to be stable.     

A flexible sigmoid function of determinate growth

A new empirical equation for the sigmoid pattern of determinate growth, 'the beta growth function', is presented. It calculates weight (w) in dependence of time, using the following three parameters: t(m), the time at which the maximum growth rate is obtained; t(e), the time at the end of growth; and w(max), the maximal value for w, which is achieved at t(e). The beta growth function was compared with four classical (logistic, Richards, Gompertz and Weibull) growth equations, and two expolinear equations. All equations described successfully the sigmoid dynamics of seed filling, plant growth and crop biomass production. However, differences were found in estimating w(max). Features of the beta function are: (1) like the Richards equation it is flexible in describing various asymmetrical sigmoid patterns (its symmetrical form is a cubic polynomial); (2) like the logistic and the Gompertz equations its parameters are numerically stable in statistical estimation; (3) like the Weibull function it predicts zero mass at time zero, but its extension to deal with various initial conditions can be easily obtained; (4) relative to the truncated expolinear equation it provides more reasonable estimates of final quantity and duration of a growth process. In addition, the new function predicts a zero growth rate at both the start and end of a precisely defined growth period. Therefore, it is unique for dealing with determinate growth, and is more suitable than other functions for embedding in process-based crop simulation models to describe the dynamics of organs as sinks to absorb assimilates. Because its parameters correspond to growth traits of interest to crop scientists, the beta growth function is suitable for characterization of environmental and genotypic influences on growth processes. However, it is not suitable for estimating maximum relative growth rate to characterize early growth that is expected to be close to exponential.     

Rat alveolar macrophages require NADPH for superoxide production in the respiratory burst. Effect of NADPH depletion by paraquat

Alveolar macrophages can be stimulated by concanavalin A to produce extracellular superoxide. Conflicting opinions exist, however, concerning the relative importance of the oxidation of either NADPH or NADH in the generation of (Formula: see text) by surface membrane-stimulated phagocytic cells. Alveolar macrophages were obtained from adult male rats by lavage with phosphate-buffered saline. Cells (approximately 10(6)/ml) were incubated in Krebs-Ringer phosphate 4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid buffer and ferricytochrome c for 15 min at 37 degrees C before addition of concanavalin A. Release of (Formula: see text) was detected as the difference in cytochrome c reduction, followed at 550 nm, in the absence and presence of superoxide dismutase. Superoxide production by concanavalin A-stimulated alveolar macrophages was markedly increased in the presence of glucose but fructose, lactate, and pyruvate were without effect. Paraquat (methylviologen), an oxidation-reduction dye, significantly reduced concanavalin A-stimulated (Formula: see text) production when incubated at 1 mM with alveolar macrophages in the absence of glucose. The effect of paraquat was reversed by glucose, but fructose, lactate, and pyruvate could not reverse paraquat inhibition. Paraquat enhanced oxidation of NADPH (but not NADH) by cell supernatant and increased pentose phosphate shunt activity in resting macrophages, but did not affect mitochondrial respiration or ATP content of alveolar macrophages. These results suggest that paraquat is able to specifically deplete NADPH in alveolar macrophages while not affecting NADH or ATP. Our conclusion is that NADPH is essential for the production of (Formula: see text) by concanavalin A-stimulated alveolar macrophages.