Mathematical modeling and optimal control of malaria

Malaria continues to infect millions and kill hundreds of thousands of people worldwide each year, despite over a century of research and attempts to control and eliminate this infectious disease. Challenges such as the development and spread of drug resistant malaria parasites, insecticide resistance to mosquitoes, climate change, the presence of individuals with subpatent malaria infections which normally are asymptomatic and behavioral plasticity in the mosquito hinder the prospects of malaria control and elimination. In this thesis, mathematical models of malaria transmission and control that address the role of drug resistance, immunity, iron supplementation and anemia, immigration and visitation, and the presence of asymptomatic carriers in malaria transmission are developed. A within-host mathematical model of severe Plasmodium falciparum malaria is also developed. First, a deterministic mathematical model for transmission of antimalarial drug resistance parasites with superinfection is developed and analyzed. The possibility of increase in the risk of superinfection due to iron supplementation and fortification in malaria endemic areas is discussed. The model results calls upon stakeholders to weigh the pros and cons of iron supplementation to individuals living in malaria endemic regions. Second, a deterministic model of transmission of drug resistant malaria parasites, including the inflow of infective immigrants, is presented and analyzed. The optimal control theory is applied to this model to study the impact of various malaria and vector control strategies, such as screening of immigrants, treatment of drug-sensitive infections, treatment of drug-resistant infections, and the use of insecticide-treated bed nets and indoor spraying of mosquitoes. The results of the model emphasize the importance of using a combination of all four controls tools for effective malaria intervention. Next, a two-age-class mathematical model for malaria transmission with asymptomatic carriers is developed and analyzed. In development of this model, four possible control measures are analyzed: the use of long-lasting treated mosquito nets, indoor residual spraying, screening and treatment of symptomatic, and screening and treatment of asymptomatic individuals. The numerical results show that a disease-free equilibrium can be attained if all four control measures are used. A common pitfall for most epidemiological models is the absence of real data; model-based conclusions have to be drawn based on uncertain parameter values. In this thesis, an approach to study the robustness of optimal control solutions under such parameter uncertainty is presented. Numerical analysis of the optimal control problem in the presence of parameter uncertainty demonstrate the robustness of the optimal control approach that: when a comprehensive control strategy is used the main conclusions of the optimal control remain unchanged, even if inevitable variability remains in the control profiles. The results provide a promising framework for the design of cost-effective strategies for disease control with multiple interventions, even under considerable uncertainty of model parameters. Finally, a separate work modeling the within-host Plasmodium falciparum infection in humans is presented. The developed model allows re-infection of already-infected red blood cells. The model hypothesizes that in severe malaria due to parasite quest for survival and rapid multiplication, the Plasmodium falciparum can be absorbed in the already-infected red blood cells which accelerates the rupture rate and consequently cause anemia. Analysis of the model and parameter identifiability using Markov chain Monte Carlo methods is presented.

alpha-Smooth muscle actin and proliferating cell nuclear antigen expression in focal segmental glomerulosclerosis: functional and structural parameters of renal disease progression

The aim of the present study was to investigate the expression of alpha-smooth muscle actin (alpha-SM-actin) and proliferating cell nuclear antigen (PCNA) in renal cortex from patients with focal segmental glomerulosclerosis (FSGS) and their correlations with parameters of renal disease progression. We analyzed renal biopsies from 41 patients with idiopathic FSGS and from 14 control individuals. The alpha-SM-actin immunoreaction was evaluated using a score that reflected the changes in the extent and intensity of staining in the glomerular or cortical area. The PCNA reaction was quantified by counting the labeled cells of the glomeruli or renal cortex. The results, reported as median ± percentile (25th; 75th), showed that the alpha-SM-actin scores in the glomeruli and tubulointerstitium from the renal cortex were 2.0 (2.0; 4.0) and 3.0 (3.0; 4.0), respectively, in patients with FSGS, and 0.5 (0.0; 1.0) and 0.0 (0.0; 0.5) in the controls. The number of PCNA-positive cells per glomerulus and graded field of tubulointerstitium from the renal cortex was 0.2 (0.0; 0.4) and 1.1 (0.3; 2.2), respectively, for patients with FSGS, and 0.0 (0.0; 0.5) and 0.0 (0.0; 0.0) for controls. The present data showed an increase of alpha-SM-actin and PCNA expression in glomeruli and renal cortex from FSGS patients. The extent of immunoreaction for alpha-SM-actin in the tubulointerstitial area was correlated with the intensity of proteinuria. However, there was no correlation between the kidney expression of these proteins and the reciprocal of plasma creatinine level or renal fibrosis. These findings suggest that the immunohistochemical alterations may be reversible.

