Spiral-Wave Dynamics in a Mathematical Model of Human Ventricular Tissue with Myocytes and Fibroblasts

Cardiac fibroblasts, when coupled functionally with myocytes, can modulate the electrophysiological properties of cardiac tissue. We present systematic numerical studies of such modulation of electrophysiological properties in mathematical models for (a) single myocyte-fibroblast (MF) units and (b) two-dimensional (2D) arrays of such units; our models build on earlier ones and allow for zero-, one-, and two-sided MF couplings. Our studies of MF units elucidate the dependence of the action-potential (AP) morphology on parameters such as E-f, the fibroblast resting-membrane potential, the fibroblast conductance G(f), and the MF gap-junctional coupling G(gap). Furthermore, we find that our MF composite can show autorhythmic and oscillatory behaviors in addition to an excitable response. Our 2D studies use (a) both homogeneous and inhomogeneous distributions of fibroblasts, (b) various ranges for parameters such as G(gap), G(f), and E-f, and (c) intercellular couplings that can be zero-sided, one-sided, and two-sided connections of fibroblasts with myocytes. We show, in particular, that the plane-wave conduction velocity CV decreases as a function of G(gap), for zero-sided and one-sided couplings; however, for two-sided coupling, CV decreases initially and then increases as a function of G(gap), and, eventually, we observe that conduction failure occurs for low values of G(gap). In our homogeneous studies, we find that the rotation speed and stability of a spiral wave can be controlled either by controlling G(gap) or E-f. Our studies with fibroblast inhomogeneities show that a spiral wave can get anchored to a local fibroblast inhomogeneity. We also study the efficacy of a low-amplitude control scheme, which has been suggested for the control of spiral-wave turbulence in mathematical models for cardiac tissue, in our MF model both with and without heterogeneities.

Monte Carlo simulation of light transport in tissue for optimizing light delivery in photoacoustic imaging of the sentinel lymph node

Noninvasive or minimally invasive identification of sentinel lymph node (SLN) is essential to reduce the surgical effects of SLN biopsy. Photoacoustic (PA) imaging of SLN in animal models has shown its promise for clinical use in the future. Here, we present a Monte Carlo simulation for light transport in the SLN for various light delivery configurations with a clinical ultrasound probe. Our simulation assumes a realistic tissue layer model and also can handle the transmission/reflectance at SLN-tissue boundary due to the mismatch of refractive index. Various light incidence angles show that for deeply situated SLNs the maximum absorption of light in the SLN is for normal incidence. We also show that if a part of the diffused reflected photons is reflected back into the skin using a reflector, the absorption of light in the SLN can be increased significantly to enhance the PA signal. (C) 2013 Society of Photo-Optical Instrumentation Engineers (SPIE)

Effect of SST Variation on ITCZ in APE Simulations

The effect of meridional variation of sea surface temperature (SST) on tropical atmospheric circulation is analyzed using Aqua-planet Experiment (APE) simulations. The meridional SST gradient around the narrow SST peak in CONTROL simulation favours a strong and single equatorial Intertropical Convergence Zone (ITCZ, defined by the maximum of zonally averaged total precipitation) in all APE models. In contrast, flat equatorial SST peak (FLAT simulation) favours split/double ITCZs flanking the SST maximum, in the majority of the APE models. Although there is reasonable agreement for SST sensitivity of ITCZ among the APE models in CONTROL, there exists disparity among them in FLAT case. Similarly, while the total and convective precipitation responses are consistent among the models, the large-scale precipitation response shows considerable inter-model variations in FLAT case. The APE intercomparison indicates that the occurrence and positioning of the ITCZ are primarily related to boundary layer moisture convergence as a response to the meridional variation of SST. Furthermore, the meridional gradient of tropospheric temperature is found to be an important factor that can influence the positioning of ITCZ. FLAT SST distribution is found to be similar to the observed distribution over the Indian region during summer season. Models that yield double ITCZs in this case simulate an easterly jet over the equatorial region (similar to 15 degrees equatorward of the ITCZ). This is analogous to the Tropical Easterly Jet (TEJ), which is a unique feature observed over the Indian region during summer monsoon season, with its core at 12 degrees N, equatorward of the seasonal convergence zone centered along 25 degrees N. In these models, positive meridional temperature gradient and the associated easterly shear in the atmosphere strengthened by moisture convergence penetrate up to the upper troposphere, with which TEJ is in thermal wind balance.

Structure-function hierarchies and von Karman-Howarth relations for turbulence in magnetohydrodynamical equations

