Viral Kinetics Suggests a Reconciliation of the Disparate Observations of the Modulation of Claudin-1 Expression on Cells Exposed to Hepatitis C Virus

The tight junction protein claudin-1 (CLDN1) is necessary for hepatitis C virus (HCV) entry into target cells. Recent studies have made disparate observations of the modulation of the expression of CLDN1 on cells following infection by HCV. In one study, the mean CLDN1 expression on cells exposed to HCV declined, whereas in another study HCV infected cells showed increased CLDN1 expression compared to uninfected cells. Consequently, the role of HCV in modulating CLDN1 expression, and hence the frequency of cellular superinfection, remains unclear. Here, we present a possible reconciliation of these disparate observations. We hypothesized that viral kinetics and not necessarily HCV-induced receptor modulation underlies these disparate observations. To test this hypothesis, we constructed a mathematical model of viral kinetics in vitro that mimicked the above experiments. Model predictions provided good fits to the observed evolution of the distribution of CLDN1 expression on cells following exposure to HCV. Cells with higher CLDN1 expression were preferentially infected and outgrown by cells with lower CLDN1 expression, resulting in a decline of the mean CLDN1 expression with time. At the same time, because the susceptibility of cells to infection increased with CLDN1 expression, infected cells tended to have higher CLDN1 expression on average than uninfected cells. Our study thus presents an explanation of the disparate observations of CLDN1 expression following HCV infection and points to the importance of considering viral kinetics in future studies of receptor expression on cells exposed to HCV.

Extracting Three Dimensional Surface Model of Human Kidney from the Visible Human Data Set using Free Software

Three dimensional digital model of a representative human kidney is needed for a surgical simulator that is capable of simulating a laparoscopic surgery involving kidney. Buying a three dimensional computer model of a representative human kidney, or reconstructing a human kidney from an image sequence using commercial software, both involve (sometimes significant amount of) money. In this paper, author has shown that one can obtain a three dimensional surface model of human kidney by making use of images from the Visible Human Data Set and a few free software packages (ImageJ, ITK-SNAP, and MeshLab in particular). Images from the Visible Human Data Set, and the software packages used here, both do not cost anything. Hence, the practice of extracting the geometry of a representative human kidney for free, as illustrated in the present work, could be a free alternative to the use of expensive commercial software or to the purchase of a digital model.

A finite population model of mlecular evolution: theory and computation

This article is concerned with the evolution of haploid organisms that reproduce asexually. In a seminal piece of work, Eigen and coauthors proposed the quasispecies model in an attempt to understand such an evolutionary process. Their work has impacted antiviral treatment and vaccine design strategies. Yet, predictions of the quasispecies model are at best viewed as a guideline, primarily because it assumes an infinite population size, whereas realistic population sizes can be quite small. In this paper we consider a population genetics-based model aimed at understanding the evolution of such organisms with finite population sizes and present a rigorous study of the convergence and computational issues that arise therein. Our first result is structural and shows that, at any time during the evolution, as the population size tends to infinity, the distribution of genomes predicted by our model converges to that predicted by the quasispecies model. This justifies the continued use of the quasispecies model to derive guidelines for intervention. While the stationary state in the quasispecies model is readily obtained, due to the explosion of the state space in our model, exact computations are prohibitive. Our second set of results are computational in nature and address this issue. We derive conditions on the parameters of evolution under which our stochastic model mixes rapidly. Further, for a class of widely used fitness landscapes we give a fast deterministic algorithm which computes the stationary distribution of our model. These computational tools are expected to serve as a framework for the modeling of strategies for the deployment of mutagenic drugs.

Improvement in Bone Properties by Using Risedronate Adsorbed Hydroxyapatite Novel Nanoparticle Based Formulation in a Rat Model of Osteoporosis

