Computational Modelling of Cardiac Electrophysiology: from Cell to Bedside

Heart diseases are the leading cause of death worldwide, both for men and women. However, the ionic mechanisms underlying many cardiac arrhythmias and genetic disorders are not completely understood, thus leading to a limited efficacy of the current available therapies and leaving many open questions for cardiac electrophysiologists. On the other hand, experimental data availability is still a great issue in this field: most of the experiments are performed in vitro and/or using animal models (e.g. rabbit, dog and mouse), even when the final aim is to better understand the electrical behaviour of in vivo human heart either in physiological or pathological conditions. Computational modelling constitutes a primary tool in cardiac electrophysiology: in silico simulations, based on the available experimental data, may help to understand the electrical properties of the heart and the ionic mechanisms underlying a specific phenomenon. Once validated, mathematical models can be used for making predictions and testing hypotheses, thus suggesting potential therapeutic targets. This PhD thesis aims to apply computational cardiac modelling of human single cell action potential (AP) to three clinical scenarios, in order to gain new insights into the ionic mechanisms involved in the electrophysiological changes observed in vitro and/or in vivo. The first context is blood electrolyte variations, which may occur in patients due to different pathologies and/or therapies. In particular, we focused on extracellular Ca2+ and its effect on the AP duration (APD). The second context is haemodialysis (HD) therapy: in addition to blood electrolyte variations, patients undergo a lot of other different changes during HD, e.g. heart rate, cell volume, pH, and sympatho-vagal balance. The third context is human hypertrophic cardiomyopathy (HCM), a genetic disorder characterised by an increased arrhythmic risk, and still lacking a specific pharmacological treatment.

Viral load monitoring of antiretroviral therapy, cohort viral load and HIV transmission in Southern Africa: A mathematical modelling analysis

In low-income settings, treatment failure is often identified using CD4 cell count monitoring. Consequently, patients remain on a failing regimen, resulting in a higher risk of transmission. We investigated the benefit of routine viral load monitoring for reducing HIV transmission.

Using scenario tree modelling for targeted herd sampling to substantiate freedom from disease

BACKGROUND: In order to optimise the cost-effectiveness of active surveillance to substantiate freedom from disease, a new approach using targeted sampling of farms was developed and applied on the example of infectious bovine rhinotracheitis (IBR) and enzootic bovine leucosis (EBL) in Switzerland. Relevant risk factors (RF) for the introduction of IBR and EBL into Swiss cattle farms were identified and their relative risks defined based on literature review and expert opinions. A quantitative model based on the scenario tree method was subsequently used to calculate the required sample size of a targeted sampling approach (TS) for a given sensitivity. We compared the sample size with that of a stratified random sample (sRS) with regard to efficiency. RESULTS: The required sample sizes to substantiate disease freedom were 1,241 farms for IBR and 1,750 farms for EBL to detect 0.2% herd prevalence with 99% sensitivity. Using conventional sRS, the required sample sizes were 2,259 farms for IBR and 2,243 for EBL. Considering the additional administrative expenses required for the planning of TS, the risk-based approach was still more cost-effective than a sRS (40% reduction on the full survey costs for IBR and 8% for EBL) due to the considerable reduction in sample size. CONCLUSIONS: As the model depends on RF selected through literature review and was parameterised with values estimated by experts, it is subject to some degree of uncertainty. Nevertheless, this approach provides the veterinary authorities with a promising tool for future cost-effective sampling designs.

Resurgence of HIV infection among men who have sex with men in Switzerland: mathematical modelling study

Background New HIV infections in men who have sex with men (MSM) have increased in Switzerland since 2000 despite combination antiretroviral therapy (cART). The objectives of this mathematical modelling study were: to describe the dynamics of the HIV epidemic in MSM in Switzerland using national data; to explore the effects of hypothetical prevention scenarios; and to conduct a multivariate sensitivity analysis. Methodology/Principal Findings The model describes HIV transmission, progression and the effects of cART using differential equations. The model was fitted to Swiss HIV and AIDS surveillance data and twelve unknown parameters were estimated. Predicted numbers of diagnosed HIV infections and AIDS cases fitted the observed data well. By the end of 2010, an estimated 13.5% (95% CI 12.5, 14.6%) of all HIV-infected MSM were undiagnosed and accounted for 81.8% (95% CI 81.1, 82.4%) of new HIV infections. The transmission rate was at its lowest from 1995–1999, with a nadir of 46 incident HIV infections in 1999, but increased from 2000. The estimated number of new infections continued to increase to more than 250 in 2010, although the reproduction number was still below the epidemic threshold. Prevention scenarios included temporary reductions in risk behaviour, annual test and treat, and reduction in risk behaviour to levels observed earlier in the epidemic. These led to predicted reductions in new infections from 2 to 26% by 2020. Parameters related to disease progression and relative infectiousness at different HIV stages had the greatest influence on estimates of the net transmission rate. Conclusions/Significance The model outputs suggest that the increase in HIV transmission amongst MSM in Switzerland is the result of continuing risky sexual behaviour, particularly by those unaware of their infection status. Long term reductions in the incidence of HIV infection in MSM in Switzerland will require increased and sustained uptake of effective interventions.

