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.

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.

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.

Dynamics and transport properties of Kondo insulators

A many-body theory of paramagnetic Kondo insulators is described, focusing specifically on single-particle dynamics, scattering rates, dc transport and optical conductivities. This is achieved by development of a non-perturbative local moment approach to the symmetric periodic Anderson model within the framework of dynamical mean-field theory. Our natural focus is the strong-coupling, Kondo lattice regime, in particular the resultant 'universal' scaling behaviour in terms of the single, exponentially small low-energy scale characteristic of the problem. Dynamics/transport on all relevant (ω, T)-scales are considered, from the gapped/activated behaviour characteristic of the low-temperature insulator through to explicit connection to single-impurity physics at high ω and/or T; and for optical conductivities emphasis is given to the nature of the optical gap, the temperature scale responsible for its destruction and the consequent clear distinction between indirect and direct gap scales. Using scaling, explicit comparison is also made to experimental results for dc transport and optical conductivities of Ce3Bi4Pt3, SmB6 and YbB12. Good agreement is found, even quantitatively; and a mutually consistent picture of transport and optics results.

Dynamics of bound dyes in a nonpolymeric aqueous gel derived from a tripodal bile salt

Rotational dynamics of polarity sensitive fluorescent dyes (ANS and DPH) in a nonpolymertic aqueous gel derived from tripodal cholamide I was studied using ultrafast time-resolved fluorescence technique. Results were compared with that of naturally occurring di- and trihydroxy bile salts. ANS in the gel showed two rotational correlation time (phi) components, 13.2 ns (bound to the hydrophobic region of the gel) and 1.0 ns (free aqueous ANS), whereas DPH showed only one component (4.8 ns). In the sol state, faster rotational motion was observed, both for ANS and DPH. Our data revealed that dyes get encapsulated more tightly in the gel network when compared to the micellar aggregates. ANS has more restrained rotation compared to DPH. This was attributed to the interaction of the sulfonate group of ANS with water molecules and hydrophilic parts of the gelator molecule. No restricted rotation was observed for DPH in the gel state unlike when it is in the gel phase of lipid bilayer.

Dynamics of pulsed flow in an elastic tube

Internal haemorrhage, often leading to cardio-vascular arrest happens to be one of the prime sources of high fatality rates in mammals. We propose a simplistic model of fluid flow in our attempt to specify the location of the haemorrhagic spot, which, if located accurately, could possibly be operated leading to an instant cure. The model we employ for the purpose is basically fluid mechanical in origin and consists of a viscous fluid, pumped by a periodic force and flowing through an elastic tube. The analogy is with that of blood, pumped from the heart and flowing through an artery or vein. Our results, aided by graphical illustrations, match reasonably well with experimental observations.