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.

Enhanced optimization-based method for the generation of patient-specific models of Purkinje networks.

Cardiac Purkinje networks are a fundamental part of the conduction system and are known to initiate a variety of cardiac arrhythmias. However, patient-specific modeling of Purkinje networks remains a challenge due to their high morphological complexity. This work presents a novel method based on optimization principles for the generation of Purkinje networks that combines geometric and activation accuracy in branch size, bifurcation angles, and Purkinje-ventricular-junction activation times. Three biventricular meshes with increasing levels of complexity are used to evaluate the performance of our approach. Purkinje-tissue coupled monodomain simulations are executed to evaluate the generated networks in a realistic scenario using the most recent Purkinje/ventricular human cellular models and physiological values for the Purkinje-ventricular-junction characteristic delay. The results demonstrate that the new method can generate patient-specific Purkinje networks with controlled morphological metrics and specified local activation times at the Purkinje-ventricular junctions.

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.

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.

Stabilized hybrid discontinuous Galerkin finite element method for the cardiac monodomain equation.

Numerical methods for solving the cardiac electrophysiology model, which describes the electrical activity in the heart, are proposed. The model problem consists of a nonlinear reaction-diffusion partial differential equation coupled to systems of ordinary differential equations that describes electrochemical reactions in cardiac cells. The proposed methods combine an operator splitting technique for the reaction-diffusion equation with primal hybrid methods for spatial discretization considering continuous or discontinuous approximations for the Lagrange multiplier. A static condensation is adopted to form a reduced global system in terms of the multiplier only. Convergence studies exhibit optimal rates of convergence and numerical experiments show that the proposed schemes can be more efficient than standard numerical techniques commonly used in this context when preconditioned iterative methods are used for the solution of linear systems.

Preconditioned augmented Lagrangian formulation for nearly incompressible cardiac mechanics.

Computational modeling of the heart is a subject of substantial medical and scientific interest, which may contribute to increase the understanding of several phenomena associated with cardiac physiological and pathological states. Modeling the mechanics of the heart have led to considerable insights, but it still represents a complex and a demanding computational problem, especially in a strongly coupled electromechanical setting. Passive cardiac tissue is commonly modeled as hyperelastic and is characterized by quasi-incompressible, orthotropic, and nonlinear material behavior. These factors are known to be very challenging for the numerical solution of the model. The near-incompressibility is known to cause numerical issues such as the well-known locking phenomenon and ill-conditioning of the stiffness matrix. In this work, the augmented Lagrangian method is used to handle the nearly incompressible condition. This approach can potentially improve computational performance by reducing the condition number of the stiffness matrix and thereby improving the convergence of iterative solvers. We also improve the performance of iterative solvers by the use of an algebraic multigrid preconditioner. Numerical results of the augmented Lagrangian method combined with a preconditioned iterative solver for a cardiac mechanics benchmark suite are presented to show its improved performance.

Ectopic beats arise from micro-reentries near infarct regions in simulations of a patient-specific heart model.

Ectopic beats are known to be involved in the initiation of a variety of cardiac arrhythmias. Although their location may vary, ectopic excitations have been found to originate from infarct areas, regions of micro-fibrosis and other heterogeneous tissues. However, the underlying mechanisms that link ectopic foci to heterogeneous tissues have yet to be fully understood. In this work, we investigate the mechanism of micro-reentry that leads to the generation of ectopic beats near infarct areas using a patient-specific heart model. The patient-specific geometrical model of the heart, including scar and peri-infarct zones, is obtained through magnetic resonance imaging (MRI). The infarct region is composed of ischemic myocytes and non-conducting cells (fibrosis, for instance). Electrophysiology is captured using an established cardiac myocyte model of the human ventricle modified to describe ischemia. The simulation results clearly reveal that ectopic beats emerge from micro-reentries that are sustained by the heterogeneous structure of the infarct regions. Because microscopic information about the heterogeneous structure of the infarct regions is not available, Monte-Carlo simulations are used to identify the probabilities of an infarct region to behave as an ectopic focus for different levels of ischemia and different percentages of non-conducting cells. From the proposed model, it is observed that ectopic beats are generated when a percentage of non-conducting cells is near a topological metric known as the percolation threshold. Although the mechanism for micro-reentries was proposed half a century ago to be a source of ectopic beats or premature ventricular contractions during myocardial infarction, the present study is the first to reproduce this mechanism in-silico using patient-specific data.

Uniformization method for solving cardiac electrophysiology models based on the Markov-chain formulation.

This paper compares different numerical methods for the solution of myocyte models of cardiac electrophysiology. In particular, it presents how the technique called uniformization method substantially increases the stability of simple first-order methods such as Euler explicit method and Rush-Larsen (RL) method, for the solution of modern electrophysiology models that are based on continuous-time Markov chains (MCs) for the description of subcellular structures, such as ion channels. The MCs are often associated with stiff ordinary differential equations that severely limit the time step used by these traditional methods. By using the uniformization method, we could significantly increase the time steps for the solution of different cardiac electrophysiology models and improve the computational performance up to 150 times compared to the performance of Euler's and RL's methods.

Automatic Parameterization Strategy for Cardiac Electrophysiology Simulations.

Driven by recent advances in medical imaging, image segmentation and numerical techniques, computer models of ventricular electrophysiology account for increasingly finer levels of anatomical and biophysical detail. However, considering the large number of model parameters involved parameterization poses a major challenge. A minimum requirement in combined experimental and modeling studies is to achieve good agreement in activation and repolarization sequences between model and experiment or patient data. In this study, we propose basic techniques which aid in determining bidomain parameters to match activation sequences. An iterative parameterization algorithm is implemented which determines appropriate bulk conductivities which yield prescribed velocities. In addition, a method is proposed for splitting the computed bulk conductivities into individual bidomain conductivities by prescribing anisotropy ratios.