Predicting perceived vulnerability for breast cancer among women with an average breast cancer risk

OBJECTIVES: The purpose of the present study was to investigate predictors of perceived vulnerability for breast cancer in women with an average risk for breast cancer. On the basis of empirical findings that suggested which variables might be associated with perceived vulnerability for breast cancer, we investigated whether knowledge of breast cancer risk factors, cancer worry, intrusions about breast cancer, optimism about not getting cancer and perceived health status have a predictive value for perceived breast cancer vulnerability. DESIGN: In a 3-step approach, we recruited 292 women from the general public in Germany who had neither a family history of breast cancer nor breast cancer themselves. After receiving an initial informational letter about study objectives, the women were interviewed by telephone and then asked to fill in a self-administered questionnaire. METHODS: We used structural equation modelling and hypothesized that each of the included variables has a direct influence on perceived vulnerability for breast cancer. RESULTS: We found a valid model with acceptable fit indices. Optimism about not getting cancer, intrusions about breast cancer and women's perceived health status explained 32% of the variance of perceived vulnerability for breast cancer. Cancer worry and knowledge about breast cancer did not influence perceived vulnerability for breast cancer. CONCLUSION: Perceived vulnerability for breast cancer is associated with health-related variables more than with knowledge about breast cancer risk factors.

Effectiveness and cost-effectiveness of traditional and new partner notification technologies for curable sexually transmitted infections: observational study, systematic reviews and mathematical modelling

BACKGROUND Partner notification is essential to the comprehensive case management of sexually transmitted infections. Systematic reviews and mathematical modelling can be used to synthesise information about the effects of new interventions to enhance the outcomes of partner notification. OBJECTIVE To study the effectiveness and cost-effectiveness of traditional and new partner notification technologies for curable sexually transmitted infections (STIs). DESIGN Secondary data analysis of clinical audit data; systematic reviews of randomised controlled trials (MEDLINE, EMBASE and Cochrane Central Register of Controlled Trials) published from 1 January 1966 to 31 August 2012 and of studies of health-related quality of life (HRQL) [MEDLINE, EMBASE, ISI Web of Knowledge, NHS Economic Evaluation Database (NHS EED), Database of Abstracts of Reviews of Effects (DARE) and Health Technology Assessment (HTA)] published from 1 January 1980 to 31 December 2011; static models of clinical effectiveness and cost-effectiveness; and dynamic modelling studies to improve parameter estimation and examine effectiveness. SETTING General population and genitourinary medicine clinic attenders. PARTICIPANTS Heterosexual women and men. INTERVENTIONS Traditional partner notification by patient or provider referral, and new partner notification by expedited partner therapy (EPT) or its UK equivalent, accelerated partner therapy (APT). MAIN OUTCOME MEASURES Population prevalence; index case reinfection; and partners treated per index case. RESULTS Enhanced partner therapy reduced reinfection in index cases with curable STIs more than simple patient referral [risk ratio (RR) 0.71; 95% confidence interval (CI) 0.56 to 0.89]. There are no randomised trials of APT. The median number of partners treated for chlamydia per index case in UK clinics was 0.60. The number of partners needed to treat to interrupt transmission of chlamydia was lower for casual than for regular partners. In dynamic model simulations, > 10% of partners are chlamydia positive with look-back periods of up to 18 months. In the presence of a chlamydia screening programme that reduces population prevalence, treatment of current partners achieves most of the additional reduction in prevalence attributable to partner notification. Dynamic model simulations show that cotesting and treatment for chlamydia and gonorrhoea reduce the prevalence of both STIs. APT has a limited additional effect on prevalence but reduces the rate of index case reinfection. Published quality-adjusted life-year (QALY) weights were of insufficient quality to be used in a cost-effectiveness study of partner notification in this project. Using an intermediate outcome of cost per infection diagnosed, doubling the efficacy of partner notification from 0.4 to 0.8 partners treated per index case was more cost-effective than increasing chlamydia screening coverage. CONCLUSIONS There is evidence to support the improved clinical effectiveness of EPT in reducing index case reinfection. In a general heterosexual population, partner notification identifies new infected cases but the impact on chlamydia prevalence is limited. Partner notification to notify casual partners might have a greater impact than for regular partners in genitourinary clinic populations. Recommendations for future research are (1) to conduct randomised controlled trials using biological outcomes of the effectiveness of APT and of methods to increase testing for human immunodeficiency virus (HIV) and STIs after APT; (2) collection of HRQL data should be a priority to determine QALYs associated with the sequelae of curable STIs; and (3) standardised parameter sets for curable STIs should be developed for mathematical models of STI transmission that are used for policy-making. FUNDING The National Institute for Health Research Health Technology Assessment programme.