A superior drug formulation capable of achieving efficient osteogenesis is in imperative demand for the treatment of osteoporosis. In the present study we investigated the potential of using novel risedronate-hydroxyapatite (HA) nanoparticle based formulation in an animal model of established osteoporosis. Nanoparticles of HA loaded with risedronate (NHLR) of various sizes (80-130 nm) were generated for bone targeted drug delivery. Three months after ovariectomy, 36 ovariectomized (OVX) rats were divided into 6 equal groups and treated with various doses of NHLR (500,350 and 250 mu g/kg intravenous single dose) and sodium risedronate (500 mu g/kg, intravenous single dose). Untreated OVX and sham OVX served as controls. One month after drug administration, the left tibia and femur were tested for bone mechanical properties and histology, respectively. In the right femur, bone density was measured by method based on Archimedes principle and bone porosity analyses were performed using X-ray imaging. NHLR (250 mu g/kg) showed a significant increase in bone density and reduced bone porosity when compared with OVX control. Moreover, NHLR (250 mu g/kg) significantly increased bone mechanical properties and bone quality when compared with OVX control. The results strongly suggest that the NHLR, which is a novel nanoparticle based formulation, has a therapeutic advantage over risedronate sodium monotherapy for the treatment of osteoporosis in a rat model of postmenopausal osteoporosis.

A finite-rate-of-innovation signal sampling perspective of the source-filter model of speech production

We consider the speech production mechanism and the asso- ciated linear source-filter model. For voiced speech sounds in particular, the source/glottal excitation is modeled as a stream of impulses and the filter as a cascade of second-order resonators. We show that the process of sampling speech signals can be modeled as filtering a stream of Dirac impulses (a model for the excitation) with a kernel function (the vocal tract response),and then sampling uniformly. We show that the problem of esti- mating the excitation is equivalent to the problem of recovering a stream of Dirac impulses from samples of a filtered version. We present associated algorithms based on the annihilating filter and also make a comparison with the classical linear prediction technique, which is well known in speech analysis. Results on synthesized as well as natural speech data are presented.

Tethering preferences of domain families co-occurring in multi-domain proteins

Genomic data of several organisms have revealed the presence of a vast repertoire of multi-domain proteins. The role played by individual domains in a multi-domain protein has a profound influence on the overall function of the protein. In the present analysis an attempt has been made to better understand the tethering preferences of domain families that occur in multi-domain proteins. The analysis has been carried out on an exhaustive dataset of 2 961 898 sequences of proteins from 930 organisms, where 741 274 proteins are comprised of at least two domain families. For every domain family, the number of other domain families with which it co-occurs within a protein in this dataset has been enumerated and is referred to as the tethering number of the domain family. It was found that, in the general dataset, the AAA ATPase family and the family of Ser/Thr kinases have the highest tethering numbers of 450 and 444 respectively. Further analysis reveals significant correlation between the number of members in a family and its tethering number. Positive correlation was also observed for the extent of a sequence and functional diversity within a family and the tethering numbers of domain families. Domain families that are present ubiquitously in diverse organisms tend to have large tethering numbers, while organism/kingdom-specific families have low tethering numbers. Thus, the analysis uncovers how domain families recombine and evolve to give rise to multi-domain proteins.

Stability analysis of aperiodic messages scheduled in the dynamic segment of flex ray protocol

In a cyber physical system like vehicles number of signals to be communicated in a network system has an increasing trend. More and more mechanical and hydraulic parts are replaced by electronic control units and infotainment and multimedia applications has increased in vehicles. Safety critical hard real time messages and aperiodic messages communicated between electronic control units have been increased in recent times. Flexray is a high bandwidth protocol consisting of static segment for supporting hard real time messages and a dynamic segment for transmitting soft and non real time messages. In this paper, a method to obtain the stability region for the random arrival of messages in each electronic control units which is scheduled in the dynamic segment of Flexray protocol is presented. Number of mini slots available in the dynamic segment of Flexray restricts the arrival rate of tasks to the micro controllers or the number of micro controllers connected to the Flexray bus. Stability region of mathematical model of the system is compared with the Flexray protocol simulation results.

Stabilization of diborane(4) by transition metal fragments and a novel metal to pi Dewar-Chatt-Duncanson model of back donation

