Adhesion failures determine the pattern of choroidal neovascularization in the eye: a computer simulation study.

Choroidal neovascularization (CNV) of the macular area of the retina is the major cause of severe vision loss in adults. In CNV, after choriocapillaries initially penetrate Bruch's membrane (BrM), invading vessels may regress or expand (CNV initiation). Next, during Early and Late CNV, the expanding vasculature usually spreads in one of three distinct patterns: in a layer between BrM and the retinal pigment epithelium (sub-RPE or Type 1 CNV), in a layer between the RPE and the photoreceptors (sub-retinal or Type 2 CNV) or in both loci simultaneously (combined pattern or Type 3 CNV). While most studies hypothesize that CNV primarily results from growth-factor effects or holes in BrM, our three-dimensional simulations of multi-cell model of the normal and pathological maculae recapitulate the three growth patterns, under the hypothesis that CNV results from combinations of impairment of: 1) RPE-RPE epithelial junctional adhesion, 2) Adhesion of the RPE basement membrane complex to BrM (RPE-BrM adhesion), and 3) Adhesion of the RPE to the photoreceptor outer segments (RPE-POS adhesion). Our key findings are that when an endothelial tip cell penetrates BrM: 1) RPE with normal epithelial junctions, basal attachment to BrM and apical attachment to POS resists CNV. 2) Small holes in BrM do not, by themselves, initiate CNV. 3) RPE with normal epithelial junctions and normal apical RPE-POS adhesion, but weak adhesion to BrM (e.g. due to lipid accumulation in BrM) results in Early sub-RPE CNV. 4) Normal adhesion of RBaM to BrM, but reduced apical RPE-POS or epithelial RPE-RPE adhesion (e.g. due to inflammation) results in Early sub-retinal CNV. 5) Simultaneous reduction in RPE-RPE epithelial binding and RPE-BrM adhesion results in either sub-RPE or sub-retinal CNV which often progresses to combined pattern CNV. These findings suggest that defects in adhesion dominate CNV initiation and progression.

A multi-cell, multi-scale model of vertebrate segmentation and somite formation.

Somitogenesis, the formation of the body's primary segmental structure common to all vertebrate development, requires coordination between biological mechanisms at several scales. Explaining how these mechanisms interact across scales and how events are coordinated in space and time is necessary for a complete understanding of somitogenesis and its evolutionary flexibility. So far, mechanisms of somitogenesis have been studied independently. To test the consistency, integrability and combined explanatory power of current prevailing hypotheses, we built an integrated clock-and-wavefront model including submodels of the intracellular segmentation clock, intercellular segmentation-clock coupling via Delta/Notch signaling, an FGF8 determination front, delayed differentiation, clock-wavefront readout, and differential-cell-cell-adhesion-driven cell sorting. We identify inconsistencies between existing submodels and gaps in the current understanding of somitogenesis mechanisms, and propose novel submodels and extensions of existing submodels where necessary. For reasonable initial conditions, 2D simulations of our model robustly generate spatially and temporally regular somites, realistic dynamic morphologies and spontaneous emergence of anterior-traveling stripes of Lfng. We show that these traveling stripes are pseudo-waves rather than true propagating waves. Our model is flexible enough to generate interspecies-like variation in somite size in response to changes in the PSM growth rate and segmentation-clock period, and in the number and width of Lfng stripes in response to changes in the PSM growth rate, segmentation-clock period and PSM length.

Front instabilities and invasiveness of simulated avascular tumors.

We study the interface morphology of a 2D simulation of an avascular tumor composed of identical cells growing in an homogeneous healthy tissue matrix (TM), in order to understand the origin of the morphological changes often observed during real tumor growth. We use the Glazier-Graner-Hogeweg model, which treats tumor cells as extended, deformable objects, to study the effects of two parameters: a dimensionless diffusion-limitation parameter defined as the ratio of the tumor consumption rate to the substrate transport rate, and the tumor-TM surface tension. We model TM as a nondiffusing field, neglecting the TM pressure and haptotactic repulsion acting on a real growing tumor; thus, our model is appropriate for studying tumors with highly motile cells, e.g., gliomas. We show that the diffusion-limitation parameter determines whether the growing tumor develops a smooth (noninvasive) or fingered (invasive) interface, and that the sensitivity of tumor morphology to tumor-TM surface tension increases with the size of the dimensionless diffusion-limitation parameter. For large diffusion-limitation parameters, we find a transition (missed in previous work) between dendritic structures, produced when tumor-TM surface tension is high, and seaweed-like structures, produced when tumor-TM surface tension is low. This observation leads to a direct analogy between the mathematics and dynamics of tumors and those observed in nonbiological directional solidification. Our results are also consistent with the biological observation that hypoxia promotes invasive growth of tumor cells by inducing higher levels of receptors for scatter factors that weaken cell-cell adhesion and increase cell motility. These findings suggest that tumor morphology may have value in predicting the efficiency of antiangiogenic therapy in individual patients.

Different ionic conditions prompt NHE2 and NHE3 translocation to the plasma membrane.

We tested whether NHE3 and NHE2 Na(+)/H(+) exchanger isoforms were recruited to the plasma membrane (PM) in response to changes in ion homeostasis. NHE2-CFP or NHE3-CFP fusion proteins were functional Na(+)/H(+) exchangers when transiently expressed in NHE-deficient PS120 fibroblasts. Confocal morphometry of cells whose PM was labeled with FM4-64 measured the fractional amount of fusion protein at the cell surface. In resting cells, 10-20% of CFP fluorescence was at PM and stable over time. A protocol commonly used to activate the Na(+)/H(+) exchange function (NH(4)-prepulse acid load sustained in Na(+)-free medium), increased PM percentages of PM NHE3-CFP and NHE2-CFP. Separation of cellular acidification from Na(+) removal revealed that only NHE3-CFP translocated when medium Na(+) was removed, and only NHE2-CFP translocated when the cell was acidified. NHE2/NHE3 chimeric proteins demonstrate that the Na(+)-removal response element resides predominantly in the NHE3 cytoplasmic tail and is distinct from the acidification response sequence of NHE2.