Cost-effectiveness of point-of-care viral load monitoring of antiretroviral therapy in resource-limited settings: mathematical modelling study

BACKGROUND Monitoring of HIV viral load in patients on combination antiretroviral therapy (ART) is not generally available in resource-limited settings. We examined the cost-effectiveness of qualitative point-of-care viral load tests (POC-VL) in sub-Saharan Africa. DESIGN Mathematical model based on longitudinal data from the Gugulethu and Khayelitsha township ART programmes in Cape Town, South Africa. METHODS Cohorts of patients on ART monitored by POC-VL, CD4 cell count or clinically were simulated. Scenario A considered the more accurate detection of treatment failure with POC-VL only, and scenario B also considered the effect on HIV transmission. Scenario C further assumed that the risk of virologic failure is halved with POC-VL due to improved adherence. We estimated the change in costs per quality-adjusted life-year gained (incremental cost-effectiveness ratios, ICERs) of POC-VL compared with CD4 and clinical monitoring. RESULTS POC-VL tests with detection limits less than 1000 copies/ml increased costs due to unnecessary switches to second-line ART, without improving survival. Assuming POC-VL unit costs between US$5 and US$20 and detection limits between 1000 and 10,000 copies/ml, the ICER of POC-VL was US$4010-US$9230 compared with clinical and US$5960-US$25540 compared with CD4 cell count monitoring. In Scenario B, the corresponding ICERs were US$2450-US$5830 and US$2230-US$10380. In Scenario C, the ICER ranged between US$960 and US$2500 compared with clinical monitoring and between cost-saving and US$2460 compared with CD4 monitoring. CONCLUSION The cost-effectiveness of POC-VL for monitoring ART is improved by a higher detection limit, by taking the reduction in new HIV infections into account and assuming that failure of first-line ART is reduced due to targeted adherence counselling.

A review of dendrogeomorphological research applied to flood risk analysis in Spain

Over the last forty years, applying dendrogeomorphology to palaeoflood analysis has improved estimates of the frequency and magnitude of past floods worldwide. This paper reviews the main results obtained by applying dendrogeomorphology to flood research in several case studies in Central Spain. These dendrogeomorphological studies focused on the following topics: (1) anatomical analysis to understand the physiological response of trees to flood damage and improve sampling efficiency; (2) compiling robust flood chronologies in ungauged mountain streams, (3) determining flow depth and estimating flood discharge using two-dimensional hydraulic modelling, and comparing them with other palaeostage indicators; (4) calibrating hydraulic model parameters (i.e. Manning roughness); and (5) implementing stochastic-based, cost–benefit analysis to select optimal mitigation measures. The progress made in these areas is presented with suggestions for further research to improve the applicability of dendrogeochronology to palaeoflood studies. Further developments will include new methods for better identification of the causes of specific types of flood damage to trees (e.g. tilted trees) or stable isotope analysis of tree rings to identify the climatic conditions associated with periods of increasing flood magnitude or frequency.

Risk adjustment in aging societies

Background In Switzerland, age is the predominant driver of solidarity transfers in risk adjustment (RA). Concerns have been voiced regarding growing imbalances in cost sharing between young and old insured due to demographic changes (larger fraction of elderly >65 years and rise in average age). Particularly young adults aged 19–25 with limited incomes have to shoulder increasing solidarity burdens. Between 1996 and 2011, monthly intergenerational solidarity payments for young adults have doubled from CHF 87 to CHF 182, which corresponds to the highest absolute transfer increase of all age groups. Results By constructing models for age-specific RA growth and for calculating the lifetime sum of RA transfers we investigated the causes and consequences of demographic changes on RA payments. The models suggest that the main driver for RA increases in the past was below average health care expenditure (HCE) growth in young adults, which was only half as high (average 2% per year) compared with older adults (average 4% per year). Shifts in age group distributions were only accountable for 2% of the CHF 95 rise in RA payments. Despite rising risk adjustment debts for young insured the balance of lifetime transfers remains positive as long as HCE growth rates are greater than the discount rate used in this model (3%). Moreover, the life-cycle model predicts that the lifetime rate of return on RA payments may even be further increased by demographic changes. Nevertheless, continued growth of RA contributions may overwhelm vulnerable age groups such as young adults. We therefore propose methods to limit the burden of social health insurance for specific age groups (e.g. young adults in Switzerland) by capping solidarity payments. Conclusions Taken together, our mathematical modelling framework helps to gain a better understanding of how demographic changes interact with risk adjustment and how redistribution of funds between age groups can be controlled without inducing further selection incentives. Those methods can help to construct more equitable systems of health financing in light of population aging.

