Computation of restoration of ligand response in the random kinetics of a prostate cancer cell signaling pathway

Mutation and/or dysfunction of signaling proteins in the mitogen activated protein kinase (MAPK) signal transduction pathway are frequently observed in various kinds of human cancer. Consistent with this fact, in the present study, we experimentally observe that the epidermal growth factor (EGF) induced activation profile of MAP kinase signaling is not straightforward dose-dependent in the PC3 prostate cancer cells. To find out what parameters and reactions in the pathway are involved in this departure from the normal dose-dependency, a model-based pathway analysis is performed. The pathway is mathematically modeled with 28 rate equations yielding those many ordinary differential equations (ODE) with kinetic rate constants that have been reported to take random values in the existing literature. This has led to us treating the ODE model of the pathways kinetics as a random differential equations (RDE) system in which the parameters are random variables. We show that our RDE model captures the uncertainty in the kinetic rate constants as seen in the behavior of the experimental data and more importantly, upon simulation, exhibits the abnormal EGF dose-dependency of the activation profile of MAP kinase signaling in PC3 prostate cancer cells. The most likely set of values of the kinetic rate constants obtained from fitting the RDE model into the experimental data is then used in a direct transcription based dynamic optimization method for computing the changes needed in these kinetic rate constant values for the restoration of the normal EGF dose response. The last computation identifies the parameters, i.e., the kinetic rate constants in the RDE model, that are the most sensitive to the change in the EGF dose response behavior in the PC3 prostate cancer cells. The reactions in which these most sensitive parameters participate emerge as candidate drug targets on the signaling pathway. (C) 2011 Elsevier Ireland Ltd. All rights reserved.

Kinetics of oxidation of boron powder

The kinetics of the oxidation of electrodeposited boron powder and the boron powder produced by the reduction process were studied using thermogravimetry (TG). The oxidation was carried out by heating boron powder in a stream of oxygen. Both isothermal and non-isothermal methods were used to study the kinetics. Model-free isoconversional method was used to derive the kinetics parameters. A two step oxidation reaction (exothermic) was observed. The oxidation reaction could not be completed due to the formation of glassy layer of boric oxide on the surface of boron powder which acts as a barrier for further diffusion of oxygen into the particle. The activation energy obtained using model-free method for electrodeposited boron is 122 +/- 7 kJ mol(-1) whereas a value of 205 +/- 9 kJ mol(-1) was obtained for boron produced by the reduction process (commercially procured boron). Mechanistic interpretation of the oxidation reaction was done using model based method. The activation energy was found to depend on the size distribution of the particles and specific surface area of the powder. (C) 2010 Elsevier B.V. All rights reserved.

Shear alignment of a disordered lamellar mesophase

The shear alignment of an initially disordered lamellar phase is examined using lattice Boltzmann simulations of a mesoscopic model based on a free-energy functional for the concentration modulation. For a small shear cell of width 8 lambda, the qualitative features of the alignment process are strongly dependent on the Schmidt number Sc = nu/D (ratio of kinematic viscosity and mass diffusion coefficient). Here, lambda is the wavelength of the concentration modulation. At low Schmidt number, it is found that there is a significant initial increase in the viscosity, coinciding with the alignment of layers along the extensional axis, followed by a decrease at long times due to the alignment along the flow direction. At high Schmidt number, alignment takes place due to the breakage and reformation of layers because diffusion is slow compared to shear deformation; this results in faster alignment. The system size has a strong effect on the alignment process; perfect alignment takes place for a small systems of width 8 lambda and 16 lambda, while a larger system of width 32 lambda does not align completely even at long times. In the larger system, there appears to be a dynamical steady state in which the layers are not perfectly aligned-where there is a balance between the annealing of defects due to shear and the creation due to an instability of the aligned lamellar phase under shear. We observe two types of defect creation mechanisms: the buckling instability under dilation, which was reported earlier, as well as a second mechanism due to layer compression.

Formation and propagation of bands in jerky flow:a coupled lattice map description

There has been revival of interest in Jerky flow from the point of view of dynamical systems. The earliest attempt in this direction was from our group. One of the predictions of the theory is that Jerky flow could be chaotic. This has been recently verified by us. We have recently extended the earlier model to account for the spatial aspect as well. Both these models are in the form of coupled set of nonlinear differential equations and hence, they are complicated in their structure. For this reason we wish to devise a model based on the results of these two theories in the form of coupled lattice map for the description of the formation and propagation of dislocation bands. We report here one such model and its results.

