A fast, high resolution, second-order central scheme for incompressible flows

A high resolution, second-order central difference method for incompressible flows is presented. The method is based on a recent second-order extension of the classic Lax–Friedrichs scheme introduced for hyperbolic conservation laws (Nessyahu H. & Tadmor E. (1990) J. Comp. Physics. 87, 408-463; Jiang G.-S. & Tadmor E. (1996) UCLA CAM Report 96-36, SIAM J. Sci. Comput., in press) and augmented by a new discrete Hodge projection. The projection is exact, yet the discrete Laplacian operator retains a compact stencil. The scheme is fast, easy to implement, and readily generalizable. Its performance was tested on the standard periodic double shear-layer problem; no spurious vorticity patterns appear when the flow is underresolved. A short discussion of numerical boundary conditions is also given, along with a numerical example.

DNA strand annealing is promoted by the yeast Rad52 protein.

The Saccharomyces cerevisiae RAD52 gene plays a pivotal role in genetic recombination. Here we demonstrate that yeast Rad52 is a DNA binding protein. To show that the interaction between Rad52 and DNA is direct and not mediated by other yeast proteins and to facilitate protein purification, a recombinant expression system was developed. The recombinant protein can bind both single- and double-stranded DNA and the addition of either Mg2+ or ATP does not enhance the binding of single-stranded DNA. Furthermore, a DNA binding domain was found in the evolutionary conserved N terminus of the protein. More importantly, we show that the protein stimulates DNA annealing even in the presence of a large excess of nonhomologous DNA. Rad52-promoted annealing follows second-order kinetics and the rate is 3500-fold faster than that of the spontaneous reaction. How this annealing activity relates to the genetic phenotype associated with rad52 mutant cells is discussed.

Modifications of retinal neurons in a mouse model of retinitis pigmentosa

Animal models of retinitis pigmentosa include the rd mouse, in which a mutation of a rod-specific phosphodiesterase leads to the rapid loss of photoreceptors during the early postnatal life. Very little is known about changes occurring in inner retinal neurons after photoreceptor loss. These changes are important in view of the possibility of restoring vision in retinas with photoreceptor degeneration by means of cell transplantation or direct stimulation of inner layers. In this paper, we show that bipolar and horizontal cells of the rd mouse retina undergo dramatic morphological modifications accompanying photoreceptor loss, demonstrating a dependence of second order neurons on these cells. While describing modifications of the rd retina, we also provide quantitative information about neurons of the wild-type mouse retina, useful for future studies on genetically altered animals.

Kinetics of photo-induced electron transfer from high-potential iron-sulfur protein to the photosynthetic reaction center of the purple phototroph Rhodoferax fermentans.

The kinetics of photo-induced electrontransfer from high-potential iron-sulfur protein (HiPIP) to the photosynthetic reaction center (RC) of the purple phototroph Rhodoferarfermentans were studied. The rapid photooxidation of heme c-556 belonging to RC is followed, in the presence of HiPIP, by a slower reduction having a second-order rate constant of 4.8 x 10(7) M(-1) x s(-1). The limiting value of kobs at high HiPIP concentration is 95 s(-1). The amplitude of this slow process decreases with increasing HiPIP concentration. The amplitude of a faster phase, observed at 556 and 425 nm and involving heme c-556 reduction, increases proportionately. The rate constant of this fast phase, determined at 425 and 556 nm, is approximately 3 x 10(5) s(-1). This value is not dependent on HiPIP concentration, indicating that it is related to a first-order process. These observations are interpreted as evidence for the formation of a HiPIP-RC complex prior to the excitation flash, having a dissociation constant of -2.5 microM. The fast phase is absent at high ionic strength, indicating that the complex involves mainly electrostatic interactions. The ionic strength dependence of kobs for the slow phase yields a second-order rate constant at infinite ionic strength of 5.4 x 10(6) M(-1) x s(-1) and an electrostatic interaction energy of -2.1 kcal/mol (1 cal = 4.184 J). We conclude that Rhodoferar fermentans HiPIP is a very effective electron donor to the photosynthetic RC.

Blind separation of auditory event-related brain responses into independent components

Averaged event-related potential (ERP) data recorded from the human scalp reveal electroencephalographic (EEG) activity that is reliably time-locked and phase-locked to experimental events. We report here the application of a method based on information theory that decomposes one or more ERPs recorded at multiple scalp sensors into a sum of components with fixed scalp distributions and sparsely activated, maximally independent time courses. Independent component analysis (ICA) decomposes ERP data into a number of components equal to the number of sensors. The derived components have distinct but not necessarily orthogonal scalp projections. Unlike dipole-fitting methods, the algorithm does not model the locations of their generators in the head. Unlike methods that remove second-order correlations, such as principal component analysis (PCA), ICA also minimizes higher-order dependencies. Applied to detected—and undetected—target ERPs from an auditory vigilance experiment, the algorithm derived ten components that decomposed each of the major response peaks into one or more ICA components with relatively simple scalp distributions. Three of these components were active only when the subject detected the targets, three other components only when the target went undetected, and one in both cases. Three additional components accounted for the steady-state brain response to a 39-Hz background click train. Major features of the decomposition proved robust across sessions and changes in sensor number and placement. This method of ERP analysis can be used to compare responses from multiple stimuli, task conditions, and subject states.