Interaction between paraventricular nucleus and septal area in the control of physiological responses induced by angiotensin II

We determined the effects of losartan (40 nmol) and PD 123319 (40 nmol) (both non-peptides and selective antagonists of the AT1 and AT2 angiotensin receptors, respectively), and [Sar¹, Ala8] angiotensin II (ANG II) (40 nmol) (a non-selective peptide antagonist of angiotensin receptors) injected into the paraventricular nucleus (PVN) on the water and salt appetite, diuresis and natriuresis and mean arterial pressure (MAP) induced by administration of 10 nmol of ANG II into the medial septal area (MSA) of male Holtzman rats weighing 250-300 g. The volume of drug solution injected was 0.5 µl over a period of 10-15 s. The responses were measured over a period of 120 min. ANG II alone injected into the MSA induced an increase in all the above parameters (8.1 ± 1.2, 1.8 ± 0.3, and 17.1 ± 1.0 ml, 217 ± 25 µEq/120 min, and 24 ± 4 mmHg, respectively, N = 10-12) compared with vehicle-treated rats (1.4 ± 0.2, 0.6 ± 0.1, and 9.3 ± 0.5 ml, 47 ± 5 µEq/120 min, and 4.1 ± 0.8 mmHg, respectively, N = 10-14). Pretreatment with losartan and [Sar¹, Ala8] ANG II completely abolished the water and sodium intake, and the pressor increase (0.5 ± 0.2, 1.1 ± 0.2, 0.5 ± 0.2, and 0.8 ± 0.2 ml, and 1.2 ± 3.9, 31 ± 4.6 mmHg, respectively, N = 9-12), whereas losartan blunted the urinary and sodium excretion induced by ANG II (13.9 ± 1.0 ml and 187 ± 10 µEq/120 min, respectively, N = 9). Pretreatment with PD 123319 and [Sar¹, Ala8] ANG II blocked the urinary and sodium excretion (10.7 ± 0.8, 9.8 ± 0.7 ml, and 67 ± 13 and 57 ± 17 µEq/120 min, respectively, N = 9), whereas pretreatment with PD 123319 partially blocked the water and sodium intake, and the MAP induced by ANG II administration (2.3 ± 0.3, 1.1 ± 0.1 ml, and 12 ± 3 mmHg, respectively, N = 9-10). These results suggest the angiotensinergic effect of the MSA on the AT1 and AT2 receptors of the PVN in terms of water and sodium homeostasis and MAP modulation.

Bone mineral density in juvenile systemic lupus erythematosus

We evaluated spine bone mineral density (BMD) in Brazilian children with juvenile systemic lupus erythematosus (JSLE) in order to detect potential predictors of reduction in bone mass. A cross-sectional study of BMD at the lumbar spine level (L2-L4) was conducted on 16 female JSLE patients aged 6-17 years. Thirty-two age-matched healthy girls were used as control. BMD at the lumbar spine was measured by dual-energy X-ray absorptiometry. Weight, height and pubertal Tanner stage were determined in patients and controls. Disease duration, mean daily steroid doses, mean cumulative steroid doses and JSLE activity measured by the systemic lupus erythematosus disease activity index (SLEDAI) were determined for all JSLE patients based on their medical charts. All parameters were used as potential determinant factors for bone loss. Lumbar BMD tended to be lower in the JSLE patients, however, this difference was not statistically significant (P = 0.10). No significant correlation was observed in JSLE girls between BMD and age, height, Tanner stage, disease duration, corticosteroid use or disease activity. We found a weak correlation between BMD and weight (r = 0.672). In the JSLE group we found no significant parameters to correlate with reduced bone mass. Disease activity and mean cumulative steroid doses were not related to BMD values. We did not observe reduced bone mass in female JSLE.

