Optimal fits of diffusion constants from single-time data points of Brownian trajectories

Experimental methods based on single particle tracking (SPT) are being increasingly employed in the physical and biological sciences, where nanoscale objects are visualized with high temporal and spatial resolution. SPT can probe interactions between a particle and its environment but the price to be paid is the absence of ensemble averaging and a consequent lack of statistics. Here we address the benchmark question of how to accurately extract the diffusion constant of one single Brownian trajectory. We analyze a class of estimators based on weighted functionals of the square displacement. For a certain choice of the weight function these functionals provide the true ensemble averaged diffusion coefficient, with a precision that increases with the trajectory resolution.

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.

Configurational diffusion down a folding funnel describes the dynamics of DNA hairpins

Elucidating the mechanism of folding of polynucleotides depends on accurate estimates of free energy surfaces and a quantitative description of the kinetics of structure formation. Here, the kinetics of hairpin formation in single-stranded DNA are measured after a laser temperature jump. The kinetics are modeled as configurational diffusion on a free energy surface obtained from a statistical mechanical description of equilibrium melting profiles. The effective diffusion coefficient is found to be strongly temperature-dependent in the nucleation step as a result of formation of misfolded loops that do not lead to subsequent zipping. This simple system exhibits many of the features predicted from theoretical studies of protein folding, including a funnel-like energy surface with many folding pathways, trapping in misfolded conformations, and non-Arrhenius folding rates.

Diffusion modeling of percutaneous absorption kinetics: 3. Variable diffusion and partition coefficients, consequences for stratum corneum depth profiles and desorption kinetics

Stratum corneum (SC) desorption experiments have yielded higher calculated steady-state fluxes than those obtained by epidermal penetration studies. A possible explanation of this result is a variable diffusion or partition coefficient across the SC. We therefore developed the diffusion model for percutaneous penetration and desorption to study the effects of either a variable diffusion coefficient or variable partition coefficient in the SC over the diffusion path length. Steady-state flux, lag time, and mean desorption time were obtained from Laplace domain solutions. Numerical inversion of the Laplace domain solutions was used for simulations of solute concentration-distance and amount penetrated (desorbed)-time profiles. Diffusion and partition coefficients heterogeneity were examined using six different models. The effect of heterogeneity on predicted flux from desorption studies was compared with that obtained in permeation studies. Partition coefficient heterogeneity had a more profound effect on predicted fluxes than diffusion coefficient heterogeneity. Concentration-distance profiles show even larger dependence on heterogeneity, which is consistent with experimental tape-stripping data reported for clobetasol propionate and other solutes. The clobetasol propionate tape-stripping data were most consistent with the partition coefficient decreasing exponentially for half the SC and then becoming a constant for the remaining SC. (C) 2004 Wiley-Liss, Inc.

Diffusion of n-decane in mesoporous MCM-41 silicas

We investigate here the diffusion of n-decane in nanoporous MCM-41 silicas with pore diameters between 3.0 and 4.3 nm, and at various temperatures and purge flow rates, by the Zero Length Column method. A complete-time-range analysis of desorption curves is used to derive the diffusion coefficient, and the effect of pore size, purge flow rate and temperature on the diffusion character is systematically studied. The results show that the calculated low-coverage diffusivity values are strongly dependent on temperature but only weakly dependent on pore size. The study reveals that transport is controlled by intracrystalline diffusion and dominated by sorbate-sorbent interaction, with the experimental isosteric heat matching the potential energy of flat-lying n-decane molecules on the surface, determined using a united atom model. The diffusion activation energy and adsorption isosteric heat at zero loading for the different pore size MCM-41 samples vary in a narrow range respectively, and their ratio is essentially constant over the pore size range studied. The study shows that the ZLC method is an effective tool to investigate the diffusion kinetics of hydrocarbons in mesoporous MCM-41 materials. (c) 2005 Elsevier Inc. All rights reserved.