We generalize the method of A. M. Polyakov, Phys. Rev. E 52, 6183 (1995)] for obtaining structure-function relations in turbulence in the stochastically forced Burgers equation, to develop structure-function hierarchies for turbulence in three models for magnetohydrodynamics (MHD). These are the Burgers analogs of MHD in one dimension Eur. Phys. J.B 9, 725 (1999)], and in three dimensions (3DMHD and 3D Hall MHD). Our study provides a convenient and unified scheme for the development of structure-function hierarchies for turbulence in a variety of coupled hydrodynamical equations. For turbulence in the three sets of MHD equations mentioned above, we obtain exact relations for third-order structure functions and their derivatives; these expressions are the analogs of the von Karman-Howarth relations for fluid turbulence. We compare our work with earlier studies of such relations in 3DMHD and 3D Hall MHD.

Critical cysteines in Akt1 regulate its activity and proteasomal degradation: implications for neurodegenerative diseases

Impaired Akt1 signaling is observed in neurodegenerative diseases, including Parkinson's disease (PD). In PD models oxidative modification of Akt1 leads to its dephosphorylation and consequent loss of its kinase activity. To explore the underlying mechanism we exposed Neuro2A cells to cadmium, a pan inhibitor of protein thiol disulfide oxidoreductases, including glutaredoxin 1 (Grx1), or downregulated Grx1, which led to dephosphorylation of Akt1, loss of its kinase activity, and also decreased Akt1 protein levels. Mutation of cysteines to serines at 296 and 310 in Akt1 did not affect its basal kinase activity but abolished cadmium- and Grx1 downregulation-induced reduction in Akt1 kinase activity, indicating their critical role in redox modulation of Akt1 function and turnover. Cadmium-induced decrease in phosphorylated Akt1 correlated with increased association of wild-type (WT) Akt1 with PP2A, which was absent in the C296-310S Akt1 mutant and was also abolished by N-acetylcysteine treatment. Further, increased proteasomal degradation of Akt1 by cadmium was not seen in the C296-310S Akt1 mutant, indicating that oxidation of cysteine residues facilitates degradation of WT Akt1. Moreover, preventing oxidative modification of Akt1 cysteines 296 and 310 by mutating them to serines increased the cell survival effects of Akt1. Thus, in neurodegenerative states such as PD, maintaining the thiol status of cysteines 296 and 310 in Akt1 would be critical for Akt1 kinase activity and for preventing its degradation by proteasomes. Preventing downregulation of Akt signaling not only has long-range consequences for cell survival but could also affect the multiple roles that Ala plays, including in the Akt-mTOR signaling cascade. (C) 2014 Elsevier Inc. All rights reserved.

Current trends in modeling hostpathogen interactions

The rapid emergence of infectious diseases calls for immediate attention to determine practical solutions for intervention strategies. To this end, it becomes necessary to obtain a holistic view of the complex hostpathogen interactome. Advances in omics and related technology have resulted in massive generation of data for the interacting systems at unprecedented levels of detail. Systems-level studies with the aid of mathematical tools contribute to a deeper understanding of biological systems, where intuitive reasoning alone does not suffice. In this review, we discuss different aspects of hostpathogen interactions (HPIs) and the available data resources and tools used to study them. We discuss in detail models of HPIs at various levels of abstraction, along with their applications and limitations. We also enlist a few case studies, which incorporate different modeling approaches, providing significant insights into disease. (c) 2013 Wiley Periodicals, Inc.

Modelling staphylococcal pneumonia in a human 3D lung tissue model system delineates toxin-mediated pathology

Staphylococcus aureus necrotizing pneumonia is recognized as a toxin-mediated disease, yet the tissue-destructive events remain elusive, partly as a result of lack of mechanistic studies in human lung tissue. In this study, a three-dimensional (3D) tissue model composed of human lung epithelial cells and fibroblasts was used to delineate the role of specific staphylococcal exotoxins in tissue pathology associated with severe pneumonia. To this end, the models were exposed to the mixture of exotoxins produced by S. aureus strains isolated from patients with varying severity of lung infection, namely necrotizing pneumonia or lung empyema, or to purified toxins. The necrotizing pneumonia strains secreted high levels of alpha-toxin and Panton-Valentine leukocidin (PVL), and triggered high cytotoxicity, inflammation, necrosis and loss of E-cadherin from the lung epithelium. In contrast, the lung empyema strain produced moderate levels of PVL, but negligible amounts of alpha-toxin, and triggered limited tissue damage. alpha-toxin had a direct damaging effect on the epithelium, as verified using toxin-deficient mutants and pure alpha-toxin. Moreover, PVL contributed to pathology through the lysis of neutrophils. A combination of alpha-toxin and PVL resulted in the most severe epithelial injury. In addition, toxin-induced release of pro-inflammatory mediators from lung tissue models resulted in enhanced neutrophil migration. Using a collection of 31 strains from patients with staphylococcal pneumonia revealed that strains producing high levels of alpha-toxin and PVL were cytotoxic and associated with fatal outcome. Also, the strains that produced the highest toxin levels induced significantly greater epithelial disruption. Of importance, toxin-mediated lung epithelium destruction could be inhibited by polyspecific intravenous immunoglobulin containing antibodies against alpha-toxin and PVL. This study introduces a novel model system for study of staphylococcal pneumonia in a human setting. The results reveal that the combination and levels of alpha-toxin and PVL correlate with tissue pathology and clinical outcome associated with pneumonia.