Alternative parameters for echocardiographic assessment of fetal diastolic function

Alternative methods to assess ventricular diastolic function in the fetus are proposed. Fetal myocardial hypertrophy in maternal diabetes was used as a model of decreased left ventricular compliance (LVC), and fetal respiratory movements as a model of increased LVC. Comparison of three groups of fetuses showed that, in 10 fetuses of diabetic mothers (FDM) with septal hypertrophy (SH), the mean excursion index of the septum primum (EISP) (ratio between the linear excursion of the flap valve and the left atrial diameter) was 0.36 ± 0.09, in 8 FDM without SH it was 0.51 ± 0.09 (P = 0.001), and in the 8 normal control fetuses (NCF) it was 0.49 ± 0.12 (P = 0.003). In another study, 28 fetuses in apnea had a mean EISP of 0.39 ± 0.05 which increased to 0.57 ± 0.07 during respiration (P < 0.001). These two studies showed that the mobility of the septum primum was reduced when LVC was decreased and was increased when LVC was enhanced. Mean pulmonary vein pulsatility was higher in 14 FDM (1.83 ± 1.21) than in 26 NCF (1.02 ± 0.31; P = 0.02). In the same fetuses, mean left atrial shortening was decreased (0.40 ± 0.11) in relation to NCF (0.51 ± 0.09; P = 0.011). These results suggest that FDM may have a higher preload than normal controls, probably as a result of increased myocardial mass and LV hypertrophy. Prenatal assessment of LV diastolic function by fetal echocardiography should include analysis of septum primum mobility, pulmonary vein pulsatility, and left atrial shortening.

The selective and non-selective cyclooxygenase inhibitors valdecoxib and piroxicam induce the same postoperative analgesia and control of trismus and swelling after lower third molar removal

We compared the clinical efficacy of orally administered valdecoxib and piroxicam for the prevention of pain, trismus and swelling after removal of horizontally and totally intrabony impacted lower third molars. Twenty-five patients were scheduled to undergo removal of symmetrically positioned lower third molars in two separate appointments. Valdecoxib (40 mg) or piroxicam (20 mg) was administered in a double-blind, randomized and crossed manner for 4 days after the surgical procedures. Objective and subjective parameters were recorded for comparison of postoperative courses. Both agents were effective for postoperative pain relief (N = 19). There was a similar mouth opening at suture removal compared with the preoperative values (86.14 ± 4.36 and 93.12 ± 3.70% of the initial measure for valdecoxib and piroxicam, respectively; ANOVA). There was no significant difference regarding the total amount of rescue medication taken by the patients treated with valdecoxib or piroxicam (173.08 ± 91.21 and 461.54 ± 199.85 mg, respectively; Wilcoxon test). There were no significant differences concerning the swelling observed on the second postoperative day compared to baseline measures (6.15 ± 1.84 and 8.46 ± 2.04 mm for valdecoxib and piroxicam, respectively; ANOVA) or on the seventh postoperative day (1.69 ± 1.61 and 2.23 ± 2.09 mm for valdecoxib and piroxicam, respectively; ANOVA). The cyclooxygenase-2 selective inhibitor valdecoxib is as effective as the non-selective cyclooxygenase inhibitor piroxicam for pain, trismus and swelling control after removal of horizontally and totally intrabony impacted lower third molars.

Autonomic control of cardiorespiratory interactions in fish, amphibians and reptiles

Control of the heart rate and cardiorespiratory interactions (CRI) is predominantly parasympathetic in all jawed vertebrates, with the sympathetic nervous system having some influence in tetrapods. Respiratory sinus arrhythmia (RSA) has been described as a solely mammalian phenomenon but respiration-related beat-to-beat control of the heart has been described in fish and reptiles. Though they are both important, the relative roles of feed-forward central control and peripheral reflexes in generating CRI vary between groups of fishes and probably between other vertebrates. CRI may relate to two locations for the vagal preganglionic neurons (VPN) and in particular cardiac VPN in the brainstem. This has been described in representatives from all vertebrate groups, though the proportion in each location is variable. Air-breathing fishes, amphibians and reptiles breathe discontinuously and the onset of a bout of breathing is characteristically accompanied by an immediate increase in heart rate plus, in the latter two groups, a left-right shunting of blood through the pulmonary circuit. Both the increase in heart rate and opening of a sphincter on the pulmonary artery are due to withdrawal of vagal tone. An increase in heart rate following a meal in snakes is related to withdrawal of vagal tone plus a non-adrenergic-non-cholinergic effect that may be due to humoral factors released by the gut. Histamine is one candidate for this role.

