=Paper=
{{Paper
|id=Vol-2086/AICS2017_paper_13
|storemode=property
|title=A Computational Lymph Tissue Model for Long Term HIV Infection Progression and Immune Fitness
|pdfUrl=https://ceur-ws.org/Vol-2086/AICS2017_paper_13.pdf
|volume=Vol-2086
|authors=Andreas Hillmann,Martin Crane,Heather Ruskin
|dblpUrl=https://dblp.org/rec/conf/aics/HillmannCR17
}}
==A Computational Lymph Tissue Model for Long Term HIV Infection Progression and Immune Fitness==
A computational lymph tissue model for long term HIV infection progression and immune fitness Andreas Hillmann1, Martin Crane1 and Heather J. Ruskin1 1 Advanced Research Computing Centre for Complex Systems Modelling (ARC-SYM), School of Computing, Dublin City University, Dublin, Ireland andreas.hillmann2@mail.dcu.ie Abstract. Given the complexity of biological systems, modelling the immune response to the spread of infectious disease is non-trivial and has been the subject of considerable computational efforts. HIV infection, if untreated, leads to a slow depletion of the immune system over many years, resulting in terminal AIDS disease and death. Despite decades of research, the biological mechanisms un- derpinning the immune system response and implicit in its impairment remain a subject of discussion. Recent biological findings indicate the importance of func- tional lymphatic tissue in maintaining homoeostasis of the immune system, which is vital for stability. HIV infection-related immune activation leads to chronic inflammation of this tissue and formation of non-functional scar tissue with detrimental effects on homeostasis. We propose a computational model based on the Cellular Automata formalism to quantify these effects in a simulated fraction of a lymphatic tissue. We also include effects of antiretroviral treatment and highlight implications for system representation through larger-scale simu- lations, as well as optimizations of treatment schedules. Keywords: HIV, computational model, Cellular Automata, antiretroviral treat- ment 1 Introduction 1.1 HIV infection Acquired Human Immunodeficiency Syndrome (AIDS) remains a global health con- cern even three decades or more after the discovery of the Human Immunodeficiency Virus as its root cause, [1]. The virus preferentially infects white blood cells of a certain type (CD4+), which are a crucial part of the immune response, leading to failure of the immune system approximately 10 years after initial infection. Potent drugs, termed An- tiretroviral Therapy (ART), have been developed, which aim at the suppression of viral activities. However, these medications do not provide a cure for the disease, due to the persistence of impermeable viral reservoirs in the body; instead they have to be admin- istered life-long to maintain sufficient viral suppression, [2]. Daily administration of multiple drugs can become a problem and interruptions are common, especially in re- source-limited settings, where drug availability is not guaranteed or social factors in- terfere with adherence, [3]. Moreover, Structured Treatment Interruptions (STIs) have been suggested to enable recovery from drug toxicities, to promote resilience of the immune system and to reduce drug costs, [4]. Though initial clinical studies on STIs provided encouraging evidence regarding safety of interrupted treatment, [5, 6], a subsequent large-scale study (SMART) indicated potential harmful effects, [7]. Despite these effects have not been characterized in detail and may have been related to study design, [8], the results led to a general loss of interest in assessing treatment interruptions in a clinical setting. Ques- tions regarding safe margins of interruption regimen and the nature of detrimental ef- fects thus remain open. 1.2 Lymphatic tissue involvement It is well-known that viral replication occurs predominantly in lymphatic tissue, where most of the susceptible cells are located, [9]. The lymphatic system is a widespread network of nodes and ducts which is crucial for immune response but also for the spread of infected cells throughout the body. Recent advances in in-vivo microscopy have enabled determination of the detailed structure of a certain kind of lymphatic tissue, the Fibroblastic Reticular Cell Network (FRCn). Its mesh-like structure has been reported to mimic a small world topology with lattice-like properties, [10]. This view is complemented by histological imaging show- ing that the Fibroblastic Reticular Cells (FRCs) are fixed to a scaffolding of collagen fibers. During inflammation, this formation of collagen fibers is greatly increased, [11], leading to collagen deposits (see Fig. 1) which block off FRCs from the chemical sig- nals (cytokines) which maintain their viability, leading to depletion of these cells and formation of non-functional collagen tissue. Antiretroviral treatment has been shown to stop the formation of the scar tissue but not to enable its reversal, [12]. The formation of FRCn has not yet been observed in vivo due to the long timescales involved (weeks to years), however. Fig. 1. Distinct stages of collagenation of lymphatic tissue for uninfected, early, and late stage HIV infection (from left). Green staining denotes functional tissue, red non-functional collagen. Image source: [12]. 1.3 Modelling approaches Various models have been proposed in the HIV context to help understanding the spread of the virus between, [13], and within hosts, [14]. The majority of these models have found to be based on systems of Ordinary Differential Equations (ODEs), [15], representing a variant of the SIR model originally used to model spread of communi- cable diseases throughout populations. In context of within-host models, each DE de- notes the change of a quantitative property related to HIV like viral load, healthy and infected cell concentrations over time. Starting from basic forms with two or three equa- tions, models have been extended to cover additional features, like effects of antiretro- viral treatment or immune response (resulting in systems of 20 or more equations). Simpler forms allow further mathematical analysis, [16], whereas more complex forms may be used for numerical simulation studies, [17]. For a more detailed review on ODE based HIV-modelling see [14] and references therein. In general, ODE-based approaches assume virus und susceptible cells being in an environment satisfying the conditions of well-mixed. That assumption fits loosely best to the bloodstream, which initially has been considered a main source of viral replica- tion, [14]. However, following the increasing availability of data on tissue involvement, the suitability of the well-mixed assumption has been questioned recently, [18]. Treatment interruptions and optimizations have also been subject to modelling stud- ies using ODEs, [15]. A popular approach, besides pure simulation studies, [19], con- siders an optimal control problem, [20–22], to obtain treatment regimen which keep biological properties, like viral load or healthy cell count within predefined bounds which are considered to be non-critical. This view of the biological system as a kind of thermostat was especially useful with only limited data on anatomical compartments other than the bloodstream available. However, as clinical studies have indicated the complexity of HIV infection progression, consideration of alternative modelling ap- proaches might be worthwhile. We therefore propose a computational model, which considers both spatial structure of lymphatic tissue and the effects of inflammation, to examine long-term development and to allow predictions to be made on effects of treatment interruptions, where clinical data is not available for reasons mentioned in section 1. Our model presents an exten- sion of the approach reported in [15] with improved biological realism obtained from drawing upon recent studies. 2 Methods - Numerical simulations 2.1 Surrounding conditions For our model we assume two effects to be responsible for the depletion of CD4+ cells, as reported by clinical studies, [11]. • HIV-induced direct killing of susceptible CD4+ cells • Impairment of survival through depletion of FRCs Both factors influence the concentration of CD4+ cells in the body, with the time-scale of the former relatively short (days to weeks) while the latter is effective at longer time- scales (months to years). However due to regenerative capabilities of CD4+ cells, the first effect is reversible under antiretroviral treatment, while only minimal regeneration has been observed, [12], for the second. We apply a bottom-up modelling approach, by using the Stochastic Cellular Autom- ata (SCA) formalism, i.e. a stochastic extension of Cellular Automata (CA), [23]. The CA formalism has successfully been applied for modelling environments involving fixed matrices, including HIV, [24]. It utilizes a regular grid or lattice structure, where each site represents a single cell and its neighborhood the cell-to-cell contacts. We fur- ther assume that each cell has a connection to its 4 immediate neighbors as suggested from analysis of in-vivo imaging data, [10]. In terms of Cellular Automata computation, this corresponds to a von Neumann neighborhood. An average lymph node measures about 1 cm in diameter and is of approximately spherical shape. From murine experiments, it is known that in FRCn, mean separation of the centers of mass of single cells is ~ 23 µm, [10]. A human lymph node, would therefore consist of around 40 Million FRCs. For our simulations, however, we con- sider a smaller fraction, namely a 2-dimensional slice through a lymph node, consisting approximately of 250,000 single cells. This slice represents less than 1% of the total volume, however this simplified view is acceptable in this case since we do not aim at reproducing an anatomically correct representation of the organ (see e.g. [25]). From clinical data, we know that collagenation mainly starts from connections in the blood circulation (High Endothelial Venules). These connections appear to be scattered randomly throughout the tissue, [12]. We assume that collagenation increase occurs predominantly through elongation of existing collagen fibers as suggested by those data available to date on this aspect. 2.2 Model rules Before simulated infection occurs, the lattice is populated with both healthy FRCs and co-located CD4+ cells. The CD4+ cells are assumed to maintain equilibrium (homeo- stasis) and their concentrations in lymphatic tissue to be stable, [26]. However, HIV infection has a two-fold effect in our model: destroying CD4+ cells through viral infec- tion and mechanisms of cellular apoptosis and pyroptosis, and through upregulation of the immune system, [26]. The upregulated immune response in turn causes increased collagenation of lymphatic tissue. The fraction of CD4+ cells which is co-located with collagenated tissue, tends to exhibit poorer survival than usual (in healthy tissue), in addition to the depletion due to HIV infection. Treatment, however reduces the latter, as well as inflammation stimula- tion, since viral load is reduced. CD4+ cells are therefore regenerated and the colla- genation processes in lymphatic tissue are slowed down. Upon reinstatement of upreg- ulation of the immune system due to a treatment interruption, however, homeostasis disruption occurs again. In terms of model rules, no direct death and/or inflammation takes place other than by viral infection. This context gives rise to the state transition model shown in Fig. 1, which contains the four possible combinations of CD4+ cells and collagen status as well as permitted changes. Fig. 2. Block diagram of state changes of the SCA with probabilities of transitions, where appli- cable (for details see Table 1.). H: healthy CD4+ cell; HC: healthy CD4+ cell co-located with collagen; D: dead CD4+ cell (empty); C: empty site co-located with collagen. As mentioned in the previous section, collagenation originates from the High Endothe- lial Venules, where CD4+ cells enter the lymph nodes from the blood circulatory sys- tem. The assumption, that the blood circulation acts as a means to spread infection throughout the body, allows us to calibrate the model using available clinical data on CD4+ cell long-term degradation under HIV infection, obtained from blood samples, [27]. Table 1 lists the model parameters together with values obtained by calibration. Table 1. Model parameter value selection Parameter Description Default value PT Treatment effectiveness 0.9 PINF Infection induced inflammation, drives collagenation 0.1 PREG CD4+ cell regeneration 0.3 PCYT CD4+ cell death due to signaling (cytokines) 0.9 PAP CD4+ cell death through infection (apoptosis) 0.1 N Collagenated neighbors 1-4 2.3 Simulation setup Using the model parametrization in this table, we aim to test the effect of different patterns of unstructured interrupted treatment in the long term on collagenation (and thus immune system degradation). To simulate interrupted treatment, we assume ‘av- erage’ adherence (the probability that the simulated patient takes the medication each week). As the treatment initiation timing has also been deemed to be important in clin- ical studies, [28], we allow for this in the simulations as well. Table 2. summarizes the parameter values chosen for the simulation runs. After the treatment period (corre- sponding to ~ 5 years), the increase in collagenation compared to that evident at the start of treatment is determined. Table 2. Simulation parameters Description Value range Time step 1 week Treatment initiation 0, 50, 100, 150, 200 time steps Treatment adherence 0%, 20%, 40%, 60%, 80% 100% Treatment period duration 250 time steps Total simulation length 500 time steps (~ 10 years) In a second experiment, we aim to investigate the effects of structured treatment inter- ruptions. Thus, we adopt the approach of CD4+ guidance, which has been pursued in recent clinical trials, [4].This approach uses a threshold value of CD4+ cell count in a patient as a surrogate marker to decide on interruption or (re-)initiation of antiretroviral treatment. For our simulations, we apply the limits used in two major clinical trials, [7, 29] and again observe the outcomes after 500 time steps. 3 Results and Discussion 3.1 Unstructured interruptions An example model run is depicted in Fig. 3, exhibiting unstructured treatment inter- ruptions due to incomplete adherence. In general, the model runs demonstrate good correspondence to clinical results, [6]. For an infection simulated to occur at time step 0, the count of CD4+ cells drops rapidly, (characteristic of the acute phase). The cell count continues to drop at a slower rate until initiation of simulated treatment, which leads to a significant recovery in number of CD4+ cells. Incomplete adherence to the treatment regime causes temporary drops in cell numbers (signified by the irregular saw tooth pattern), as well as an overall decline in the long term. This decline accelerates once the treatment period concludes. Collagenation, however, exhibits a constant in- crease, albeit with less-pronounced slope during treatment. In Fig. 4 we show the loss of CD4+ cells and the increase in collagenation during treatment periods. The results illustrate the influence of treatment initiation timing and of adherence to the regimen. We observe that the latter has a fundamental impact on long-term progression of collagenation and, consequently, on immune system fitness. Collagen CD4+ cells 100% 90% 80% tissue area covered 70% 60% 50% 40% 30% 20% 10% 0% 0 100 200 300 400 500 time step [weeks] Fig. 3. Time course of collagenation (blue line) and population of CD4+ cells (red line) during a period of interrupted treatment (100 to 350 time steps; 0.8 adherence) 0.0 0.2 0.4 0.6 0.8 1.0 80% 70% collogenation increase 60% 50% 40% 30% 20% 10% 0% 0 50 100 150 200 time step of treatment initiation [weeks] Fig. 4. Effect of different adherence levels (lines) and treatment initiation on collagenation in- crease at end of treatment period. Treatment initiation timing has a less-pronounced effect. Results show that, if treatment is introduced during the early phases of infection (within the first 2 years), initiation time has greater impact than for later stages. However, with high adherence (80% and above), collagenation can be kept within acceptable bounds even for late treatment in- itiation. As an acceptable outcome, we take a total collagenated area of less than 75% at the end of the simulated 10 years. This level has been determined as a threshold from clinical data, [12], above which AIDS disease is imminent. 3.2 Structured interruptions Example results of our second experiment concerning structured treatment interrup- tions are shown in Fig. 5. Using the same thresholds for CD4+ guidance as two large- scale clinical studies (SMART, LOTTI, see Table 3.), simulated data is in good corre- spondence with clinical data. The two sets of clinical data indicate that counts of CD4+ cells are constantly decreasing during the course of these studies, reflected by the sim- ulation data. Since the duration of the clinical trials was limited due to practical reasons, we extended simulation times to allow for an estimation of long-term effects of the interruption thresholds used in the studies. LOTTI SMART 1200 CD4+ count [cells/µl] 1000 800 600 400 200 0 10 20 30 40 50 study duration [months] Fig. 5. Comparison of clinical data from two CD4+ guided STI studies (LOTTI, SMART; dots) with model output (lines) from corresponding simulation experiments. On this behalf, absolute values of CD4+ count limits from clinical literature were nor- malized for use with the model by assuming an initial cell count of 1000 cells/µl.and results were again obtained after 500 time steps (~ 10 years). The results (see Table 3.) show a relatively small difference both in terms of collagenation and CD4+ counts. These results are remarkable in a way since the SMART study strongly indicated ad- verse effects of interrupted treatment, whereas the LOTTI study did not find significant differences in this regard compared to continuous treatment. The different outcomes of our simulations in terms of collagenation and CD4+ count may therefore account for some threshold between non-critical progression and adverse effects. Table 3. Experimental results for CD4+ guided structured treatment interruptions CD4+ count limits [cells/µl] Simulation results (~ 10 years) Clinical study Interruption Resumption Collagenation CD4+ count SMART [7] 350 250 97% 260 LOTTI [29] 700 350 81% 380 4 Conclusion We have proposed a simple Stochastic Cellular Automata model, which describes the effects of HIV infection and antiretroviral therapy on both CD4+ cell counts and colla- gen deposition in lymphatic tissue. Additionally, the model can predict possible out- comes of incomplete antiretroviral treatment, both unstructured and structured (difficult to obtain from controlled clinical studies). Simulated adherence patterns indicated the importance for long-term immunological fitness, of overall adherence to medication, and of treatment initiation timing during the early phases of infection, with regard to improved management. These experiments on antiretroviral treatment timing and ad- herence explore the implications for outcomes of Structured Treatment Interruptions, formerly discounted due to lack of long term predictions, despite some evident patient benefits. Further steps might include an extension of our approach to evaluate strategies for an improved monitoring of antiretroviral therapy by helping to identify patterns, like certain concentration thresholds, which may indicate a future loss of viral control given poor adherence by employing AI techniques. Applicability in real-life settings may fur- therly be improved by incorporation of realistic sampling intervals to obtain biological properties from blood or tissue with least patient stress. Regarding structured interruptions, numerical simulation in conjunction with appro- priate optimization techniques could be utilized for identifying optimal treatment regi- men with improved biological realism compared to ODE-based approaches. A conceiv- able way would be the use of evolutionary algorithms to optimize treatment on/off pat- terns. In addition, the proposed model offers further development potential, not least as it represents to date only a small part of human lymphatic tissue. A realistic simulation size even for one lymph node requires around 100 times more lattice sites than the current model. This problem can be alleviated using state-of-the-art parallelization ap- proaches, [30], to simulate different regions and this forms part of ongoing investiga- tion. Unfortunately, clinical data for validation of model performance are sparse at present due to ethical considerations as well as measurement limitations, but recent advances in biological imaging and in-vivo microscopy may well provide additional insights and enable refinement of key parameters. In summary, computational modelling of lymphatic tissue features involving cell- to-cell connections offers considerable potential for further investigation and optimiza- tion of long-term HIV treatment. Promising foundations using bottom-up modelling have already been laid for investigating infectious diseases. We think it worthwhile to pursue further, given growing availability of appropriate experimental and clinical data. References 1. Barré-Sinoussi, F., Chermann, J.C., Rey, F., Nugeyre, M.T., Chamaret, S., Gruest, J., Dauguet, C., Axler-Blin, C., Vézinet-Brun, F., Rouzioux, C., Rozenbaum, W., Montagnier, L.: Isolation of a T-lymphotropic retrovirus from a patient at risk for acquired immune deficiency syndrome (AIDS). Science. 220, 868–871 (1983). 2. WHO: Consolidated guidelines on the use of antiretroviral drugs for treating and preventing HIV infection, http://www.who.int/hiv/pub/guidelines/arv2013/download/en/. 3. Mann, M., Lurie, M.N., Kimaiyo, S., Kantor, R.: Effects of political conflict-induced treatment interruptions on HIV drug resistance. AIDS Rev. 15, 15–24 (2013). 4. Benson, C.A.: Structured treatment interruptions--new findings. Top. HIV Med. 14, 107–11 (2006). 5. Rosenberg, E.S., Altfeld, M., Poon, S.H., Phillips, M.N., Wilkes, B.M., Eldridge, R.L., Robbins, G.K., D’Aquila, R.T., Goulder, P.J., Walker, B.D.: Immune control of HIV-1 after early treatment of acute infection. Nature. 407, 523–526 (2000). 6. Dybul, M., Chun, T.W., Yoder, C., Hidalgo, B., Belson, M., Hertogs, K., Larder, B., Dewar, R.L., Fox, C.H., Hallahan, C.W., Justement, J.S., Migueles, S.A., Metcalf, J.A., Davey, R.T., Daucher, M., Pandya, P., Baseler, M., Ward, D.J., Fauci, A.S.: Short-cycle structured intermittent treatment of chronic HIV infection with highly active antiretroviral therapy: effects on virologic, immunologic, and toxicity parameters. Proc. Natl. Acad. Sci. U. S. A. 98, 15161–6 (2001). 7. Siegel, L., El-Sadr, W.: New perspectives in HIV treatment interruption: The SMART study. PRN Noteb. 11, 8–9 (2006). 8. Hirschel, B., Flanigan, T.: Is it smart to continue to study treatment interruptions? AIDS. 23, 757–759 (2009). 9. Pantaleo, G., Graziosi, C., Demarest, J.F., Butini, L., Montroni, M., Fox, C.H., Orenstein, J.M., Kotler, D.P., Fauci, A.S.: HIV infection is active and progressive in lymphoid tissue during the clinically latent stage of disease. Nature. 362, 355–358 (1993). 10. Novkovic, M., Onder, L., Cupovic, J., Abe, J., Bomze, D., Cremasco, V., Scandella, E., Stein, J. V., Bocharov, G., Turley, S.J., Ludewig, B.: Topological Small-World Organization of the Fibroblastic Reticular Cell Network Determines Lymph Node Functionality. PLOS Biol. 14, e1002515 (2016). 11. Schacker, T.W., Reilly, C., Beilman, G.J., Taylor, J., Skarda, D., Krason, D., Larson, M., Haase, A.T.: Amount of lymphatic tissue fibrosis in HIV infection predicts magnitude of HAART-associated change in peripheral CD4 cell count. AIDS. 19, 2169– 2171 (2005). 12. Zeng, M., Southern, P.J., Reilly, C.S., Beilman, G.J., Chipman, J.G., Schacker, T.W., Haase, A.T.: Lymphoid Tissue Damage in HIV-1 Infection Depletes Naïve T Cells and Limits T Cell Reconstitution after Antiretroviral Therapy. PLoS Pathog. 8, e1002437 (2012). 13. May, R.M., Anderson, R.M.: Transmission dynamics of HIV infection. Nature. 326, 137–142 (1987). 14. Perelson, A.S., Nelson, P.W.: Mathematical Analysis of HIV-1 Dynamics in Vivo. SIAM Rev. 41, 3–44 (1999). 15. Hillmann, A., Crane, M., Ruskin, H.J.: HIV models for treatment interruption: Adaptation and comparison. Phys. A Stat. Mech. its Appl. 483, 44–56 (2017). 16. Pawelek, K. a., Liu, S., Pahlevani, F., Rong, L.: A model of HIV-1 infection with two time delays: Mathematical analysis and comparison with patient data. Math. Biosci. 235, 98–109 (2012). 17. Sun, Q., Min, L.: Dynamics Analysis and Simulation of a Modified HIV Infection Model with a Saturated Infection Rate. Comput. Math. Methods Med. 2014, 1–14 (2014). 18. Cohen, J.: Tissue Says Blood Is Misleading, Confusing HIV Cure Efforts. Science (80- . ). 334, 1614–1614 (2011). 19. Ferreira, J., Hernandez-Vargas, E. a., Middleton, R.H.: Computer simulation of structured treatment interruption for HIV infection. Comput. Methods Programs Biomed. 104, 50–61 (2011). 20. Adams, B.M., Banks, H.T., Kwon, H.-D., Tran, H.T.: Dynamic multidrug therapies for hiv: optimal and sti control approaches. Math. Biosci. Eng. 1, 223–241 (2004). 21. Krakovska, O., Wahl, L.M.: Optimal drug treatment regimens for HIV depend on adherence. J. Theor. Biol. 246, 499–509 (2007). 22. Rivadeneira, P.S., Moog, C.H., Stan, G.-B., Brunet, C., Raffi, F., Ferré, V., Costanza, V., Mhawej, M.J., Biafore, F., Ouattara, D. a., Ernst, D., Fonteneau, R., Xia, X.: Mathematical Modeling of HIV Dynamics After Antiretroviral Therapy Initiation: A Review. Biores. Open Access. 3, 233–241 (2014). 23. von Neumann, J.: Theory of self-reproducing automata. Inf. Storage Retr. 5, 151 (1969). 24. Wodarz, D., Levy, D.N.: Effect of different modes of viral spread on the dynamics of multiply infected cells in human immunodeficiency virus infection. J. R. Soc. Interface. 8, 289–300 (2011). 25. Kislitsyn, A., Savinkov, R., Novkovic, M., Onder, L., Bocharov, G.: Computational Approach to 3D Modeling of the Lymph Node Geometry. Computation. 3, 222–234 (2015). 26. Fletcher, A.L., Acton, S.E., Knoblich, K.: Lymph node fibroblastic reticular cells in health and disease. Nat. Rev. Immunol. 15, 350–361 (2015). 27. Fauci, A.S., Pantaleo, G., Stanley, S., Weissman, D.: Immunopathogenic mechanisms of HIV infection. Ann. Intern. Med. 124, 654–663 (1996). 28. Lundgren, J.D., Babiker, A., Gordin, F., Emery, S., Grund, B., Sharma, S., Avihingsanon, A., Cooper, D.A., Fätkenheuer, G., Llibre, J.M., Molina, J.-M., Munderi, P., Schechter, M., Wood, R., Klingman, K.L., Collins, S., Lane, H.C., Phillips, A.N., Neaton, J.D.: Initiation of Antiretroviral Therapy in Early Asymptomatic HIV Infection. N. Engl. J. Med. 373, 795–807 (2015). 29. Maggiolo, F., Airoldi, M., Callegaro, A., Martinelli, C., Dolara, A., Bini, T., Gregis, G., Quinzan, G., Ripamonti, D., Ravasio, V., Suter, F.: CD4 cell-guided scheduled treatment interruptions in HIV-infected patients with sustained immunologic response to HAART. AIDS. 23, 799–807 (2009). 30. Perrin, D., Bezbradica, M., Crane, M., Ruskin, H.J., Duhamel, C.: High-Performance Computing for Data Analytics. 2012 IEEE/ACM 16th Int. Symp. Distrib. Simul. Real Time Appl. 234–242 (2012).