Mass uptake study of the diffusion of water and SBF into poly(2-hydroxyethyl methacrylate-co-tetrahydrofurfuryI methacrylate) containing aspirin or vitamin B-12

The ingress of water and Kokubo simulated body fluid (SBF) into poly (2-hydroxyethyl methacrylate) (PHEMA), and its co-polymers with tetrahydrofurduryl methacrylate (THFMA), loaded with either one of two model drugs, vitamin 1312 or aspirin, was studied by mass uptake over the temperature range 298-318 K. The polymers were studied as cylinders and were loaded with either 5 wt% or 10 wt% of the drugs. From DSC studies it was observed that vitamin B-12 behaved as a physical cross-linker restricting chain segmental mobility, and so had a small anti-plasticisation effect on PHEMA and the co-polymers rich in HEMA, but almost no effect on the T-g of co-polymers rich in THFMA. On the other hand, aspirin exhibited a plasticising effect on PHEMA and the copolymers. All of the polymers were found to absorb water and SBF according to a Fickian diffusion mechanism. The polymers were all found to swell to a greater extent in SBF than in water, which was attributed to the presence of Tris buffer in the SBF. The sorptions of the two penetrants were found to follow Fickian kinetics in all cases and the diffusion coefficients at 310 K for SBF were found to be smaller than those for water, except for the polymers containing aspirin where the diffusion coefficients were higher than for the other systems. For example, for sorption into PHEMA the diffusion coefficient for water was 1.41 X 10(-11) m(2)/s and for SBF was 0.79 x 10-11 m(2)/s, but in the presence of 5 wt% aspirin the corresponding values were 1.27 x 10(-1)1 m(2)/s and 1.25 x 10(-11) m(2)/s, respectively. The corresponding values for PHEMA loaded with 5 wt% B-12 were 1.25 x 10(-11) m(2)/s and 0.74 x 10(-11) m(2)/s, respectively.

Salt diffusion and distribution in meat studied by Na-23 nuclear magnetic resonance imaging and relaxometry

This study introduces the use of combined Na-23 magnetic resonance imaging (MRI) and Na-23 NMR relaxometry for the study of meat curing. The diffusion of sodium ions into the meat was measured using Na-23 MRI on a 1 kg meat sample brined in 10% w/w NaCl for 3-100 h. Calculations revealed a diffusion coefficient of 1 x 10(-5) cm(2)/s after 3 h of curing and subsequently decreasing to 8 x 10(-6) cm(2)/s at longer curing times, suggesting that changes occur in the microscopic structure of the meat during curing. The microscopic mobility and distribution of sodium was measured using Na-23 relaxometry. Two sodium populations were observed, and with increasing length of curing time the relaxation times of these changed, reflecting a salt-induced swelling and increase in myofibrillar pore sizes. Accordingly, the present study demonstrated that pore size and thereby salt-induced swelling in meat can be assessed using Na-23 relaxometry.

Stress corrosion cracking and hydrogen diffusion in magnesium

Evaluation of recent data for hydrogen (H) diffusion in magnesium (Mg) yielded a new equation for the diffusion coefficient of H in Mg. This indicates that there can be significant H transport ahead of a stress corrosion crack in Mg at ambient temperature and that H may be involved in the mechanism of stress corrosion cracking in Mg.

NMR microscopy: A tool for measuring monomer diffusion in supercritical CO2

The diffusion of styrene into linear low density polyethylene in a solution of supercritical CO2 was investigated using NMR microimaging. For both pure styrene and styrene dissolved in supercritical CO2, the diffusion was found to follow Fickian kinetics. Supercritical CO2 was found to enhance the rate and extent of diffusion of styrene into the substrate by up to three times under the conditions of this investigation, compared to pure styrene. NMR imaging was used to measure the concentration profiles of the styrene penetrants in real time, and the results were fitted to a Fickian model for diffusion. At a CO2 pressure of 150 bar and temperature of 40 degrees C, the diffusion coefficient of a 30 wt-% solution of styrene into LLDPE was calculated to be 1 X 10(-11) m(2). s(-1). This is significantly faster than the diffusion coefficient measured for pure styrene diffusion at 40 degrees C (3 x 10(-12) m(2). s(-1)). The diffusion coefficients determined by gravimetric analysis were slightly higher than those determined by the imaging method. This was probably due to residual styrene and/or polystyrene adhering to the surface of the films in the gravimetric technique.

