The modulatory effect of cell­cell contact on the tumourigenic potential of pre-malignant epithelial cells: a computational exploration.

Malignant development cannot be attributed alone to genetic changes in a single cell, but occurs as a result of the complex interplay between the failure of cellular regulation mechanisms and the presence of a permissive microenvironment. Although E-cadherin is classified as a 'metastasis suppressor' owing to its role in intercellular adhesion, the observation that it may be downregulated at a premalignant stage is indicative of additional roles in neoplastic development. We have used an agent-based computational model to explore the emergent behaviour resulting from the interaction of single and subpopulations of E-cadherin-compromised cells with unaffected normal epithelial cells within a monolayer environment. We have extended this to investigate the importance of local tissue perturbations in the form of scratch-wounding, or ablation of randomly-dispersed normal cells, on the growth of a single cell exhibiting E-cadherin loss. Our results suggest that the microenvironment with respect to localized cell density and normal/E-cadherin-compromised neighbours is crucial in determining whether an abnormal individual cell proliferates or remains dormant within the monolayer. These predictions raise important questions relating to the propensity for individual mutations to give rise to disease, and future experimental exploration of these will enhance our understanding of a complex, multifactorial pathological process.

Multi-scale interaction of particulate flow and the artery wall.

We discuss, from the perspective of basic science, the physical and biological processes which underlie atherosclerotic (plaque) initiation at the vascular endothelium, identifying the widely separated spatial and temporal scales which participate. We draw on current, related models of vessel wall evolution, paying particular attention to the role of particulate flow (blood is not a continuum fluid), and proceed to propose, then validate all the key components in a multiply-coupled, multi-scale modeling strategy (in qualitative terms only, note). Eventually, this strategy should lead to a quantitative, patient-specific understanding of the coupling between particulate flow and the endothelial state.

Anti-social cells: predicting the influence of E-cadherin loss on the growth of epithelial cell populations.

The characteristics of biological tissues are determined by the interactions of large numbers of autonomous cells. These interactions can be mediated remotely by diffusive biochemical factors, or by direct cell-cell contact. E-cadherin is a protein expressed on the surface of normal epithelial cells that plays a key role in mediating intercellular adhesion via calcium-dependent homotypic interactions. E-cadherin is a metastasis-suppressor protein and its loss of function is associated with malignant progression. The purpose of this study was to apply an agent-based simulation paradigm in order to examine the emergent growth properties of mixed populations consisting of normal and E-cadherin defective cells in monolayer cell culture. Specifically, we have investigated the dynamics of normal cell:cell interactions in terms of intercellular adhesion and migration, and have used a simplified rule to represent the concepts of juxtacrine epidermal growth factor receptor (EGFR) activation and subsequent effect on cell proliferation. This cellular level control determines the overall population growth in a simulated experiment. Our approach provides a tool for modelling the development of defined biological abnormalities in epithelial and other biological tissues, raising novel predictions for future experimental testing. The results predict that even a relatively small number of abnormal ('anti-social') cells can modify the rate of the total population expansion, but the magnitude of this effect also depends on the extrinsic (culture) environment. In addition to directly influencing population dynamics, 'anti-social' cells can also disrupt the behaviour of the normal cells around them.

Safety evaluation of polydextrose in infant formula using a suckling piglet model.

Oligosaccharides, the third largest component in human milk, are virtually absent from cow's milk and most infant formula. Prebiotic carbohydrates like polydextrose (PDX) have been proposed as surrogates for human milk oligosaccharides. Safety assessments of novel infant formula ingredients include dose-response experiments in appropriate neonatal animal models such as the suckling pig. To further substantiate the safety of the ingredient, one-day old pigs were fed a cow's milk-based formula supplemented with PDX (1.7, 4.3, 8.5 or 17 g/L) for 18 days (n=13/dose) and compared to appropriate control (unsupplemented formula; n=13) and reference groups (day 0 pigs, and sow-reared pigs; n=13). Growth rate, formula intake, stool consistency, behavior score, blood chemistry and hematology, relative organ weights (% of body weight), tissue morphology (i.e. liver, kidney and pancreas) and pancreas biochemistry did not differ among formula-fed pigs (P>0.1). Polydextrose mimicked other prebiotic carbohydrates and had no adverse effect at the highest tested level 17.0 g PDX/L, equivalent to a dose of 8.35 g/kg of body weight per day.

