Li+ counterion self-diffusion in ordered DNA

The anisotropic self-diffusion coefficient of ⁷Li⁺ (I = 3/2) counterions has been studied in hydrated, macroscopically oriented Li-(B)DNA fibers at relatively high water contents, corresponding to approximate DNA-DNA helix axis distances of 22–35 Å, using the pulsed field gradient hmr spin-echo method. Self-diffusion coefficients parallel (D_∥) and perpendicular (D_?) to the DNA helix axis increase with increasing salt content and with increasing DNA-DNA helix axis distance. The observed anisotropy D_∥/D_? decreases from 1.6 to 1.2 with the DNA-DNA separation increasing from 22 to 35 Å in the salt-free sample. This result can be understood by the obstruction effect caused by the DNA molecules themselves. The values of the Li⁺ self-diffusion coefficients in the most water-rich system with no added salt (corresponding to an approximate distance of 35 Å between the DNA helix axes) were D_∥ ～ 1.15 × 10^?10 m² s^?1 and D_? ～ 0.98 × 10^?10 m² s^?1, compared to 9.14 × 10^?10 m² s^?1 for the diffusion of Li⁺ in an aqueous solution of LiCl (～ 2.1M). The possible occurrence of restriction effects in the DNA fibers have also been studied by determining the self-diffusion coefficient at different effective diffusion times. The self-diffusion coefficient of Li⁺ in the sample with the largest DNA-DNA helix axis distance seems to be independent of the effective diffusion time, which indicates that the lithium ions are not trapped within impermeable barriers. The possibility of diffusion through permeable barriers has also been investigated, and is discussed. © 1994 John Wiley & Sons, Inc.     

How to predict diffusion of medium-sized molecules in polymer matrices. From atomistic to coarse grain simulations

The normal diffusion regime of many small and medium-sized molecules occurs on a time scale that is too long to be studied by atomistic simulations. Coarse-grained (CG) molecular simulations allow to investigate length and time scales that are orders of magnitude larger compared to classical molecular dynamics simulations, hence providing a valuable approach to span time and length scales where normal diffusion occurs. Here we develop a novel multi-scale method for the prediction of diffusivity in polymer matrices which combines classical and CG molecular simulations. We applied an atomistic-based method in order to parameterize the CG MARTINI force field, providing an extension for the study of diffusion behavior of penetrant molecules in polymer matrices. As a case study, we found the parameters for benzene (as medium sized penetrant molecule whose diffusivity cannot be determined through atomistic models) and Poly (vinyl alcohol) (PVA) as polymer matrix. We validated our extended MARTINI force field determining the self diffusion coefficient of benzene (2.27·10⁻⁹ m² s⁻¹) and the diffusion coefficient of benzene in PVA (0.263·10⁻¹² m² s⁻¹). The obtained diffusion coefficients are in remarkable agreement with experimental data (2.20·10⁻⁹ m² s⁻¹ and 0.25·10⁻¹² m² s⁻¹, respectively). We believe that this method can extend the application range of computational modeling, providing modeling tools to study the diffusion of larger molecules and complex polymeric materials.     

Brownian Dynamics Simulations of Colloidal Liquids: Hydrodynamics and Stress Relaxation

Abstract

We describe the statistical mechanics background and additional algorithmic features of a recently proposed simple mean-field Brownian Dynamics algorithm formulated to include many-body hydrodynamics, using a local density approximation for the friction coefficient. We show that the equations of motion satisfy the incompressibility of phase space. We make further developments to the model, computing the hydrodynamic effects on the shear stress relaxation function. We show that stress relaxation takes place over two well-defined regimes, in both cases with and without mean field hydrodynamics, MFH. At short times ta ²/D ₀ < 10^?3, where a is the radius of the colloidal particle and D ₀ is the self-diffusion coefficient at infinite dilution, decay of the stress autocorrelation function, C_s(t) is essentially independent of volume fraction and does not fit to a simple analytic form. At longer times than ta ²/D ₀ < 10^?2 the decay has the fractional exponential form ～exp(-t ^β) with β ? 1. The transition between these two regimes coincides with a rapid fall in the time-dependent diffusion coefficient from the so-called short-time to long-time values. We do not find any evidence for power law decay in the C_s(t) as predicted by recent mode-coupling based analytical expansions.     

