RESUMO
Malaria remains a global health problem despite the many attempts to control and eradicate it. There is an urgent need to understand the current transmission dynamics of malaria and to determine the interventions necessary to control malaria. In this paper, we seek to develop a fit-for-purpose mathematical model to assess the interventions needed to control malaria in an endemic setting. To achieve this, we formulate a malaria transmission model to analyse the spread of malaria in the presence of interventions. A sensitivity analysis of the model is performed to determine the relative impact of the model parameters on disease transmission. We explore how existing variations in the recruitment and management of intervention strategies affect malaria transmission. Results obtained from the study imply that the discontinuation of existing interventions has a significant effect on malaria prevalence. Thus, the maintenance of interventions is imperative for malaria elimination and eradication. In a scenario study aimed at assessing the impact of long-lasting insecticidal nets (LLINs), indoor residual spraying (IRS), and localized individual measures, our findings indicate that increased LLINs utilization and extended IRS coverage (with longer-lasting insecticides) cause a more pronounced reduction in symptomatic malaria prevalence compared to a reduced LLINs utilization and shorter IRS coverage. Additionally, our study demonstrates the impact of localized preventive measures in mitigating the spread of malaria when compared to the absence of interventions.
Assuntos
Mosquiteiros Tratados com Inseticida , Inseticidas , Malária , Conceitos Matemáticos , Modelos Biológicos , Controle de Mosquitos , Humanos , Malária/prevenção & controle , Malária/epidemiologia , Malária/transmissão , Controle de Mosquitos/métodos , Controle de Mosquitos/estatística & dados numéricos , Mosquiteiros Tratados com Inseticida/estatística & dados numéricos , Animais , Mosquitos Vetores/parasitologia , Prevalência , Simulação por Computador , Anopheles/parasitologia , Doenças Endêmicas/prevenção & controle , Doenças Endêmicas/estatística & dados numéricosRESUMO
Many baleen whales are renowned for their acoustic communication. Under pristine conditions, this communication can plausibly occur across hundreds of kilometres. Frequent vocalisations may allow a dispersed migrating group to maintain contact, and therefore benefit from improved navigation via the "wisdom of the crowd". Human activities have considerably inflated ocean noise levels. Here we develop a data-driven mathematical model to investigate how ambient noise levels may inhibit whale migration. Mathematical models allow us to simultaneously simulate collective whale migration behaviour, auditory cue detection, and noise propagation. Rising ambient noise levels are hypothesised to influence navigation through three mechanisms: (i) diminished communication space; (ii) reduced ability to hear external sound cues and; (iii) triggering noise avoidance behaviour. Comparing pristine and current soundscapes, we observe navigation impairment that ranges from mild (increased journey time) to extreme (failed navigation). Notably, the three mechanisms induce qualitatively different impacts on migration behaviour. We demonstrate the model's potential predictive power, exploring the extent to which migration may be altered under future shipping and construction scenarios.
RESUMO
Background: Mechanosensation is an important trigger of physiological processes in the gastrointestinal tract. Aberrant responses to mechanical input are associated with digestive disorders, including visceral hypersensitivity. Transient Receptor Potential Vanilloid 4 (TRPV4) is a mechanosensory ion channel with proposed roles in visceral afferent signaling, intestinal inflammation, and gut motility. While TRPV4 is a potential therapeutic target for digestive disease, current mechanistic understanding of how TRPV4 may influence gut function is limited by inconsistent reports of TRPV4 expression and distribution. Methods: In this study we profiled functional expression of TRPV4 using Ca2+ imaging of wholemount preparations of the mouse, monkey, and human intestine in combination with immunofluorescent labeling for established cellular markers. The involvement of TRPV4 in colonic motility was assessed in vitro using videomapping and contraction assays. Results: The TRPV4 agonist GSK1016790A evoked Ca2+ signaling in muscularis macrophages, enteric glia, and endothelial cells. TRPV4 specificity was confirmed using TRPV4 KO mouse tissue or antagonist pre-treatment. Calcium responses were not detected in other cell types required for neuromuscular signaling including enteric neurons, interstitial cells of Cajal, PDGFRα+ cells, and intestinal smooth muscle. TRPV4 activation led to rapid Ca2+ responses by a subpopulation of glial cells, followed by sustained Ca2+ signaling throughout the enteric glial network. Propagation of these waves was suppressed by inhibition of gap junctions or Ca2+ release from intracellular stores. Coordinated glial signaling in response to GSK1016790A was also disrupted in acute TNBS colitis. The involvement of TRPV4 in the initiation and propagation of colonic motility patterns was examined in vitro. Conclusions: We reveal a previously unappreciated role for TRPV4 in the initiation of distension-evoked colonic motility. These observations provide new insights into the functional role of TRPV4 activation in the gut, with important implications for how TRPV4 may influence critical processes including inflammatory signaling and motility.
RESUMO
Coordination of cell behaviour is key to a myriad of biological processes including tissue morphogenesis, wound healing, and tumour growth. As such, individual-based computational models, which explicitly describe inter-cellular interactions, are commonly used to model collective cell dynamics. However, when using individual-based models, it is unclear how descriptions of cell boundaries affect overall population dynamics. In order to investigate this we define three cell boundary descriptions of varying complexities for each of three widely used off-lattice individual-based models: overlapping spheres, Voronoi tessellation, and vertex models. We apply our models to multiple biological scenarios to investigate how cell boundary description can influence tissue-scale behaviour. We find that the Voronoi tessellation model is most sensitive to changes in the cell boundary description with basic models being inappropriate in many cases. The timescale of tissue evolution when using an overlapping spheres model is coupled to the boundary description. The vertex model is demonstrated to be the most stable to changes in boundary description, though still exhibits timescale sensitivity. When using individual-based computational models one should carefully consider how cell boundaries are defined. To inform future work, we provide an exploration of common individual-based models and cell boundary descriptions in frequently studied biological scenarios and discuss their benefits and disadvantages.
Assuntos
Conceitos Matemáticos , Modelos Biológicos , Software , Comunicação Celular , MorfogêneseRESUMO
Understanding the intestinal transport of particles is critical in several fields ranging from optimizing drug delivery systems to capturing health risks from the increased presence of nano- and micro-sized particles in human environment. While Caco-2 cell monolayers grown on permeable supports are the traditional in vitro model used to probe intestinal absorption of dissolved molecules, they fail to recapitulate the transcytotic activity of polarized enterocytes. Here, an intestine-on-chip model is combined with in silico modeling to demonstrate that the rate of particle transcytosis is ≈350× higher across Caco-2 cell monolayers exposed to fluid shear stress compared to Caco-2 cells in standard "static" configuration. This relates to profound phenotypical alterations and highly polarized state of cells grown under mechanical stimulation and it is shown that transcytosis in the microphysiological model is energy-dependent and involves both clathrin and macropinocytosis mediated endocytic pathways. Finally, it is demonstrated that the increased rate of transcytosis through cells exposed to flow is explained by a higher rate of internal particle transport (i.e., vesicular cellular trafficking and basolateral exocytosis), rather than a change in apical uptake (i.e., binding and endocytosis). Taken together, the findings have important implications for addressing research questions concerning intestinal transport of engineered and environmental particles.
Assuntos
Endocitose , Transcitose , Humanos , Células CACO-2 , Endocitose/fisiologia , Intestinos , Transporte BiológicoRESUMO
Calcium (Ca2ï¼) plays a critical role in the excitation contraction coupling (ECC) process that mediates the contraction of cardiomyocytes during each heartbeat. While ryanodine receptors (RyRs) are the primary Ca2ï¼ channels responsible for generating the cell-wide Ca2ï¼ transients during ECC, Ca2ï¼ release, via inositol 1,4,5-trisphosphate (IP3) receptors (IP3Rs) are also reported in cardiomyocytes to elicit ECC-modulating effects. Recent studies suggest that the localization of IP3Rs at dyads grant their ability to modify the occurrence of Ca2ï¼ sparks (elementary Ca2ï¼ release events that constitute cell wide Ca2ï¼ releases associated with ECC) which may underlie their modulatory influence on ECC. Here, we aim to uncover the mechanism by which dyad-localized IP3Rs influence Ca2ï¼ spark dynamics. To this end, we developed a mathematical model of the dyad that incorporates the behaviour of IP3Rs, in addition to RyRs, to reveal the impact of their activity on local Ca2ï¼ handling and consequent Ca2ï¼ spark occurrence and its properties. Consistent with published experimental data, our model predicts that the propensity for Ca2ï¼ spark formation increases in the presence of IP3R activity. Our simulations support the hypothesis that IP3Rs elevate Ca2ï¼ in the dyad, sensitizing proximal RyRs towards activation and hence Ca2ï¼ spark formation. The stochasticity of IP3R gating is an important aspect of this mechanism. However, dyadic IP3R activity lowers the Ca2ï¼ available in the junctional sarcoplasmic reticulum (JSR) for release, thus resulting in Ca2ï¼ sparks with similar durations but lower amplitudes.
Assuntos
Sinalização do Cálcio , Miócitos Cardíacos , Sinalização do Cálcio/fisiologia , Miócitos Cardíacos/metabolismo , Canal de Liberação de Cálcio do Receptor de Rianodina/metabolismo , Retículo Sarcoplasmático/metabolismo , Modelos Teóricos , Cálcio/metabolismoRESUMO
Malaria remains a global health problem despite the many attempts to control and eradicate it. There is an urgent need to understand the current transmission dynamics of malaria and to determine the interventions necessary to control malaria. In this paper, we seek to develop a fit-for-purpose mathematical model to assess the interventions needed to control malaria in an endemic setting. To achieve this, we formulate a malaria transmission model to analyse the spread of malaria in the presence of interventions. A sensitivity analysis of the model is performed to determine the relative impact of the model parameters on disease transmission. We explore how existing variations in the recruitment and management of intervention strategies affect malaria transmission. Results obtained from the study imply that the discontinuation of existing interventions has a significant effect on malaria prevalence. Thus, the maintenance of interventions is imperative for malaria elimination and eradication. In a scenario study aimed at assessing the impact of long-lasting insecticidal nets (LLINs), indoor residual spraying (IRS), and localized individual measures, our findings indicate that increased LLINs utilization and extended IRS coverage (with longer-lasting insecticides) cause a more pronounced reduction in symptomatic malaria prevalence compared to a reduced LLINs utilization and shorter IRS coverage. Additionally, our study demonstrates the impact of localized preventive measures in mitigating the spread of malaria when compared to the absence of interventions.
RESUMO
Nanoparticles are increasingly employed as a vehicle for the targeted delivery of therapeutics to specific cell types. However, much remains to be discovered about the fundamental biology that dictates the interactions between nanoparticles and cells. Accordingly, few nanoparticle-based targeted therapeutics have succeeded in clinical trials. One element that hinders our understanding of nanoparticle-cell interactions is the presence of heterogeneity in nanoparticle dosage data obtained from standard experiments. It is difficult to distinguish between heterogeneity that arises from stochasticity in nanoparticle-cell interactions, and that which arises from heterogeneity in the cell population. Mathematical investigations have revealed that both sources of heterogeneity contribute meaningfully to the heterogeneity in nanoparticle dosage. However, these investigations have relied on simplified models of nanoparticle internalisation. Here we present a stochastic mathematical model of nanoparticle internalisation that incorporates a suite of relevant biological phenomena such as multistage internalisation, cell division, asymmetric nanoparticle inheritance and nanoparticle saturation. Critically, our model provides information about nanoparticle dosage at an individual cell level. We perform model simulations to examine the influence of specific biological phenomena on the heterogeneity in nanoparticle dosage in the absence of heterogeneity in the cell population. Under certain modelling assumptions, we derive analytic approximations of the nanoparticle dosage distribution. We demonstrate that the analytic approximations are accurate, and show that nanoparticle dosage can be described by a Poisson mixture distribution with rate parameters that are a function of Beta-distributed random variables. We discuss the implications of the analytic results with respect to parameter estimation and model identifiability from standard experimental data. Finally, we highlight extensions and directions for future research.
Assuntos
Nanopartículas , Modelos Teóricos , Distribuição de Poisson , Divisão CelularRESUMO
Designing nano-engineered particles capable of the delivery of therapeutic and diagnostic agents to a specific target remains a significant challenge. Understanding how interactions between particles and cells are impacted by the physicochemical properties of the particle will help inform rational design choices. Mathematical and computational techniques allow for details regarding particle-cell interactions to be isolated from the interwoven set of biological, chemical, and physical phenomena involved in the particle delivery process. Here we present a machine learning framework capable of elucidating particle-cell interactions from experimental data. This framework employs a data-driven modelling approach, augmented by established biological knowledge. Crucially, the model of particle-cell interactions learned by the framework can be interpreted and analysed, in contrast to the 'black box' models inherent to other machine learning approaches. We apply the framework to association data for thirty different particle-cell pairs. This library of data contains both adherent and suspension cell lines, as well as a diverse collection of particles. We consider hyperbranched polymer and poly(methacrylic acid) particles, from 6 nm to 1032 nm in diameter, with small molecule, monoclonal antibody, and peptide surface functionalisations. Despite the diverse nature of the experiments, the learned models of particle-cell interactions for each particle-cell pair are remarkably consistent: out of 2048 potential models, only four unique models are learned. The models reveal that nonlinear saturation effects are a key feature governing particle-cell interactions. Further, the framework provides robust estimates of particle performance, which facilitates quantitative evaluation of particle design choices.
Assuntos
Aprendizado de Máquina , Polímeros , PeptídeosRESUMO
The maintenance of tissue and organ structures during dynamic homeostasis is often not well understood. In order for a system to be stable, cell renewal, cell migration and cell death must be finely balanced. Moreover, a tissue's shape must remain relatively unchanged. Simple epithelial tissues occur in various structures throughout the body, such as the endothelium, mesothelium, linings of the lungs, saliva and thyroid glands, and gastrointestinal tract. Despite the prevalence of simple epithelial tissues, there are few models which accurately describe how these tissues maintain a stable structure. Here, we present a novel, 3D, deformable, multilayer, cell-centre model of a simple epithelium. Cell movement is governed by the minimisation of a bending potential across the epithelium, cell-cell adhesion, and viscous effects. We show that the model is capable of maintaining a consistent tissue structure while undergoing self renewal. We also demonstrate the model's robustness under tissue renewal, cell migration and cell removal. The model presented here is a valuable advancement towards the modelling of tissues and organs with complex and generalised structures.
Assuntos
Epitélio , Adesão Celular , Morte Celular , Movimento Celular , HomeostaseRESUMO
Nanoparticles hold great preclinical promise in cancer therapy but continue to suffer attrition through clinical trials. Advanced, three dimensional (3D) cellular models such as tumor spheroids can recapitulate elements of the tumor environment and are considered the superior model to evaluate nanoparticle designs. However, there is an important need to better understand nanoparticle penetration kinetics and determine how different cell characteristics may influence this nanoparticle uptake. A key challenge with current approaches for measuring nanoparticle accumulation in spheroids is that they are often static, losing spatial and temporal information which may be necessary for effective nanoparticle evaluation in 3D cell models. To overcome this challenge, we developed an analysis platform, termed the Determination of Nanoparticle Uptake in Tumor Spheroids (DONUTS), which retains spatial and temporal information during quantification, enabling evaluation of nanoparticle uptake in 3D tumor spheroids. Outperforming linear profiling methods, DONUTS was able to measure silica nanoparticle uptake to 10 µm accuracy in both isotropic and irregularly shaped cancer cell spheroids. This was then extended to determine penetration kinetics, first by a forward-in-time, center-in-space model, and then by mathematical modelling, which enabled the direct evaluation of nanoparticle penetration kinetics in different spheroid models. Nanoparticle uptake was shown to inversely relate to particle size and varied depending on the cell type, cell stiffness and density of the spheroid model. The automated analysis method we have developed can be applied to live spheroids in situ, for the advanced evaluation of nanoparticles as delivery agents in cancer therapy.
Assuntos
Nanopartículas , Neoplasias , Humanos , Tamanho da Partícula , Análise Espaço-Temporal , Esferoides CelularesRESUMO
The question of whether biological populations survive or are eventually driven to extinction has long been examined using mathematical models. In this work, we study population survival or extinction using a stochastic, discrete lattice-based random walk model where individuals undergo movement, birth and death events. The discrete model is defined on a two-dimensional hexagonal lattice with periodic boundary conditions. A key feature of the discrete model is that crowding effects are introduced by specifying two different crowding functions that govern how local agent density influences movement events and birth/death events. The continuum limit description of the discrete model is a nonlinear reaction-diffusion equation, and we focus on crowding functions that lead to linear diffusion and a bistable source term that is often associated with the strong Allee effect. Using both the discrete and continuum modelling tools, we explore the complicated relationship between the long-term survival or extinction of the population and the initial spatial arrangement of the population. In particular, we study different spatial arrangements of initial distributions: (i) a well-mixed initial distribution where the initial density is independent of position in the domain; (ii) a vertical strip initial distribution where the initial density is independent of vertical position in the domain; and, (iii) several forms of two-dimensional initial distributions where the initial population is distributed in regions with different shapes. Our results indicate that the shape of the initial spatial distribution of the population affects extinction of bistable populations. All software required to solve the discrete and continuum models used in this work are available on GitHub .
Assuntos
Conceitos Matemáticos , Modelos Biológicos , Difusão , Extinção Biológica , Humanos , Densidade Demográfica , Dinâmica PopulacionalRESUMO
Flours were made from the sprouted seeds of the low- and high-tannin faba bean cultivars Fabelle, FB9-4, Snowbird, and Snowdrop. Headspace measurements on sprouted flours found the most favourable aroma profiles following 48 h sprouting and 24 h drying at 60 °C. Lipoxygenase activity, and the tannin, protein, and moisture contents were determined for unsprouted and sprouted faba bean flours. Lipoxygenase activity was higher in sprouted seeds before drying. Protein content increased after sprouting, whereas the tannin content decreased, especially for high-tannin varieties. Key volatile flavour compounds of faba bean flours included pentanal, hexanal, heptanal, octanal, nonanal, decanal, 1-hexanol, 1-octen-3-ol, 3-methylbutanal, phenyl acetaldehyde, 3-methylbutyric acid, d-limonene, ß-linalool, menthol, and estragole; these include oxidative degradation products of oleic, linoleic, and some amino acids. An overall flavour improvement was achieved after germination, as indicated by a decrease in bitter compounds (tannins) and beany flavours (hexanal, nonanal, 2-heptanone, and 2-pentylfuran).
Assuntos
Fabaceae , Vicia faba , Farinha , Taninos , PaladarRESUMO
Understanding the interactions between nano-engineered particles and cells is necessary for the rational design of particles for therapeutic, diagnostic and imaging purposes. In particular, the informed design of particles relies on the quantification of the relationship between the physicochemical properties of the particles and the rate at which cells interact with, and subsequently internalise, particles. Quantitative models, both mathematical and computational, provide a powerful tool for elucidating this relationship, as well as for understanding the mechanisms governing the intertwined processes of interaction and internalisation. Here we review the different types of mathematical and computational models that have been used to examine particle-cell interactions and particle internalisation. We detail the mathematical methodology for each type of model, the benefits and limitations associated with the different types of models, and highlight the advances in understanding gleaned from the application of these models to experimental observations of particle internalisation. We discuss the recent proposal and ongoing community adoption of standardised experimental reporting, and how this adoption is an important step toward unlocking the full potential of modelling approaches. Finally, we consider future directions in quantitative models of particle-cell interactions and highlight the need for hybrid experimental and theoretical investigations to address hitherto unanswered questions.
RESUMO
We present a solid theoretical foundation for interpreting the origin of Allee effects by providing the missing link in understanding how local individual-based mechanisms translate to global population dynamics. Allee effects were originally proposed to describe population dynamics that cannot be explained by exponential and logistic growth models. However, standard methods often calibrate Allee effect models to match observed global population dynamics without providing any mechanistic insight. By introducing a stochastic individual-based model, with proliferation, death and motility rates that depend on local density, we present a modelling framework that translates particular global Allee effects to specific individual-based mechanisms. Using data from ecology and cell biology, we unpack individual-level mechanisms implicit in an Allee effect model and provide simulation tools for others to repeat this analysis.
RESUMO
The question of whether a population will persist or go extinct is of key interest throughout ecology and biology. Various mathematical techniques allow us to generate knowledge regarding individual behaviour, which can be analysed to obtain predictions about the ultimate survival or extinction of the population. A common model employed to describe population dynamics is the lattice-based random walk model with crowding (exclusion). This model can incorporate behaviour such as birth, death and movement, while including natural phenomena such as finite size effects. Performing sufficiently many realizations of the random walk model to extract representative population behaviour is computationally intensive. Therefore, continuum approximations of random walk models are routinely employed. However, standard continuum approximations are notoriously incapable of making accurate predictions about population extinction. Here, we develop a new continuum approximation, the state-space diffusion approximation, which explicitly accounts for population extinction. Predictions from our approximation faithfully capture the behaviour in the random walk model, and provides additional information compared to standard approximations. We examine the influence of the number of lattice sites and initial number of individuals on the long-term population behaviour, and demonstrate the reduction in computation time between the random walk model and our approximation.
RESUMO
Nano-engineered particles have the potential to enhance therapeutic success and reduce toxicity-based treatment side effects via the targeted delivery of drugs to cells. This delivery relies on complex interactions between numerous biological, chemical and physical processes. The intertwined nature of these processes has thus far hindered attempts to understand their individual impact. Variation in experimental data, such as the number of particles inside each cell, further inhibits understanding. Here, we present a mathematical framework that is capable of examining the impact of individual processes during particle delivery. We demonstrate that variation in experimental particle uptake data can be explained by three factors: random particle motion; variation in particle-cell interactions; and variation in the maximum particle uptake per cell. Without all three factors, the experimental data cannot be explained. This work provides insight into biological mechanisms that cause heterogeneous responses to treatment, and enables precise identification of treatment-resistant cell subpopulations.
Assuntos
Comunicação Celular , Nanopartículas , Tamanho da PartículaRESUMO
Nanoengineering has the potential to revolutionize medicine by designing drug delivery systems that are both efficacious and highly selective. Determination of the affinity between cell lines and nanoparticles is thus of central importance, both to enable comparison of particles and to facilitate prediction of in vivo response. Attempts to compare particle performance can be dominated by experimental artifacts (including settling effects) or variability in experimental protocol. Instead, qualitative methods are generally used, limiting the reusability of many studies. Herein, we introduce a mathematical model-based approach to quantify the affinity between a cell-particle pairing, independent of the aforementioned confounding artifacts. The analysis presented can serve as a quantitative metric of the stealth, fouling, and targeting performance of nanoengineered particles in vitro. We validate this approach using a newly created in vitro dataset, consisting of seven different disulfide-stabilized poly(methacrylic acid) particles ranging from ~100 to 1000â¯nm in diameter that were incubated with three different cell lines (HeLa, THP-1, and RAW 264.7). We further expanded this dataset through the inclusion of previously published data and use it to determine which of five mathematical models best describe cell-particle association. We subsequently use this model to perform a quantitative comparison of cell-particle association for cell-particle pairings in our dataset. This analysis reveals a more complex cell-particle association relationship than a simplistic interpretation of the data, which erroneously assigns high affinity for all cell lines examined to large particles. Finally, we provide an online tool (http://bionano.xyz/estimator), which allows other researchers to easily apply this modeling approach to their experimental results.
Assuntos
Modelos Teóricos , Nanopartículas/administração & dosagem , Animais , Dissulfetos/administração & dosagem , Dissulfetos/química , Ouro/administração & dosagem , Ouro/química , Células HeLa , Humanos , Camundongos , Nanopartículas/química , Tamanho da Partícula , Ácidos Polimetacrílicos/administração & dosagem , Ácidos Polimetacrílicos/química , Células RAW 264.7 , Dióxido de Silício/administração & dosagem , Dióxido de Silício/química , Células THP-1RESUMO
Selective self-assembly in multicomponent mixtures offers a method for isolating desired components from complex systems for the rapid production of functional materials. Developing approaches capable of selective assembly of "target" components into intended three-dimensional structures is challenging because of the intrinsically high complexity of multicomponent systems. Herein, we report the selective coordination-driven self-assembly of metal-phenolic networks (MPNs) from a series of complex multicomponent systems (including crude plant extracts) into thin films via metal chelation with phenolic ligands. The metal (FeIII) selectively assembles low abundant phenolic components (e.g., myricetrin and quercetrin) from plant extracts into thin films. This selective metal-phenolic assembly is independent of the substrate properties (e.g., size, surface charge, and shape). Moreover, the high selectivity is consistent across different target phenolic ligands in model mixtures, even though each individual component can form thin films from single-component systems. A computational simulation of film formation suggests that the driving force for the selective behavior stems from differences in the number of chelating sites in the phenolic structures. The MPN films are shown to demonstrate improved antioxidant properties compared with the corresponding phenolic compounds in their free form, therefore exhibiting potential as free-standing antioxidant films.
RESUMO
Environments with immobile obstacles or void regions that inhibit and alter the motion of individuals within that environment are ubiquitous. Correlation in the location of individuals within such environments arises as a combination of the mechanisms governing individual behavior and the heterogeneous structure of the environment. Measures of spatial structure and correlation have been successfully implemented to elucidate the roles of the mechanisms underpinning the behavior of individuals. In particular, the pair correlation function has been used across biology, ecology, and physics to obtain quantitative insight into a variety of processes. However, naively applying standard pair correlation functions in the presence of obstacles may fail to detect correlation, or suggest false correlations, due to a reliance on a distance metric that does not account for obstacles. To overcome this problem, here we present an analytic expression for calculating a corrected pair correlation function for lattice-based domains containing obstacles. We demonstrate that this obstacle pair correlation function is necessary for isolating the correlation associated with the behavior of individuals, rather than the structure of the environment. Using simulations that mimic cell migration and proliferation we demonstrate that the obstacle pair correlation function recovers the short-range correlation known to be present in this process, independent of the heterogeneous structure of the environment. Further, we show that the analytic calculation of the obstacle pair correlation function derived here is significantly faster to implement than the corresponding numerical approach.