From pathway to population--a multiscale model of juxtacrine EGFR-MAPK signalling.

BACKGROUND: Most mathematical models of biochemical pathways consider either signalling events that take place within a single cell in isolation, or an 'average' cell which is considered to be representative of a cell population. Likewise, experimental measurements are often averaged over populations consisting of hundreds of thousands of cells. This approach ignores the fact that even within a genetically-homogeneous population, local conditions may influence cell signalling and result in phenotypic heterogeneity. We have developed a multi-scale computational model that accounts for emergent heterogeneity arising from the influences of intercellular signalling on individual cells within a population. Our approach was to develop an ODE model of juxtacrine EGFR-ligand activation of the MAPK intracellular pathway and to couple this to an agent-based representation of individual cells in an expanding epithelial cell culture population. This multi-scale, multi-paradigm approach has enabled us to simulate Extracellular signal-regulated kinase (Erk) activation in a population of cells and to examine the consequences of interpretation at a single cell or population-based level using virtual assays. RESULTS: A model consisting of a single pair of interacting agents predicted very different Erk activation (phosphorylation) profiles, depending on the formation rate and stability of intercellular contacts, with the slow formation of stable contacts resulting in low but sustained activation of Erk, and transient contacts resulting in a transient Erk signal. Extension of this model to a population consisting of hundreds to thousands of interacting virtual cells revealed that the activated Erk profile measured across the entire cell population was very different and may appear to contradict individual cell findings, reflecting heterogeneity in population density across the culture. This prediction was supported by immunolabelling of an epithelial cell population grown in vitro, which confirmed heterogeneity of Erk activation. CONCLUSION: These results illustrate that mean experimental data obtained from analysing entire cell populations is an oversimplification, and should not be extrapolated to deduce the signal:response paradigm of individual cells. This multi-scale, multi-paradigm approach to biological simulation provides an important conceptual tool in addressing how information may be integrated over multiple scales to predict the behaviour of a biological system.

Modelling the interaction of haemodynamics and the artery wall: current status and future prospects.

Clinical research has historically focused on the two main strategies of in vivo and in vitro experimentation. The concept of applying scientific theory to direct clinical applications is relatively recent. In this paper we focus on the interaction of wall shear stress with the endothelium and discuss how 'state of the art' computer modelling techniques can provide valuable data to aid understanding. Such data may be used to inform experiment and further, may help identify the key features of this complex system. Current emphasis is on coupling haemodynamics with models of biological phenomena to test hypotheses or predict the likely outcome of a disease or an intervention. New technologies to enable the integration of models of different types, levels of complexity and scales, are being developed. As will be discussed, the ultimate goal is the translation of this technology to the clinical arena.

Electrical bioimpedance readings increase with higher pressure applied to the measuring probe.

Electrical bioimpedance spectroscopy (EBIS) is a technique that uses a probe to calculate the transfer impedance from tissues. This transfer impedance can give information about the normal or pathological condition of the tissue. To take readings, pressure has to be applied to the probe in order to get a good contact between the electrodes and the tissue. We have been using EBIS to investigate the early diagnosis of dysplasia and cancer in the human cervix, oesophagus and bladder. We have found that, with increasing pressure (range used here was approximately 1 kPa to approximately 50 kPa), the resistivity readings increase in a consistent way up to 80%. In this paper, we show how this is a case in three different tissue types (oesophageal, gastric and vesical samples). These increases can be higher than those associated with the pathological changes that we are investigating (non-inflamed columnar tissue, for instance, shows values 50% higher than dysplastic columnar tissue). Finite-element modelling was also used to investigate the effect of volume reduction in the connective tissue or stroma. This simulation suggests no strong correlation between reduction of this structure and increase in resistivity. We hypothesize therefore that these changes may be mainly associated with the squeezing of water from the extracellular space. Finally, as pressure is difficult to control by hand, we raise the issue of the necessity of considering this variable when making EIS measurements.