Adhesion between cells, diffusion of growth factors, and elasticity of the AER produce the paddle shape of the chick limb.

A central question in developmental biology is how cells interact to organize into tissues? In this paper, we study the role of mesenchyme-ectoderm interaction in the growing chick limb bud using Glazier and Graner's cellular Potts model, a grid-based stochastic framework designed to simulate cell interactions and movement. We simulate cellular mechanisms including cell adhesion, growth, and division and diffusion of morphogens, to show that differential adhesion between the cells, diffusion of growth factors through the extracellular matrix, and the elastic properties of the apical ectodermal ridge together can produce the proper shape of the limb bud.

Membrane insertion of betaine/GABA transporter during hypertonic stress correlates with nuclear accumulation of TonEBP.

MDCK cells stably transfected with betaine/GABA transporter tagged with EGFP (EGFP-BGT) were used to study plasma membrane insertion of EGFP-BGT. Adaptive response to hypertonicity requires nuclear migration of TonEBP. Confocal microscopy showed that after 6 h hypertonicity, the nuclear/cytoplasmic ratio of TonEBP fluorescence was increased to 2.4 compared to 1.4 in isotonic controls (P<0.001). The ratio in hypertonic cells was reduced by the proteasome inhibitor MG-132 in a dose-dependent way. Inhibition was 50% at 3 microM. After 6 h, hypertonicity expressed EGFP-BGT was localized in the plasma membrane, but there was no change in total EGFP-BGT abundance compared to isotonic controls. In contrast, EGFP-BGT remained mostly intracellular when 3 microM MG-132 was included in the hypertonic medium. The transport function of EGFP-BGT was studied as Na(+)-dependent uptake of [(3)H]GABA. This was not changed by MG-132 in isotonic controls, but MG-132 produced dose-dependent inhibition of hypertonic upregulation of Na(+)/GABA cotransport. Inhibition was 80% at 3 muM MG-132. Transport likely reflects membrane insertion of EGFP-BGT and there was a positive correlation (P<0.05) between Na(+)/GABA cotransport and the N/C ratio of TonEBP. Results are consistent with a role for TonEBP-mediated transcription in synthesis of additional proteins required for membrane insertion of EGFP-BGT protein.

Front instabilities and invasiveness of simulated 3D avascular tumors.

We use the Glazier-Graner-Hogeweg model to simulate three-dimensional (3D), single-phenotype, avascular tumors growing in an homogeneous tissue matrix (TM) supplying a single limiting nutrient. We study the effects of two parameters on tumor morphology: a diffusion-limitation parameter defined as the ratio of the tumor-substrate consumption rate to the substrate-transport rate, and the tumor-TM surface tension. This initial model omits necrosis and oxidative/hypoxic metabolism effects, which can further influence tumor morphology, but our simplified model still shows significant parameter dependencies. The diffusion-limitation parameter determines whether the growing solid tumor develops a smooth (noninvasive) or fingered (invasive) interface, as in our earlier two-dimensional (2D) simulations. The sensitivity of 3D tumor morphology to tumor-TM surface tension increases with the size of the diffusion-limitation parameter, as in 2D. The 3D results are unexpectedly close to those in 2D. Our results therefore may justify using simpler 2D simulations of tumor growth, instead of more realistic but more computationally expensive 3D simulations. While geometrical artifacts mean that 2D sections of connected 3D tumors may be disconnected, the morphologies of 3D simulated tumors nevertheless correlate with the morphologies of their 2D sections, especially for low-surface-tension tumors, allowing the use of 2D sections to partially reconstruct medically-important 3D-tumor structures.

3D multi-cell simulation of tumor growth and angiogenesis.

We present a 3D multi-cell simulation of a generic simplification of vascular tumor growth which can be easily extended and adapted to describe more specific vascular tumor types and host tissues. Initially, tumor cells proliferate as they take up the oxygen which the pre-existing vasculature supplies. The tumor grows exponentially. When the oxygen level drops below a threshold, the tumor cells become hypoxic and start secreting pro-angiogenic factors. At this stage, the tumor reaches a maximum diameter characteristic of an avascular tumor spheroid. The endothelial cells in the pre-existing vasculature respond to the pro-angiogenic factors both by chemotaxing towards higher concentrations of pro-angiogenic factors and by forming new blood vessels via angiogenesis. The tumor-induced vasculature increases the growth rate of the resulting vascularized solid tumor compared to an avascular tumor, allowing the tumor to grow beyond the spheroid in these linear-growth phases. First, in the linear-spherical phase of growth, the tumor remains spherical while its volume increases. Second, in the linear-cylindrical phase of growth the tumor elongates into a cylinder. Finally, in the linear-sheet phase of growth, tumor growth accelerates as the tumor changes from cylindrical to paddle-shaped. Substantial periods during which the tumor grows slowly or not at all separate the exponential from the linear-spherical and the linear-spherical from the linear-cylindrical growth phases. In contrast to other simulations in which avascular tumors remain spherical, our simulated avascular tumors form cylinders following the blood vessels, leading to a different distribution of hypoxic cells within the tumor. Our simulations cover time periods which are long enough to produce a range of biologically reasonable complex morphologies, allowing us to study how tumor-induced angiogenesis affects the growth rate, size and morphology of simulated tumors.