Applying Mathematical Models to Study the Role of the Immune System in Chronic Myelogenous Leukemia

Although tyrosine kinase inhibitors (TKIs) such as imatinib have transformed chronic myelogenous leukemia (CML) into a chronic condition, these therapies are not curative in the majority of cases. Most patients must continue TKI therapy indefinitely, a requirement that is both expensive and that compromises a patient's quality of life. While TKIs are known to reduce leukemic cells' proliferative capacity and to induce apoptosis, their effects on leukemic stem cells, the immune system, and the microenvironment are not fully understood. A more complete understanding of their global therapeutic effects would help us to identify any limitations of TKI monotherapy and to address these issues through novel combination therapies. Mathematical models are a complementary tool to experimental and clinical data that can provide valuable insights into the underlying mechanisms of TKI therapy. Previous modeling efforts have focused on CML patients who show biphasic and triphasic exponential declines in BCR-ABL ratio during therapy. However, our patient data indicates that many patients treated with TKIs show fluctuations in BCR-ABL ratio yet are able to achieve durable remissions. To investigate these fluctuations, we construct a mathematical model that integrates CML with a patient's autologous immune response to the disease. In our model, we define an immune window, which is an intermediate range of leukemic concentrations that lead to an effective immune response against CML. While small leukemic concentrations provide insufficient stimulus, large leukemic concentrations actively suppress a patient's immune system, thus limiting it's ability to respond. Our patient data and modeling results suggest that at diagnosis, a patient's high leukemic concentration is able to suppress their immune system. TKI therapy drives the leukemic population into the immune window, allowing the patient's immune cells to expand and eventually mount an efficient response against the residual CML. This response drives the leukemic population below the immune window, causing the immune population to contract and allowing the leukemia to partially recover. The leukemia eventually reenters the immune window, thus stimulating a sequence of weaker immune responses as the two populations approach equilibrium. We hypothesize that a patient's autologous immune response to CML may explain the fluctuations in BCR-ABL ratio that are regularly seen during TKI therapy. These fluctuations may serve as a signature of a patient's individual immune response to CML. By applying our modeling framework to patient data, we are able to construct an immune profile that can then be used to propose patient-specific combination therapies aimed at further reducing a patient's leukemic burden. Our characterization of a patient's anti-leukemia immune response may be especially valuable in the study of drug resistance, treatment cessation, and combination therapy.

A Numerical study of fluid flow in peristaltic pumps and polymeric elastic pipes

This thesis work deals with a mathematical description of flow in polymeric pipe and in a specific peristaltic pump. This study involves fluid-structure interaction analysis in presence of complex-turbulent flows treated in an arbitrary Lagrangian-Eulerian (ALE) framework. The flow simulations are performed in COMSOL 4.4, as 2D axial symmetric model, and ABAQUS 6.14.1, as 3D model with symmetric boundary conditions. In COMSOL, the fluid and structure problems are coupled by monolithic algorithm, while ABAQUS code links ABAQUS CFD and ABAQUS Standard solvers with single block-iterative partitioned algorithm. For the turbulent features of the flow, the fluid model in both codes is described by RNG k-ϵ. The structural model is described, on the basis of the pipe material, by Elastic models or Hyperelastic Neo-Hookean models with Rayleigh damping properties. In order to describe the pulsatile fluid flow after the pumping process, the available data are often defective for the fluid problem. Engineering measurements are normally able to provide average pressure or velocity at a cross-section. This problem has been analyzed by McDonald's and Womersley's work for average pressure at fixed cross section by Fourier analysis since '50, while nowadays sophisticated techniques including Finite Elements and Finite Volumes exist to study the flow. Finally, we set up peristaltic pipe simulations in ABAQUS code, by using the same model previously tested for the fl uid and the structure.

Numerical and analytical study of an atherosclerosis inflammatory disease model