Modelling the electrical properties of bladder tissue--quantifying impedance changes due to inflammation and oedema.

Electrical impedance spectroscopy has been developed as a potential method for the diagnosis of carcinoma in epithelial tissues. An understanding of the influence of structural changes in the tissue on the properties measured using this technique is essential for interpreting measured data and optimization of probe design. In contrast to other tissue types, carcinoma in situ of the bladder gives rise to an increase in electrical impedance over the kHz-MHz frequency range in comparison to normal tissue. Finite element models of the urothelium and the underlying superficial lamina propria have been constructed and solved in order to ascertain the influence of structural changes associated with malignancy, oedema and inflammation on the measured electrical properties of the tissue. Sensitivity analysis of results from a composite tissue model suggests that the increase in lymphocyte density in the lamina propria associated with an inflammatory response to the infiltration of urine into the tissue may explain these unusual electrical properties.

Development and validation of computational models of cellular interaction.

In this paper we take the view that computational models of biological systems should satisfy two conditions - they should be able to predict function at a systems biology level, and robust techniques of validation against biological models must be available. A modelling paradigm for developing a predictive computational model of cellular interaction is described, and methods of providing robust validation against biological models are explored, followed by a consideration of software issues.

Agent-based computational modeling of wounded epithelial cell monolayers.

Computational modeling of biological systems, or in silico biology, is an emerging tool for understanding structure and order in biological tissues. Computational models of the behavior of epithelial cells in monolayer cell culture have been developed and used to predict the healing characteristics of scratch wounds made to urothelial cell cultures maintained in low- and physiological [Ca2+] environments. Both computational models and in vitro experiments demonstrated that in low exogenous [Ca2+], the closure of 500-microm scratch wounds was achieved primarily by cell migration into the denuded area. The wound healing rate in low (0.09 mM) [Ca2+] was approximately twice as rapid as in physiological (2 mM) [Ca2+]. Computational modeling predicted that in cell cultures that are actively proliferating, no increase in the fraction of cells in the S-phase would be expected, and this conclusion was supported experimentally in vitro by bromodeoxyuridine incorporation assay. We have demonstrated that a simple rule-based model of cell behavior, incorporating rules relating to contact inhibition of proliferation and migration, is sufficient to qualitatively predict the calcium-dependent pattern of wound closure observed in vitro. Differences between the in vitro and in silico models suggest a role for wound-induced signaling events in urothelial cell cultures.

The epitheliome: agent-based modelling of the social behaviour of cells.

We have developed a new computational modelling paradigm for predicting the emergent behaviour resulting from the interaction of cells in epithelial tissue. As proof-of-concept, an agent-based model, in which there is a one-to-one correspondence between biological cells and software agents, has been coupled to a simple physical model. Behaviour of the computational model is compared with the growth characteristics of epithelial cells in monolayer culture, using growth media with low and physiological calcium concentrations. Results show a qualitative fit between the growth characteristics produced by the simulation and the in vitro cell models.

Modelling of epithelial tissue impedance measured using three different designs of probe.