The feasibility of using transition metal fragments to stabilize B2H4 in planar configuration by donating 2 electrons to the boron moiety is investigated. Building upon the existing theoretical and experimental data and aided by the isolobal analogy, the model transition metal complexes Cr(CO)(4)B2H4 (6), Mn(CO)-CpB2H4 (7), Fe(CO)(3)B2H4 (8) and CoCpB2H4 (9) are chosen to illustrate this unique bonding feature bond strengthening with pi-back donation. Other possible types of complexes with B2H4 and the metal fragment are also explored and the energies are compared. One of the low energy isomers wherein the planar B2H4 interacts with the metal fragment in an in-plane fashion represents a unique case study for the Dewar-Chatt-Duncanson model. In this complex the back-donation from the metal fills the p bonding orbital between the two boron atoms thus forming a B=B double bond.

EFFECT OF CORE-SHELL STRUCTURE OF HYDROGEL BEADS ON THE THRESHOLD CONCENTRATION OF WATER FOR SWELLING AND ITS PH SENSITIVITY

In this paper we investigate the effect of core-shell structure of Sodium Alginate based hydrogel beads and their size on certain activation threshold concentration of water for applications in swelling and pH sensing. This type of hydrogel experiences diffusive pressure due to transport of certain free charges across its interface with a solvent or electrolyte. This process is essentially a dynamic equilibrium of the electric force field, stress in the polymeric network with cage like structure and molecular diffusion including phase transformation due to pressure imbalance between the hydrogel and its surroundings. The effect of pH of the solvant on the swelling rate of these beads has been studied experimentally. A mathematical model of the swelling process has been developed by considering Nernst-Planck equation representing the migration of mobile ions and Er ions, Poisson equation representing the equilibrium of the electric field and mechanical field equation representing swelling of the gel. An attempt has been made to predict the experimentally observed phenomena using these numerical simulations. It is observed experimentally that certain minimum concentration called activation threshold concentration of the water molecules must be present in the hydrogel in order to activate the swelling process. For the required activation threshold concentration of water in the beads, the pH induced change in the rate of swelling is also investigated. This effect is analyzed for various different core-shell structures of the beads.

Interfacial growth as a model of tube-width heterogeneities in concentrated solutions of stiff polymers

Recent experimental measurements of the distribution P(w) of transverse chain fluctuations w in concentrated solutions of F-actin filaments B. Wang, J Guan, S. M. Anthony, S. C. Bae, K. S. Schweizer, and S. Granick, Phys. Rev. Lett. 104, 118301 (2010); J. Glaser, D. Chakraborty, K. Kroy, I. Lauter, M. Degawa, N. Kirchgessner, B. Hoffmann, R. Merkel, and M. Giesen, Phys. Rev. Lett. 105, 037801 (2010)] are shown to be well-fit to an expression derived from a model of the conformations of a single harmonically confined weakly bendable rod. The calculation of P(w) is carried out essentially exactly within a path integral approach that was originally applied to the study of one-dimensional randomly growing interfaces. Our results are generally as successful in reproducing experimental trends as earlier approximate results obtained from more elaborate many-chain treatments of the confining tube potential.

A model of coupled thermal, mechanical, and electrostatic field effects in III-N thin film heterostructures

A one-dimensional coupled multi-physics based model has been developed to accurately compute the effects of electrostatic, mechanical, and thermal field interactions on the electronic energy band structure in group III-nitrides thin film heterostructures. Earlier models reported in published literature assumes electro-mechanical field with uniform temperature thus neglecting self-heating. Also, the effects of diffused interface on the energy band structure were not studied. We include these effects in a self-consistent manner wherein the transport equation is introduced along with the electro-mechanical models, and the lattice structural variation as observed in experiments are introduced at the interface. Due to these effects, the electrostatic potential distribution in the heterostructure is altered. The electron and hole ground state energies decrease by 5% and 9%, respectively, at a relative temperature of 700 K, when compared with the results obtained from the previously reported electro-mechanical model assuming constant and uniform temperature distribution. A diffused interface decreases the ground state energy of electrons and holes by about 11% and 9%, respectively, at a relative temperature of 700 K when compared with the predictions based on uniform temperature based electro-mechanical model. (C) 2013 AIP Publishing LLC.

Interfacial growth as a model of tube-width heterogeneities in concentrated solutions of stiff polymers

