A 3D Approach Using a Control Algorithm to Minimize the Effects on the Healthy Tissue in the Hyperthermia for Cancer Treatment.

According to the World Health Organization, cancer is a worldwide health problem. Its high mortality rate motivates scientists to study new treatments. One of these new treatments is hyperthermia using magnetic nanoparticles. This treatment consists in submitting the target region with a low-frequency magnetic field to increase its temperature over 43 °C, as the threshold for tissue damage and leading the cells to necrosis. This paper uses an in silico three-dimensional Pennes' model described by a set of partial differential equations (PDEs) to estimate the percentage of tissue damage due to hyperthermia. Differential evolution, an optimization method, suggests the best locations to inject the nanoparticles to maximize tumor cell death and minimize damage to healthy tissue. Three different scenarios were performed to evaluate the suggestions obtained by the optimization method. The results indicate the positive impact of the proposed technique: a reduction in the percentage of healthy tissue damage and the complete damage of the tumors were observed. In the best scenario, the optimization method was responsible for decreasing the healthy tissue damage by 59% when the nanoparticles injection sites were located in the non-intuitive points indicated by the optimization method. The numerical solution of the PDEs is computationally expensive. This work also describes the implemented parallel strategy based on CUDA to reduce the computational costs involved in the PDEs resolution. Compared to the sequential version executed on the CPU, the proposed parallel implementation was able to speed the execution time up to 84.4 times.

A Poroelastic Approach for Modelling Myocardial Oedema in Acute Myocarditis.

Myocarditis is a general set of mechanisms that manifest themselves into the inflammation of the heart muscle. In 2017, more than 3 million people were affected by this disease worldwide, causing about 47,000 deaths. Many aspects of the origin of this disease are well known, but several important questions regarding the disease remain open. One of them is why some patients develop a significantly localised inflammation while others develop a much more diffuse inflammation, reaching across large portions of the heart. Furthermore, the specific role of the pathogenic agent that causes inflammation as well as the interaction with the immune system in the progression of the disease are still under discussion. Providing answers to these crucial questions can have an important impact on patient treatment. In this scenario, computational methods can aid specialists to understand better the relationships between pathogens and the immune system and elucidate why some patients develop diffuse myocarditis. This paper alters a recently developed model to study the myocardial oedema formation in acute infectious myocarditis. The model describes the finite deformation regime using partial differential equations to represent tissue displacement, fluid pressure, fluid phase, and the concentrations of pathogens and leukocytes. A sensitivity analysis was performed to understand better the influence of the most relevant model parameters on the disease dynamics. The results showed that the poroelastic model could reproduce local and diffuse myocarditis dynamics in simplified and complex geometrical domains.

Timing the race of vaccination, new variants, and relaxing restrictions during COVID-19 pandemic.

Late in 2019, China identified a new type of coronavirus, SARS-CoV-2, and due to its fast spread, the World Health Organisation (WHO) declared a pandemic named COVID-19. Some variants of this virus were detected, including the Delta, which caused new waves of infections. This work uses an extended version of a SIRD model that includes vaccination effects to measure the impact of the Delta variant in three countries: Germany, Israel and Brazil. The calibrated models were able to reproduce the dynamics of the above countries. In addition, hypothetical scenarios were simulated to quantify the impact of vaccination and mitigation policies during the Delta wave. The results showed that the model could reproduce the complex dynamics observed in the different countries. The estimated increase of transmission rate due to the Delta variant was highest in Israel (7.9), followed by Germany (2.7) and Brazil (1.5). These values may support the hypothesis that people immunised against COVID-19 may lose their defensive antibodies with time since Israel, Germany, and Brazil fully vaccinated half of the population in March, July, and October. The scenario to study the impact of vaccination revealed relative reductions in the total number of deaths between 30% and 250%; an absolute reduction of 300 thousand deaths in Brazil due to vaccination during the Delta wave. The second hypothetical scenario revealed that mitigation policies saved up to 300 thousand Brazilians; relative reductions in the total number of deaths between 24% and 120% in the three analysed countries. Therefore, the results suggest that both vaccination and mitigation policies were crucial in decreasing the spread and the number of deaths during the Delta wave.

A Validated Mathematical Model of the Cytokine Release Syndrome in Severe COVID-19.

By June 2021, a new contagious disease, the Coronavirus disease 2019 (COVID-19), has infected more than 172 million people worldwide, causing more than 3.7 million deaths. Many aspects related to the interactions of the disease's causative agent, SAR2-CoV-2, and the immune response are not well understood: the multiscale interactions among the various components of the human immune system and the pathogen are very complex. Mathematical and computational tools can help researchers to answer these open questions about the disease. In this work, we present a system of fifteen ordinary differential equations that models the immune response to SARS-CoV-2. The model is used to investigate the hypothesis that the SARS-CoV-2 infects immune cells and, for this reason, induces high-level productions of inflammatory cytokines. Simulation results support this hypothesis and further explain why survivors have lower levels of cytokines levels than non-survivors.