Impedance measurement is a promising technique for detecting pre-malignant changes in epithelial tissue. This paper considers how the design of the impedance probe affects the ability to discriminate between tissue types. To do this, finite element models of the electrical properties of squamous and glandular columnar epithelia have been used. The glandular tissue model is described here for the first time. Glandular mucosa is found in many regions of the gastrointestinal tract, such as the stomach and intestine, and has a large effective surface area. Firstly, the electrical properties of a small section of gland, with epithelial cells and supportive tissue, are determined. These properties are then used to build up a three-dimensional model of a whole section of mucosa containing many thousands of glands. Measurements using different types of impedance probe were simulated by applying different boundary conditions to the models. Transepithelial impedance, and tetrapolar measurement with a probe placed on the tissue surface have been modelled. In the latter case, the impedance can be affected by conductive fluid, such as mucus, on the tissue surface. This effect has been investigated, and a new design of probe, which uses a guard electrode to counteract this potential source of variability, is proposed.

A study of the morphological parameters of cervical squamous epithelium.

Electrical impedance spectroscopy is a technique that has been investigated as a potential method for the diagnosis of epithelial carcinomas. Finite element modelling can provide an insight into the patterns of current flow in normal and pathological epithelium and hence aid in the process of probe design optimization. In order to develop a finite element model of the structure of normal and precancerous cervical squamous epithelium, it was first necessary to obtain the mean values and ranges of a number of morphological tissue parameters. The most important parameters in discriminating normal from neoplastic tissue were identified as being cell size and shape distribution, nuclear-to-cytoplasmic volume ratio and volume of extracellular space. A survey of the literature revealed an absence of reliable quantitative data for these parameters. We therefore present the results of our own basic image analysis on normal and pathological tissue sections, which we hope will be of use to other workers wishing to model cervical squamous epithelium, or other similar tissue structures.

Modelled current distribution in cervical squamous tissue.

The electrical properties of cervical squamous epithelium have been modelled in the frequency range 100 Hz to 10 MHz. The hierarchical modelling process comprises a cellular level stage, which includes detailed models of cells typical of different depths within the epithelium and a tissue model, which utilizes electrical properties obtained from the cellular models. The fit between the modelled and measured impedance spectra and the distribution of current with depth depends on the macroscopic model structure. Both the properties of the basement membrane and the presence of a surface mucus layer are shown to have a significant effect. The best fit with measured data is obtained when a 10 microm thick, high-conductivity surface layer is included in the tissue model.

Dermal fibroblast morphology is affected by stretching and not by C48/80.

Both stretching and C48/80 have been hypothesized to cause disruption of cell-matrix adhesions and thereby affect the dynamics of fluid balance in tissues. We investigated the effect of sinusoidal stretching and/or C48/80 on the morphology of fibroblasts in skin excised from the backs of Wistar-Möller rats in order to assess how these stimuli affect cellular interactions in tissues. Tissue samples were either soaked in Krebs' buffer with and without C48/80, or sinusoidally stretched (20% strain) in buffer with and without C48/80. Control skin was fixed immediately after excision. All tissues were processed for transmission electron microscopy. Morphometric analyses demonstrated that sinusoidal stretching of the skin results in the retraction or disruption of fibroblast cytoplasmic extensions, rounding up of the cell bodies and subsequently in increased tissue water content. C48/80 had no apparent effect on fibroblast morphology and adherence in tissues.

Privacy considerations in the context of an Australian observational database.

Observational databases are increasingly acknowledged for their value in clinical investigation. Australian general practice in particular presents an exciting opportunity to examine treatment in a natural setting. The paper explores issues such as privacy and confidentiality--foremost considerations when conducting this form of pharmacoepidemiological research. Australian legislation is currently addressing these exact issues in order to establish clear directives regarding ethical concerns. The development of a pharmacoepidemiological database arising from the integration of computerized Australian general practice records is described in addition, to the challenges associated with creating a database which considers patient privacy. The database known as 'Medic-GP', presently contains more than 950,000 clinical notes (including consultations, pathology, diagnostic imaging and adverse reactions) over a 5-year time period and relates to 55,000 patients. The paper then details a retrospective study which utilized the database to examine the interaction between antibiotic prescribing and patient outcomes from a community perspective, following a policy intervention. This study illustrates the application of computerized general practice records in research.