Diurnal variation of atmospheric CO2 concentration and delta C-13 in an urban atmosphere during winter-role of the Nocturnal Boundary Layer

The paper presents the importance of the Nocturnal Boundary Layer in driving the diurnal variability of the atmospheric CO2 mixing ratio and the carbon isotope ratio at ground level from an urban station in India. Our observations are the first of their kind from this region. The atmospheric CO2 mixing ratio and the carbon isotopic ratio were measured for both the morning (05:30-07:30 IST) and afternoon time (16:00-18:00 IST) air samples at 5 m above ground level in Bangalore city, Karnataka State (12 degrees 58' N, 77 degrees 38' E, masl = 920 m) for a 10 day period during the winter of 2008. We observed a change of similar to 7% the in CO2 mixing ratio between the morning and afternoon time air samples. A stable isotope analysis of CO2 from morning samples showed a depletion in the carbon isotope ratio by similar to 2 parts per thousand compared to the afternoon samples. Along with the ground-based measurement of air samples, data of radiosonde measurements were also obtained from the Indian Meteorological Department to identify the vertical atmospheric structure at different time in a day. We proposed the presence or absence of the NBL as a controlling factor for the observed variability in the mixing ratio as well as its isotopic composition. Here we used the Keeling model approach to find out the carbon isotope ratio for the local sources. The local sources have further been characterized as anthropogenic and biological respiration (in %) using a two-component mixing model. We also used a vertical mixing model based on the concept of the mixing of isotopically depleted (carbon isotope) ``polluted air'' (PA) with isotopically enriched ``free atmospheric air'' (FA) above. Using this modeling approach, the contribution of FA at ground level is being estimated for both the morning and afternoon time air samples.

Electrical transport studies in alkali borovanadate glasses

The ac conductivity and dielectric behaviors of sodium borovanadate glasses have been studied over wide ranges of composition and frequency. The de activation energies calculated from the complex impedance plots decrease linearly with the Na2O concentration, indicating that ionic conductivity dominates in these glasses. The possible origin of low-temperature departures of conductivity curves (from linearity) of vanadium-rich glasses in log sigma versus 1/T plots is discussed. The ac conductivities have been fitted to the Almond-West type power law expression with use of a single value of s. It is found that in most of the glasses s exhibits a temperature-dependent minimum. The dielectric data are converted into moduli (M*) and are analyzed using the Kohlrausch-William-Watts stretched exponential function, The activation barriers, W, calculated from the temperature-dependent dielectric loss peaks compare well with the activation barriers calculated from the de conductivity plots. The stretching exponent beta is found to be temperature independent and is not likely to be related as in the equation beta = 1 - s, An attempt is made to elucidate the origin of the stretching phenomena. It appears that either a model of the increased contribution of polarization energy (caused by the increased modifier concentration) and hence the increased monopole-induced dipole interactions or a model based on increased intercationic interactions can explain the slowing down of the primitive relaxation in ionically conducting glasses.

Orientational relaxation in dipolar systems: How much do we understand the role of correlations

Dipolar systems, both liquids and solids, constitute a class of naturally abundant systems that are important in all branches of natural science. The study of orientational relaxation provides a powerful method to understand the microscopic properties of these systems and, fortunately, there are many experimental tools to study orientational relaxation in the condensed phases. However, even after many years of intense research, our understanding of orientational relaxation in dipolar systems has remained largely imperfect. A major hurdle towards achieving a comprehensive understanding is the long range and complex nature of dipolar interactions which also made reliable theoretical study extremely difficult. These difficulties have led to the development of continuum model based theories, which although they provide simple, elegant expressions for quantities of interest, are mostly unsatisfactory as they totally neglect the molecularity of inter-molecular interactions. The situation has improved in recent years because of renewed studies, led by computer simulations. In this review, we shall address some of the recent advances, with emphasis on the work done in our laboratory at Bangalore. The reasons for the failure of the continuum model, as revealed by the recent Brownian dynamics simulations of the dipolar lattice, are discussed. The main reason is that the continuum model predicts too fast a decay of the torque-torque correlation function. On the other hand, a perturbative calculation, based on Zwanzig's projection operator technique, provides a fairly satisfactory description of the single particle orientational dynamics for not too strongly polar dipolar systems. A recently developed molecular hydrodynamic theory that properly includes the effects of intermolecular orientational pair correlations provides an even better description of the single-particle orientational dynamics. We also discuss the rank dependence of the dielectric friction. The other topics reviewed here includes dielectric relaxation and solvation dynamics, as they are intimately connected with orientational relaxation. Recent molecular dynamics simulations of the dipolar lattice are also discussed. The main theme of the present review is to understand the effects of intermolecular interactions on orientational relaxation. The presence of strong orientational pair correlation leads to a strong coupling between the single particle and the collective dynamics. This coupling can lead to rich dynamical properties, some of which are detailed here, while a major part remains yet unexplored.