Immunological and biochemical parameters of patients with metabolic syndrome and the participation of oxidative and nitroactive stress

Metabolic syndrome (MS) is a multifactorial disease involving inflammatory activity and endothelial dysfunction. The aim of the present study was to evaluate the relationship between the changes in lipoperoxidation, in immunological and biochemical parameters and nitric oxide metabolite (NOx) levels in MS patients. Fifty patients with MS (4 males/46 females) and 50 controls (3 males/47 females) were studied. Compared to control (Mann-Whitney test), MS patients presented higher serum levels (P < 0.05) of fibrinogen: 314 (185-489) vs 262 (188-314) mg/dL, C-reactive protein (CRP): 7.80 (1.10-46.50) vs 0.70 (0.16-5.20) mg/dL, interleukin-6: 3.96 (3.04-28.18) vs 3.33 (2.55-9.63) pg/mL, uric acid: 5.45 (3.15-9.65) vs 3.81 (2.70-5.90) mg/dL, and hydroperoxides: 20,689 (19,076-67,182) vs 18,636 (15,926-19,731) cpm. In contrast, they presented lower (P < 0.05) adiponectin: 7.11 (3.19-18.22) vs 12.31 (9.11-27.27) µg/mL, and NOx levels: 5.69 (2.36-8.18) vs 6.72 (5.14-12.43) µM. NOx was inversely associated (Spearman’s rank correlation) with body mass index (r = -0.2858, P = 0.0191), insulin resistance determined by the homeostasis model assessment (r = -0.2530, P = 0.0315), CRP (r = -0.2843, P = 0.0171) and fibrinogen (r = -0.2464, P = 0.0413), and positively correlated with hydroperoxides (r = 0.2506, P = 0.0408). In conclusion, NOx levels are associated with obesity, insulin resistance, oxidative stress, and inflammatory markers. The high uric acid levels together with reactive oxygen species generation may be responsible for the reduced NO levels, which in turn lead to endothelial dysfunction. The elevated plasma chemiluminescence reflecting both increased plasma oxidation and reduced antioxidant capacity may play a role in the MS mechanism.

Influence of eNOS gene polymorphism on cardiometabolic parameters in response to physical training in postmenopausal women

The health-promoting effects of exercise training (ET) are related to nitric oxide (NO) production and/or its bioavailability. The objective of this study was to determine whether single nucleotide polymorphism of the endothelial NO synthase (eNOS) gene at positions -786T>C, G894T (Glu298Asp) and at the variable number of tandem repeat (VNTR) Intron 4b/a would interfere with the cardiometabolic responses of postmenopausal women submitted to physical training. Forty-nine postmenopausal women were trained in sessions of 30-40 min, 3 days a week for 8 weeks. Genotypes, oxidative stress status and cardiometabolic parameters were then evaluated in a double-blind design. Both systolic and diastolic blood pressure values were significantly reduced after ET, which was genotype-independent. However, women without eNOS gene polymorphism at position -786T>C (TT genotype) and Intron 4b/a (bb genotype) presented a better reduction of total cholesterol levels (-786T>C: before = 213 ± 12.1, after = 159.8 ± 14.4, Δ = -24.9% and Intron 4b/a: before = 211.8 ± 7.4, after = 180.12 ± 6.4 mg/dL, Δ = -15%), and LDL cholesterol (-786T>C: before = 146.1 ± 13.3, after = 82.8 ± 9.2, Δ = -43.3% and Intron 4b/a: before = 143.2 ± 8, after = 102.7 ± 5.8 mg/dL, Δ = -28.3%) in response to ET compared to those who carried the mutant allele. Superoxide dismutase activity was significantly increased in trained women whereas no changes were observed in malondialdehyde levels. Women without eNOS gene polymorphism at position -786T>C and Intron 4b/a showed a greater reduction of plasma cholesterol levels in response to ET. Furthermore, no genotype influence was observed on arterial blood pressure or oxidative stress status in this population.