Analysis of tight junctions during neutrophil transendothelial migration.

Intercellular junctions have long been considered the main sites through which adherent neutrophils (PMNs) penetrate the endothelium. Tight junctions (TJs; zonula occludens) are the most apical component of the intercellular cleft and they form circumferential belt-like regions of intimate contact between adjacent endothelial cells. Whether PMN transmigration involves disruption of the TJ complex is unknown. We report here that endothelial TJs appear to remain intact during PMN adhesion and transmigration. Human umbilical vein endothelial cell (HUVEC) monolayers, a commonly used model for studying leukocyte trafficking, were cultured in astrocyte-conditioned medium to enhance TJ expression. Immunofluorescence microscopy and immunoblot analysis showed that activated PMN adhesion to resting monolayers or PMN migration across interleukin-1-treated monolayers does not result in widespread proteolytic loss of TJ proteins (ZO-1, ZO-2, and occludin) from endothelial borders. Ultrastructurally, TJs appear intact during and immediately following PMN transendothelial migration. Similarly, transendothelial electrical resistance is unaffected by PMN adhesion and migration. Previously, we showed that TJs are inherently discontinuous at tricellular corners where the borders of three endothelial cells meet and PMNs migrate preferentially at tricellular corners. Collectively, these results suggest that PMN migration at tricellular corners preserves the barrier properties of the endothelium and does not involve widespread disruption of endothelial TJs.

The groESL chaperone operon of Lactobacillus johnsonii.

The Lactobacillus johnsonii VPI 11088 groESL operon was localized on the chromosome near the insertion element IS1223. The operon was initially cloned as a series of three overlapping PCR fragments, which were sequenced and used to design primers to amplify the entire operon. The amplified fragment was used as a probe to recover the chromosomal copy of the groESL operon from a partial library of L. johnsonii VPI 11088 (NCK88) DNA, cloned in the shuttle vector pTRKH2. The 2,253-bp groESL fragment contained three putative open reading frames, two of which encoded the ubiquitous GroES and GroEL chaperone proteins. Analysis of the groESL promoter region revealed three transcription initiation sites, as well as three sets of inverted repeats (IR) positioned between the transcription and translation start sites. Two of the three IR sets bore significant homology to the CIRCE elements, implicated in negative regulation of the heat shock response in many bacteria. Northern analysis and primer extension revealed that multiple temperature-sensitive promoters preceded the groESL chaperone operon, suggesting that stress protein production in L. johnsonii is strongly regulated. Maximum groESL transcription activity was observed following a shift to 55 degrees C, and a 15 to 30-min exposure of log-phase cells to this temperature increased the recovery of freeze-thawed L. johnsonii VPI 11088. These results suggest that a brief, preconditioning heat shock can be used to trigger increased chaperone production and provide significant cross-protection from the stresses imposed during the production of frozen culture concentrates.

P-selectin mediates neutrophil adhesion to endothelial cell borders.

During an acute inflammatory response, endothelial P-selectin (CD62P) can mediate the initial capture of neutrophils from the free flowing bloodstream. P-selectin is stored in secretory granules (Weibel-Palade bodies) and is rapidly expressed on the endothelial surface after stimulation with histamine or thrombin. Because neutrophil transmigration occurs preferentially at endothelial borders, we wished to determine whether P-selectin-dependent neutrophil capture (adhesion) occurs at endothelial cell borders. Under static or hydrodynamic flow (2 dyn/cm2) conditions, histamine (10(-4) M) or thrombin (0.2 U/mL) treatment induced preferential (> or = 75%) neutrophil adhesion to the cell borders of endothelial monolayers. Blocking antibody studies established that neutrophil adhesion was completely P-selectin dependent. P-selectin surface expression increased significantly after histamine treatment and P-selectin immunostaining was concentrated along endothelial borders. We conclude that preferential P-selectin expression along endothelial borders may be an important mechanism for targeting neutrophil migration at endothelial borders.