We study a reaction–diffusion mathematical model for the evolution of atherosclerosis as an inflammation process by combining analytical tools with computer-intensive numerical calculations. The computational work involved the calculation of more than sixty thousand solutions of the full reaction–diffusion system and lead to the complete characterisation of the ωω-limit for every initial condition. Qualitative properties of the solution are rigorously proved, some of them hinted at by the numerical study

Chebanov law and Vakar formula in mathematical models of complex systems

Ecological models written in a mathematical language L(M) or model language, with a given style or methodology can be considered as a text. It is possible to apply statistical linguistic laws and the experimental results demonstrate that the behaviour of a mathematical model is the same of any literary text of any natural language. A text has the following characteristics: (a) the variables, its transformed functions and parameters are the lexic units or LUN of ecological models; (b) the syllables are constituted by a LUN, or a chain of them, separated by operating or ordering LUNs; (c) the flow equations are words; and (d) the distribution of words (LUM and CLUN) according to their lengths is based on a Poisson distribution, the Chebanov's law. It is founded on Vakar's formula, that is calculated likewise the linguistic entropy for L(M). We will apply these ideas over practical examples using MARIOLA model. In this paper it will be studied the problem of the lengths of the simple lexic units composed lexic units and words of text models, expressing these lengths in number of the primitive symbols, and syllables. The use of these linguistic laws renders it possible to indicate the degree of information given by an ecological model.

Estudo comparativo do desempenho de controladores robustos aplicados a um sistema de tanques acoplados.

Currently the uncertain system has attracted much academic community from the standpoint of scientific research and also practical applications. A series of mathematical approaches emerge in order to troubleshoot the uncertainties of real physical systems. In this context, the work presented here focuses on the application of control theory in a nonlinear dynamical system with parametric variations in order and robustness. We used as the practical application of this work, a system of tanks Quanser associates, in a configuration, whose mathematical model is represented by a second order system with input and output (SISO). The control system is performed by PID controllers, designed by various techniques, aiming to achieve robust performance and stability when subjected to parameter variations. Other controllers are designed with the intention of comparing the performance and robust stability of such systems. The results are obtained and compared from simulations in Matlab-simulink.

Modeling and Model Predictive Control of an Industrial Hydraulic System

This thesis aims to illustrate the construction of a mathematical model of a hydraulic system, oriented to the design of a model predictive control (MPC) algorithm. The modeling procedure starts with the basic formulation of a piston-servovalve system. The latter is a complex non linear system with some unknown and not measurable effects that constitute a challenging problem for the modeling procedure. The first level of approximation for system parameters is obtained basing on datasheet informations, provided workbench tests and other data from the company. Then, to validate and refine the model, open-loop simulations have been made for data matching with the characteristics obtained from real acquisitions. The final developed set of ODEs captures all the main peculiarities of the system despite some characteristics due to highly varying and unknown hydraulic effects, like the unmodeled resistive elements of the pipes. After an accurate analysis, since the model presents many internal complexities, a simplified version is presented. The latter is used to linearize and discretize correctly the non linear model. Basing on that, a MPC algorithm for reference tracking with linear constraints is implemented. The results obtained show the potential of MPC in this kind of industrial applications, thus a high quality tracking performances while satisfying state and input constraints. The increased robustness and flexibility are evident with respect to the standard control techniques, such as PID controllers, adopted for these systems. The simulations for model validation and the controlled system have been carried out in a Python code environment.

Anaerobic And Aerobic Performances In Elite Basketball Players.

