Genome-wide identification of replication fork stalling/pausing sites and the interplay between RNA Pol II transcription and DNA replication progression.

BACKGROUND: DNA replication progression can be affected by the presence of physical barriers like the RNA polymerases, leading to replication stress and DNA damage. Nonetheless, we do not know how transcription influences overall DNA replication progression. RESULTS: To characterize sites where DNA replication forks stall and pause, we establish a genome-wide approach to identify them. This approach uses multiple timepoints during S-phase to identify replication fork/stalling hotspots as replication progresses through the genome. These sites are typically associated with increased DNA damage, overlapped with fragile sites and with breakpoints of rearrangements identified in cancers but do not overlap with replication origins. Overlaying these sites with a genome-wide analysis of RNA polymerase II transcription, we find that replication fork stalling/pausing sites inside genes are directly related to transcription progression and activity. Indeed, we find that slowing down transcription elongation slows down directly replication progression through genes. This indicates that transcription and replication can coexist over the same regions. Importantly, rearrangements found in cancers overlapping transcription-replication collision sites are detected in non-transformed cells and increase following treatment with ATM and ATR inhibitors. At the same time, we find instances where transcription activity favors replication progression because it reduces histone density. CONCLUSIONS: Altogether, our findings highlight how transcription and replication overlap during S-phase, with both positive and negative consequences for replication fork progression and genome stability by the coexistence of these two processes.

Discretised Flux Balance Analysis for Reaction-Diffusion Simulation of Single-Cell Metabolism.

Metabolites have to diffuse within the sub-cellular compartments they occupy to specific locations where enzymes are, so reactions could occur. Conventional flux balance analysis (FBA), a method based on linear programming that is commonly used to model metabolism, implicitly assumes that all enzymatic reactions are not diffusion-limited though that may not always be the case. In this work, we have developed a spatial method that implements FBA on a grid-based system, to enable the exploration of diffusion effects on metabolism. Specifically, the method discretises a living cell into a two-dimensional grid, represents the metabolic reactions in each grid element as well as the diffusion of metabolites to and from neighbouring elements, and simulates the system as a single linear programming problem. We varied the number of rows and columns in the grid to simulate different cell shapes, and the method was able to capture diffusion effects at different shapes. We then used the method to simulate heterogeneous enzyme distribution, which suggested a theoretical effect on variability at the population level. We propose the use of this method, and its future extensions, to explore how spatiotemporal organisation of sub-cellular compartments and the molecules within could affect cell behaviour.

A mechanical modeling framework to study endothelial permeability.

The inner lining of blood vessels, the endothelium, is made up of endothelial cells. Vascular endothelial (VE)-cadherin protein forms a bond with VE-cadherin from neighboring cells to determine the size of gaps between the cells and thereby regulate the size of particles that can cross the endothelium. Chemical cues such as thrombin, along with mechanical properties of the cell and extracellular matrix are known to affect the permeability of endothelial cells. Abnormal permeability is found in patients suffering from diseases including cardiovascular diseases, cancer, and COVID-19. Even though some of the regulatory mechanisms affecting endothelial permeability are well studied, details of how several mechanical and chemical stimuli acting simultaneously affect endothelial permeability are not yet understood. In this article, we present a continuum-level mechanical modeling framework to study the highly dynamic nature of the VE-cadherin bonds. Taking inspiration from the catch-slip behavior that VE-cadherin complexes are known to exhibit, we model the VE-cadherin homophilic bond as cohesive contact with damage following a traction-separation law. We explicitly model the actin cytoskeleton and substrate to study their role in permeability. Our studies show that mechanochemical coupling is necessary to simulate the influence of the mechanical properties of the substrate on permeability. Simulations show that shear between cells is responsible for the variation in permeability between bicellular and tricellular junctions, explaining the phenotypic differences observed in experiments. An increase in the magnitude of traction force due to disturbed flow that endothelial cells experience results in increased permeability, and it is found that the effect is higher on stiffer extracellular matrix. Finally, we show that the cylindrical monolayer exhibits higher permeability than the planar monolayer under unconstrained cases. Thus, we present a contact mechanics-based mechanochemical model to investigate the variation in the permeability of endothelial monolayer due to multiple loads acting simultaneously.

Mathematical reconstruction of the metabolic network in an in-vitro multiple myeloma model.