A molecular-dynamics simulation study of diffusion of a single model carbonic chain on a graphite (001) surface

Molecular-dynamics simulations have been used to study the diffusion of a short single model carbonic chain on the graphite (001) surface. The calculated diffusion coefficient (D) first increases, then decreases with increasing chain length (N). This abnormal behavior is similar to polymer lateral diffusion at the solid–liquid interface. Furthermore, we have studied the relation between the mean-square gyration radius and N. Figure Log–log plot of the self-diffusion coefficient D versus the chain length N. The error bars are the standard deviation measured in three repeated simulations     

A facile synthesis of ceramides.

Lateral diffusion of phosphatide molecules in liquid crystalline bilayers has been analysed as a case of co-operative lattice diffusion. The potential energy of interaction between two molecules is assumed to arise from Van der Waals interactions of the hydrocarbon chains, and to have the form suggested by Salem [6]. From the observed values of the self-diffusion constant (of the order of 10^?8 cm² sec^?1) the depth of the potential “well” for two molecules at the equilibrium separation was estimated to have a lower limit of 1.95 kcal per mole, and the energy barrier to lateral motion was estimated to have an upper limit of 7.21 kcal per mole.     

Modelling Diffusion in Zeolites by Molecular Dynamics Simulations

Abstract

The diffusion of molecules sorbed in zeolites is of growing interest for understanding the mechanisms of chemical processes with regard to selectivity and reactivity [1].

MD simulations give insight into physical systems on the molecular level allowing to study and visualize the motion of molecules even beyond the possibilities of experiments [2,3]. Single system parameters can easily be varied to study their influence, also those parameters that are fixed in reality (e.g., the size of particles). We present a cross section of our recent work to illustrate the capabilities of MD: The self diffusion coefficients (D) of a mixture of methane and xenon in silicalite show remarkable deviations from those of the pure species. This is shown and confirmed by PFG NMR experiments [4].

Simulating ethane in zeolite A the mechanism of diffusion has been studied. The effects of rotation on the diffusion lead to cases where D decreases with growing temperature [5].

The independence of self diffusion on lattice vibrations is proven even for zeolites with windows of guest particle size comparing simulations with rigid and vibrating zeolite lattice [6].     

Analysis of diffusion and relaxation behavior of water in apple parenchymal cells

It has been demonstrated by an example of apple parenchymal cells that NMR spectroscopy can be used to analyze the relaxation and diffusion of water molecules in plant cells. With small diffusion times, three relaxation components have been distinguished, which correspond to water in a vacuole, in the cytoplasm, and in intercellular liquid. The coefficient of self-diffusion corresponding to these components have been determined. With large diffusion times, it is possible to distinguish two components. For the slowly relaxing component (which corresponds to water in a vacuole), the regime of restricted diffusion was observed. For a quickly relaxing component, an anomalous increase in the coefficient of self-diffusion with the time of diffusion took place.     

Intracellular Transport Dynamics of Endosomes Containing DNA Polyplexes along the Microtubule Network

We have explored the transport of DNA polyplexes enclosed in endosomes within the cellular environment by multiple particle tracking (MPT). The polyplex-loaded endosomes demonstrate enhanced diffusion at short timescales (t < 7 s) with their mean-square displacement (MSD) 〈Δx(t)²〉 scaling as t^1.25. For longer time intervals they exhibit subdiffusive transport and have an MSD scaling as t^0.7. This crossover from an enhanced diffusion to a subdiffusive regime can be explained by considering the action of motor proteins that actively transport these endosomes along the cellular microtubule network and the thermal bending modes of the microtubule network itself.     