THE PURPOSE OF THIS STUDY WAS TO PROPOSE A SPECIFIC LACTATE MINIMUM TEST FOR ELITE BASKETBALL PLAYERS CONSIDERING THE: Running Anaerobic Sprint Test (RAST) as a hyperlactatemia inductor, short distances (specific distance, 20 m) during progressive intensity and mathematical analysis to interpret aerobic and anaerobic variables. The basketball players were assigned to four groups: All positions (n=26), Guard (n= 7), Forward (n=11) and Center (n=8). The hyperlactatemia elevation (RAST) method consisted of 6 maximum sprints over 35 m separated by 10 s of recovery. The progressive phase of the lactate minimum test consisted of 5 stages controlled by an electronic metronome (8.0, 9.0, 10.0, 11.0 and 12.0 km/h) over a 20 m distance. The RAST variables and the lactate values were analyzed using visual and mathematical models. The intensity of the lactate minimum test, determined by a visual method, reduced in relation to polynomial fits (2nd degree) for the Small Forward positions and General groups. The Power and Fatigue Index values, determined by both methods, visual and 3rd degree polynomial, were not significantly different between the groups. In conclusion, the RAST is an excellent hyperlactatemia inductor and the progressive intensity of lactate minimum test using short distances (20 m) can be specifically used to evaluate the aerobic capacity of basketball players. In addition, no differences were observed between the visual and polynomial methods for RAST variables, but lactate minimum intensity was influenced by the method of analysis.

Modelagem matemática e simulação numérica do resfriamento rápido de morango com ar forçado

This work approaches the forced air cooling of strawberry by numerical simulation. The mathematical model that was used describes the process of heat transfer, based on the Fourier's law, in spherical coordinates and simplified to describe the one-dimensional process. For the resolution of the equation expressed for the mathematical model, an algorithm was developed based on the explicit scheme of the numerical method of the finite differences and implemented in the scientific computation program MATLAB 6.1. The validation of the mathematical model was made by the comparison between theoretical and experimental data, where strawberries had been cooled with forced air. The results showed to be possible the determination of the convective heat transfer coefficient by fitting the numerical and experimental data. The methodology of the numerical simulations was showed like a promising tool in the support of the decision to use or to develop equipment in the area of cooling process with forced air of spherical fruits.

Identificação eletrônica para avaliação do comportamento dos suínos na fase de gestação

The efficiency of swine production performance depends on the herd administration, such as good nutrition, sanitary control, facilities and appropriate environmental conditions. The concept of this production model is directly related with the reduction of selective losses and the process control. Each production segment is controlled to reach the optimization in the system totality, it is necessary to apply animals handling concepts, environmental control implementation, diseases control, nutrition control, information concerning in guaranteeing the animal welfare and individual identification. The present work presents as objective the development of the mathematical model to evaluate interactions among the internal atmosphere of the installation and the thermal animals preference, in the expectation of detecting a relationship among the frequency access to the drinking fountain and the atmosphere conditions - temperature, black globe temperature and relative humidity, using as tool the electronic identification. The results obtained by the mathematical model, allowed to conclude accurately the evaluation of the swine thermal preference correlating with the climatic variables in the pregnancy stage.

Determinação da difusividade efetiva de raiz de chicória

Inulin is a fructooligosacharide found in diverse agricultural products, amongst them garlic, banana, Jerusalem artichoke and chicory root. Inulin generally is used in developed countries, as a substitute of sugar and/or fat due to its characteristics of fitting as functional and dietary food. Chicory root is usually used as source and raw material for commercial extration of inulin. The experiments consisted on drying sliced chicory roots based on a factorial experimental design in a convective dryer whose alows the air to pass perpendicularly through the tray. Effective diffusivity (dependent variable) has been determined for each experimental combination of independent variables (air temperature and velocity). The data curves have been fitted by the solution of the second Fick law and Page's model. Effective difusivity varied from 3.51 x 10-10 m² s-1 to 1.036 x 10-10 m² s-1. It is concluded that, for the range of studied values, air temperature is the only statistically significant variable. So, a first order mathematical model was obtained, representing effective diffusivity behavior as function of air temperature. The best drying condition was correspondent to the trial using the highest drying air temperature.

Modelos de transferência de metal pesado na cana-de-açúcar adubada com composto de lixo urbano

