Impact of implementation choices on quantitative predictions of cell-based computational models

‘Cell-based’ models provide a powerful computational tool for studying the mechanisms underlying the growth and dynamics of biological tissues in health and disease. An increasing amount of quantitative data with cellular resolution has paved the way for the quantitative parameterisation and validation of such models. However, the numerical implementation of cell-based models remains challenging, and little work has been done to understand to what extent implementation choices may influence model predictions. Here, we consider the numerical implementation of a popular class of cell-based models called vertex models, which are often used to study epithelial tissues. In two-dimensional vertex models, a tissue is approximated as a tessellation of polygons and the vertices of these polygons move due to mechanical forces originating from the cells. Such models have been used extensively to study the mechanical regulation of tissue topology in the literature. Here, we analyse how the model predictions may be affected by numerical parameters, such as the size of the time step, and non-physical model parameters, such as length thresholds for cell rearrangement. We find that vertex positions and summary statistics are sensitive to several of these implementation parameters. For example, the predicted tissue size decreases with decreasing cell cycle durations, and cell rearrangement may be suppressed by large time steps. These findings are counter-intuitive and illustrate that model predictions need to be thoroughly analysed and implementation details carefully considered when applying cell-based computational models in a quantitative setting.Publisher PDFPeer reviewe

Robust cell tracking in epithelial tissues through identification of maximum common subgraphs

Tracking of cells in live-imaging microscopy videos of epithelial sheets is a powerful tool for investigating fundamental processes in embryonic development. Characterizing cell growth, proliferation, intercalation and apoptosis in epithelia helps us to understand how morphogenetic processes such as tissue invagination and extension are locally regulated and controlled. Accurate cell tracking requires correctly resolving cells entering or leaving the field of view between frames, cell neighbour exchanges, cell removals and cell divisions. However, current tracking methods for epithelial sheets are not robust to large morphogenetic deformations and require significant manual interventions. Here, we present a novel algorithm for epithelial cell tracking, exploiting the graph-theoretic concept of a ‘maximum common subgraph’ to track cells between frames of a video. Our algorithm does not require the adjustment of tissue-specific parameters, and scales in sub-quadratic time with tissue size. It does not rely on precise positional information, permitting large cell movements between frames and enabling tracking in datasets acquired at low temporal resolution due to experimental constraints such as phototoxicity. To demonstrate the method, we perform tracking on the Drosophila embryonic epidermis and compare cell–cell rearrangements to previous studies in other tissues. Our implementation is open source and generally applicable to epithelial tissues

Multifractals of Normalized First Passage Time in Sierpinski Gasket