Alkaline Polymer Membrane‐Based Ultrathin,Flexible, and High‐Performance Solid‐State Zn‐Air Battery

The increasing demand for portable and wearable electronics requires lightweight, thin, and highly flexible power sources, for example, flexible zinc‐air batteries (ZABs). The so‐far reported flexible ZAB devices mostly remain bulky, with a design consisting of two relatively thick substrates (e.g., carbon cloths and/or metal foams) and a gel electrolyte‐coated separator in between. Herein, an ultrathin (≈0.2 mm) solid‐state ZAB with high flexibility and performance is introduced by directly forming self‐standing active layers on each surface of an alkaline polymer membrane through an ink‐casting/hot‐pressing approach. A Fe/N‐doped 3D carbon with hierarchic pores and an interconnected network structure is used as cathode electrocatalyst, so that the backing gas‐diffusion layer (e.g., carbon cloth) can be abandoned. What is further, a microstructure‐modulating method to significantly increase the FeN₄ active sites for oxygen reduction reaction is developed, thus significantly boosting the performance of the ZAB. The assembled solid‐state ZAB manifests remarkable peak power density of 250 mW cm^?3 and high capacity of 150.4 mAh cm^?3 at 8.3 mA cm^?3, as well as excellent flexibility. The new design should provide valuable opportunity to the portable and wearable electronics.     

Molecular dynamics investigation of nitric oxide (II) interaction with a model biological membrane

Nitric oxide (II) diffusion through a model two-component (phosphatidylcholine and phosphatidylethanolamine molecules) biological membrane is investigated using the molecular dynamics method. It is shown that NO molecules are rotating in the process of diffusion into the phospholipid bilayer. The calculated diffusion coefficient D[NO] = 0.35 (±0.23) × 10^?5 (cm²/s) is in a good agreement with literature data. This testifies that free diffusion of NO molecules may be a plausible mechanism of the NO permeation through the membranes.     

Molecular Dynamics Study of Water in Hydrogels

Abstract

Molecular dynamics simulations are performed for aqueous solutions of polymers: Poly (vinyl alcohol) (PVA), Poly (vinyl methylether) (PVME), and Poly (N-isopropyl acrylamide) (PNiPAM). The distributions and dynamics of hydrogen-bonds, the translational diffusion of water, and the orientational relaxation of water are analyzed to investigate the properties of water which is highly influenced by the surrounding polymer chains. The water molecules around the polymer chains are highly hindered by the chains.     

Static and dynamic properties of mid-size liquid n- alkanes,C12∼C400: a molecular dynamics simulation study

ABSTRACT

In this paper, we have extended our previous study of the static and dynamic properties (self-diffusion coefficient D_self and friction coefficient ζ) of liquid n-alkane systems up C₄₀₀ at several temperatures (～2300?K) using molecular dynamics (MD) simulations in the canonical ensembles. For the small n-alkanes with n?≤?120 (n: the chain length), the chains are clearly ?R² _ee?/6?R² _g? ≥ 1 (1.06 ～ 1.44), which leads to the conclusion that the liquid n-alkanes are far away from the ideal chain regime. But for the n-alkanes of n?≥?160, the chains are ?R² _ee?/6?R² _g? ≈ 1, indicating that they are Gaussian. It is found that the long chains of these n-alkanes at high temperatures show abnormalities in density and friction coefficient. We observed a clear transition in the power law dependence of n-alkane self-diffusion coefficient on the molecular weight (M) of n-alkane, D_self ～ M^?γ, occurs in the range C₁₂₀～C₁₆₀ at temperatures of 318, and 618?K, corresponding to a crossover from the ‘oligomer’ to the ‘Rouse’ regime. The entanglement lengths (N_e) are calculated by the Z1 code and discussed shortly.     

