Synthesis, characterization, thermal degradation, and comparative chain dynamics studies of weak-link polysulfide polymers

Synthesis, spectroscopic and thermal characterization of two new classes of polysulfide polymers: poly[1(phenoxymethyl) ethylene polysulfide] (PPMEP), and poly [1-(phenoxy) ethylene polysulfide] (PPEP) is presented. The direct pyrolysis mass spectrometry (DP-MS) technique, used to study the thermal degradation behavior of these polysulfide polymers, indicated that the polymers underwent degradation through the weak-links scission. The thermal stability of the polysulfide polymers decreased as the ``rank'' (number of sulfur atoms in the polysulfide linkage; n=1, 2, 4) increased. The main-chain flexibility of these polysulfide polymers in terms of their C-13 NMR spinlattice relaxation time (T-1) measurements on the backbone methine (-CH-) and methylene (-CH2-) carbons are reported here for the first time. A comparative study of the solution chain dynamics indicated that it increased as ``rank'' of the polysulfide linkages decreased as well as by introducing side chain spacers such as, ether (-O-) or methyleneoxy (-CH2O-) groups.

A molecular dynamics study based post facto free energy analysis of the binding of bovine angiogenin with UMP and CMP ligands

Angiogenin is a protein belonging to the superfamily of RNase A. The RNase activity of this protein is essential for its angiogenic activity. Although members of the RNase A family carry out RNase activity, they differ markedly in their strength and specificity. In this paper, we address the problem of higher specificity of angiogenin towards cytosine against uracil in the first base binding position. We have carried out extensive nano-second level molecular dynamics(MD) computer simulations on the native bovine angiogenin and on the CMP and UMP complexes of this protein in aqueous medium with explicit molecular solvent. The structures thus generated were subjected to a rigorous free energy component analysis to arrive at a plausible molecular thermodynamic explanation for the substrate specificity of angiogenin.

Lysis dynamics and membrane oligomerization pathways for Cytolysin A (ClyA) pore-forming toxin

Pore-forming toxins are known for their ability to efficiently form transmembrane pores which eventually leads to cell lysis. The dynamics of lysis and underlying self-assembly or oligomerization pathways leading to pore formation are incompletely understood. In this manuscript the pore-forming kinetics and lysis dynamics of Cytolysin-A (ClyA) toxins on red blood cells (RBCs) are quantified and compared with experimental lysis data. Lysis experiments are carried out on a fixed mass of RBCs, under isotonic conditions in phosphate-buffered saline, for different initial toxin concentrations ranging from 2.94-14.7 nM. Kinetic models which account for monomer binding, conformation and oligomerization to form the dodecameric ClyA pore complex are developed and lysis is assumed to occur when the number of pores per RBC (n(p)) exceeds a critical number, n(pc). By analysing the model in a sublytic regime (n(p) < n(pc)) the number of pores per RBC to initiate lysis is found to lie between 392 and 768 for the sequential oligomerization mechanism and between 5300 and 6300 for the non-sequential mechanism. Rupture rates which are first order in the number of RBCs are seen to provide the best agreement with the lysis experiments. The time constants for pore formation are estimated to lie between 1 and 20 s and monomer conformation time scales were found to be 2-4 times greater than the oligomerization times. Cell rupture takes places in 100s of seconds, and occurs predominantly with a steady number of pores ranging from 515 to 11 000 on the RBC surface for the sequential mechanism. Both the sequential irreversible and non-sequential kinetics provide similar predictions of the hemoglobin release dynamics, however the hemoglobin released as a function of the toxin concentration was accurately captured only with the sequential model. Each mechanism develops a distinct distribution of mers on the surface, providing a unique experimentally observable fingerprint to identify the underlying oligomerization pathways. Our study offers a method to quantify the extent and dynamics of lysis which is an important aspect of developing novel drug and gene delivery strategies based on pore-forming toxins.

Self-organized criticality in dynamics without branching