Differential diffusion: Often a finite-mising length effect

Many instances of differential diffusion, i e, different species having different turbulent diffusion coefficients in the same flow, can be explained as a finite mixing length effect. That is, in a simple mixing length scenario, the turbulent diffusion coefficient has the form 1 ( m )2 m m c l K w l OL   =  +    where, wm is the mixing velocity, lm the mixing length and Lc the overall distribution scale for a particular species. The first term represents the familiar gradient diffusion while the second term becomes important when lm/Lc is finite. This second term shows that different species will have different diffusion coefficients if they have different overall distribution scales. Such different Lcs may come about due to different boundary conditions and different intrinsic properties (molecular diffusivity, settling velocity etc) for different species. For momentum transfer in turbulent oscillatory boundary layers the second term is imaginary and explains observed phase leads of shear stresses ahead of velocity gradients.

Quantitative interpretation of diffusion-ordered NMR spectra:can we rationalize small molecule diffusion coefficients?

Appealingly simple: A new method is described that allows the diffusion coefficient of a small molecule to be estimated given only the molecular weight and the viscosity of the solvent used. This method makes possible the quantitative interpretation of the diffusion domain of diffusion-ordered NMR spectra (see picture). © 2013 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

Inverse Problem for Fractional Diffusion Equation

MSC 2010: 26A33, 33E12, 34K29, 34L15, 35K57, 35R30

Radiotherapy Assessment Using Diffusion Weighted MRI

Purpose: There are two goals of this study. The first goal of this study is to investigate the feasibility of using classic textural feature extraction in radiotherapy response assessment among a unique cohort of early stage breast cancer patients who received the single-dose preoperative radiotherapy. The second goal of this study is to investigate the clinical feasibility of using classic texture features as potential biomarkers which are supplementary to regional apparent diffusion coefficient in gynecological cancer radiotherapy response assessment.

Methods and Materials: For the breast cancer study, 15 patients with early stage breast cancer were enrolled in this retrospective study. Each patient received a single-fraction radiation treatment, and DWI and DCE-MRI scans were conducted before and after the radiotherapy. DWI scans were acquired using a spin-echo EPI sequence with diffusion weighting factors of b = 0 and b = 500 mm2/s, and the apparent diffusion coefficient (ADC) maps were calculated. DCE-MRI scans were acquired using a T1-weighted 3D SPGR sequence with a temporal resolution of about 1 minute. The contrast agent (CA) was intravenously injected with a 0.1 mmol/kg bodyweight dose at 2 ml/s. Two parameters, volume transfer constant (Ktrans) and kep were analyzed using the two-compartment Tofts pharmacokinetic model. For pharmacokinetic parametric maps and ADC maps, 33 textural features were generated from the clinical target volume (CTV) in a 3D fashion using the classic gray level co-occurrence matrix (GLCOM) and gray level run length matrix (GLRLM). Wilcoxon signed-rank test was used to determine the significance of each texture feature’s change after the radiotherapy. The significance was set to 0.05 with Bonferroni correction.