Bacterial denitrifying nitric oxide reductases and aerobic respiratory terminal oxidases use similar delivery pathways for their molecular substrates

The superfamily of heme?copper oxidoreductases (HCOs) include both NO and O₂ reductases. Nitric oxide reductases (NORs) are bacterial membrane enzymes that catalyze an intermediate step of denitrification by reducing nitric oxide (NO) to nitrous oxide (N₂O). They are structurally similar to heme?copper oxygen reductases (HCOs), which reduce O₂ to water. The experimentally observed apparent bimolecular rate constant of NO delivery to the deeply buried catalytic site of NORs was previously reported to approach the diffusion-controlled limit (10⁸–10⁹?M^?1?s^?1). Using the crystal structure of cytochrome-c dependent NOR (cNOR) from Pseudomonas aeruginosa, we employed several protocols of molecular dynamics (MD) simulation, which include flooding simulations of NO molecules, implicit ligand sampling and umbrella sampling simulations, to elucidate how NO in solution accesses the catalytic site of this cNOR. The results show that NO partitions into the membrane, enters the enzyme from the lipid bilayer and diffuses to the catalytic site via a hydrophobic tunnel that is resolved in the crystal structures. This is similar to what has been found for O₂ diffusion through the closely related O₂ reductases. The apparent second order rate constant approximated using the simulation data is ~5?×?10⁸?M^?1?s^?1, which is optimized by the dynamics of the amino acid side chains lining in the tunnel. It is concluded that both NO and O₂ reductases utilize well defined hydrophobic tunnels to assure that substrate diffusion to the buried catalytic sites is not rate limiting under physiological conditions.     

Molecular Dynamics Simulation of Self Diffusion in Microclusters

Abstract

Molecular dynamics simulation is carried out to study the mechanism of self diffusion which is characteristic of solid-like microclusters. A two-dimensional system with the Lennard-Jones potential is employed and the temperatures near the triple point of the two-dimensional bulk system are adopted for the simulation. The results show that: a) microclusters consist of two regions, i.e., solid-like cores and liquid-like surface regions, b) the size dependence of the diffusion coefficient for microclusters is weak, and the value of the solid-like core region is not much different from that of the bulk liquid, c) the activation energy of diffusion for microclusters is twenty to thirty times larger than that for the bulk liquid, d) the diffusion mechanism in the solid-like region involves the collective motion of small domains containing ten to twenty atoms which results in the formation of low density regions, sometimes even vacancy clusters, between them, and atoms in the low density regions change their positions to cause diffusion.     

Triplet excited states of polycyclic aromatic compounds as probes of their microenvironment in serum albumin complexes.

Using flash photolysis techniques, the triplet excited states of benzo(a)pyrene, pyrene, benz(a)anthracene and other aromatic hydrocarbons have been detected in complexes of bovine (and human) serum albumin dissolved in aqueous solutions at room temperature. The triplet lifetimes can be adjusted to any value within the microsecond-millisecond time domains by varying the partial pressure of oxygen from zero to one atmosphere, thus providing a useful probe on these time scales. Local oxygen concentrations as low as ～ 2 × 10^?7M can be detected. In air saturated solutions, the triplet lifetimes are sensitive to pH dependent conformational changes of the host bovine serum albumin molecules.     

Flow extremes and benthic organic matter shape the metabolism of a headwater Mediterranean stream