We demonstrate the phenomenon of self-organized criticality (SOC) in a simple random walk model described by a random walk of a myopic ant, i.e., a walker who can see only nearest neighbors. The ant acts on the underlying lattice aiming at uniform digging, i.e., reduction of the height profile of the surface but is unaffected by the underlying lattice. In one, two, and three dimensions we have explored this model and have obtained power laws in the time intervals between consecutive events of "digging." Being a simple random walk, the power laws in space translate to power laws in time. We also study the finite size scaling of asymptotic scale invariant process as well as dynamic scaling in this system. This model differs qualitatively from the cascade models of SOC.

Shock propagation in gas dynamics: explicit form of higher order compatibility conditions

With the use of tensor analysis and the method of singular surfaces, an infinite system of equations can be derived to study the propagation of curved shocks of arbitrary strength in gas dynamics. The first three of these have been explicitly given here. This system is further reduced to one involving scalars only. The choice of dependent variables in the infinite system is quite important, it leads to coefficients free from singularities for all values of the shock strength.

Experiments on a plume with off-source heating: Implications for cloud fluid dynamics

We report here on a series of laboratory experiments on plumes, undertaken with the object of simulating the effect of the heat release that occurs in clouds on condensation of water vapor. The experimental technique used for this purpose relies on ohmic heating generated in an electrically conducting plume fluid subjected to a suitable alternating voltage across specified axial stations in the plume flow [Bhat et al., 1989]. The present series of experiments achieves a value of the Richardson number that is toward the lower end of the range that characteristics cumulus clouds. It is found that the buoyancy enhancement due to heating disrupts the eddy structures in the flow and reduces the dilution owing to entrainment of ambient fluid that would otherwise have occurred in the central region of the plume. Heating also reduces the spread rate of the plume, but as it accelerates the flow as well, the overall specific mass flux in the plume does not show a very significant change at the heat input employed in the experiment. However, there is some indication that the entrainment rate (proportional to the streamwise derivative of the mass flux) is slightly higher immediately after heat injection and slightly lower farther downstream. The measurements support a previous proposal for a cloud scenario [Bhat and Narasimha, 1996] and demonstrate how fresh insights into certain aspects of the fluid dynamics of clouds may be derived from the experimental techniques employed here.

Cell Dynamics of Model Proteins

We design rapidly folding sequences by assigning the strongest couplings to the contacts present in a target native state in a two dimensional model of heteropolymers. The pathways to folding and their dependence on the temperature are illustrated via a mapping of the dynamics into motion within the space of the maximally compact cells.

Dynamics of ribonuclease A and ribonuclease S: Computational and experimental studies

RNase S is a complex consisting of two proteolytic fragments of RNase A: the S peptide (residues 1-20) and S protein (residues 21-124). RNase S and RNase A have very similar X-ray structures and enzymatic activities. previous experiments have shown increased rates of hydrogen exchange and greater sensitivity to tryptic cleavage for RNase S relative to RNase A. It has therefore been asserted that the RNase S complex is considerably more dynamically flexible than RNase A. In the present study we examine the differences in the dynamics of RNaseS and RNase A computationally, by MD simulations, and experimentally, using trypsin cleavage as a probe of dynamics. The fluctuations around the average solution structure during the simulation were analyzed by measuring the RMS deviation in coordinates. No significant differences between RNase S and RNase A dynamics were observed in the simulations. We were able to account for the apparent discrepancy between simulation and experiment by a simple model, According to this model, the experimentally observed differences in dynamics can be quantitatively explained by the small amounts of free S peptide and S protein that are present in equilibrium with the RNase S complex. Thus, folded RNase A and the RNase S complex have identical dynamic behavior, despite the presence of a break in polypeptide chain between residues 20 and 21 in the latter molecule. This is in contrast to what has been widely believed for over 30 years about this important fragment complementation system.

Self-Absorption and Focusing Effects in Carrier Relaxation Dynamics under Depletion Conditions

