Particle dynamics in the channel flow of a turbulent particle-gas suspension at high Stokes number. Part 1. DNS and fluctuating force model

The fluctuating force model is developed and applied to the turbulent flow of a gas-particle suspension in a channel in the limit of high Stokes number, where the particle relaxation time is large compared to the fluid correlation time, and low particle Reynolds number where the Stokes drag law can be used to describe the interaction between the particles and fluid. In contrast to the Couette flow, the fluid velocity variances in the different directions in the channel are highly non-homogeneous, and they exhibit significant variation across the channel. First, we analyse the fluctuating particle velocity and acceleration distributions at different locations across the channel. The distributions are found to be non-Gaussian near the centre of the channel, and they exhibit significant skewness and flatness. However, acceleration distributions are closer to Gaussian at locations away from the channel centre, especially in regions where the variances of the fluid velocity fluctuations are at a maximum. The time correlations for the fluid velocity fluctuations and particle acceleration fluctuations are evaluated, and it is found that the time correlation of the particle acceleration fluctuations is close to the time correlations of the fluid velocity in a `moving Eulerian' reference, moving with the mean fluid velocity. The variances of the fluctuating force distributions in the Langevin simulations are determined from the time correlations of the fluid velocity fluctuations and the results are compared with direct numerical simulations. Quantitative agreement between the two simulations are obtained provided the particle viscous relaxation time is at least five times larger than the fluid integral time.

Mathematical Model of Viral Kinetics In Vitro Estimates the Number of E2-CD81 Complexes Necessary for Hepatitis C Virus Entry

Interaction between the hepatitis C virus (HCV) envelope protein E2 and the host receptor CD81 is essential for HCV entry into target cells. The number of E2-CD81 complexes necessary for HCV entry has remained difficult to estimate experimentally. Using the recently developed cell culture systems that allow persistent HCV infection in vitro, the dependence of HCV entry and kinetics on CD81 expression has been measured. We reasoned that analysis of the latter experiments using a mathematical model of viral kinetics may yield estimates of the number of E2-CD81 complexes necessary for HCV entry. Here, we constructed a mathematical model of HCV viral kinetics in vitro, in which we accounted explicitly for the dependence of HCV entry on CD81 expression. Model predictions of viral kinetics are in quantitative agreement with experimental observations. Specifically, our model predicts triphasic viral kinetics in vitro, where the first phase is characterized by cell proliferation, the second by the infection of susceptible cells and the third by the growth of cells refractory to infection. By fitting model predictions to the above data, we were able to estimate the threshold number of E2-CD81 complexes necessary for HCV entry into human hepatoma-derived cells. We found that depending on the E2-CD81 binding affinity, between 1 and 13 E2-CD81 complexes are necessary for HCV entry. With this estimate, our model captured data from independent experiments that employed different HCV clones and cells with distinct CD81 expression levels, indicating that the estimate is robust. Our study thus quantifies the molecular requirements of HCV entry and suggests guidelines for intervention strategies that target the E2-CD81 interaction. Further, our model presents a framework for quantitative analyses of cell culture studies now extensively employed to investigate HCV infection.

Validation of a Modified Eddy Dissipation Concept Model for Stationary and Nonstationary Turbulent Diffusion Flames

A transient flame simulation tool based on unsteady Reynolds average Navier Stokes (RANS) is characterized for stationary and nonstationary flame applications with a motivation of performing computationally affordable flame stability studies. Specifically, the KIVA-3V code is utilized with incorporation of a recently proposed modified eddy dissipation concept for simulating turbulence-chemistry interaction along with a model for radiation loss. Detailed comparison of velocities, turbulent kinetic energies, temperature, and species are made with the experimental data of the turbulent, non-premixed DLR_A CH4/H-2/N-2 jet flame. The comparison shows that the model is able to predict flame structure very well. The effect of some of the modeling assumptions is assessed, and strategies to model a stationary diffusion flame are recommended. Unsteady flame simulation capabilities of the numerical model are assessed by simulating an acoustically excited, experimental, oscillatory H-2-air diffusion flame. Comparisons are made with oscillatory velocity field and OH plots, and the numerical code is observed to predict transient flame structure well.

Influence of the Methane-Zeolite A Interaction Potential on the Concentration Dependence of Self-diffusivity

Studies on the diffusion of methane in a zeolite structure type LTA (as per IZA nomenclature) have indicated that different types of methane zeolite potentials exist in the literature in which methane is treated within the united-atom model. One set of potentials, referred to as model A, has a methane oxygen diameter of 3.14 angstrom, while another set of potential parameters, model B, employs a larger value of 3.46 angstrom. Fritzsche and co-workers (1993) have shown that these two potentials lead to two distinctly different energetic barriers for the passage of methane through the eight-ring window in the cation-free form of zeolite A. Here, we compute the variation of the self-diffusivity (D) with loading (c) for these two types of potentials and show that this slight variation in the diameter changes the concentration dependence qualitatively: thus, D decreases monotonically with c for model A, while D increases and goes through a maximum before finally decreasing for model B. This effect and the surprising congruence of the diffusion coefficients for both models at high loadings is examined in detail at the molecular level. Simulations for different temperatures reveal the Arrhenius behaviour of the self-diffusion coefficient. The apparent activation energy is found to vary with the loading. We conclude that beside the cage-to-cage jumps, which are essential for the migration of the guest molecules, at high concentrations migration within the cage and guest guest interactions with other molecules become increasingly dominant influences on the diffusion coefficient and make the guest zeolite interaction less important for both model A and model B.

beta-Turn Analogues in Model a ss-Hybrid Peptides: Structural Characterization of Peptides Containing ss 2,2Ac6c and ss 3,3Ac6c Residues

The effect of gem-dialkyl substituents on the backbone conformations of beta-amino acid residues in peptides has been investigated by using four model peptides: Boc-Xxx-beta 2,2Ac6c(1-aminomethylcyclohexanecarboxylic acid)-NHMe (Xxx=Leu (1), Phe (2); Boc=tert-butyloxycarbonyl) and Boc-Xxx-beta 3,3Ac6c(1-aminocyclohexaneacetic acid)-NHMe (Xxx=Leu (3), Phe (4)). Tetrasubstituted carbon atoms restrict the ranges of stereochemically allowed conformations about flanking single bonds. The crystal structure of Boc-Leu-beta 2,2Ac6c-NHMe (1) established a C11 hydrogen-bonded turn in the a beta-hybrid sequence. The observed torsion angles (a(similar to-60 degrees, similar to-30 degrees), beta(similar to-90 degrees, similar to 60 degrees, similar to-90 degrees)) corresponded to a C11 helical turn, which was a backbone-expanded analogue of the type III beta turn in aa sequences. The crystal structure of the peptide Boc-Phe-beta 3,3Ac6c-NHMe (4) established a C11 hydrogen-bonded turn with distinctly different backbone torsion angles (a(similar to-60 degrees, similar to 120 degrees), beta(similar to 60 degrees, ?60 degrees, similar to-60 degrees)), which corresponded to a backbone-expanded analogue of the type II beta turn observed in aa sequences. In peptide 4, the two molecules in the asymmetric unit adopted backbone torsion angles of opposite signs. In one of the molecules, the Phe residue adopted an unfavorable backbone conformation, with the energetic penalty being offset by a favorable aromatic interaction between proximal molecules in the crystal. NMR spectroscopy studies provided evidence for the maintenance of folded structures in solution in these a beta-hybrid sequences.

An association and Wilson activity coefficient model for solubilities of aromatic solid pollutants in supercritical carbon dioxide

The solubilities of various solid pollutants in supercritical carbon dioxide were investigated. The intermolecular interactions play a significant role in determining the solubilities of solids in supercritical carbon dioxide. A new model equation was derived by using the concepts of association and activity coefficient model to correlate the solubilities of solids. The model equation combines the association and Wilson activity coefficient models and includes the interaction potentials between the molecules, which are useful in understanding the behavior of the solid solutes in SCCO2. The new model equation involves five adjustable parameters to correlate the solubilities of solids by incorporating the interactions between the molecules. The equation correlated 75 solid systems with an average AARD of around 9%, which was better than the correlations obtained from standard models such as Mendez Santiago-Teja (MT) model and association model. (C) 2012 Elsevier B.V. All rights reserved.

Enhanced J-protein interaction and compromised protein stability of mtHsp70 variants lead to mitochondrial dysfunction in Parkinsons disease

Parkinsons disease (PD) is the second most prevalent progressive neurological disorder commonly associated with impaired mitochondrial function in dopaminergic neurons. Although familial PD is multifactorial in nature, a recent genetic screen involving PD patients identified two mitochondrial Hsp70 variants (P509S and R126W) that are suggested in PD pathogenesis. However, molecular mechanisms underlying how mtHsp70 PD variants are centrally involved in PD progression is totally elusive. In this article, we provide mechanistic insights into the mitochondrial dysfunction associated with human mtHsp70 PD variants. Biochemically, the R126W variant showed severely compromised protein stability and was found highly susceptible to aggregation at physiological conditions. Strikingly, on the other hand, the P509S variant exhibits significantly enhanced interaction with J-protein cochaperones involved in folding and import machinery, thus altering the overall regulation of chaperone-mediated folding cycle and protein homeostasis. To assess the impact of mtHsp70 PD mutations at the cellular level, we developed yeast as a model system by making analogous mutations in Ssc1 ortholog. Interestingly, PD mutations in yeast (R103W and P486S) exhibit multiple in vivo phenotypes, which are associated with omitochondrial dysfunction', including compromised growth, impairment in protein translocation, reduced functional mitochondrial mass, mitochondrial DNA loss, respiratory incompetency and increased susceptibility to oxidative stress. In addition to that, R103W protein is prone to aggregate in vivo due to reduced stability, whereas P486S showed enhanced interaction with J-proteins, thus remarkably recapitulating the cellular defects that are observed in human PD variants. Taken together, our findings provide evidence in favor of direct involvement of mtHsp70 as a susceptibility factor in PD.

Causality analysis of groundwater dynamics based on a Vector Autoregressive model in the semi-arid basin of Gundal (South India)

Causal relationships existing between observed levels of groundwater in a semi-arid sub-basin of the Kabini River basin (Karnataka state, India) are investigated in this study. A Vector Auto Regressive model is used for this purpose. Its structure is built on an upstream/downstream interaction network based on observed hydro-physical properties. Exogenous climatic forcing is used as an input based on cumulated rainfall departure. Optimal models are obtained thanks to a trial approach and are used as a proxy of the dynamics to derive causal networks. It appears to be an interesting tool for analysing the causal relationships existing inside the basin. The causal network reveals 3 main regions: the Northeastern part of the Gundal basin is closely coupled to the outlet dynamics. The Northwestern part is mainly controlled by the climatic forcing and only marginally linked to the outlet dynamic. Finally, the upper part of the basin plays as a forcing rather than a coupling with the lower part of the basin allowing for a separate analysis of this local behaviour. The analysis also reveals differential time scales at work inside the basin when comparing upstream oriented with downstream oriented causalities. In the upper part of the basin, time delays are close to 2 months in the upward direction and lower than 1 month in the downward direction. These time scales are likely to be good indicators of the hydraulic response time of the basin which is a parameter usually difficult to estimate practically. This suggests that, at the sub-basin scale, intra-annual time scales would be more relevant scales for analysing or modelling tropical basin dynamics in hard rock (granitic and gneissic) aquifers ubiquitous in south India. (c) 2012 Elsevier B.V. All rights reserved.

Interaction of nucleic acids with carbon nanotubes and dendrimers

Nucleic acid interaction with nanoscale objects like carbon nanotubes (CNTs) and dendrimers is of fundamental interest because of their potential application in CNT separation, gene therapy and antisense therapy. Combining nucleic acids with CNTs and dendrimers also opens the door towards controllable self-assembly to generate various supra-molecular and nano-structures with desired morphologies. The interaction between these nanoscale objects also serve as a model system for studying DNA compaction, which is a fundamental process in chromatin organization. By using fully atomistic simulations, here we report various aspects of the interactions and binding modes of DNA and small interfering RNA (siRNA) with CNTs, graphene and dendrimers. Our results give a microscopic picture and mechanism of the adsorption of single- and double-strand DNA (ssDNA and dsDNA) on CNT and graphene. The nucleic acid-CNT interaction is dominated by the dispersive van der Waals (vdW) interaction. In contrast, the complexation of DNA (both ssDNA and dsDNA) and siRNA with various generations of poly-amido-amine (PAMAM) dendrimers is governed by electrostatic interactions. Our results reveal that both the DNA and siRNA form stable complex with the PAMAM dendrimer at a physiological pH when the dendrimer is positively charged due to the protonation of the primary amines. The size and binding energy of the complex increase with increase in dendrimer generation. We also give a summary of the current status in these fields and discuss future prospects.

Dendrimer building toolkit: Model building and characterization of various dendrimer architectures

We have developed a graphical user interface based dendrimer builder toolkit (DBT) which can be used to generate the dendrimer configuration of desired generation for various dendrimer architectures. The validation of structures generated by this tool was carried out by studying the structural properties of two well known classes of dendrimers: ethylenediamine cored poly(amidoamine) (PAMAM) dendrimer, diaminobutyl cored poly(propylene imine) (PPI) dendrimer. Using full atomistic molecular dynamics (MD) simulation we have calculated the radius of gyration, shape tensor and monomer density distribution for PAMAM and PPI dendrimer at neutral and high pH. A good agreement between the available simulation and experimental (small angle X-ray and neutron scattering; SAXS, SANS) results and calculated radius of gyration was observed. With this validation we have used DBT to build another new class of nitrogen cored poly(propyl ether imine) dendrimer and study it's structural features using all atomistic MD simulation. DBT is a versatile tool and can be easily used to generate other dendrimer structures with different chemistry and topology. The use of general amber force field to describe the intra-molecular interactions allows us to integrate this tool easily with the widely used molecular dynamics software AMBER. This makes our tool a very useful utility which can help to facilitate the study of dendrimer interaction with nucleic acids, protein and lipid bilayer for various biological applications. © 2012 Wiley Periodicals, Inc.

Dendrimer building toolkit: Model building and characterization of various dendrimer architectures

We have developed a graphical user interface based dendrimer builder toolkit (DBT) which can be used to generate the dendrimer configuration of desired generation for various dendrimer architectures. The validation of structures generated by this tool was carried out by studying the structural properties of two well known classes of dendrimers: ethylenediamine cored poly(amidoamine) (PAMAM) dendrimer, diaminobutyl cored poly(propylene imine) (PPI) dendrimer. Using full atomistic molecular dynamics (MD) simulation we have calculated the radius of gyration, shape tensor and monomer density distribution for PAMAM and PPI dendrimer at neutral and high pH. A good agreement between the available simulation and experimental (small angle X-ray and neutron scattering; SAXS, SANS) results and calculated radius of gyration was observed. With this validation we have used DBT to build another new class of nitrogen cored poly(propyl ether imine) dendrimer and study it's structural features using all atomistic MD simulation. DBT is a versatile tool and can be easily used to generate other dendrimer structures with different chemistry and topology. The use of general amber force field to describe the intra-molecular interactions allows us to integrate this tool easily with the widely used molecular dynamics software AMBER. This makes our tool a very useful utility which can help to facilitate the study of dendrimer interaction with nucleic acids, protein and lipid bilayer for various biological applications. (c) 2012 Wiley Periodicals, Inc.

Monsoon sensitivity to aerosol direct radiative forcing in the community atmosphere model

Aerosol forcing remains a dominant uncertainty in climate studies. The impact of aerosol direct radiative forcing on Indian monsoon is extremely complex and is strongly dependent on the model, aerosol distribution and characteristics specified in the model, modelling strategy employed as well as on spatial and temporal scales. The present study investigates (i) the aerosol direct radiative forcing impact on mean Indian summer monsoon when a combination of quasi-realistic mean annual cycles of scattering and absorbing aerosols derived from an aerosol transport model constrained with satellite observed Aerosol Optical Depth (AOD) is prescribed, (ii) the dominant feedback mechanism behind the simulated impact of all-aerosol direct radiative forcing on monsoon and (iii) the relative impacts of absorbing and scattering aerosols on mean Indian summer monsoon. We have used CAM3, an atmospheric GCM (AGCM) that has a comprehensive treatment of the aerosol-radiation interaction. This AGCM has been used to perform climate simulations with three different representations of aerosol direct radiative forcing due to the total, scattering aerosols and black carbon aerosols. We have also conducted experiments without any aerosol forcing. Aerosol direct impact due to scattering aerosols causes significant reduction in summer monsoon precipitation over India with a tendency for southward shift of Tropical Convergence Zones (TCZs) over the Indian region. Aerosol forcing reduces surface solar absorption over the primary rainbelt region of India and reduces the surface and lower tropospheric temperatures. Concurrent warming of the lower atmosphere over the warm oceanic region in the south reduces the land-ocean temperature contrast and weakens the monsoon overturning circulation and the advection of moisture into the landmass. This increases atmospheric convective stability, and decreases convection, clouds, precipitation and associated latent heat release. Our analysis reveals a defining negative moisture-advection feedback that acts as an internal damping mechanism spinning down the regional hydrological cycle and leading to significant circulation changes in response to external radiative forcing perturbations. When total aerosol loading (both absorbing and scattering aerosols) is prescribed, dust and black carbon aerosols are found to cause significant atmospheric heating over the monsoon region but the aerosol-induced weakening of meridional lower tropospheric temperature gradient (leading to weaker summer monsoon rainfall) more than offsets the increase in summer-time rainfall resulting from the atmospheric heating effect of absorbing aerosols, leading to a net decrease of summer monsoon rainfall. Further, we have carried out climate simulations with globally constant AODs and also with the constant AODs over the extended Indian region replaced by realistic AODs. Regional aerosol radiative forcing perturbations over the Indian region is found to have impact not only over the region of loading but over remote tropical regions as well. This warrants the need to prescribe realistic aerosol properties in strategic regions such as India in order to accurately assess the aerosol impact.

Salmonella enterica serovars Typhimurium and Typhi as model organisms Revealing paradigm of host-pathogen interactions

The lifestyle of intracellular pathogens has always questioned the skill of a microbiologist in the context of finding the permanent cure to the diseases caused by them. The best tool utilized by these pathogens is their ability to reside inside the host cell, which enables them to easily bypass the humoral immunity of the host, such as the complement system. They further escape from the intracellular immunity, such as lysosome and inflammasome, mostly by forming a protective vacuole-bound niche derived from the host itself. Some of the most dreadful diseases are caused by these vacuolar pathogens, for example, tuberculosis by Mycobacterium or typhoid fever by Salmonella. To deal with such successful pathogens therapeutically, the knowledge of a host-pathogen interaction system becomes primarily essential, which further depends on the use of a model system. A well characterized pathogen, namely Salmonella, suits the role of a model for this purpose, which can infect a wide array of hosts causing a variety of diseases. This review focuses on various such aspects of research on Salmonella which are useful for studying the pathogenesis of other intracellular pathogens.

A Comparative Kinetic and Thermodynamic Perspective of the sigma-Competition Model in Escherichia coli

Transcription is the most fundamental step in gene expression in any living organism. Various environmental cues help in the maturation of core RNA polymerase (RNAP; alpha(2)beta beta'omega) with different sigma-factors, leading to the directed recruitment of RNAP to different promoter DNA sequences. Thus it is essential to determine the sigma-factors that affect the preferential partitioning of core RNAP among various a-actors, and the role of sigma-switching in transcriptional gene regulation. Further, the macromolecular assembly of holo RNAP takes place in an extremely crowded environment within a cell, and thus far the kinetics and thermodynamics of this molecular recognition process have not been well addressed. In this study we used a site-directed bioaffinity immobilization method to evaluate the relative binding affinities of three different Escherichia coli sigma-factors to the same core RNAP with variations in temperature and ionic strength while emulating the crowded cellular milieu. Our data indicate that the interaction of core RNAP-sigma is susceptible to changes in external stimuli such as osmolytic and thermal stress, and the degree of susceptibility varies among different sigma-factors. This allows for a reversible sigma-switching from housekeeping factors to alternate sigma-factors when the organism senses a change in its physiological conditions.

Recapitulating tumour microenvironment in chitosan-gelatin three-dimensional scaffolds: an improved in vitro tumour model

Owing to the reduced co-relationship between conventional flat Petri dish culture (two-dimensional) and the tumour microenvironment, there has been a shift towards three-dimensional culture systems that show an improved analogy to the same. In this work, an extracellular matrix (ECM)-mimicking three-dimensional scaffold based on chitosan and gelatin was fabricated and explored for its potential as a tumour model for lung cancer. It was demonstrated that the chitosan-gelatin (CG) scaffolds supported the formation of tumoroids that were similar to tumours grown in vivo for factors involved in tumour-cell-ECM interaction, invasion and metastasis, and response to anti-cancer drugs. On the other hand, the two-dimensional Petri dish surfaces did not demonstrate gene-expression profiles similar to tumours grown in vivo. Further, the three-dimensional CG scaffolds supported the formation of tumoroids, using other types of cancer cells such as breast, cervix and bone, indicating a possible wider potential for in vitro tumoroid generation. Overall, the results demonstrated that CG scaffolds can be an improved in vitro tool to study cancer progression and drug screening for solid tumours.