1. Single‐station diel oxygen curves were used to monitor the oxygen metabolism of an intermittent, forested third‐order stream (Fuirosos) in the Mediterranean area, over a period of 22 months. Ecosystem respiration (ER) and gross primary production (GPP) were estimated and related to organic matter inputs and photosynthetically active radiation (PAR) in order to understand the effect of the riparian forest on stream metabolism. 2. Annual ER was 1690 g O₂ m^?2 year^?1 and annual GPP was 275 g O₂ m^?2 year^?1. Fuirosos was therefore a heterotrophic stream, with P : R ratios averaging 0.16. 3. GPP rates were relatively low, ranging from 0.05 to 1.9 g O₂ m^?2 day^?1. The maximum values of GPP occurred during a few weeks in spring, and ended when the riparian canopy was fully closed. The phenology of the riparian vegetation was an important determinant of light availability, and consequently, of GPP. 4. On a daily scale, light and temperature were the most important factors governing the shape of photosynthesis–irradiance (P–I) curves. Several patterns could be generalised in the P–I relationships. Hysteresis‐type curves were characteristic of late autumn and winter. Light saturation responses (that occurred at irradiances higher than 90 μE m^?2 s^?1) were characteristic of early spring. Linear responses occurred during late spring, summer and early autumn when there was no evidence of light saturation. 5. Rates of ER were high when compared with analogous streams, ranging from 0.4 to 32 g O₂ m^?2 day^?1. ER was highest in autumn 2001, when organic matter accumulations on the streambed were extremely high. By contrast, the higher discharge in autumn 2002 prevented these accumulations and caused lower ER. The Mediterranean climate, and in its effect the hydrological regime, were mainly responsible for the temporal variation in benthic organic matter, and consequently of ER.     

Dynamic light scattering from collagen solutions. I. Translational diffusion coefficient and aggregation effects

Solutions of native collagen extracted from rat tail tendons in neutral salt solution have been studied by dynamic light scattering. The spectra obtained are consistent with the presence in solution of both single rod-shaped collagen molecules and aggregates of molecules. No contribution to the spectrum has been detected at any scattering angle from rotational diffusion of single molecules, although a measurable broadening effect is expected at high angles. The translational diffusion coefficient D of single molecules, calculated from the broader spectral component, shows an anomalous dependence on collagen concentration with a maximum value of D_20,w = 8.6 ± 0.2 × 10^?12 m²/sec near the concentration 0.04% by weight. Above 0.05% D falls linearly with increasing concentration and takes the value D _20,w = 8.1 ± 0.2 × 10^?12 m²/sec at 0.064% collagen.     

Molecular dynamics simulation of liquid water under the influence of an external electric field

Molecular dynamics simulations of liquid water were performed at 258K and a density of 1.0?g/cm³ under various applied external electric field, ranging 0～10¹⁰?V/m. The influence of external field on structural and dynamical properties of water was investigated. The simple point charge (SPC) model is used for water molecules. An enhancement of the water hydrogen bond structure with increasing strength of the electric field has been deduced from the radial distribution functions and the analysis of hydrogen bonds structure. With increasing field strength, water system has a more perfect structure, which is similar to ice structure. However, the electrofreezing phenomenon of liquid water has not been detected since the self-diffusion coefficient was very large. The self-diffusion coefficient decreases remarkably with increasing strength of electric field and the self-diffusion coefficient is anisotropic.     

Molecular Dynamics Simulation Studies of the Density and Temperature Dependence of Self-diffusion in a Cylindrical Micropore

We report Molecular Dynamics calculations of radial density profiles and self-diffusion coefficients of Lennard-Jones fluids in a cylindrical pore of radius 2σ, for a wide range of temperatures and densities. At n _p σ³ = 0.825 the self-diffusion coefficient parallel to the pore walls D _∥*. follows a monotonic (nearly linear) increase with kT/ε and is very similar to that of the bulk self-diffusion coefficient D _b *. At n _p σ³ = 0.4 and kT/ε ≤ 1.0 the curve of D _∥* vs. kT/ε shows a distinct inflection in the region 0.7 ≤ kT/ε ≤ 0.9 and values of D _∥* are much less than D _b * decreasing to near solid state values at very low temperatures. At the highest temperature studied, kT/ε = 2.98, D _∥* is almost inversely proportional to density and in a fairly close agreement with that of D _b *. At KT/ε = 0.49, D _∥* is much smaller than D _b *. The motion of adsorbate particles normal to the walls is also discussed.