The effect of a one-dimensional field (1) on the self-absorption characteristics and (2) when we have a finite numerical aperture for the objective lens that focuses the laser beam on the solid are considered here. Self-absorption, in particular its manifestation as an inner filter for the emitted signal, has been observed in luminescence experiments. Models for this effect exist and have been analyzed, but only in the absence of space charge. Using our previous results on minority carrier relaxation in the presence of a field, we obtain expressions incorporating inner filter effects. Focusing of a light beam on the sample, by an objective lens, results in a three-dimensional source and consequently a three-dimensional continuity equation to be solved for the minority carrier concentration. Assuming a one-dimensional electric field and employing Fourier-Bessel transforms, we recast the problem of carrier relaxation and solve the same via an identity that relates it to solutions obtained in the absence of focusing effects. The inner filter effect as well as focusing introduces new time scales in the problem of carrier relaxation. The interplay between the electric field and the parameters which characterize these effects and the consequent modulation of the intensity and time scales of carrier decay signals are analyzed and discussed.

Vortex dynamics at rf frequencies in Bi2Sr2CaCu2O8 single crystals

We have studied the low magnetic field high temperature region of the H-T phase diagram of Bi2Sr2CaCu2O8 single crystals using the technique of non-resonant rf response at a frequency of 20 MHz. With H(rf)parallel to a, H parallel to c, the isothermal magnetic field scans below T-c show that the frequency f(H) of the tank circuit decreases continuously with increase in H before saturating at H similar to H-D(T). Such a decrease in f(H) reflects increasing rf penetration into the weakly screened region between CuO bilayers. The saturation of f(H) at its lowest value for H similar to H-D(T) indicates complete rf penetration land hence the disappearance of field dependence) due to the vanishing of the screening rf currents I-rf(c) in those regions or equivalently when the phase coherence between adjacent superconducting layers vanishes. Therefore H,(T) represents the decoupling of the adjacent superconducting bilayers, and hence also a 3D to 2D decoupling transition of the vortex structure. Simultaneous monitoring of the field dependent rf power dissipation P(H) shows a maximum in dP/dH at H-D(T). The observed H-D(T) line in many crystals is in excellent agreement with the (l/t-1) behavior proposed for decoupling.

Dynamics of a dilute sheared inelastic fluid. I. Hydrodynamic modes and velocity correlation functions

The hydrodynamic modes and the velocity autocorrelation functions for a dilute sheared inelastic fluid are analyzed using an expansion in the parameter epsilon=(1-e)(1/2), where e is the coefficient of restitution. It is shown that the hydrodynamic modes for a sheared inelastic fluid are very different from those for an elastic fluid in the long-wave limit, since energy is not a conserved variable when the wavelength of perturbations is larger than the ``conduction length.'' In an inelastic fluid under shear, there are three coupled modes, the mass and the momenta in the plane of shear, which have a decay rate proportional to k(2/3) in the limit k -> 0, if the wave vector has a component along the flow direction. When the wave vector is aligned along the gradient-vorticity plane, we find that the scaling of the growth rate is similar to that for an elastic fluid. The Fourier transforms of the velocity autocorrelation functions are calculated for a steady shear flow correct to leading order in an expansion in epsilon. The time dependence of the autocorrelation function in the long-time limit is obtained by estimating the integral of the Fourier transform over wave number space. It is found that the autocorrelation functions for the velocity in the flow and gradient directions decay proportional to t(-5/2) in two dimensions and t(-15/4) in three dimensions. In the vorticity direction, the decay of the autocorrelation function is proportional to t(-3) in two dimensions and t(-7/2) in three dimensions.

Dynamics of a dilute sheared inelastic fluid. II. The effect of correlations

The effect of correlations on the viscosity of a dilute sheared inelastic fluid is analyzed using the ring-kinetic equation for the two-particle correlation function. The leading-order contribution to the stress in an expansion in epsilon=(1-e)(1/2) is calculated, and it is shown that the leading-order viscosity is identical to that obtained from the Green-Kubo formula, provided the stress autocorrelation function in a sheared steady state is used in the Green-Kubo formula. A systemmatic extension of this to higher orders is also formulated, and the higher-order contributions to the stress from the ring-kinetic equation are determined in terms of the terms in the Chapman-Enskog solution for the Boltzmann equation. The series is resummed analytically to obtain a renormalized stress equation. The most dominant contributions to the two-particle correlation function are products of the eigenvectors of the conserved hydrodynamic modes of the two correlated particles. In Part I, it was shown that the long-time tails of the velocity autocorrelation function are not present in a sheared fluid. Using those results, we show that correlations do not cause a divergence in the transport coefficients; the viscosity is not divergent in two dimensions, and the Burnett coefficients are not divergent in three dimensions. The equations for three-particle and higher correlations are analyzed diagrammatically. It is found that the contributions due to the three-particle and higher correlation functions to the renormalized viscosity are smaller than those due to the two-particle distribution function in the limit epsilon -> 0. This implies that the most dominant correlation effects are due to the two-particle correlations.

