“Coarse” stability and bifurcation analysis using time-steppers: A reaction-diffusion example

Evolutionary, pattern forming partial differential equations (PDEs) are often derived as limiting descriptions of microscopic, kinetic theory-based models of molecular processes (e.g., reaction and diffusion). The PDE dynamic behavior can be probed through direct simulation (time integration) or, more systematically, through stability/bifurcation calculations; time-stepper-based approaches, like the Recursive Projection Method [Shroff, G. M. & Keller, H. B. (1993) SIAM J. Numer. Anal. 30, 1099–1120] provide an attractive framework for the latter. We demonstrate an adaptation of this approach that allows for a direct, effective (“coarse”) bifurcation analysis of microscopic, kinetic-based models; this is illustrated through a comparative study of the FitzHugh-Nagumo PDE and of a corresponding Lattice–Boltzmann model.

Ecological chaos in the wake of invasion.

Irregularities in observed population densities have traditionally been attributed to discretization of the underlying dynamics. We propose an alternative explanation by demonstrating the evolution of spatiotemporal chaos in reaction-diffusion models for predator-prey interactions. The chaos is generated naturally in the wake of invasive waves of predators. We discuss in detail the mechanism by which the chaos is generated. By considering a mathematical caricature of the predator-prey models, we go on to explain the dynamical origin of the irregular behavior and to justify our assertion that the behavior we present is a genuine example of spatiotemporal chaos.

Twist grain boundaries in three-dimensional lamellar Turing structures

Steady spatial self-organization of three-dimensional chemical reaction-diffusion systems is discussed with the emphasis put on the possible defects that may alter the Turing patterns. It is shown that one of the stable defects of a three-dimensional lamellar Turing structure is a twist grain boundary embedding a Scherk minimal surface.

Microtubule self-organization is gravity-dependent

Although weightlessness is known to affect living cells, the manner by which this occurs is unknown. Some reaction-diffusion processes have been theoretically predicted as being gravity-dependent. Microtubules, a major constituent of the cellular cytoskeleton, self-organize in vitro by way of reaction-diffusion processes. To investigate how self-organization depends on gravity, microtubules were assembled under low gravity conditions produced during space flight. Contrary to the samples formed on an in-flight 1 × g centrifuge, the samples prepared in microgravity showed almost no self-organization and were locally disordered.

Proton transfer from the bulk to the bound ubiquinone QB of the reaction center in chromatophores of Rhodobacter sphaeroides: Retarded conveyance by neutral water

The mechanism of proton transfer from the bulk into the membrane protein interior was studied. The light-induced reduction of a bound ubiquinone molecule QB by the photosynthetic reaction center is accompanied by proton trapping. We used kinetic spectroscopy to measure (i) the electron transfer to QB (at 450 nm), (ii) the electrogenic proton delivery from the surface to the QB site (by electrochromic carotenoid response at 524 nm), and (iii) the disappearance of protons from the bulk solution (by pH indicators). The electron transfer to QB− and the proton-related electrogenesis proceeded with the same time constant of ≈100 μs (at pH 6.2), whereas the alkalinization in the bulk was distinctly delayed (τ ≈ 400 μs). We investigated the latter reaction as a function of the pH indicator concentration, the added pH buffers, and the temperature. The results led us to the following conclusions: (i) proton transfer from the surface-located acidic groups into the QB site followed the reduction of QB without measurable delay; (ii) the reprotonation of these surface groups by pH indicators and hydronium ions was impeded, supposedly, because of their slow diffusion in the surface water layer; and (iii) as a result, the protons were slowly donated by neutral water to refill the proton vacancies at the surface. It is conceivable that the same mechanism accounts for the delayed relaxation of the surface pH changes into the bulk observed previously with bacteriorhodopsin membranes and thylakoids. Concerning the coupling between proton pumps in bioenergetic membranes, our results imply a tendency for the transient confinement of protons at the membrane surface.

Fluorescence correlation analysis of probe diffusion simplifies quantitative pathogen detection by PCR

A sensitive, labor-saving, and easily automatable nonradioactive procedure named APEX-FCS (amplified probe extension detected by fluorescence correlation spectroscopy) has been established to detect specific in vitro amplification of pathogen genomic sequences. As an example, Mycobacterium tuberculosis genomic DNA was subjected to PCR amplification with the Stoffel fragment of Thermus aquaticus DNA polymerase in the presence of nanomolar concentrations of a rhodamine-labeled probe (third primer), binding to the target in between the micromolar amplification primers. The probe becomes extended only when specific amplification occurs. Its low concentration avoids false-positives due to unspecific hybridization under PCR conditions. With increasing portion of extended probe molecules, the probe’s average translational diffusion properties gradually change over the course of the reaction, reflecting amplification kinetics. Following PCR, this change from a stage of high to a stage of low mobility can directly be monitored during a 30-s measurement using a fluorescence correlation spectroscopy device. Quantitation down to 10 target molecules in a background of 2.5 μg unspecific DNA without post-PCR probe manipulations could be achieved with different primer/probe combinations. The assay holds the promise to concurrently perform amplification, probe hybridization, and specific detection without opening the reaction chamber, if sealable foils are used.

