Nonlinear dynamics of eucaryotic pyruvate dehydrogenase multienzyme complex: Decarboxylation rate, oscillations, and multiplicity

Pyruvate conversion to acetyl-CoA by the pyruvate dehydrogenase (PDH) multienzyme complex is known as a key node in affecting the metabolic fluxes of animal cell culture. However, its possible role in causing possible nonlinear dynamic behavior such as oscillations and multiplicity of animal cells has received little attention. In this work, the kinetic and dynamic behavior of PDH of eucaryotic cells has been analyzed by using both in vitro and simplified in vivo models. With the in vitro model the overall reaction rate (v(1)) of PDH is shown to be a nonlinear function of pyruvate concentration, leading to oscillations under certain conditions. All enzyme components affect v, and the nonlinearity of PDH significantly, the protein X and the core enzyme dihydrolipoamide acyltransferase (E2) being mostly predominant. By considering the synthesis rates of pyruvate and PDH components the in vitro model is expanded to emulate in vivo conditions. Analysis using the in vivo model reveals another interesting kinetic feature of the PDH system, namely, multiple steady states. Depending on the pyruvate and enzyme levels or the operation mode, either a steady state with high pyruvate decarboxylation rate or a steady state with significantly lower decarboxylation rate can be achieved under otherwise identical conditions. In general, the more efficient steady state is associated with a lower pyruvate concentration. A possible time delay in the substrate supply and enzyme synthesis can also affect the steady state to be achieved and lead's to oscillations under certain conditions. Overall, the predictions of multiplicity for the PDH system agree qualitatively well with recent experimental observations in animal cell cultures. The model analysis gives some hints for improving pyruavte metabolism in animal cell culture.

Photodissociation dynamics of CH2ICl at 222, 236, 266, 280, and-304 nm

Dynamics of I*(P-2(1/2)) formation from CH2ICl dissociation has-been investigated at five different ultraviolet excitation wavelengths, e.g., 222, 236, 266, 280, and similar to304 nm. The quantum yield of I*((2)p(1/2)) production, phi*, has been measured by monitoring nascent I(P-2(3/2)) and I* concentrations using a resonance enhanced multiphoton ionization detection scheme. The measured quantum yield as a function of excitation energy follows the same trend as that of methyl iodide except at 236 run. The photodissociation dynamics of CH2ICl also involves three upper states similar to methyl iodide, and a qualitative correlation diagram has been constructed to account for the observed quantum yield. From the difference in behavior at 236 nm, it appears that the crossing region between the two excited states ((3)Q(0) and (1)Q(1)) is located near the exit valley away from the Franck Condon excitation region. The B- and C-band transitions do not participate in the dynamics, and the perturbation of the methyl iodide states due to Cl-I interaction is relatively weak at the photolysis wavelengths employed in this investigation.

Role of length scales on microwave thawing dynamics in 2D cylinders

Microwave (MW) thawing of 2D frozen cylinders exposed to uniform plane waves from one face, is modeled using the effective heat capacity formulation with the MW power obtained from the electric field equations. Computations are illustrated for tylose (23% methyl cellulose gel) which melts over a range of temperatures giving rise to a mushy zone. Within the mushy region the dielectric properties are functions of the liquid volume fraction. The resulting coupled, time dependent non-linear equations are solved using the Galerkin finite element method with a fixed mesh. Our method efficiently captures the multiple connected thawed domains that arise due to the penetration of MWs in the sample. For a cylinder of diameter D, the two length scales that control the thawing dynamics are D/D-p and D/lambda(m), where D-p and lambda(m) are the penetration depth and wavelength of radiation in the sample respectively. For D/D-p, D/lambda(m) much less than 1 power absorption is uniform and thawing occurs almost simultaneously across the sample (Regime I). For D/D-p much greater than 1 thawing is seen to occur from the incident face, since the power decays exponentially into the sample (Regime III). At intermediate values, 0.2 < D/D-p, D/lambda(m) < 2.0 (Regime II) thawing occurs from the unexposed face at smaller diameters, from both faces at intermediate diameters and from the exposed and central regions at larger diameters. Average power absorption during thawing indicates a monotonic rise in Regime I and a monotonic decrease in Regime III. Local maxima in the average power observed for samples in Regime II are due to internal resonances within the sample. Thawing time increases monotonically with sample diameter and temperature gradients in the sample generally increase from Regime I to Regime III. (C) 2002 Elsevier Science Ltd. All rights reserved.