Anomalous thermal dynamics of Bragg gratings inscribed in germanosilicate optical fiber

An interesting, periodic appearance of a new peak has been observed in the reflected spectrum of a Fiber Bragg Grating (FBG) inscribed in a germanosilicate fiber during thermal treatment. The new peak occurs on the longer wavelength side of the spectrum during heating and on the shorter wavelength side during cooling, following an identical reverse dynamics. Comparison with a commercial grating with 99.9% reflectivity shows a similar decay dynamics. It is proposed that the distortion due to simultaneous erasure and thermal expansion of the index modulation profile may be responsible for the observed anomaly. The reported results help us in understanding the thermal behavior of FBGs and provide additional insights into the mechanisms responsible for the photosensitivity in germanosilicate fibers.

Intermittent Dynamics, Stochastic Resonance and Dynamical Heterogeneity in Supercooled Liquid Water

Here we find through computer simulations and theoretical analysis that the low temperature thermodynamic anomalies of liquid water arises from the intermittent fluctuation between its high density and low density forms, consisting largely of 5-coordinated and 4-coordinated water molecules, respectively. The fluctuations exhibit strong dynamic heterogeneity (defined by the four point time correlation function), accompanied by a divergence like growth of the dynamic correlation length, of the type encountered in fragile supercooled liquids. The intermittency has been explained by invoking a two state model often employed to understand stochastic resonance, with the relevant periodic perturbation provided here by the fluctuation of the total volume of the system.

Thermal Degradation Processes in Poly(xylylene sulfides) Investigated by Comparative Direct Pyrolysis MS and Flash Pyrolysis GC/MS Experiments

The thermal degradation processes of two sulfur polymers, poly(xylylene sulfide) (PXM) and poly(xylylene disulfide) (PXD), were investigated in parallel by direct pyrolysis mass spectrometry (DPMS) and flash pyrolysis GC/MS (Py-GC/MS). Thermogravimetric data showed that these polymers decompose with two separate steps in the temperature ranges of 250-280 and 600-650 degrees C, leaving a high amount of residue (about 50% at 800 degrees C). The pyrolysis products detected by DPMS in the first degradation step of PXM and PXD were terminated by three types of end groups, -CH3, -CH2SH, and -CH=S, originating from thermal cleavage reactions involving a series of homolytic chain scissions followed by hydrogen transfer reactions, generating several oligomers containing some intact xylylene sulfide repeating units. The presence of pyrolysis compounds containing some stilbene-like units in the first degradation step has also been observed. Their formation has been accounted for with a parallel cleavage involving the elimination of H2S from the PXM main chains. These unsaturated units can undergo cross-linking at higher temperatures, producing the high amount of char residue observed. The thermal degradation compounds detected by DPMS in the second decomposition step at about 600-650 degrees C were constituted of condensed aromatic molecules containing dihydrofenanthrene and fenanthrene units. These compounds might be generated from the polymer chains containing stilbene units, by isomerization and dehydrogenation reactions. The pyrolysis products obtained in the Py-GC/MS of PXM and PXD at 610 degrees C are almost identical. The relative abundance in the pyrolysate and the spectral properties of the main pyrolysis products were found to be in generally good agreement with those obtained by DPMS. Polycyclic aromatic hydrocarbons (PAHs) were also detected by Py-GC/MS but in minor amounts with respect to DPMS. This apparent discrepancy was due to the simultaneous detection of PAHs together with all pyrolysis products in the Py-GC/MS, whereas in DPMS they were detected in the second thermal degradation step without the greatest part of pyrolysis compounds generated in the first degradation step. The results obtained by DPMS and PSI-GC/MS experiments showed complementary data for the degradation of PXM and PXD and, therefore, allowed the unequivocal formulation of the thermal degradation mechanism for these sulfur-containing polymers.