Two amyloid precursor protein transgenic mouse models with Alzheimer disease-like pathology

Mutations in the amyloid precursor protein (APP) gene cause early-onset familial Alzheimer disease (AD) by affecting the formation of the amyloid β (Aβ) peptide, the major constituent of AD plaques. We expressed human APP751 containing these mutations in the brains of transgenic mice. Two transgenic mouse lines develop pathological features reminiscent of AD. The degree of pathology depends on expression levels and specific mutations. A 2-fold overexpression of human APP with the Swedish double mutation at positions 670/671 combined with the V717I mutation causes Aβ deposition in neocortex and hippocampus of 18-month-old transgenic mice. The deposits are mostly of the diffuse type; however, some congophilic plaques can be detected. In mice with 7-fold overexpression of human APP harboring the Swedish mutation alone, typical plaques appear at 6 months, which increase with age and are Congo Red-positive at first detection. These congophilic plaques are accompanied by neuritic changes and dystrophic cholinergic fibers. Furthermore, inflammatory processes indicated by a massive glial reaction are apparent. Most notably, plaques are immunoreactive for hyperphosphorylated tau, reminiscent of early tau pathology. The immunoreactivity is exclusively found in congophilic senile plaques of both lines. In the higher expressing line, elevated tau phosphorylation can be demonstrated biochemically in 6-month-old animals and increases with age. These mice resemble major features of AD pathology and suggest a central role of Aβ in the pathogenesis of the disease.

Limited internal friction in the rate-limiting step of a two-state protein folding reaction

Small, single-domain proteins typically fold via a compact transition-state ensemble in a process well fitted by a simple, two-state model. To characterize the rate-limiting conformational changes that underlie two-state folding, we have investigated experimentally the effects of changing solvent viscosity on the refolding of the IgG binding domain of protein L. In conjunction with numerical simulations, our results indicate that the rate-limiting conformational changes of the folding of this domain are strongly coupled to solvent viscosity and lack any significant “internal friction” arising from intrachain collisions. When compared with the previously determined solvent viscosity dependencies of other, more restricted conformational changes, our results suggest that the rate-limiting folding transition involves conformational fluctuations that displace considerable amounts of solvent. Reconciling evidence that the folding transition state ensemble is comprised of highly collapsed species with these and similar, previously reported results should provide a significant constraint for theoretical models of the folding process.

Single-molecule protein folding: Diffusion fluorescence resonance energy transfer studies of the denaturation of chymotrypsin inhibitor 2

We report single-molecule folding studies of a small, single-domain protein, chymotrypsin inhibitor 2 (CI2). CI2 is an excellent model system for protein folding studies and has been extensively studied, both experimentally (at the ensemble level) and theoretically. Conformationally assisted ligation methodology was used to synthesize the proteins and site-specifically label them with donor and acceptor dyes. Folded and denatured subpopulations were observed by fluorescence resonance energy transfer (FRET) measurements on freely diffusing single protein molecules. Properties of these subpopulations were directly monitored as a function of guanidinium chloride concentration. It is shown that new information about different aspects of the protein folding reaction can be extracted from such subpopulation properties. Shifts in the mean transfer efficiencies are discussed, FRET efficiency distributions are translated into potentials, and denaturation curves are directly plotted from the areas of the FRET peaks. Changes in stability caused by mutation also are measured by comparing pseudo wild-type CI2 with a destabilized mutant (K17G). Current limitations and future possibilities and prospects for single-pair FRET protein folding investigations are discussed.

Stretching lattice models of protein folding

A new class of experiments that probe folding of individual protein domains uses mechanical stretching to cause the transition. We show how stretching forces can be incorporated in lattice models of folding. For fast folding proteins, the analysis suggests a complex relation between the force dependence and the reaction coordinate for folding.

Turbulent diffusion phase transition is due to singular energy spectrum.

The phase transition for turbulent diffusion, reported by Avellaneda and Majda [Avellaneda, M. & Majda, A. J. (1994) Philos. Trans. R. Soc. London A 346, 205-233, and several earlier papers], is traced to a modeling assumption in which the energy spectrum of the turbulent fluid is singularly dependent on the viscosity in the inertial range. Phenomenological models of turbulence and intermittency, by contrast, require that the energy spectrum be independent of the viscosity in the inertial range. When the energy spectrum is assumed to be consistent with the phenomenological models, there is no phase transition for turbulent diffusion.