A transient expansion of the native state precedes aggregation of recombinant human interferon-γ

Aggregation of proteins, even under conditions favoring the native state, is a ubiquitous problem in biotechnology and biomedical engineering. Providing a mechanistic basis for the pathways that lead to aggregation should allow development of rational approaches for its prevention. We have chosen recombinant human interferon-γ (rhIFN-γ) as a model protein for a mechanistic study of aggregation. In the presence of 0.9 M guanidinium hydrochloride, rhIFN-γ aggregates with first order kinetics, a process that is inhibited by addition of sucrose. We describe a pathway that accounts for both the observed first-order aggregation of rhIFN-γ and the effect of sucrose. In this pathway, aggregation proceeds through a transient expansion of the native state. Sucrose shifts the equilibrium within the ensemble of rhIFN-γ native conformations to favor the most compact native species over more expanded ones, thus stabilizing rhIFN-γ against aggregation. This phenomenon is attributed to the preferential exclusion of sucrose from the protein surface. In addition, kinetic analysis combined with solution thermodynamics shows that only a small (9%) expansion surface area is needed to form the transient native state that precedes aggregation. The approaches used here link thermodynamics and aggregation kinetics to provide a powerful tool for understanding both the pathway of protein aggregation and the rational use of excipients to inhibit the process.

The role of London forces in defining noncentrosymmetric order of high dipole moment–high hyperpolarizability chromophores in electrically poled polymeric thin films

Graphs of second harmonic generation coefficients and electro-optic coefficients (measured by ellipsometry, attenuated total reflection, and two-slit interference modulation) as a function of chromophore number density (chromophore loading) are experimentally observed to exhibit maxima for polymers containing chromophores characterized by large dipole moments and polarizabilities. Modified London theory is used to demonstrated that this behavior can be attributed to the competition of chromophore-applied electric field and chromophore–chromophore electrostatic interactions. The comparison of theoretical and experimental data explains why the promise of exceptional macroscopic second-order optical nonlinearity predicted for organic materials has not been realized and suggests routes for circumventing current limitations to large optical nonlinearity. The results also suggest extensions of measurement and theoretical methods to achieve an improved understanding of intermolecular interactions in condensed phase materials including materials prepared by sequential synthesis and block copolymer methods.

Linear and nonlinear pathways of spectral information transmission in the cochlear nucleus

At the level of the cochlear nucleus (CN), the auditory pathway divides into several parallel circuits, each of which provides a different representation of the acoustic signal. Here, the representation of the power spectrum of an acoustic signal is analyzed for two CN principal cells—chopper neurons of the ventral CN and type IV neurons of the dorsal CN. The analysis is based on a weighting function model that relates the discharge rate of a neuron to first- and second-order transformations of the power spectrum. In chopper neurons, the transformation of spectral level into rate is a linear (i.e., first-order) or nearly linear function. This transformation is a predominantly excitatory process involving multiple frequency components, centered in a narrow frequency range about best frequency, that usually are processed independently of each other. In contrast, type IV neurons encode spectral information linearly only near threshold. At higher stimulus levels, these neurons are strongly inhibited by spectral notches, a behavior that cannot be explained by level transformations of first- or second-order. Type IV weighting functions reveal complex excitatory and inhibitory interactions that involve frequency components spanning a wider range than that seen in choppers. These findings suggest that chopper and type IV neurons form parallel pathways of spectral information transmission that are governed by two different mechanisms. Although choppers use a predominantly linear mechanism to transmit tonotopic representations of spectra, type IV neurons use highly nonlinear processes to signal the presence of wide-band spectral features.

Kinetics and mechanism of polyamide ("peptide") nucleic acid binding to duplex DNA.

To elucidate the mechanism of recognition of double-stranded DNA (dsDNA) by homopyrimidine polyamide ("peptide") nucleic acid (PNA) leading to the strand-displacement, the kinetics of the sequence-specific PNA/DNA binding have been studied. The binding was monitored with time by the gel retardation and nuclease S1 cleavage assays. The experimental kinetic curves obey pseudo-first-order kinetics and the dependence of the pseudo-first-order rate constant, kps, on PNA concentration, P, obeys a power law kps approximately P gamma with 2 < gamma < 3. The kps values for binding of decamer PNA to dsDNA target sites with one mismatch are hundreds of times slower than for the correct site. A detailed kinetic scheme for PNA/DNA binding is proposed that includes two major steps of the reaction of strand invasion: (i) a transient partial opening of the PNA binding site on dsDNA and incorporation of one PNA molecule with the formation of an intermediate PNA/DNA duplex and (ii) formation of a very stable PNA2/DNA triplex. A simple theoretical treatment of the proposed kinetic scheme is performed. The interpretation of our experimental data in the framework of the proposed kinetic scheme leads to the following conclusions. The sequence specificity of the recognition is essentially provided at the "search" step of the process, which consists in the highly reversible transient formation of duplex between one PNA molecule and the complementary strand of duplex DNA while the other DNA strand is displaced. This search step is followed by virtually irreversible "locking" step via PNA2/DNA triplex formation. The proposed mechanism explains how the binding of homopyrimidine PNA to dsDNA meets two apparently mutually contradictory features: high sequence specificity of binding and remarkable stability of both correct and mismatched PNA/DNA complexes.