Origin of carrier-type reversal in Pb-Ge-Se glasses: A detailed thermal, electrical, and structural study

A long-standing and important problem in glass science has been carrier-type reversal (CTR) in semiconducting glasses. This phenomenon is exhibited by Pb-Ge-Se glasses also. It has been addressed here by carrying out detailed electrical, thermal, and spectroscopic investigations. PbxGe42-xSe58 (x = 0-20) glasses were prepared by a two stage melt-quenching process and characterized using x-ray diffraction, high-resolution electron microscropy, and energy dispersive analysis of x-rays. Thermoelectric power and high-pressure electrical resistivity have been measured. IR, Raman, and X-ray adsorption near edge structure spectroscopies have been used for examining the glass structures as well as differential scanning calorimetry (DSC) for studying the thermal properties. A structural model based on the chemical nature of the constituents has been proposed to account for the observed properties of these glasses. Effect of Pb incorporation on local structures and qualitative consequences on the energy band structures of Ge-Se glasses has been considered. The p -->n transition has been attributed to the energetic disposition of the sp(3)d(2) band of Pb atoms, which is located closely above the lone pair band of selenium. This feature makes Pb unique in the context of p -->n transition of chalcogenide glasses. The model can be extended successfully to account for the CTR behavior observed in Bi containing chalcogenide glasses also.

Relative performances of textural models in estimating soil moisture characteristic

The soil moisture characteristic (SMC) forms an important input to mathematical models of water and solute transport in the unsaturated-soil zone. Owing to their simplicity and ease of use, texture-based regression models are commonly used to estimate the SMC from basic soil properties. In this study, the performances of six such regression models were evaluated on three soils. Moisture characteristics generated by the regression models were statistically compared with the characteristics developed independently from laboratory and in-situ retention data of the soil profiles. Results of the statistical performance evaluation, while providing useful information on the errors involved in estimating the SMC, also highlighted the importance of the nature of the data set underlying the regression models. Among the models evaluated, the one possessing an underlying data set of in-situ measurements was found to be the best estimator of the in-situ SMC for all the soils. Considerable errors arose when a textural model based on laboratory data was used to estimate the field retention characteristics of unsaturated soils.

Molybdenum tetracarbonyl complexes of (aryloxy) lambda(3)-cyclotriphosphazanes; X-ray crystal structure of [Mo(CO)(4){EtNP(OC6H4Br-4)}(3)]

Treatment of the lambda(3)-cyclotriphosphazanes, cis-{EtNP(OR)}(3) [R = C6H4Br-4 (L-1) or C6H4Br-2 (L-2)] with [Mo(CO)(4)(NBD)] (NBD = norbornadiene) yields the mononuclear complexes [Mo(CO)(4)L-1] (1) and [Mo(CO)(4)L-2] (2). which have been characterised by IR, NMR (P-31 and H-1) and FAB mass spectral data. The structure of 1 has been confirmed by single crystal X-ray analysis. The structural and conformational changes brought about by complexation are discussed in terms of a bonding model based on "negative hyperconjugation". (C) 1998 Elsevier Science Ltd. All rights reserved.

Stochastic dynamics modeling of the protein sequence length distribution in genomes: implications for microbial evolution

In this paper, we report an analysis of the protein sequence length distribution for 13 bacteria, four archaea and one eukaryote whose genomes have been completely sequenced, The frequency distribution of protein sequence length for all the 18 organisms are remarkably similar, independent of genome size and can be described in terms of a lognormal probability distribution function. A simple stochastic model based on multiplicative processes has been proposed to explain the sequence length distribution. The stochastic model supports the random-origin hypothesis of protein sequences in genomes. Distributions of large proteins deviate from the overall lognormal behavior. Their cumulative distribution follows a power-law analogous to Pareto's law used to describe the income distribution of the wealthy. The protein sequence length distribution in genomes of organisms has important implications for microbial evolution and applications. (C) 1999 Elsevier Science B.V. All rights reserved.

Identification of delamination in composite beams using spectral estimation and a genetic algorithm