It is increasingly apparent that cancer cells, in addition to remodelling their metabolism to survive and proliferate, adapt and manipulate the metabolism of other cells. This property may be a telling sign that pre-clinical tumour metabolism studies exclusively utilising in-vitro mono-culture models could prove to be limited for uncovering novel metabolic targets able to translate into clinical therapies. Although this is increasingly recognised, and work towards addressing the issue is becoming routinary much remains poorly understood. For instance, knowledge regarding the biochemical mechanisms through which cancer cells manipulate non-cancerous cell metabolism, and the subsequent impact on their survival and proliferation remains limited. Additionally, the variations in these processes across different cancer types and progression stages, and their implications for therapy, also remain largely unexplored. This study employs an interdisciplinary approach that leverages the predictive power of mathematical modelling to enrich experimental findings. We develop a functional multicellular in-silico model that facilitates the qualitative and quantitative analysis of the metabolic network spawned by an in-vitro co-culture model of bone marrow mesenchymal stem- and myeloma cell lines. To procure this model, we devised a bespoke human genome constraint-based reconstruction workflow that combines aspects from the legacy mCADRE & Metabotools algorithms, the novel redHuman algorithm, along with 13C-metabolic flux analysis. Our workflow transforms the latest human metabolic network matrix (Recon3D) into two cell-specific models coupled with a metabolic network spanning a shared growth medium. When cross-validating our in-silico model against the in-vitro model, we found that the in-silico model successfully reproduces vital metabolic behaviours of its in-vitro counterpart; results include cell growth predictions, respiration rates, as well as support for observations which suggest cross-shuttling of redox-active metabolites between cells.

Interplay of p53 and XIAP protein dynamics orchestrates cell fate in response to chemotherapy.

Chemotherapeutic drugs are used to treat almost all types of cancer, but the intended response, i.e., elimination, is often incomplete, with a subset of cancer cells resisting treatment. Two critical factors play a role in chemoresistance: the p53 tumour suppressor gene and the X-linked inhibitor of apoptosis (XIAP). These proteins have been shown to act synergistically to elicit cellular responses upon DNA damage induced by chemotherapy, yet, the mechanism is poorly understood. This study introduces a mathematical model characterising the apoptosis pathway activation by p53 before and after mitochondrial outer membrane permeabilisation upon treatment with the chemotherapy Doxorubicin (Dox). "In-silico" simulations show that the p53 dynamics change dose-dependently. Under medium to high doses of Dox, p53 concentration ultimately stabilises to a high level regardless of XIAP concentrations. However, caspase-3 activation may be triggered or not depending on the XIAP induction rate, ultimately determining whether the cell will perish or resist. Consequently, the model predicts that failure to activate apoptosis in some cancer cells expressing wild-type p53 might be due to heterogeneity between cells in upregulating the XIAP protein, rather than due to the p53 protein concentration. Our model suggests that the interplay of the p53 dynamics and the XIAP induction rate is critical to determine the cancer cells' therapeutic response.

Endothelium and Subendothelial Matrix Mechanics Modulate Cancer Cell Transendothelial Migration.

Cancer cell extravasation, a key step in the metastatic cascade, involves cancer cell arrest on the endothelium, transendothelial migration (TEM), followed by the invasion into the subendothelial extracellular matrix (ECM) of distant tissues. While cancer research has mostly focused on the biomechanical interactions between tumor cells (TCs) and ECM, particularly at the primary tumor site, very little is known about the mechanical properties of endothelial cells and the subendothelial ECM and how they contribute to the extravasation process. Here, an integrated experimental and theoretical framework is developed to investigate the mechanical crosstalk between TCs, endothelium and subendothelial ECM during in vitro cancer cell extravasation. It is found that cancer cell actin-rich protrusions generate complex push-pull forces to initiate and drive TEM, while transmigration success also relies on the forces generated by the endothelium. Consequently, mechanical properties of the subendothelial ECM and endothelial actomyosin contractility that mediate the endothelial forces also impact the endothelium's resistance to cancer cell transmigration. These results indicate that mechanical features of distant tissues, including force interactions between the endothelium and the subendothelial ECM, are key determinants of metastatic organotropism.

Feedback between mechanosensitive signaling and active forces governs endothelial junction integrity.

The formation and recovery of gaps in the vascular endothelium governs a wide range of physiological and pathological phenomena, from angiogenesis to tumor cell extravasation. However, the interplay between the mechanical and signaling processes that drive dynamic behavior in vascular endothelial cells is not well understood. In this study, we propose a chemo-mechanical model to investigate the regulation of endothelial junctions as dependent on the feedback between actomyosin contractility, VE-cadherin bond turnover, and actin polymerization, which mediate the forces exerted on the cell-cell interface. Simulations reveal that active cell tension can stabilize cadherin bonds, but excessive RhoA signaling can drive bond dissociation and junction failure. While actin polymerization aids gap closure, high levels of Rac1 can induce junction weakening. Combining the modeling framework with experiments, our model predicts the influence of pharmacological treatments on the junction state and identifies that a critical balance between RhoA and Rac1 expression is required to maintain junction stability. Our proposed framework can help guide the development of therapeutics that target the Rho family of GTPases and downstream active mechanical processes.

Effects of altered cellular ultrastructure on energy metabolism in diabetic cardiomyopathy: an in silico study.

Diabetic cardiomyopathy is a leading cause of heart failure in diabetes. At the cellular level, diabetic cardiomyopathy leads to altered mitochondrial energy metabolism and cardiomyocyte ultrastructure. We combined electron microscopy (EM) and computational modelling to understand the impact of diabetes-induced ultrastructural changes on cardiac bioenergetics. We collected transverse micrographs of multiple control and type I diabetic rat cardiomyocytes using EM. Micrographs were converted to finite-element meshes, and bioenergetics was simulated over them using a biophysical model. The simulations also incorporated depressed mitochondrial capacity for oxidative phosphorylation (OXPHOS) and creatine kinase (CK) reactions to simulate diabetes-induced mitochondrial dysfunction. Analysis of micrographs revealed a 14% decline in mitochondrial area fraction in diabetic cardiomyocytes, and an irregular arrangement of mitochondria and myofibrils. Simulations predicted that this irregular arrangement, coupled with the depressed activity of mitochondrial CK enzymes, leads to large spatial variation in adenosine diphosphate (ADP)/adenosine triphosphate (ATP) ratio profile of diabetic cardiomyocytes. However, when spatially averaged, myofibrillar ADP/ATP ratios of a cardiomyocyte do not change with diabetes. Instead, average concentration of inorganic phosphate rises by 40% owing to lower mitochondrial area fraction and dysfunction in OXPHOS. These simulations indicate that a disorganized cellular ultrastructure negatively impacts metabolite transport in diabetic cardiomyopathy. This article is part of the theme issue 'The cardiomyocyte: new revelations on the interplay between architecture and function in growth, health, and disease'.

Integrated in silico and 3D in vitro model of macrophage migration in response to physical and chemical factors in the tumor microenvironment.

Macrophages are abundant in the tumor microenvironment (TME), serving as accomplices to cancer cells for their invasion. Studies have explored the biochemical mechanisms that drive pro-tumor macrophage functions; however the role of TME interstitial flow (IF) is often disregarded. Therefore, we developed a three-dimensional microfluidic-based model with tumor cells and macrophages to study how IF affects macrophage migration and its potential contribution to cancer invasion. The presence of either tumor cells or IF individually increased macrophage migration directedness and speed. Interestingly, there was no additive effect on macrophage migration directedness and speed under the simultaneous presence of tumor cells and IF. Further, we present an in silico model that couples chemokine-mediated signaling with mechanosensing networks to explain our in vitro observations. In our model design, we propose IL-8, CCL2, and ß-integrin as key pathways that commonly regulate various Rho GTPases. In agreement, in vitro macrophage migration remained elevated when exposed to a saturating concentration of recombinant IL-8 or CCL2 or to the co-addition of a sub-saturating concentration of both cytokines. Moreover, antibody blockade against IL-8 and/or CCL2 inhibited migration that could be restored by IF, indicating cytokine-independent mechanisms of migration induction. Importantly, we demonstrate the utility of an integrated in silico and 3D in vitro approach to aid the design of tumor-associated macrophage-based immunotherapeutic strategies.

Are the Effects of Independent Biophysical Factors Linearly Additive? A 3D Tumor Migration Model.

Interstitial fluid flow plays a critical role in tumor cell invasion, yet this role has not been explored extensively in combination with other microenvironmental factors. Here, we establish a novel computational model of three-dimensional breast cancer cell migration to unveil the effect of interstitial fluid flow in the dependence of various extracellular matrix (ECM) physical properties. Our model integrates several principal factors: fluid dynamics, autologous chemotaxis, collagen fiber network structure, ECM stiffness, and cell-fiber and cell-flow interaction. First, independently with an aligned collagen fiber network and interstitial fluid flow, this model is validated by successfully reproducing the cell migration patterns. In the model, the interstitial fluid flow leads to directional symmetry breaking of chemotactic migration and synergizes with the ECM orientation to regulate cell migration. This synergy is universal in both the mesenchymal and the amoeboid migration modes, despite the fact that the cell-ECM interaction are different. Consequently, we construct a cell displacement function depending on these factors. Our cell migration model enables three-dimensional cancer migration prediction, mechanism exploration, and inhibition treatment design in a complex tumor microenvironment.

Balance of mechanical forces drives endothelial gap formation and may facilitate cancer and immune-cell extravasation.

The formation of gaps in the endothelium is a crucial process underlying both cancer and immune cell extravasation, contributing to the functioning of the immune system during infection, the unfavorable development of chronic inflammation and tumor metastasis. Here, we present a stochastic-mechanical multiscale model of an endothelial cell monolayer and show that the dynamic nature of the endothelium leads to spontaneous gap formation, even without intervention from the transmigrating cells. These gaps preferentially appear at the vertices between three endothelial cells, as opposed to the border between two cells. We quantify the frequency and lifetime of these gaps, and validate our predictions experimentally. Interestingly, we find experimentally that cancer cells also preferentially extravasate at vertices, even when they first arrest on borders. This suggests that extravasating cells, rather than initially signaling to the endothelium, might exploit the autonomously forming gaps in the endothelium to initiate transmigration.

Computational model of wound healing: EGF secreted by fibroblasts promotes delayed re-epithelialization of epithelial keratinocytes.

It is widely agreed that keratinocyte migration plays a crucial role in wound re-epithelialization. Defects in this function contribute to wound reoccurrence causing significant clinical problems. Several in vitro studies have shown that the speed of migrating keratinocytes can be regulated by epidermal growth factor (EGF) which affects keratinocyte's integrin expression. The relationship between integrin expression (through cell-matrix adhesion) stimulated by EGF and keratinocyte migration speed is not linear since increased adhesion, due to increased integrin expression, has been experimentally shown to slow down cell migration due to the biphasic dependence of cell speed on adhesion. In our previous work we showed that keratinocytes that were co-cultured with EGF-enhanced fibroblasts formed an asymmetric migration pattern, where, the cumulative distances of keratinocytes migrating toward fibroblasts were smaller than those migrating away from fibroblasts. This asymmetric pattern is thought to be provoked by high EGF concentration secreted by fibroblasts. The EGF stimulates the expression of integrin receptors on the surface of keratinocytes migrating toward fibroblasts via paracrine signaling. In this paper, we present a computational model of keratinocyte migration that is controlled by EGF secreted by fibroblasts using the Cellular Potts Model (CPM). Our computational simulation results confirm the asymmetric pattern observed in experiments. These results provide a deeper insight into our understanding of the complexity of keratinocyte migration in the presence of growth factor gradients and may explain re-epithelialization failure in impaired wound healing.

Mechanical and Systems Biology of Cancer.

Mechanics and biochemical signaling are both often deregulated in cancer, leading toincreased cell invasiveness, proliferation, and survival. The dynamics and interactions of cytoskeletal components control basic mechanical properties, such as cell tension, stiffness, and engagement with the extracellular environment, which can lead to extracellular matrix remodeling. Intracellular mechanics can alter signaling and transcription factors, impacting cell decision making. Additionally, signaling from soluble and mechanical factors in the extracellular environment, such as substrate stiffness and ligand density, can modulate cytoskeletal dynamics. Computational models closely integrated with experimental support, incorporating cancer-specific parameters, can provide quantitative assessments and serve as predictive tools toward dissecting the feedback between signaling and mechanics and across multiple scales and domains in tumor progression.

Solid dissolution in a fluid solvent is characterized by the interplay of surface area-dependent diffusion and physical fragmentation.

The processes of dissolution and fragmentation have high relevance in pharmaceutical research, medicine, digestive physiology, and engineering design. Experimentally, dissolution and fragmentation are observed to occur simultaneously, yet little is known about the relative importance of each of these processes and their impact on the dissolution process as a whole. Thus, in order to better explain these phenomena and the manner in which they interact, we have developed a novel mathematical model of dissolution, based on partial differential equations, taking into consideration the two constituent processes of surface area-dependent diffusive mass removal and physical fragmentation of the solid particles, and the basic physical laws governing these processes. With this model, we have been able to quantify the effects of the interplay between these two processes and determine the optimal conditions for rapid solid dissolution in liquid solvents. We were able to reproduce experimentally observed phenomena and simulate dissolution under a wide range of experimentally occurring conditions to give new perspectives into the kinetics of this common, yet complex process. Finally, we demonstrated the utility of this model to aid in experiment and device design as an optimisation tool.

Integrated Analysis of Intracellular Dynamics of MenaINV Cancer Cells in a 3D Matrix.

The intracellular environment is composed of a filamentous network that exhibits dynamic turnover of cytoskeletal components and internal force generation from molecular motors. Particle tracking microrheology enables a means to probe the internal mechanics and dynamics. Here, we develop an analytical model to capture the basic features of the active intracellular mechanical environment, including both thermal and motor-driven effects, and show consistency with a diverse range of experimental microrheology data. We further perform microrheology experiments, integrated with Brownian dynamics simulations of the active cytoskeleton, on metastatic breast cancer cells embedded in a three-dimensional collagen matrix with and without the presence of epidermal growth factor to probe the intracellular mechanical response in a physiologically mimicking scenario. Our results demonstrate that EGF stimulation can alter intracellular stiffness and power output from molecular motor-driven fluctuations in cells overexpressing an invasive isoform of the actin-associated protein Mena.

Dynamic interplay between tumour, stroma and immune system can drive or prevent tumour progression.

In the tumour microenvironment, cancer cells directly interact with both the immune system and the stroma. It is firmly established that the immune system, historically believed to be a major part of the body's defence against tumour progression, can be reprogrammed by tumour cells to be ineffective, inactivated, or even acquire tumour promoting phenotypes. Likewise, stromal cells and extracellular matrix can also have pro-and anti-tumour properties. However, there is strong evidence that the stroma and immune system also directly interact, therefore creating a tripartite interaction that exists between cancer cells, immune cells and tumour stroma. This interaction contributes to the maintenance of a chronically inflamed tumour microenvironment with pro-tumorigenic immune phenotypes and facilitated metastatic dissemination. A comprehensive understanding of cancer in the context of dynamical interactions of the immune system and the tumour stroma is therefore required to truly understand the progression toward and past malignancy.

Controlling uncertainty in aptamer selection.

The search for high-affinity aptamers for targets such as proteins, small molecules, or cancer cells remains a formidable endeavor. Systematic Evolution of Ligands by EXponential Enrichment (SELEX) offers an iterative process to discover these aptamers through evolutionary selection of high-affinity candidates from a highly diverse random pool. This randomness dictates an unknown population distribution of fitness parameters, encoded by the binding affinities, toward SELEX targets. Adding to this uncertainty, repeating SELEX under identical conditions may lead to variable outcomes. These uncertainties pose a challenge when tuning selection pressures to isolate high-affinity ligands. Here, we present a stochastic hybrid model that describes the evolutionary selection of aptamers to explore the impact of these unknowns. To our surprise, we find that even single copies of high-affinity ligands in a pool of billions can strongly influence population dynamics, yet their survival is highly dependent on chance. We perform Monte Carlo simulations to explore the impact of environmental parameters, such as the target concentration, on selection efficiency in SELEX and identify strategies to control these uncertainties to ultimately improve the outcome and speed of this time- and resource-intensive process.

Stochastic multi-scale models of competition within heterogeneous cellular populations: Simulation methods and mean-field analysis.

We propose a modelling framework to analyse the stochastic behaviour of heterogeneous, multi-scale cellular populations. We illustrate our methodology with a particular example in which we study a population with an oxygen-regulated proliferation rate. Our formulation is based on an age-dependent stochastic process. Cells within the population are characterised by their age (i.e. time elapsed since they were born). The age-dependent (oxygen-regulated) birth rate is given by a stochastic model of oxygen-dependent cell cycle progression. Once the birth rate is determined, we formulate an age-dependent birth-and-death process, which dictates the time evolution of the cell population. The population is under a feedback loop which controls its steady state size (carrying capacity): cells consume oxygen which in turn fuels cell proliferation. We show that our stochastic model of cell cycle progression allows for heterogeneity within the cell population induced by stochastic effects. Such heterogeneous behaviour is reflected in variations in the proliferation rate. Within this set-up, we have established three main results. First, we have shown that the age to the G1/S transition, which essentially determines the birth rate, exhibits a remarkably simple scaling behaviour. Besides the fact that this simple behaviour emerges from a rather complex model, this allows for a huge simplification of our numerical methodology. A further result is the observation that heterogeneous populations undergo an internal process of quasi-neutral competition. Finally, we investigated the effects of cell-cycle-phase dependent therapies (such as radiation therapy) on heterogeneous populations. In particular, we have studied the case in which the population contains a quiescent sub-population. Our mean-field analysis and numerical simulations confirm that, if the survival fraction of the therapy is too high, rescue of the quiescent population occurs. This gives rise to emergence of resistance to therapy since the rescued population is less sensitive to therapy.