Many-Body Localization in the Presence of a Single-Particle Mobility Edge

In one dimension, noninteracting particles can undergo a localization-delocalization transition in a quasiperiodic potential. Recent studies have suggested that this transition transforms into a many-body localization (MBL) transition upon the introduction of interactions. It has also been shown that mobility edges can appear in the single particle spectrum for certain types of quasiperiodic potentials. Here, we investigate the effect of interactions in two models with such mobility edges. Employing the technique of exact diagonalization for finite-sized systems, we calculate the level spacing distribution, time evolution of entanglement entropy, optical conductivity, and return probability to detect MBL. We find that MBL does indeed occur in one of the two models we study, but the entanglement appears to grow faster than logarithmically with time unlike in other MBL systems.

Inter-comparison and performance evaluation of chemistry transport models over Indian region

Aerosol loading over the South Asian region has the potential to affect the monsoon rainfall, Himalayan glaciers and regional air-quality, with implications for the billions in this region. While field campaigns and network observations provide primary data, they tend to be location/season specific. Numerical models are useful to regionalize such location-specific data. Studies have shown that numerical models underestimate the aerosol scenario over the Indian region, mainly due to shortcomings related to meteorology and the emission inventories used. In this context, we have evaluated the performance of two such chemistry-transport models: WRF-Chem and SPRINTARS over an India-centric domain. The models differ in many aspects including physical domain, horizontal resolution, meteorological forcing and so on etc. Despite these differences, both the models simulated similar spatial patterns of Black Carbon (BC) mass concentration, (with a spatial correlation of 0.9 with each other), and a reasonable estimates of its concentration, though both of them under-estimated vis-a-vis the observations. While the emissions are lower (higher) in SPRINTARS (WRF-Chem), overestimation of wind parameters in WRF-Chem caused the concentration to be similar in both models. Additionally, we quantified the under-estimations of anthropogenic BC emissions in the inventories used these two models and three other widely used emission inventories. Our analysis indicates that all these emission inventories underestimate the emissions of BC over India by a factor that ranges from 1.5 to 2.9. We have also studied the model simulations of aerosol optical depth over the Indian region. The models differ significantly in simulations of AOD, with WRF-Chem having a better agreement with satellite observations of AOD as far as the spatial pattern is concerned. It is important to note that in addition to BC, dust can also contribute significantly to AOD. The models differ in simulations of the spatial pattern of mineral dust over the Indian region. We find that both meteorological forcing and emission formulation contribute to these differences. Since AOD is column integrated parameter, description of vertical profiles in both models, especially since elevated aerosol layers are often observed over Indian region, could be also a contributing factor. Additionally, differences in the prescription of the optical properties of BC between the models appear to affect the AOD simulations. We also compared simulation of sea-salt concentration in the two models and found that WRF-Chem underestimated its concentration vis-a-vis SPRINTARS. The differences in near-surface oceanic wind speeds appear to be the main source of this difference. In-spite of these differences, we note that there are similarities in their simulation of spatial patterns of various aerosol species (with each other and with observations) and hence models could be valuable tools for aerosol-related studies over the Indian region. Better estimation of emission inventories could improve aerosol-related simulations. (C) 2015 Elsevier Ltd. All rights reserved.

A dilation-driven vortex flow in sheared granular materials explains a rheometric anomaly

Granular flows occur widely in nature and industry, yet a continuum description that captures their important features is yet not at hand. Recent experiments on granular materials sheared in a cylindrical Couette device revealed a puzzling anomaly, wherein all components of the stress rise nearly exponentially with depth. Here we show, using particle dynamics simulations and imaging experiments, that the stress anomaly arises from a remarkable vortex flow. For the entire range of fill heights explored, we observe a single toroidal vortex that spans the entire Couette cell and whose sense is opposite to the uppermost Taylor vortex in a fluid. We show that the vortex is driven by a combination of shear-induced dilation, a phenomenon that has no analogue in fluids, and gravity flow. Dilatancy is an important feature of granular mechanics, but not adequately incorporated in existing models.