A physical approach on flood risk vulnerability of buildings

The design of efficient hydrological risk mitigation strategies and their subsequent implementation relies on a careful vulnerability analysis of the elements exposed. Recently, extensive research efforts were undertaken to develop and refine empirical relationships linking the structural vulnerability of buildings to the impact forces of the hazard processes. These empirical vulnerability functions allow estimating the expected direct losses as a result of the hazard scenario based on spatially explicit representation of the process patterns and the elements at risk classified into defined typological categories. However, due to the underlying empiricism of such vulnerability functions, the physics of the damage-generating mechanisms for a well-defined element at risk with its peculiar geometry and structural characteristics remain unveiled, and, as such, the applicability of the empirical approach for planning hazard-proof residential buildings is limited. Therefore, we propose a conceptual assessment scheme to close this gap. This assessment scheme encompasses distinct analytical steps: modelling (a) the process intensity, (b) the impact on the element at risk exposed and (c) the physical response of the building envelope. Furthermore, these results provide the input data for the subsequent damage evaluation and economic damage valuation. This dynamic assessment supports all relevant planning activities with respect to a minimisation of losses, and can be implemented in the operational risk assessment procedure.

Viral load versus CD4⁺ monitoring and 5-year outcomes of antiretroviral therapy in HIV-positive children in Southern Africa: a cohort-based modelling study.

OBJECTIVES Many paediatric antiretroviral therapy (ART) programmes in Southern Africa rely on CD4⁺ to monitor ART. We assessed the benefit of replacing CD4⁺ by viral load monitoring. DESIGN A mathematical modelling study. METHODS A simulation model of HIV progression over 5 years in children on ART, parameterized by data from seven South African cohorts. We simulated treatment programmes with 6-monthly CD4⁺ or 6- or 12-monthly viral load monitoring. We compared mortality, second-line ART use, immunological failure and time spent on failing ART. In further analyses, we varied the rate of virological failure, and assumed that the rate is higher with CD4⁺ than with viral load monitoring. RESULTS About 7% of children were predicted to die within 5 years, independent of the monitoring strategy. Compared with CD4⁺ monitoring, 12-monthly viral load monitoring reduced the 5-year risk of immunological failure from 1.6 to 1.0% and the mean time spent on failing ART from 6.6 to 3.6 months; 1% of children with CD4⁺ compared with 12% with viral load monitoring switched to second-line ART. Differences became larger when assuming higher rates of virological failure. When assuming higher virological failure rates with CD4⁺ than with viral load monitoring, up to 4.2% of children with CD4⁺ compared with 1.5% with viral load monitoring experienced immunological failure; the mean time spent on failing ART was 27.3 months with CD4⁺ monitoring and 6.0 months with viral load monitoring. Conclusion: Viral load monitoring did not affect 5-year mortality, but reduced time on failing ART, improved immunological response and increased switching to second-line ART.

Modelling susceptibility of grassland soil to macropore flow

Investigating preferential flow, including macropore flow, is crucial to predicting and preventing point sources of contamination in soil, for example in the vicinity of pumping wells. With a view to advancing groundwater protection, this study aimed (i) to quantify the strength of macropore flow in four representative natural grassland soils on the Swiss plateau, and (ii) to define the parameters that significantly control macropore flow in grassland soil. For each soil type we selected three measurement points on which three successive irrigation experiments were carried out, resulting in a total of 36 irrigations. The strength of macropore flow, parameterized as the cumulated water volume flowing from macropores at a depth of 1 m in response to an irrigation of 60 mm h−1 intensity and 1 h duration, was simulated using the dual-permeability MACRO model. The model calibration was based on the key soil parameters and fine measurements of water content at different depths. Modelling results indicate high performance of macropore flow in all investigated soil types except in gleysols. The volume of water that flowed from macropores and was hence expected to reach groundwater varied between 81% and 94% in brown soils, 59% and 67% in para-brown soils, 43% and 56% in acid brown soils, and 22% and 35% in gleysols. These results show that spreading pesticides and herbicides in pumping well protection zones poses a high risk of contamination and must be strictly prohibited. We also found that organic carbon content was not correlated with the strength of macropore flow, probably due to its very weak variation in our study, while saturated water content showed a negative correlation with macropore flow. The correlation between saturated hydraulic conductivity (Ks) and macropore flow was negative as well, but weak. Macropore flow appears to be controlled by the interaction between the bulk density of the uppermost topsoil layer (0–0.10 m) and the macroporosity of the soil below. This interaction also affects the variations in Ks and saturated water content. Further investigations are needed to better understand the combined effect of all these processes including the exchange between micropore and macropore domains.