For the gynecological cancer study, 12 female patients with gynecologic cancer treated with fractionated external beam radiotherapy (EBRT) combined with high dose rate (HDR) intracavitary brachytherapy were studied. Each patient first received EBRT treatment followed by five fractions of HDR treatment. Before EBRT and before each fraction of brachytherapy, Diffusion Weighted MRI (DWI-MRI) and CT scans were acquired. DWI scans were acquired in sagittal plane utilizing a spin-echo echo-planar imaging sequence with weighting factors of b = 500 s/mm2 and b = 1000 s/mm2, one set of images of b = 0 s/mm2 were also acquired. ADC maps were calculated using linear least-square fitting method. Distributed diffusion coefficient (DDC) maps and stretching parameter α were calculated. For ADC and DDC maps, 33 classic texture features were generated utilizing the classic gray level run length matrix (GLRLM) and gray level co-occurrence matrix (GLCOM) from high-risk clinical target volume (HR-CTV). Wilcoxon signed-rank statistics test was applied to determine the significance of each feature’s numerical value change after radiotherapy. Significance level was set to 0.05 with multi-comparison correction if applicable.

Results: For the breast cancer study, regarding ADC maps calculated from DWI-MRI, 24 out of 33 CTV features changed significantly after the radiotherapy. For DCE-MRI pharmacokinetic parameters, all 33 CTV features of Ktrans and 33 features of kep changed significantly.

For the gynecological cancer study, regarding ADC maps, 28 out of 33 HR-CTV texture features showed significant changes after the EBRT treatment. 28 out of 33 HR-CTV texture features indicated significant changes after HDR treatments. The texture features that indicated significant changes after HDR treatments are the same as those after EBRT treatment. 28 out of 33 HR-CTV texture features showed significant changes after whole radiotherapy treatment process. The texture features that indicated significant changes for the whole treatment process are the same as those after HDR treatments.

Conclusion: Initial results indicate that certain classic texture features are sensitive to radiation-induced changes. Classic texture features with significant numerical changes can be used in monitoring radiotherapy effect. This might suggest that certain texture features might be used as biomarkers which are supplementary to ADC and DDC for assessment of radiotherapy response in breast cancer and gynecological cancer.

Diagnostic value of diffusion weighted MRI and ADC in differential diagnosis of cavernous hemangioma of the liver.

Aims: To investigate the use of diffusion weighted magnetic resonance imaging (DWI) and the apparent diffusion coefficient (ADC) values in the diagnosis of hemangioma. Materials and methods: The study population consisted of 72 patients with liver masses larger than 1 cm (72 focal lesions). DWI examination with a b value of 600 s/mm2 was carried out for all patients. After DWI examination, an ADC map was created and ADC values were measured for 72 liver masses and normal liver tissue (control group). The average ADC values of normal liver tissue and focal liver lesions, the “cut-off” ADC values, and the diagnostic sensitivity and specificity of the ADC map in diagnosing hemangioma, benign and malignant lesions were researched. Results: Of the 72 liver masses, 51 were benign and 21 were malignant. Benign lesions comprised 38 hemangiomas and 13 simple cysts. Malignant lesions comprised 9 hepatocellular carcinomas, and 12 metastases. The highest ADC values were measured for cysts (3.782±0.53×10-3 mm2/s) and hemangiomas (2.705±0.63×10-3 mm2/s). The average ADC value of hemangiomas was significantly higher than malignant lesions and the normal control group (p<0.001). The average ADC value of cysts were significantly higher when compared to hemangiomas and normal control group (p<0.001). To distinguish hemangiomas from malignant liver lesions, the “cut-off” ADC value of 1.800×10-3 mm2/s had a sensitivity of 97.4% and a specificity of 90.9%. To distinguish hemangioma from normal liver parenchyma the “cut-off” value of 1.858×10-3 mm2/s had a sensitivity of 97.4% and a specificity of 95.7%. To distinguish benign liver lesions from malignant liver lesions the “cut-off” value of 1.800×10-3 mm2/s had a sensitivity of 96.1% and a specificity of 90.0%. Conclusion: DWI and quantitative measurement of ADC values can be used in differential diagnosis of benign and malignant liver lesions and also in the diagnosis and differentiation of hemangiomas. When dynamic examination cannot distinguish cases with vascular metastasis and lesions from hemangioma, DWI and ADC values can be useful in the primary diagnosis and differential diagnosis. The technique does not require contrast material, so it can safely be used in patients with renal failure. Keywords: