Generalized flows, intrinsic stochasticity, and turbulent transport

The study of passive scalar transport in a turbulent velocity field leads naturally to the notion of generalized flows, which are families of probability distributions on the space of solutions to the associated ordinary differential equations which no longer satisfy the uniqueness theorem for ordinary differential equations. Two most natural regularizations of this problem, namely the regularization via adding small molecular diffusion and the regularization via smoothing out the velocity field, are considered. White-in-time random velocity fields are used as an example to examine the variety of phenomena that take place when the velocity field is not spatially regular. Three different regimes, characterized by their degrees of compressibility, are isolated in the parameter space. In the regime of intermediate compressibility, the two different regularizations give rise to two different scaling behaviors for the structure functions of the passive scalar. Physically, this means that the scaling depends on Prandtl number. In the other two regimes, the two different regularizations give rise to the same generalized flows even though the sense of convergence can be very different. The “one force, one solution” principle is established for the scalar field in the weakly compressible regime, and for the difference of the scalar in the strongly compressible regime, which is the regime of inverse cascade. Existence and uniqueness of an invariant measure are also proved in these regimes when the transport equation is suitably forced. Finally incomplete self similarity in the sense of Barenblatt and Chorin is established.

Single-molecule spectroscopy of the β2 adrenergic receptor: Observation of conformational substates in a membrane protein

Single-molecule studies of the conformations of the intact β2 adrenergic receptor were performed in solution. Photon bursts from the fluorescently tagged adrenergic receptor in a micelle were recorded. A photon-burst algorithm and a Poisson time filter were implemented to characterize single molecules diffusing across the probe volume of a confocal microscope. The effects of molecular diffusion and photon number fluctuations were deconvoluted by assuming that Poisson distributions characterize the molecular occupation and photon numbers. Photon-burst size histograms were constructed, from which the source intensity distributions were extracted. Different conformations of the β2 adrenergic receptor cause quenching of the bound fluorophore to different extents and hence produce different photon-burst sizes. An analysis of the photon-burst histograms shows that there are at least two distinct substates for the native adrenergic membrane receptor. This behavior is in contrast to one peak observed for the dye molecule, rhodamine 6G. We test the reliability and robustness of the substate number determination by investigating the application of different binning criteria. Conformational changes associated with agonist binding result in a marked change in the distribution of photon-burst sizes. These studies provide insight into the conformational heterogeneity of G protein-coupled receptors in the presence and absence of a bound agonist.

Sorting by diffusion: An asymmetric obstacle course for continuous molecular separation

A separation technique employing a microfabricated sieve has been demonstrated by observing the motion of DNA molecules of different size. The sieve consists of a two-dimensional lattice of obstacles whose asymmetric disposition rectifies the Brownian motion of molecules driven through the device, causing them to follow paths that depend on their diffusion coefficient. A nominal 6% resolution by length of DNA molecules in the size range 15–30 kbp may be achieved in a 4-inch (10-cm) silicon wafer. The advantage of this method is that samples can be loaded and sorted continuously, in contrast to the batch mode commonly used in gel electrophoresis.

Intranuclear diffusion and hybridization state of oligonucleotides measured by fluorescence correlation spectroscopy in living cells

Fluorescein-labeled oligodeoxynucleotides (oligos) were introduced into cultured rat myoblasts, and their molecular movements inside the nucleus were studied by fluorescence correlation spectroscopy (FCS) and fluorescence recovery after photobleaching (FRAP). FCS revealed that a large fraction of both intranuclear oligo(dT) (43%) and oligo(dA) (77%) moves rapidly with a diffusion coefficient of 4 × 10−7 cm2/s. Interestingly, this rate of intranuclear oligo movement is similar to their diffusion rates measured in aqueous solution. In addition, we detected a large fraction (45%) of the intranuclear oligo(dT), but not oligo(dA), diffusing at slower rates (≤1 × 10−7 cm2/s). The amount of this slower-moving oligo(dT) was greatly reduced if the oligo(dT) was prehybridized in solution with (unlabeled) oligo(dA) prior to introduction to cells, presumably because the oligo(dT) was then unavailable for subsequent hybridization to endogenous poly(A) RNA. The FCS-measured diffusion rate for much of the slower oligo(dT) population approximated the diffusion rate in aqueous solution of oligo(dT) hybridized to a large polyadenylated RNA (1.0 × 10−7 cm2/s). Moreover, this intranuclear movement rate falls within the range of calculated diffusion rates for an average-sized heterogeneous nuclear ribonucleoprotein particle in aqueous solution. A subfraction of oligo(dT) (15%) moved over 10-fold more slowly, suggesting it was bound to very large macromolecular complexes. Average diffusion coefficients obtained from FRAP experiments were in agreement with the FCS data. These results demonstrate that oligos can move about within the nucleus at rates comparable to those in aqueous solution and further suggest that this is true for large ribonucleoprotein complexes as well.

“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.

Role of tumor–host interactions in interstitial diffusion of macromolecules: Cranial vs. subcutaneous tumors

The large size of many novel therapeutics impairs their transport through the tumor extracellular matrix and thus limits their therapeutic effectiveness. We propose that extracellular matrix composition, structure, and distribution determine the transport properties in tumors. Furthermore, because the characteristics of the extracellular matrix largely depend on the tumor–host interactions, we postulate that diffusion of macromolecules will vary with tumor type as well as anatomical location. Diffusion coefficients of macromolecules and liposomes in tumors growing in cranial windows (CWs) and dorsal chambers (DCs) were measured by fluorescence recovery after photobleaching. For the same tumor types, diffusion of large molecules was significantly faster in CW than in DC tumors. The greater diffusional hindrance in DC tumors was correlated with higher levels of collagen type I and its organization into fibrils. For molecules with diameters comparable to the interfibrillar space the diffusion was 5- to 10-fold slower in DC than in CW tumors. The slower diffusion in DC tumors was associated with a higher density of host stromal cells that synthesize and organize collagen type I. Our results point to the necessity of developing site-specific drug carriers to improve the delivery of molecular medicine to solid tumors.