Effects of rosiglitazone on serum paraoxonase activity and metabolic parameters in patients with type 2 diabetes mellitus

Human serum paraoxonase contributes to the anti-atherogenic effect of high-density lipoprotein cholesterol (HDL-C) and has been shown to protect both low-density lipoprotein cholesterol (LDL-C) and HDL-C against lipid peroxidation. We investigated the effects of rosiglitazone on paraoxonase activity and metabolic parameters in patients with type 2 diabetes mellitus [50 patients (30 males, 20 females); mean±SD age: 58.7±9.2 years, body mass index: 28.2±4.1'kg/m2], in whom glucose control could not be achieved despite treatment with metformin, sulphonylurea, and/or insulin. The patients were given 4'mg/day rosiglitazone for 3 months in addition to their usual treatment. Serum paraoxonase activity, malondialdehyde, homocysteine, and lipid profile were measured at the time of initiation and at the end of therapy with rosiglitazone. After rosiglitazone therapy, serum levels of HDL-C, apolipoprotein A-1, and paraoxonase activity increased significantly (P<0.05) and malondialdehyde, homocysteine, lipoprotein(a), and glucose levels decreased significantly (P<0.05), but no significant changes in levels of total cholesterol and apolipoprotein B were observed. Triglyceride levels also increased significantly (P<0.05). Rosiglitazone treatment led to an improvement in glycemic control and to an increase in paraoxonase activity and HDL-C levels. Although rosiglitazone showed favorable effects on oxidant/antioxidant balance and lipid profile, further studies are needed to determine the effect of rosiglitazone on cardiovascular risk factors and cardiovascular morbidity and mortality.

Limb immobilization alters functional electrophysiological parameters of sciatic nerve

Immobilization, used in clinical practice to treat traumatologic problems, causes changes in muscle, but it is not known whether changes also occur in nerves. We investigated the effects of immobilization on excitability and compound action potential (CAP) and the ultrastructure of the rat sciatic nerve. Fourteen days after immobilization of the right leg of adult male Wistar rats (n=34), animals were killed and the right sciatic nerve was dissected and mounted in a moist chamber. Nerves were stimulated at a baseline frequency of 0.2 Hz and tested for 2 min at 20, 50, and 100 Hz. Immobilization altered nerve excitability. Rheobase and chronaxy changed from 3.13±0.05 V and 52.31±1.95 µs (control group, n=13) to 2.84±0.06 V and 59.71±2.79 µs (immobilized group, n=15), respectively. Immobilization altered the amplitude of CAP waves and decreased the conduction velocity of the first CAP wave (from 93.63±7.49 to 79.14±5.59 m/s) but not of the second wave. Transmission electron microscopy showed fragmentation of the myelin sheath of the sciatic nerve of immobilized limbs and degeneration of the axon. In conclusion, we demonstrated that long-lasting leg immobilization can induce alterations in nerve function.

Examining the side effects of sucrose for pain relief in preterm infants: a case-control study

Sucrose solution is recommended as relevant pain relief management in neonates during acute painful procedures; however, only a few studies have analyzed the potentially adverse effects of sucrose administration to preterm neonates. The goal of this study was to examine the potential side effects of sucrose for pain relief in preterm infants, assessing feeding and weight gain during hospitalization and their feeding patterns postdischarge. The study sample consisted of 43 preterm neonates divided into two groups: a sucrose group (SG, n=18) and a control group (CG, n=25) in which no sucrose was administered. The SG received 0.5 mL/kg 25% oral sucrose for 2 min prior to all acute painful procedures during three consecutive days. A prospective review of medical charts was performed for all samples. The study was done prior to implementation of the institutional sucrose guidelines as a routine service, and followed all ethical requirements. There were no statistically significant differences between groups in terms of weight gain, length of stay with orogastric tubes, and parenteral feeding. Postdischarge, infant nutritional intake included feeding human milk to 67% of the SG and 74% of the CG. There were no statistically significant differences between groups regarding human milk feeding patterns postdischarge. Neonate feeding patterns and weight gain were unaffected following the short-term use of sucrose for pain relief.

Novel control, haptic and calibration methods for teleoperated electrohydraulic servo systems

The aim of this thesis is to propose a novel control method for teleoperated electrohydraulic servo systems that implements a reliable haptic sense between the human and manipulator interaction, and an ideal position control between the manipulator and the task environment interaction. The proposed method has the characteristics of a universal technique independent of the actual control algorithm and it can be applied with other suitable control methods as a real-time control strategy. The motivation to develop this control method is the necessity for a reliable real-time controller for teleoperated electrohydraulic servo systems that provides highly accurate position control based on joystick inputs with haptic capabilities. The contribution of the research is that the proposed control method combines a directed random search method and a real-time simulation to develop an intelligent controller in which each generation of parameters is tested on-line by the real-time simulator before being applied to the real process. The controller was evaluated on a hydraulic position servo system. The simulator of the hydraulic system was built based on Markov chain Monte Carlo (MCMC) method. A Particle Swarm Optimization algorithm combined with the foraging behavior of E. coli bacteria was utilized as the directed random search engine. The control strategy allows the operator to be plugged into the work environment dynamically and kinetically. This helps to ensure the system has haptic sense with high stability, without abstracting away the dynamics of the hydraulic system. The new control algorithm provides asymptotically exact tracking of both, the position and the contact force. In addition, this research proposes a novel method for re-calibration of multi-axis force/torque sensors. The method makes several improvements to traditional methods. It can be used without dismantling the sensor from its application and it requires smaller number of standard loads for calibration. It is also more cost efficient and faster in comparison to traditional calibration methods. The proposed method was developed in response to re-calibration issues with the force sensors utilized in teleoperated systems. The new approach aimed to avoid dismantling of the sensors from their applications for applying calibration. A major complication with many manipulators is the difficulty accessing them when they operate inside a non-accessible environment; especially if those environments are harsh; such as in radioactive areas. The proposed technique is based on design of experiment methodology. It has been successfully applied to different force/torque sensors and this research presents experimental validation of use of the calibration method with one of the force sensors which method has been applied to.

Doubly fed wind turbine performance in variable grid conditions

Wind turbines based on doubly fed induction generators (DFIG) become the most popular solution in high power wind generation industry. While this topology provides great performance with the reduced power rating of power converter, it has more complicated structure in comparison with full-rated topologies, and therefore leads to complexity of control algorithms and electromechanical processes in the system. The purpose of presented study is to present a proper vector control scheme for the DFIG and overall control for the WT to investigate its behavior at different wind speeds and in different grid voltage conditions: voltage sags, magnitude and frequency variations. The key principles of variable-speed wind turbine were implemented in simulation model and demonstrated during the study. Then, based on developed control scheme and mathematical model, the set of simulation is made to analyze reactive power capabilities of the DFIG wind turbine. Further, the rating of rotor-side converter is modified to not only generate active rated active power, but also to fulfill Grid Codes. Results of modelling and analyzing of the DFIG WT behavior under different speeds and different voltage conditions are presented in the work.

Design of automatic control system for the precipitation of bromelain from the extract of pineapple wastes

In this work, bromelain was recovered from ground pineapple stem and rind by means of precipitation with alcohol at low temperature. Bromelain is the name of a group of powerful protein-digesting, or proteolytic, enzymes that are particularly useful for reducing muscle and tissue inflammation and as a digestive aid. Temperature control is crucial to avoid irreversible protein denaturation and consequently to improve the quality of the enzyme recovered. The process was carried out alternatively in two fed-batch pilot tanks: a glass tank and a stainless steel tank. Aliquots containing 100 mL of pineapple aqueous extract were fed into the tank. Inside the jacketed tank, the protein was exposed to unsteady operating conditions during the addition of the precipitating agent (ethanol 99.5%) because the dilution ratio "aqueous extract to ethanol" and heat transfer area changed. The coolant flow rate was manipulated through a variable speed pump. Fine tuned conventional and adaptive PID controllers were on-line implemented using a fieldbus digital control system. The processing performance efficiency was enhanced and so was the quality (enzyme activity) of the product.