An efficient strategy for identification of delamination in composite beams and connected structures is presented. A spectral finite-element model consisting of a damaged spectral element is used for model-based prediction of the damaged structural response in the frequency domain. A genetic algorithm (GA) specially tailored for damage identification is derived and is integrated with finite-element code for automation. For best application of the GA, sensitivities of various objective functions with respect to delamination parameters are studied and important conclusions are presented. Model-based simulations of increasing complexity illustrate some of the attractive features of the strategy in terms of accuracy as well as computational cost. This shows the possibility of using such strategies for the development of smart structural health monitoring softwares and systems.

Polarization switching studies in ferroelectric glycine phosphite single crystals

Glycine Phosphite [NH3CH2COOH3PO3], abbreviated as GPI, undergoes a para-ferroelectric phase transition from the monoclinic symmetry P2(1)/a to P2(1) at 224.7 K. We report here a systematic study of the polarization switching process in this crystal. Growth of these crystals from aqueous solution has been undertaken employing both solvent evaporation and slow cooling methods. Hysteresis loop measurements along the polar b-axis yielded a spontaneous polarization value of 0.5 muC/cm(2) and a coercive field of 2.5 kV/cm. Conventional Merz technique was employed for polarization switching studies, wherein bipolar square pulses were applied to the sample to induce domain reversal. The transient switching pulse that flows through the sample on application of the field was recorded. The maximum switching time required for domain switching was measured both as a function of electric field and temperature. The experimentally observed switching curves were fitted with the model based on the Pulvari-Kuebler theory of nucleation and growth of domains. From the experimental data, the values of mobility and activation field were obtained. It was observed that switching process in this crystal is predominantly governed by the forward growth of domain walls in the high field region. However, switching process in GPI crystal was found to be slower than that found in other glycine based ferroelectric crystals.

Steady compressible flow of cohesionless granular materials through a wedge-shaped bunker

A continuum model based on the critical-state theory of soil mechanics is used to generate stress, density, and velocity profiles, and to compute discharge rates for the flow of granular material in a mass flow bunker. The bin–hopper transition region is idealized as a shock across which all the variables change discontinuously. Comparison with the work of Michalowski (1987) shows that his experimentally determined rupture layer lies between his prediction and that of the present theory. However, it resembles the former more closely. The conventional condition involving a traction-free surface at the hopper exit is abandoned in favour of an exit shock below which the material falls vertically with zero frictional stress. The basic equations, which are not classifiable under any of the standard types, require excessive computational time. This problem is alleviated by the introduction of the Mohr–Coulomb approximation (MCA). The stress, density, and velocity profiles obtained by integration of the MCA converge to asymptotic fields on moving down the hopper. Expressions for these fields are derived by a perturbation method. Computational difficulties are encountered for bunkers with wall angles θw [gt-or-equal, slanted] 15° these are overcome by altering the initial conditions. Predicted discharge rates lie significantly below the measured values of Nguyen et al. (1980), ranging from 38% at θw = 15° to 59% at θw = 32°. The poor prediction appears to be largely due to the exit condition used here. Paradoxically, incompressible discharge rates lie closer to the measured values. An approximate semi-analytical expression for the discharge rate is obtained, which predicts values within 9% of the exact (numerical) ones in the compressible case, and 11% in the incompressible case. The approximate analysis also suggests that inclusion of density variation decreases the discharge rate. This is borne out by the exact (numerical) results – for the parameter values investigated, the compressible discharge rate is about 10% lower than the incompressible value. A preliminary comparison of the predicted density profiles with the measurements of Fickie et al. (1989) shows that the material within the hopper dilates more strongly than predicted. Surprisingly, just below the exit slot, there is good agreement between theory and experiment.

Full Sphere Axisymmetric Simulations of the Solar Dynamo

We explore a full sphere (2D axisymmetric) kinematic solar dynamo model based on the Babcock-Leighton idea that the poloidal field is generated in the surface layers from the decay of tilted bipolar solar active regions. This model incorporates the helioseismically deduced solar rotation profile and an algorithm for buoyancy motivated from simulations of flux tube dynamics. A prescribed deep meridional circulation plays an important role in the advection of magnetic flux. We specifically address the parity issue and show that – contrary to some recent claims – the Babcock-Leighton dynamo can reproduce solar-like dipolar parity if certain reasonable conditions are satisfied in the solar interior, the most important requirement being that the poloidal field of the two hemispheres be efficiently coupled across the equator.