Slow orientational dynamics of water molecules at a micellar surface

The dynamics of water molecules near an aqueous micellar interface is studied in an atomistic molecular dynamics simulation of cesium pentadecafluorooctanoate (CsPFO) in water. The dipolar orientational time correlation function (tcf) and the translational diffusion of the water molecules are investigated. Results show that both the reorientational and the translational motion of water molecules near the micelle are restricted. In particular, the orientational tcf exhibits a very slow component in the long time which is slower than its bulk value by 2 orders of magnitude. This slow decay seems to be related to the slow decay often observed in experiments. The origin of the slow decay is analyzed.

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.

Effect of orientational motion of mobile chromophores on the dynamics of Forster energy transfer in polymers

Brownian dynamics (BD) simulations have been carried out to explore the effects of the orientational motion of the donor-acceptor (D-A) chromophore pair on the Forster energy transfer between the D-A pair embedded in a polymer chain in solution. It is found that the usually employed orientational averaging (that is, replacing the orientational factor, kappa, by kappa (2) = 2/3) may lead to an error in the estimation of the rate of the reaction by about 20%. In the limit of slow orientational relaxation, the preaveraging of the orientational factor leads to an overestimation of the rate, while in the opposite limit of very fast orientational relaxation, the usual scheme underestimates the rate. The latter results from an interesting interplay between reaction and diffusion. On the other hand, when one of the chromophores is fixed, the preaveraged rate is found to be fairly reliable if the rotational relaxation of the chromophore is sufficiently fast. The present study also reveals a power law dependence of the FRET rate on the chain length (rate proportional to N- alpha, with alpha approximate to 2.6).

Nonlinear Dynamics of Beams with Stochastic Parameter Variations

The aim of this paper is to investigate the steady state response of beams under the action of random support motions. The study is of relevance in the context of earthquake response of extended land based structures such as pipelines and long span bridges, and, secondary systems such as piping networks in nuclear power plant installations. The following complicating features are accounted for in the response analysis: (a) differential support motions: this is characterized in terms of cross power spectral density functions associated with distinct support motions, (b) nonlinear support conditions, and (c) stochastically inhomogeneous stiffness and mass variations of the beam structure; questions on non-Gaussian models for these variations are considered. The method of stochastic finite elements is combined with equivalent linearization technique and Monte Carlo simulations to obtain response moments.

Synthesis, structural characterization, thermal studies and chain dynamics of poly(methacrylonitrile peroxide) by NMR spectroscopy

Poly(methacrylonitrile peroxide) (PMNP) has been synthesized from methacrylonitrile by free radical initiated oxidative polymerization and characterized by different spectroscopic methods. NMR spectroscopy confirmed the alternating copolymer structure with labile peroxy bonds in the main chain. The extreme instability of PMNP was noted from FTIR spectroscopy. Thermal degradation studies by using differential scanning calorimetry and thermogravimetry have revealed that PMNP degrades highly exothermically and the heat of degradation, 42.5 kcal mol−1, is of the same order as that reported for other vinyl polyperoxides. Mass spectral fragmentation pattern under electron impact (EI) condition has also been investigated. The mechanism of the primary exothermic degradation has been substantiated by thermochemical calculations. The chain dynamics of the polyperoxide chain has been studied by means of 13C spin–lattice relaxation times (T1) of the main chain as well as the side chain carbons. The temperature dependence of the spin–lattice relaxation times shows that the PMNP is more flexible compared to the analogous poly(styrene peroxide).