The research approaches recycling of urban waste compost (UWC) as an alternative fertilizer for sugarcane crop and as a social and environmental solution to the solids residuals growth in urban centers. A mathematical model was used in order to know the metal dynamics as decision support tool, aiming to establish of criteria and procedures for UWC's safe use, limited by the amount of heavy metal. A compartmental model was developed from experimental data in controlled conditions and partially checked with field data. This model described the heavy metal transference in the system soil-root-aerial portion of sugarcane plants and concluded that nickel was metal to be concern, since it takes approximately three years to be attenuated in the soil, reaching the aerial portions of the plant at high concentrations. Regarding factors such as clay content, oxide level and soil pH, it was observed that for soil with higher buffering capacity, the transfer of the majority of the metals was slower. This model may become an important tool for the attainment of laws regarding the UWC use, aiming to reduce environment contamination the waste accumulation and production costs.

Otimização do descascamento químico de raízes do yacon (Polymnia sonchifolia Poepp.)

A full factorial design 2³ was used to evaluate the effect of process variables in chemical peeling of yacon roots, cultivated in Curitiba, State of Paraná. Eleven treatments, with three central points, were done in which they had been evaluated at three different levels of sodium hydroxide solution, % (g/100 mL) [6, 10, 14], temperature of the same solution, °C [70, 80, 90], and residence time in the sodium hydroxide solution, minutes [2, 4, 6]. All the studied variables had affected significantly (p<0.05) the yield of yacon roots subjected to chemical peeling. The variable that most affected the yield was the time of permanence in the sodium hydroxide solution. The mathematical model obtained for the yield (%) was good with R² aj = 0.8497, and non significant lack of fit (p=0.9312).Therefore, the model can be used for predictive purposes. In the central point a satisfactory yield (84% to 87%) with a high percentage of removed peel was obtained (96% to 98%) indicating that the treatment with 10% of sodium hydroxide solution, temperature of 80º C per 4 minutes can be used in the chemical peeling of yacon roots.

Limiar de fadiga neuromuscular determinado por diferentes periodos de analise do sinal eletromiografico

Universidade Estadual de Campinas . Faculdade de Educação Física

A hybrid heuristic for the multi-plant capacitated lot sizing problem with setup carry-over

This paper addresses the capacitated lot sizing problem (CLSP) with a single stage composed of multiple plants, items and periods with setup carry-over among the periods. The CLSP is well studied and many heuristics have been proposed to solve it. Nevertheless, few researches explored the multi-plant capacitated lot sizing problem (MPCLSP), which means that few solution methods were proposed to solve it. Furthermore, to our knowledge, no study of the MPCLSP with setup carry-over was found in the literature. This paper presents a mathematical model and a GRASP (Greedy Randomized Adaptive Search Procedure) with path relinking to the MPCLSP with setup carry-over. This solution method is an extension and adaptation of a previously adopted methodology without the setup carry-over. Computational tests showed that the improvement of the setup carry-over is significant in terms of the solution value with a low increase in computational time.

Modelagem matemática do controle de brucelose bovina por vacinação

As fêmeas bovinas, por sua importância na transmissão e na manutenção da brucelose, constituíram o alvo dos inquéritos do Programa Nacional de Controle e Erradicação da Brucelose e da Tuberculose Animal. Com base em informações obtidas em unidades federativas onde foram realizados inquéritos sorológicos e observadas prevalências de animais acima de 2%, elaborou-se um modelo para simular a dinâmica da brucelose em rebanhos bovinos formados exclusivamente por fêmeas, analisando o efeito de estratégias de vacinação. Para baixa cobertura vacinal, da ordem de 30%, o tempo para reduzir a prevalência a 2%, valor adotado como referência, pode ser longo, aproximando-se do dobro do tempo necessário para uma cobertura mais alta, de 90%. De acordo com o modelo, o tempo para reduzir a prevalência a 1% ou 2%, que permitam passar à fase de erradicação, pode chegar a uma década. Recomenda-se a intensificação do esforço para a vacinação de fêmeas, procurando atingir alta cobertura vacinal.