Recent experimental measurements of the distribution P(w) of transverse chain fluctuations w in concentrated solutions of F-actin filaments B. Wang, J Guan, S. M. Anthony, S. C. Bae, K. S. Schweizer, and S. Granick, Phys. Rev. Lett. 104, 118301 (2010); J. Glaser, D. Chakraborty, K. Kroy, I. Lauter, M. Degawa, N. Kirchgessner, B. Hoffmann, R. Merkel, and M. Giesen, Phys. Rev. Lett. 105, 037801 (2010)] are shown to be well-fit to an expression derived from a model of the conformations of a single harmonically confined weakly bendable rod. The calculation of P(w) is carried out essentially exactly within a path integral approach that was originally applied to the study of one-dimensional randomly growing interfaces. Our results are generally as successful in reproducing experimental trends as earlier approximate results obtained from more elaborate many-chain treatments of the confining tube potential. (C) 2013 AIP Publishing LLC.

An enthalpy-based model of dendritic growth in a convecting binary alloy melt

Purpose-In the present work, a numerical method, based on the well established enthalpy technique, is developed to simulate the growth of binary alloy equiaxed dendrites in presence of melt convection. The paper aims to discuss these issues. Design/methodology/approach-The principle of volume-averaging is used to formulate the governing equations (mass, momentum, energy and species conservation) which are solved using a coupled explicit-implicit method. The velocity and pressure fields are obtained using a fully implicit finite volume approach whereas the energy and species conservation equations are solved explicitly to obtain the enthalpy and solute concentration fields. As a model problem, simulation of the growth of a single crystal in a two-dimensional cavity filled with an undercooled melt is performed. Findings-Comparison of the simulation results with available solutions obtained using level set method and the phase field method shows good agreement. The effects of melt flow on dendrite growth rate and solute distribution along the solid-liquid interface are studied. A faster growth rate of the upstream dendrite arm in case of binary alloys is observed, which can be attributed to the enhanced heat transfer due to convection as well as lower solute pile-up at the solid-liquid interface. Subsequently, the influence of thermal and solutal Peclet number and undercooling on the dendrite tip velocity is investigated. Originality/value-As the present enthalpy based microscopic solidification model with melt convection is based on a framework similar to popularly used enthalpy models at the macroscopic scale, it lays the foundation to develop effective multiscale solidification.

Understanding the role of domain-domain linkers in the spatial orientation of domains in multi-domain proteins

Inter-domain linkers (IDLs)' bridge flanking domains and support inter-domain communication in multi-domain proteins. Their sequence and conformational preferences enable them to carry out varied functions. They also provide sufficient flexibility to facilitate domain motions and, in conjunction with the interacting interfaces, they also regulate the inter-domain geometry (IDG). In spite of the basic intuitive understanding of the inter-domain orientations with respect to linker conformations and interfaces, we still do not entirely understand the precise relationship among the three. We show that IDG is evolutionarily well conserved and is constrained by the domain-domain interface interactions. The IDLs modulate the interactions by varying their lengths, conformations and local structure, thereby affecting the overall IDG. Results of our analysis provide guidelines in modelling of multi-domain proteins from the tertiary structures of constituent domain components.

Reaction dynamics under confinement: an exact path integral treatment of a two-stage model of stochastic gene expression

Gene expression in living systems is inherently stochastic, and tends to produce varying numbers of proteins over repeated cycles of transcription and translation. In this paper, an expression is derived for the steady-state protein number distribution starting from a two-stage kinetic model of the gene expression process involving p proteins and r mRNAs. The derivation is based on an exact path integral evaluation of the joint distribution, P(p, r, t), of p and r at time t, which can be expressed in terms of the coupled Langevin equations for p and r that represent the two-stage model in continuum form. The steady-state distribution of p alone, P(p), is obtained from P(p, r, t) (a bivariate Gaussian) by integrating out the r degrees of freedom and taking the limit t -> infinity. P(p) is found to be proportional to the product of a Gaussian and a complementary error function. It provides a generally satisfactory fit to simulation data on the same two-stage process when the translational efficiency (a measure of intrinsic noise levels in the system) is relatively low; it is less successful as a model of the data when the translational efficiency (and noise levels) are high.