Understanding glassy dynamics from the free-energy landscape

Many interesting features of the dynamics of simple liquids near the glass transition may be understood in terms of properties of the free-energy landscape obtained from numerical studies of a model free-energy functional. Main results obtained from this approach are summarized and a list of references to relevant publications is provided. (C) 2002 Elsevier Science B.V. All rights reserved.

A nonlinear spectral finite element model for analysis of wave propagation in solid with internal friction and dissipation

A geometrically non-linear Spectral Finite Flement Model (SFEM) including hysteresis, internal friction and viscous dissipation in the material is developed and is used to study non-linear dissipative wave propagation in elementary rod under high amplitude pulse loading. The solution to non-linear dispersive dissipative equation constitutes one of the most difficult problems in contemporary mathematical physics. Although intensive research towards analytical developments are on, a general purpose cumputational discretization technique for complex applications, such as finite element, but with all the features of travelling wave (TW) solutions is not available. The present effort is aimed towards development of such computational framework. Fast Fourier Transform (FFT) is used for transformation between temporal and frequency domain. SFEM for the associated linear system is used as initial state for vector iteration. General purpose procedure involving matrix computation and frequency domain convolution operators are used and implemented in a finite element code. Convergnence of the spectral residual force vector ensures the solution accuracy. Important conclusions are drawn from the numerical simulations. Future course of developments are highlighted.

On the dynamics of activation of mammalian X chromosomes

We have investigated a mathematical model of the process of activation of the X chromosomes in eutherian mammals. The model assumes that the activation is brought about over some definite time interval T by the complete saturation of N receptor sites on an X chromosome by M activating molecules (or multiples of M). The probability λ of a first hit on the receptor site is considered to be very much lower than that of subsequent hits; that is, we assume strong co-operative binding. Assuming further that an incomplete saturation of receptor sites is malfunctional, we can show that for proper activation of X chromosomes in normal diploid males and females, we must have λMT ≥ 3 and 0·96 ≤ N/M ≤ 1. An extension of this analysis for the triploid cases shows that under these conditions, we cannot explain the activation of two X's if the number of activating molecules is fixed at M. This suggests that there must be two classes of triploid embryos differing from each other in a step-wise manner in the number of activating molecules. In other words, triploids with two active X chromosomes would require 2M activating molecules as opposed to M molecules in triploids with a single active X. This interpretation of the two classes of triploids would be consistent with differing imprinting histories of the parental contributions to the triploid zygote.

Particle dynamics in a turbulent particle–gas suspension at high Stokes number. Part 2. The fluctuating-force model

A fluctuating-force model is developed for representing the effect of the turbulent fluid velocity fluctuations on the particle phase in a turbulent gas–solid suspension in the limit of high Stokes number, where the particle relaxation time is large compared with the correlation time for the fluid velocity fluctuations. In the model, a fluctuating force is incorporated in the equation of motion for the particles, and the force distribution is assumed to be an anisotropic Gaussian white noise. It is shown that this is equivalent to incorporating a diffusion term in the Boltzmann equation for the particle velocity distribution functions. The variance of the force distribution, or equivalently the diffusion coefficient in the Boltzmann equation, is related to the time correlation functions for the fluid velocity fluctuations. The fluctuating-force model is applied to the specific case of a Couette flow of a turbulent particle–gas suspension, for which both the fluid and particle velocity distributions were evaluated using direct numerical simulations by Goswami & Kumaran (2010). It is found that the fluctuating-force simulation is able to quantitatively predict the concentration, mean velocity profiles and the mean square velocities, both at relatively low volume fractions, where the viscous relaxation time is small compared with the time between collisions, and at higher volume fractions, where the time between collisions is small compared with the viscous relaxation time. The simulations are also able to predict the velocity distributions in the centre of the Couette, even in cases in which the velocity distribution is very different from a Gaussian distribution.