The multifractal behavior of the normalized first passage time is investigated on the two dimensional Sierpinski gasket with both absorbing and reflecting barriers. The normalized first passage time for Sinai model and the logistic model to arrive at the absorbing barrier after starting from an arbitrary site, especially obtained by the calculation via the Monte Carlo simulation, is discussed numerically. The generalized dimension and the spectrum are also estimated from the distribution of the normalized first passage time, and compared with the results on the finitely square lattice.Comment: 10 pages, Latex, with 3 figures and 1 table. to be published in J. Phys. Soc. Jpn. Vol.67(1998

Inhibition of sodium–glucose cotransporter-2 preserves cardiac function during regional myocardial ischemia independent of alterations in myocardial substrate utilization

The goal of the present study was to evaluate the effects of SGLT2i on cardiac contractile function, substrate utilization, and efficiency before and during regional myocardial ischemia/reperfusion injury in normal, metabolically healthy swine. Lean swine received placebo or canagliflozin (300 mg PO) 24 h prior to and the morning of an invasive physiologic study protocol. Hemodynamic and cardiac function measurements were obtained at baseline, during a 30-min complete occlusion of the circumflex coronary artery, and during a 2-h reperfusion period. Blood pressure, heart rate, coronary flow, and myocardial oxygen consumption were unaffected by canagliflozin treatment. Ventricular volumes remained unchanged in controls throughout the protocol. At the onset of ischemia, canagliflozin produced acute large increases in left ventricular end-diastolic and systolic volumes which returned to baseline with reperfusion. Canagliflozin-mediated increases in end-diastolic volume were directly associated with increases in stroke volume and stroke work relative to controls during ischemia. Canagliflozin also increased cardiac work efficiency during ischemia relative to control swine. No differences in myocardial uptake of glucose, lactate, free fatty acids or ketones, were noted between treatment groups at any time. In separate experiments using a longer 60 min coronary occlusion followed by 2 h of reperfusion, canagliflozin increased end-diastolic volume and stroke volume and significantly diminished myocardial infarct size relative to control swine. These data demonstrate that SGLT2i with canagliflozin preserves cardiac contractile function and efficiency during regional myocardial ischemia and provides ischemia protection independent of alterations in myocardial substrate utilization

Risk-shifting Through Issuer Liability and Corporate Monitoring

This article explores how issuer liability re-allocates fraud risk and how risk allocation may reduce the incidence of fraud. In the US, the apparent absence of individual liability of officeholders and insufficient monitoring by insurers under-mine the potential deterrent effect of securities litigation. The underlying reasons why both mechanisms remain ineffective are collective action problems under the prevailing dispersed ownership structure, which eliminates the incentives to moni-tor set by issuer liability. This article suggests that issuer liability could potentially have a stronger deterrent effect when it shifts risk to individuals or entities holding a larger financial stake. Thus, it would enlist large shareholders in monitoring in much of Europe. The same risk-shifting effect also has implications for the debate about the relationship between securities litigation and creditor interests. Credi-tors’ claims should not be given precedence over claims of defrauded investors (e.g., because of the capital maintenance principle), since bearing some of the fraud risk will more strongly incentivise large creditors, such as banks, to monitor the firm in jurisdictions where corporate debt is relatively concentrated

Polymer-tethered glyconanoparticle colourimetric biosensors for lectin binding : structural and experimental parameters to ensure a robust output

Glycan–lectin interactions play essential roles in biology; as the site of attachment for pathogens, cell–cell communication, and as crucial players in the immune system. Identifying if a new glycan (natural or unnatural) binds a protein partner, or if a new protein (or mutant) binds a glycan remains a non-trivial problem, with few accessible or low-cost tools available. Micro-arrays allow for the interrogation of 100's of glycans but are not widely available in individual laboratories. Biophysical techniques such as isothermal titration calorimetry, surface plasmon resonance spectrometry, biolayer interferometry and nuclear magnetic resonance spectroscopy all provide detailed understanding of glycan binding but are relatively expensive. Glycosylated plasmonic nanoparticles based on gold cores with polymeric tethers have emerged as biosensors to detect glycan–protein binding, based on colourimetric (red to blue) outputs which can be easily interpreted by a simple UV-visible spectrometer or by eye. Despite the large number of reports there are no standard protocols for each system or recommended start points, to allow a new user to deploy this technology. Here we explore the key parameters of nanoparticle size, polymeric tether length and gold concentration to provide some guidelines for how polymer-tethered glycosylated gold nanoparticles can be used to probe a new glycan/protein interactions, with minimal optimisation barriers. This work aimed to remove the need to explore chemical and nanoparticle space and hence remove a barrier for other users when deploying this system. We show that the concentration of the gold core is crucial to balance strong responses versus false positives and recommend a gold core size and polymer tether length which balances sufficient colloidal stability and output. Whilst subtle differences between glycans/lectins will impact the outcomes, these parameters should enable a lab user to quickly evaluate binding using minimal quantities of the glycan and lectin, to select candidates for further study

Resummation of the Divergent Perturbation Series for a Hydrogen Atom in an Electric Field

We consider the resummation of the perturbation series describing the energy displacement of a hydrogenic bound state in an electric field (known as the Stark effect or the LoSurdo-Stark effect), which constitutes a divergent formal power series in the electric field strength. The perturbation series exhibits a rich singularity structure in the Borel plane. Resummation methods are presented which appear to lead to consistent results even in problematic cases where isolated singularities or branch cuts are present on the positive and negative real axis in the Borel plane. Two resummation prescriptions are compared: (i) a variant of the Borel-Pade resummation method, with an additional improvement due to utilization of the leading renormalon poles (for a comprehensive discussion of renormalons see [M. Beneke, Phys. Rep. vol. 317, p. 1 (1999)]), and (ii) a contour-improved combination of the Borel method with an analytic continuation by conformal mapping, and Pade approximations in the conformal variable. The singularity structure in the case of the LoSurdo-Stark effect in the complex Borel plane is shown to be similar to (divergent) perturbative expansions in quantum chromodynamics.Comment: 14 pages, RevTeX, 3 tables, 1 figure; numerical accuracy of results enhanced; one section and one appendix added and some minor changes and additions; to appear in phys. rev.

Comparison of Source Images for protons, π−\pi^-'s and Λ\Lambda's in 6 AGeV Au+Au collisions

Source images are extracted from two-particle correlations constructed from strange and non-strange hadrons produced in 6 AGeV Au + Au collisions. Very different source images result from pp vs p$\Lambda$ vs $\pi^-\pi^-$ correlations. These observations suggest important differences in the space-time emission histories for protons, pions and neutral strange baryons produced in the same events