Ribavirin-Induced Anemia in Hepatitis C Virus Patients Undergoing Combination Therapy

The current standard of care for hepatitis C virus (HCV) infection - combination therapy with pegylated interferon and ribavirin - elicits sustained responses in only similar to 50% of the patients treated. No alternatives exist for patients who do not respond to combination therapy. Addition of ribavirin substantially improves response rates to interferon and lowers relapse rates following the cessation of therapy, suggesting that increasing ribavirin exposure may further improve treatment response. A key limitation, however, is the toxic side-effect of ribavirin, hemolytic anemia, which often necessitates a reduction of ribavirin dosage and compromises treatment response. Maximizing treatment response thus requires striking a balance between the antiviral and hemolytic activities of ribavirin. Current models of viral kinetics describe the enhancement of treatment response due to ribavirin. Ribavirin-induced anemia, however, remains poorly understood and precludes rational optimization of combination therapy. Here, we develop a new mathematical model of the population dynamics of erythrocytes that quantitatively describes ribavirin-induced anemia in HCV patients. Based on the assumption that ribavirin accumulation decreases erythrocyte lifespan in a dose-dependent manner, model predictions capture several independent experimental observations of the accumulation of ribavirin in erythrocytes and the resulting decline of hemoglobin in HCV patients undergoing combination therapy, estimate the reduced erythrocyte lifespan during therapy, and describe inter-patient variations in the severity of ribavirin-induced anemia. Further, model predictions estimate the threshold ribavirin exposure beyond which anemia becomes intolerable and suggest guidelines for the usage of growth hormones, such as erythropoietin, that stimulate erythrocyte production and avert the reduction of ribavirin dosage, thereby improving treatment response. Our model thus facilitates, in conjunction with models of viral kinetics, the rational identification of treatment protocols that maximize treatment response while curtailing side effects.

Heuristic algorithms for scheduling of a batch processor in automobile gear manufacturing