The Quixotic Task of Forecasting Peaks of COVID-19: Rather Focus on Forward and Backward Projections.

Over the last months, mathematical models have been extensively used to help control the COVID-19 pandemic worldwide. Although extremely useful in many tasks, most models have performed poorly in forecasting the pandemic peaks. We investigate this common pitfall by forecasting four countries' pandemic peak: Austria, Germany, Italy, and South Korea. Far from the peaks, our models can forecast the pandemic dynamics 20 days ahead. Nevertheless, when calibrating our models close to the day of the pandemic peak, all forecasts fail. Uncertainty quantification and sensitivity analysis revealed the main obstacle: the misestimation of the transmission rate. Inverse uncertainty quantification has shown that significant changes in transmission rate commonly precede a peak. These changes are a key factor in forecasting the pandemic peak. Long forecasts of the pandemic peak are therefore undermined by the lack of models that can forecast changes in the transmission rate, i.e., how a particular society behaves, changes of mitigation policies, or how society chooses to respond to them. In addition, our studies revealed that even short forecasts of the pandemic peak are challenging. Backward projections have shown us that the correct estimation of any temporal change in the transmission rate is only possible many days ahead. Our results suggest that the distance between a change in the transmission rate and its correct identification in the curve of active infected cases can be as long as 15 days. This is intrinsic to the phenomenon and how it affects epidemic data: a new case is usually only reported after an incubation period followed by a delay associated with the test. In summary, our results suggest the phenomenon itself challenges the task of forecasting the peak of the COVID-19 pandemic when only epidemic data is available. Nevertheless, we show that exciting results can be obtained when using the same models to project different scenarios of reduced transmission rates. Therefore, our results highlight that mathematical modeling can help control COVID-19 pandemic by backward projections that characterize the phenomena' essential features and forward projections when different scenarios and strategies can be tested and used for decision-making.

Characterization of the COVID-19 pandemic and the impact of uncertainties, mitigation strategies, and underreporting of cases in South Korea, Italy, and Brazil.

By April 7th, 2020, the Coronavirus disease 2019 (COVID-19) has infected one and a half million people worldwide, accounting for over 80 thousand of deaths in 209 countries and territories around the world. The new and fast dynamics of the pandemic are challenging the health systems of different countries. In the absence of vaccines or effective treatments, mitigation policies, such as social isolation and lock-down of cities, have been adopted, but the results vary among different countries. Some countries were able to control the disease at the moment, as is the case of South Korea. Others, like Italy, are now experiencing the peak of the pandemic. Finally, countries with emerging economies and social issues, like Brazil, are in the initial phase of the pandemic. In this work, we use mathematical models with time-dependent coefficients, techniques of inverse and forward uncertainty quantification, and sensitivity analysis to characterize essential aspects of the COVID-19 in the three countries mentioned above. The model parameters estimated for South Korea revealed effective social distancing and isolation policies, border control, and a high number in the percentage of reported cases. In contrast, underreporting of cases was estimated to be very high in Brazil and Italy. In addition, the model estimated a poor isolation policy at the moment in Brazil, with a reduction of contact around 40%, whereas Italy and South Korea estimated numbers for contact reduction are at 75% and 90%, respectively. This characterization of the COVID-19, in these different countries under different scenarios and phases of the pandemic, supports the importance of mitigation policies, such as social distancing. In addition, it raises serious concerns for socially and economically fragile countries, where underreporting poses additional challenges to the management of the COVID-19 pandemic by significantly increasing the uncertainties regarding its dynamics.

A personalized computational model of edema formation in myocarditis based on long-axis biventricular MRI images.

BACKGROUND: Myocarditis is defined as the inflammation of the myocardium, i.e. the cardiac muscle. Among the reasons that lead to this disease, we may include infections caused by a virus, bacteria, protozoa, fungus, and others. One of the signs of the inflammation is the formation of edema, which may be a consequence of the interaction between interstitial fluid dynamics and immune response. This complex physiological process was mathematically modeled using a nonlinear system of partial differential equations (PDE) based on porous media approach. By combing a model based on Biot's poroelasticity theory with a model for the immune response we developed a new hydro-mechanical model for inflammatory edema. To verify this new computational model, T2 parametric mapping obtained by Magnetic Resonance (MR) imaging was used to identify the region of edema in a patient diagnosed with unspecific myocarditis. RESULTS: A patient-specific geometrical model was created using MRI images from the patient with myocarditis. With this model, edema formation was simulated using the proposed hydro-mechanical mathematical model in a two-dimensional domain. The computer simulations allowed us to correlate spatiotemporal dynamics of representative cells of the immune systems, such as leucocytes and the pathogen, with fluid accumulation and cardiac tissue deformation. CONCLUSIONS: This study demonstrates that the proposed mathematical model is a very promising tool to better understand edema formation in myocarditis. Simulations obtained from a patient-specific model reproduced important aspects related to the formation of cardiac edema, its area, position, and shape, and how these features are related to immune response.