The collective dynamics of self-propelled particles

We propose a method for the dynamic simulation of a collection of self-propelled particles in a viscous Newtonian fluid. We restrict attention to particles whose size and velocity are small enough that the fluid motion is in the creeping flow regime. We propose a simple model for a self-propelled particle, and extended the Stokesian Dynamics method to conduct dynamic simulations of a collection of such particles. In our description, each particle is treated as a sphere with an orientation vector p, whose locomotion is driven by the action of a force dipole Sp of constant magnitude S0 at a point slightly displaced from its centre. To simplify the calculation, we place the dipole at the centre of the particle, and introduce a virtual propulsion force Fp to effect propulsion. The magnitude F0 of this force is proportional to S0. The directions of Sp and Fp are determined by p. In isolation, a self-propelled particle moves at a constant velocity u0 p, with the speed u0 determined by S0. When it coexists with many such particles, its hydrodynamic interaction with the other particles alters its velocity and, more importantly, its orientation. As a result, the motion of the particle is chaotic. Our simulations are not restricted to low particle concentration, as we implement the full hydrodynamic interactions between the particles, but we restrict the motion of particles to two dimensions to reduce computation. We have studied the statistical properties of a suspension of self-propelled particles for a range of the particle concentration, quantified by the area fraction φa. We find several interesting features in the microstructure and statistics. We find that particles tend to swim in clusters wherein they are in close proximity. Consequently, incorporating the finite size of the particles and the near-field hydrodynamic interactions is of the essence. There is a continuous process of breakage and formation of the clusters. We find that the distributions of particle velocity at low and high φa are qualitatively different; it is close to the normal distribution at high φa, in agreement with experimental measurements. The motion of the particles is diffusive at long time, and the self-diffusivity decreases with increasing φa. The pair correlation function shows a large anisotropic build-up near contact, which decays rapidly with separation. There is also an anisotropic orientation correlation near contact, which decays more slowly with separation. Movies are available with the online version of the paper.

Particle dynamics in a turbulent particle–gas suspension at high Stokes number. Part 1. Velocity and acceleration distributions

The effect of fluid velocity fluctuations on the dynamics of the particles in a turbulent gas–solid suspension is analysed in the low-Reynolds-number and high Stokes number limits, where the particle relaxation time is long compared with the correlation time for the fluid velocity fluctuations, and the drag force on the particles due to the fluid can be expressed by the modified Stokes law. The direct numerical simulation procedure is used for solving the Navier–Stokes equations for the fluid, the particles are modelled as hard spheres which undergo elastic collisions and a one-way coupling algorithm is used where the force exerted by the fluid on the particles is incorporated, but not the reverse force exerted by the particles on the fluid. The particle mean and root-mean-square (RMS) fluctuating velocities, as well as the probability distribution function for the particle velocity fluctuations and the distribution of acceleration of the particles in the central region of the Couette (where the velocity profile is linear and the RMS velocities are nearly constant), are examined. It is found that the distribution of particle velocities is very different from a Gaussian, especially in the spanwise and wall-normal directions. However, the distribution of the acceleration fluctuation on the particles is found to be close to a Gaussian, though the distribution is highly anisotropic and there is a correlation between the fluctuations in the flow and gradient directions. The non-Gaussian nature of the particle velocity fluctuations is found to be due to inter-particle collisions induced by the large particle velocity fluctuations in the flow direction. It is also found that the acceleration distribution on the particles is in very good agreement with the distribution that is calculated from the velocity fluctuations in the fluid, using the Stokes drag law, indicating that there is very little correlation between the fluid velocity fluctuations and the particle velocity fluctuations in the presence of one-way coupling. All of these results indicate that the effect of the turbulent fluid velocity fluctuations can be accurately represented by an anisotropic Gaussian white noise.