In this paper we address a scheduling problem for minimising total weighted tardiness. The motivation for the paper comes from the automobile gear manufacturing process. We consider the bottleneck operation of heat treatment stage of gear manufacturing. Real life scenarios like unequal release times, incompatible job families, non-identical job sizes and allowance for job splitting have been considered. A mathematical model taking into account dynamic starting conditions has been developed. Due to the NP-hard nature of the problem, a few heuristic algorithms have been proposed. The performance of the proposed heuristic algorithms is evaluated: (a) in comparison with optimal solution for small size problem instances, and (b) in comparison with `estimated optimal solution' for large size problem instances. Extensive computational analyses reveal that the proposed heuristic algorithms are capable of consistently obtaining near-optimal solutions (that is, statistically estimated one) in very reasonable computational time.

Effective drug dosage design for treatment of infectious diseases using dynamic inversion

A generic nonlinear mathematical model describing the human immunological dynamics is used to design an effective automatic drug administration scheme. Even though the model describes the effects of various drugs on the dynamic system, this work is confined to the drugs that kill the invading pathogen and heal the affected organ. From a system theoretic point of view, the drug inputs can be interpreted as control inputs, which can be designed based on control theoretic concepts. The controller is designed based on the principle of dynamic inversion and is found to be effective in curing the �nominal model patient� by killing the invading microbes and healing the damaged organ. A major advantage of this technique is that it leads to a closed-form state feedback form of control. It is also proved from a rigorous mathematical analysis that the internal dynamics of the system remains stable when the proposed controller is applied. A robustness study is also carried out for testing the effectiveness of the drug administration scheme for parameter uncertainties. It is observed from simulation studies that the technique has adequate robustness for many �realistic model patients� having off-nominal parameter values as well.

Optimal blood glucose regulation of diabetic patients using single network adaptive critics

Diabetes is a long-term disease during which the body's production and use of insulin are impaired, causing glucose concentration level to increase in the bloodstream. Regulating blood glucose levels as close to normal as possible leads to a substantial decrease in long-term complications of diabetes. In this paper, an intelligent online feedback-treatment strategy is presented for the control of blood glucose levels in diabetic patients using single network adaptive critic (SNAC) neural networks (which is based on nonlinear optimal control theory). A recently developed mathematical model of the nonlinear dynamics of glucose and insulin interaction in the blood system has been revised and considered for synthesizing the neural network for feedback control. The idea is to replicate the function of pancreatic insulin, i.e. to have a fairly continuous measurement of blood glucose and a situation-dependent insulin injection to the body using an external device. Detailed studies are carried out to analyze the effectiveness of this adaptive critic-based feedback medication strategy. A comparison study with linear quadratic regulator (LQR) theory shows that the proposed nonlinear approach offers some important advantages such as quicker response, avoidance of hypoglycemia problems, etc. Robustness of the proposed approach is also demonstrated from a large number of simulations considering random initial conditions and parametric uncertainties. Copyright (C) 2009 John Wiley & Sons, Ltd.

Evaluation of relative growth limitation due to depletion of glucose and oxygen during fungal growth on a spherical solid particle

A mathematical model for glucose and oxygen consumption, and cell growth during fungal growth on a single solid particle is developed. A moving biofilm is assumed to be present on the surface of the solid particle. Initially only glucose is assumed to be growth limiting and later oxygen transferred from the gas phase on to the biofilm is also assumed to be growth limiting. Glucose is found to be severely growth limiting when assumed to be the only growth limiting factor and its limiting levels far less severe when oxygen limitation is also included. The objective of the model is to gain a better understanding of the mass transfer and relative growth limiting characteristics of glucose and oxygen in fungal growth systems. The results obtained from the model proposed here will be the subject of future work.

Estimating the Threshold Surface Density of Gp120-CCR5 Complexes Necessary for HIV-1 Envelope-Mediated Cell-Cell Fusion

Reduced expression of CCR5 on target CD4(+) cells lowers their susceptibility to infection by R5-tropic HIV-1, potentially preventing transmission of infection and delaying disease progression. Binding of the HIV-1 envelope (Env) protein gp120 with CCR5 is essential for the entry of R5 viruses into target cells. The threshold surface density of gp120-CCR5 complexes that enables HIV-1 entry remains poorly estimated. We constructed a mathematical model that mimics Env-mediated cell-cell fusion assays, where target CD4(+)CCR5(+) cells are exposed to effector cells expressing Env in the presence of a coreceptor antagonist and the fraction of target cells fused with effector cells is measured. Our model employs a reaction network-based approach to describe protein interactions that precede viral entry coupled with the ternary complex model to quantify the allosteric interactions of the coreceptor antagonist and predicts the fraction of target cells fused. By fitting model predictions to published data of cell-cell fusion in the presence of the CCR5 antagonist vicriviroc, we estimated the threshold surface density of gp120-CCR5 complexes for cell-cell fusion as similar to 20 mu m(-2). Model predictions with this threshold captured data from independent cell-cell fusion assays in the presence of vicriviroc and rapamycin, a drug that modulates CCR5 expression, as well as assays in the presence of maraviroc, another CCR5 antagonist, using sixteen different Env clones derived from transmitted or early founder viruses. Our estimate of the threshold surface density of gp120-CCR5 complexes necessary for HIV-1 entry thus appears robust and may have implications for optimizing treatment with coreceptor antagonists, understanding the non-pathogenic infection of non-human primates, and designing vaccines that suppress the availability of target CD4(+)CCR5(+) cells.

Simulation and modelling of orbital debris environment by equivalent breakups

It is well known that the increasing space activities pose a serious threat to future missions. This is mainly due to the presence of spent stages, rockets spacecraft and fragments which can lead to collisions. The calculation of the collision probability of future space vehicles with the orbital debris is necessary for estimating the risk. There is lack of adequately catalogued and openly available detailed information on the explosion characteristics of trackable and untrackable debris data. Such a situation compels one to develop suitable mathematical modelling of the explosion and the resultant debris environment. Based on a study of the available information regarding the fragmentation, subsequent evolution and observation, it turns out to be possible to develop such a mathematical model connecting the dynamical features of the fragmentation with the geometrical/orbital characteristics of the debris and representing the environment through the idea of equivalent breakup. (C) 1997 COSPAR.

Transient analysis of mammalian cell growth in hollow fibre bioreactor

A mathematical model describing the dynamics of mammalian cell growth in hollow fibre bioreactor operated in closed shell mode is developed. Mammalian cells are assumed to grow as an expanding biofilm in the extra-capillary space surrounding the fibre. Diffusion is assumed to be the dominant process in the radial direction while axial convection dominates in the lumen of the bioreactor. The transient simulation results show that steep gradients in the cell number are possible under the condition of substrate limitation. The precise conditions which result in nonuniform growth of cells along the length of the bioreactor are delineated. The effect of various operating conditions, such as substrate feed rate, length of the bioreactor and diffusivity of substrate in different regions of the bioreactor, on the bioreactor performance are evaluated in terms of time required to attain the steady-state. The rime of growth is introduced as a measure of effectiveness factor for the bioreactor and is found to be dependent on two parameters, a modified Peclet number and a Thiele modulus. Diffusion, reaction and/or convection control regimes are identified based on these two parameters. The model is further extended to include dual substrate growth limitations, and the relative growth limiting characteristics of two substrates are evaluated. (C) 1997 Elsevier Science Ltd.

Modelling, simulation, and graphical user interface for industrial gas carburising process

A mathematical model has been developed for the gas carburising (diffusion) process using finite volume method. The computer simulation has been carried out for an industrial gas carburising process. The model's predictions are in good agreement with industrial experimental data and with data collected from the literature. A study of various mass transfer and diffusion coefficients has been carried out in order to suggest which correlations should be used for the gas carburising process. The model has been interfaced in a Windows environment using a graphical user interface. In this way, the model is extremely user friendly. The sensitivity analysis of various parameters such as initial carbon concentration in the specimen, carbon potential of the atmosphere, temperature of the process, etc. has been carried out using the model.

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.

Three-dimensional double-diffusive convection and macrosegregation during non-equilibrium solidification of binary mixtures

A three-dimensional transient mathematical model (following a fixed-grid enthalpy-based continuum formulation) is used to study the interaction of double-diffusive natural convection and non-equilibrium solidification of a binary mixture in a cubic enclosure cooled from a side. Investigations are carried out for two separate test systems, one corresponding to a typical model "metal-alloy analogue" system and other corresponding to a real metal-alloy system. Due to stronger effects of solutal buoyancy in actual metal-alloy systems than in corresponding analogues, the convective transport mechanisms for the two cases are quite different. However, in both cases, similar elements of three-dimensionality are observed in the curvature and spacing of the projected streamlines. As a result of three-dimensional convective flow patterns, a significant solute macrosegregation is observed across the transverse sections of the cavity, which cannot be captured by two-dimensional simulations. (C) 2003 Elsevier Science Ltd. All rights reserved.

The Behaviour of Two-Phase Flow of DNAPL and Water through a Fractured Rock under Confining Pressure

This study presents the characterization of DNAPL and water flow in a fracture under confining pressure. A comprehensive mathematical model and the conditions under which DNAPL will enter an initially water-saturated deforming rock fracture are discussed. A numerical model with which to predict the quantity of each phase in terms of their saturations in deforming rock joint is developed. The effect of varying confining stresses on the traverse time of DNAPL across a fractured aquitard is studied. The sensitivity analysis for physical and hydraulic properties like initial fracture apertures, fracture dips, equivalent fracture aperture and confining pressures are performed and discussed.

Issues in Achievement of Haptic Realism

This paper deals with haptic realism related to Kinematic capabilities of the devices used in manipulation of virtual objects in virtual assembly environments and its effect on achieving haptic realism. Haptic realism implies realistic touch sensation. In virtual world all the operations are to be performed in the same way and with same level of accuracy as in the real world .In order to achieve realism there should be a complete mapping of real and virtual world dimensions. Experiments are conducted to know the kinematic capabilities of the device by comparing the dimensions of the object in the real and virtual world. Registered dimensions in the virtual world are found to be approximately 1.5 times that of the real world. Dimensional variations observed were discrepancy due to exoskeleton and discrepancy due to real and virtual hands. Experiments are conducted to know the discrepancy due to exoskeleton and this discrepancy can be taken care of by either at the hardware or software level. A Mathematical model is proposed to know the discrepancy between real and virtual hands. This could not give a fixed value and can not be taken care of by calibration. Experiments are conducted to figure out how much compensation can be given to achieve haptic realism.

On the fundamentally asynchronous nature of interference in cooperative base station systems

Cooperative transmission by base stations can significantly improve the spectral efficiency of multiuser, multi-cell multiple input multiple output systems. We show that in such systems the multiuser interference is asynchronous by nature, even when perfect timing-advance mechanisms ensure that the desired signal components arrive synchronously. We establish an accurate mathematical model for the asynchronism, and use it to show that the asynchronism leads to a significant performance degradation of existing linear preceding designs that assumed synchronous interference. We consider three different previously proposed precoding designs, and show how to modify them to effectively mitigate asynchronous interference.

Sequences for Phase-Encoded Optical CDMA

In phase-encoded optical CDMA (OCDMA) spreading is achieved by encoding the phase of signal spectrum. Here, a mathematical model for the output signal of a phase-encoded OCDMA system is first derived. This is shown to lead to a performance metric for the design of spreading sequences for asynchronous transmission. Generalized bent functions are used to construct a family of efficient phase-encoding sequences. It is shown how M-ary modulation of these spreading sequences is possible. The problem of designing efficient phaseencoded sequences is then related to the problem of minimizing PMEPR (peak-to-mean envelope power ratio) in an OFDM communication system.