Stochastic simulation for spatial modelling of dynamic process in a living cell

One of the fundamental motivations underlying computational cell biology is to gain insight into the complicated dynamical processes taking place, for example, on the plasma membrane or in the cytosol of a cell. These processes are often so complicated that purely temporal mathematical models cannot adequately capture the complex chemical kinetics and transport processes of, for example, proteins or vesicles. On the other hand, spatial models such as Monte Carlo approaches can have very large computational overheads. This chapter gives an overview of the state of the art in the development of stochastic simulation techniques for the spatial modelling of dynamic processes in a living cell.

Linear models for endocytic transformations from live cell imaging

Endocytosis is the process by which cells internalise molecules including nutrient proteins from the extracellular media. In one form, macropinocytosis, the membrane at the cell surface ruffles and folds over to give rise to an internalised vesicle. Negatively charged phospholipids within the membrane called phosphoinositides then undergo a series of transformations that are critical for the correct trafficking of the vesicle within the cell, and which are often pirated by pathogens such as Salmonella. Advanced fluorescent video microscopy imaging now allows the detailed observation and quantification of these events in live cells over time. Here we use these observations as a basis for building differential equation models of the transformations. An initial investigation of these interactions was modelled with reaction rates proportional to the sum of the concentrations of the individual constituents. A first order linear system for the concentrations results. The structure of the system enables analytical expressions to be obtained and the problem becomes one of determining the reaction rates which generate the observed data plots. We present results with reaction rates which capture the general behaviour of the reactions so that we now have a complete mathematical model of phosphoinositide transformations that fits the experimental observations. Some excellent fits are obtained with modulated exponential functions; however, these are not solutions of the linear system. The question arises as to how the model may be modified to obtain a system whose solution provides a more accurate fit.

Structural Systematics of the Anhydrous 1:1 Proton-Transfer Compounds of 3,5-Dinitrosalicylic Acid with Aniline and Monosubstituted Anilines

The crystal structures of the proton-transfer compounds of 3,5-dinitrosalicylic acid (DNSA) with a series of aniline-type Lewis bases [aniline, 2-hydroxyaniline, 2-methoxyaniline, 3-methoxyaniline, 4-fluoroaniline, 4-chloroaniline and 2-aminoaniline] have been determined and their hydrogen-bonding systems analysed. All are anhydrous 1:1 salts: [(C6H8N)+(C7H3N2O7)-], (1), [(C6H8NO)+(C7H3N2O7)-], (2), [(C7H10NO)+(C7H3N2O7)-], (3), [(C7H10NO)+(C7H3N2O7)-], (4), [(C6H7FN)+(C7H3N2O7)-], (5), [(C6H7ClN)+(C7H3N2O7)-], (6), and [(C6H9N2)+(C7H3N2O7)-], (7) respectively. Crystals of 1 and 6 are triclinic, space group P-1 while the remainder are monoclinic with space group either P21/n (2, 4, 5 and 7) or P21 (3). Unit cell dimensions and contents are: for 1, a = 7.2027(17), b = 7.5699(17), c = 12.9615(16) Å, α = 84.464(14), β = 86.387(15), γ = 75.580(14)o, Z = 2; for 2, a = 7.407(3), b = 6.987(3), c = 27.653(11) Å, β = 94.906(7)o, Z = 4; for 3, a = 8.2816(18), b = 23.151(6), c = 3.9338(10), β = 95.255(19)o, Z = 2; for 4, a = 11.209(2), b = 8.7858(19), c = 15.171(3) Å, β = 93.717(4)o, Z = 4; for 5, a = 26.377(3), b = 10.1602(12), c = 5.1384(10) Å, β = 91.996(13)o, Z = 4; for 6, a = 11.217(3), b = 14.156(5), c = 4.860(3) Å, α = 99.10(4), β = 96.99(4), γ = 76.35(2)o, Z = 2; for 7, a = 12.830(4), b = 8.145(3), c = 14.302(4) Å, β = 102.631(6)o, Z = 4. In all compounds at least one primary linear intermolecular N+-H…O(carboxyl) hydrogen-bonding interaction is present which, together with secondary hydrogen bonding results in the formation of mostly two-dimensional network structures, exceptions being with compounds 4 and 5 (one-dimensional) and compound 6 (three-dimensional). In only two cases [compounds 1 and 4], are weak cation-anion or cation-cation π-π interactions found while weak aromatic C-H…O interactions are insignificant. The study shows that all compounds fit the previously formulated classification scheme for primary and secondary interactive modes for proton-transfer compounds of 3,5-dinitrosalicylic acid but there are some unusual variants.

Active regulation of the epidermal calcium profile

A distinct calcium profile is strongly implicated in regulating the multi-layered structure of the epidermis. However, the mechanisms that govern the regulation of this calcium profile are currently unclear. It clearly depends on the relatively impermeable barrier of the stratum corneum (passive regulation) but may also depend on calcium exchanges between keratinocytes and extracellular fluid (active regulation). Using a mathematical model that treats the viable sublayers of unwounded human and murine epidermis as porous media and assumes that their calcium profiles are passively regulated, we demonstrate that these profiles are also actively regulated. To obtain this result, we found that diffusion governs extracellular calcium motion in the viable epidermis and hence intracellular calcium is the main source of the epidermal calcium profile. Then, by comparison with experimental calcium profiles and combination with a hypothesised cell velocity distribution in the viable epidermis, we found that the net influx of calcium ions into keratinocytes from extracellular fluid may be constant and positive throughout the stratum basale and stratum spinosum, and that there is a net outflux of these ions in the stratum granulosum. Hence the calcium exchange between keratinocytes and extracellular fluid differs distinctly between the stratum granulosum and the underlying sublayers, and these differences actively regulate the epidermal calcium profile. Our results also indicate that plasma membrane dysfunction may be an early event during keratinocyte disintegration in the stratum granulosum.

VAMP3 regulates podosome organisation in macrophages and together with Stx4/SNAP23 mediates adhesion, cell spreading and persistent migration

The ability of cells to adhere, spread and migrate is essential to many physiological processes, particularly in the immune system where cells must traffic to sites of inflammation and injury. By altering the levels of individual components of the VAMP3/Stx4/SNAP23 complex we show here that this SNARE complex regulates efficient macrophage adhesion, spreading and migration on fibronectin. During cell spreading this complex mediates the polarised exocytosis of VAMP3- positive recycling endosome membrane into areas of membrane expansion, where VAMP3's surface partner Q-SNARE complex Stx4/SNAP23 was found to accumulate. Lowering the levels of VAMP3 in spreading cells resulted in a more rounded cell morphology and most cells were found to be devoid of the typical ring-like podosome superstructures seen normally in spreading cells. In migrating cells lowering VAMP3 levels disrupted the polarised localisation of podosome clusters. The reduced trafficking of recycling endosome membrane to sites of cell spreading and the disorganised podosome localisation in migrating macrophages greatly reduced their ability to persistently migrate on fibronectin. Thus, this important SNARE complex facilitates macrophage adhesion, spreading, and persistent macrophage migration on fibronectin through the delivery of VAMP3-positive membrane with its cargo to expand the plasma membrane and to participate in organising adhesive podosome structures.

Recycling endosome membrane incorporation into the leading edge regulates lamellipodia formation and macrophage migration

In comparison to our knowledge of the recycling of adhesion receptors and actin assembly, exactly how the cell controls its surface membrane to form a lamellipodium during migration is poorly understood. Here, we show the recycling endosome membrane is incorporated into the leading edge of a migrating cell to expand lamellipodia membrane. We have identified the SNARE complex that is necessary for fusion of the recycling endosome with the cell surface, as consisting of the R-SNARE VAMP3 on the recycling endosome partnering with the surface Q-SNARE Stx4/SNAP23, which was found to translocate and accumulate on the leading edge of migrating cells. Increasing VAMP3-mediated fusion of the recycling endosome with the surface increased membrane ruffling, while inhibition of VAMP3-mediated fusion showed that incorporation of the recycling endosome is necessary for efficient lamellipodia formation. At the same time, insertion of this recycling endosome membrane also delivers its cargo integrin α5β1 to the cell surface. The loss of this extra membrane for lamellipodia expansion and delivery of cargo in cells resulted in macrophages with a diminished capacity to effectively migrate. Thus, the recycling endosome membrane is incorporated into the leading edge and this aids expansion of the lamellipodia and simultaneously delivers integrins necessary for efficient cell migration.

A biphasic scaffold design combined with cell sheet technology for simultaneous regeneration of alveolar bone/periodontal ligament complex

This study describes the design of a biphasic scaffold composed of a Fused Deposition Modeling scaffold (bone compartment) and an electrospun membrane (periodontal compartment) for periodontal regeneration. In order to achieve simultaneous alveolar bone and periodontal ligament regeneration a cell-based strategy was carried out by combining osteoblast culture in the bone compartment and placement of multiple periodontal ligament (PDL) cell sheets on the electrospun membrane. In vitro data showed that the osteoblasts formed mineralized matrix in the bone compartment after 21 days in culture and that the PDL cell sheet harvesting did not induce significant cell death. The cell-seeded biphasic scaffolds were placed onto a dentin block and implanted for 8 weeks in an athymic rat subcutaneous model. The scaffolds were analyzed by μCT, immunohistochemistry and histology. In the bone compartment, a more intense ALP staining was obtained following seeding with osteoblasts, confirming the μCT results which showed higher mineralization density for these scaffolds. A thin mineralized cementum-like tissue was deposited on the dentin surface for the scaffolds incorporating the multiple PDL cell sheets, as observed by H&E and Azan staining. These scaffolds also demonstrated better attachment onto the dentin surface compared to no attachment when no cell sheets were used. In addition, immunohistochemistry revealed the presence of CEMP1 protein at the interface with the dentine. These results demonstrated that the combination of multiple PDL cell sheets and a biphasic scaffold allows the simultaneous delivery of the cells necessary for in vivo regeneration of alveolar bone, periodontal ligament and cementum. © 2012 Elsevier Ltd.

Incipient neurovulnerability and neuroprotection in early psychosis : a proton spectroscopy study of the anterior hippocampus at 3Tesla

Programmed cell death (PCD) and progenitor cell generation (of glial and in some brain areas also neuronal fate) in the CNS is an active process throughout life and is generally not associated with gliosis which means that PCD can be pathologically silent. The striking discovery that progenitor cell generation (of glial and in some brain areas neuronal fate) is widespread in the adult CNS (especially the hippocampus) suggest a much more dynamic scenario than previously thought and transcends the dichotomy between neurodevelopmental and neurodegenerative models of schizophrenia and related disorders. We suggest that the regulatory processes that control the regulation of PCD and the generation of progenitor cells may be disturbed in the early phase of psychotic disorders underpinning a disconnectivity syndrom at the onset of clinically overt disorders. An ongoing 1H-MRS study of the anterior hippocampus at 3 Tesla in mostly drug-naive first-episode psychosis patients suggests no change in NAA, but significant increases in myo-inositol and lactate. The data suggests that neuronal integrity in the anterior hippocampus is still intact at the early stage of illness or mainly only functionally impaired. However the increase in lactate and myo-inositol may reflect a potential disturbance of generation and PCD of progenitor cells (of glial and in selected brain areas also neuronal fate) at the onset of psychosis. If true the use of neuroprotective agents such as lithium or eicosapentaenoic acid (which inhibit PCD and support cell generation)in the early phase of psychotic disorders may be a potent treatment avenue to explore.

Effect Of changes in sodium intake on cell transport parameters

Changing sodium intake from 70-200 mmol/day elevates blood pressure in normotensive volunteers by 6/4 mmHg. Older people, people with reduced renal function on a low sodium diet and people with a family history of hypertension are more likely to show this effect. The rise in blood pressure was associated with a fall in plasma volume suggesting that plasma volume changes do not initiate hypertension. In normotensive individuals the most common abnormality in membrane sodium transport induced by an extra sodium load was an increased permeability of the red cell to sodium. Some normotensive individuals also had an increase in the level of a plasma inhibitor that inhibited Na-K ATPase. These individuals also appeared to have a rise in blood pressure. Sodium intake and blood pressure are related. The relationship differs in different people and is probably controlled by the genetically inherited capacity of systems involved in membrane sodium transport.

Human Kallikrein 4 signal peptide induces cytotoxic T cell responses in healthy donors and prostate cancer patients.

Immunotherapy is a promising new treatment for patients with advanced prostate and ovarian cancer, but its application is limited by the lack of suitable target antigens that are recognized by CD8+ cytotoxic T lymphocytes (CTL). Human kallikrein 4 (KLK4) is a member of the kallikrein family of serine proteases that is significantly overexpressed in malignant versus healthy prostate and ovarian tissue, making it an attractive target for immunotherapy. We identified a naturally processed, HLA-A*0201-restricted peptide epitope within the signal sequence region of KLK4 that induced CTL responses in vitro in most healthy donors and prostate cancer patients tested. These CTL lysed HLA-A*0201+ KLK4 + cell lines and KLK4 mRNA-transfected monocyte-derived dendritic cells. CTL specific for the HLA-A*0201-restricted KLK4 peptide were more readily expanded to a higher frequency in vitro compared to the known HLA-A*0201-restricted epitopes from prostate cancer antigens; prostate-specific antigen (PSA), prostate-specific membrane antigen (PSMA) and prostatic acid phosphatase (PAP). These data demonstrate that KLK4 is an immunogenic molecule capable of inducing CTL responses and identify it as an attractive target for prostate and ovarian cancer immunotherapy.

Stratum basale keratinocyte expression of the cell surface glycoprotein CDCP1 during epidermogenesis and its role in keratinocyte migration

Background: Epidermogenesis and epidermal wound healing are tightly regulated processes during which keratinocytes must migrate, proliferate and differentiate. Cell to cell adhesion is crucial to the initiation and regulation of these processes. CUB domain containing protein 1 (CDCP1) is a transmembrane glycoprotein that is differentially tyrosine phosphorylated during changes in cell adhesion and survival signalling and is expressed by keratinocytes in native human skin, as well as in primary cultures. Objectives: To investigate the expression of CDCP1 during epidermogenesis and its role in keratinocyte migration. Methods: We examined both human skin tissue and an in vitro three-dimensional human skin equivalent model to examine the expression of CDCP1 during epidermogenesis. To examine the role of CDCP1 in keratinocyte migration we used a function blocking anti-CDCP1 antibody and a real-time Transwell™ cell migration assay. Results: Immunohistochemical analysis indicated that in native human skin CDCP1 is expressed in the stratum basale and stratum spinosum. In contrast, during epidermogenesis in a 3-dimensional human skin equivalent model CDCP1 was expressed only in the stratum basale, with localization restricted to the cell-cell membrane. No expression was detected in basal keratinocytes that were in contact with the basement membrane. Further, an anti-CDCP1 function blocking antibody was shown to disrupt keratinocyte chemotactic migration in vitro. Conclusions: These findings delineate the expression of CDCP1 in human epidermal keratinocytes during epidermogenesis and demonstrate that CDCP1 is involved in keratinocyte migration.

Serine protease inhibition and mitochondrial dysfunction associated with cisplatin resistance in human tumor cell lines : targets for therapy

Indicators of mitochondrial function were studied in two different cell culture models of cis-diamminedichloroplatinum-II (CDDP) resistance: the intrinsically resistant human ovarian cancer cell line CI-80-13S, and resistant clones (HeLa-S1a and HeLa-S1b) generated by stable expression of the serine protease inhibitor—plasminogen activator inhibitor type-2 (PAI-2), in the human cervical cancer cell line HeLa. In both models, CDDP resistance was associated with sensitivity to killing by adriamycin, etoposide, auranofin, bis[1,2-bis(diphenylphosphino)ethane]gold(I) chloride {[Au(DPPE)2]Cl}, CdCl2 and the mitochondrial inhibitors rhodamine-123 (Rhl23), dequalinium chloride (DeCH), tetraphenylphosphonium (TPP), and ethidium bromide (EtBr) and with lower constitutive levels of ATP. Unlike the HeLa clones, CI-80-13S cells were additionally sensitive to chloramphenicol, 1-methyl-4-phenylpyridinium ion (MPP+), rotenone, thenoyltrifluoroacetone (TTFA), and antimycin A, and showed poor reduction of 1-[4,5-dimethylthiazol-2-yl]-2,5-diphenyltetrazolium bromide (MTT), suggesting a deficiency in NADH dehydrogenase and/or succinate dehydrogenase activities. Total platinum uptake and DNA-bound platinum were slightly lower in CI-80-13S than in sensitive cells. The HeLa-S1a and HeLa-S1b clones, on the other hand, showed poor reduction of triphenyltetrazolium chloride (TTC), indicative of low cytochrome c oxidase activity. Total platinum uptake by HeLa-S1a was similar to HeLa, but DNA-bound platinum was much lower than for the parent cell line. The mitochondria of CI-80-13S and HeLa-S1a showed altered morphology and were fewer in number than those of JAM and HeLa. In both models, CDDP resistance was associated with less platinum accumulation and with mitochondrial and membrane defects, brought about one case with expression of a protease inhibitor which is implicated in tumor progression. Such markers may identify tumors suitable for treatment with gold phosphine complexes or other mitochondrial inhibitors.

Polarized fluid movement and not cell death, creates luminal spaces in adult prostate epithelium

There are two predominant theories for lumen formation in tissue morphogenesis: cavitation driven by cell death, and membrane separation driven by epithelial polarity. To define the mechanism of lumen formation in prostate acini, we examined both theories in several cell lines grown in three-dimensional (3D) Matrigel culture. Lumen formation occurred early in culture and preceded the expression of cell death markers for apoptosis (active caspase 3) and autophagy (LC-3). Active caspase 3 was expressed by very few cells and inhibition of apoptosis did not suppress lumen formation. Despite LC-3 expression in all cells within a spheroid, this was not associated with cell death. However, expression of a prostate-secretory protein coincided with lumen formation and subsequent disruption of polarized fluid movement led to significant inhibition of lumen formation. This work indicates that lumen formation is driven by the polarized movement of fluids and proteins in 3D prostate epithelial models and not by cavitation.

An in-vitro 3-D cell culture model for studying pathomechanisms in AMD

Purpose To develop a novel 3-D cell culture model with the view to studying the pathomechanisms underlying the development of age-related macular degeneration (AMD). Our central hypothesis is that the silk structural protein fibroin used in conjunction with cultured human cells can be used to mimic the structural relationships between the RPE and choriocapillaris in health and disease. Methods Co-cultures of human RPE cells (ARPE-19 cells grown in Miller’s medium) and microvascular endothelial cells (HMEC-1 cells grown in endothelial culture medium) were established on opposing sides of a synthetic Bruch’s membrane (3 microns thick) constructed from B mori silk fibroin. Cell attachment was facilitated by pre-coating the fibroin membrane with vitronectin (for ARPE-19 cells) and gelatin (for HMEC-1 cells) respectively. The effects of tropoelastin on attachment of ARPE-19 cells was also examined. Barrier function was examined by measurement of trans-epithelial resistance (TER) using a voltohmmeter (EVOM-2). The phagocytic activity of the synthetic RPE was tested using vitronectin-coated microspheres (2 micron diameter FluoSpheres). In some cultures, membrane defects were created by puncturing within a 24 G needle. The architecture of the synthetic tissue before and after wounding was examined by confocal microscopy after staining for ZO-1 and F-actin. Results The RPE layer of the 3D model developed a cobblestoned morphology (validated by staining for ZO-1 and F-actin), displayed barrier function (validated by measurement of TER) and demonstrated cytoplasmic uptake of vitronectin-coated microspheres. Attachment of ARPE-19 cells to fibroin was unaffected by tropoelastin. Microvascular endothelial cells attached well to the gelatin-coated surface of the fibroin membrane and remained physically separated from the overlaying RPE layer. The fibroin membranes were amenable to puncturing without collapse thus providing the opportunity to study transmembrane migration of the endothelial cells. Conclusions Synthetic Bruch’s membranes constructed from silk fibroin, vitronectin and gelatin, support the co-cultivation of RPE cells and microvascular endothelial cells. The resulting RPE layer displays functions similar to that of native RPE and the entire tri-layered structure displays potential to be used as an in vitro model of choroidal neovascularization.

Influence of fatty acid structure on fuel properties of algae derived biodiesel

Physical and chemical properties of biofuel are influenced by structural features of fatty acid such as chain length, degree of unsaturation and branching of the chain. A simple and reliable calculation method to estimate fuel property is therefore needed to avoid experimental testing which is difficult, costly and time consuming. Typically in commercial biodiesel production such testing is done for every batch of fuel produced. In this study 9 different algae species were selected that were likely to be suitable for subtropical climates. The fatty acid methyl esters (FAMEs) of all algae species were analysed and the fuel properties like cetane number (CN), cold filter plugging point (CFPP), kinematic viscosity (KV), density and higher heating value (HHV) were determined. The relation of each fatty acid with particular fuel property is analysed using multivariate and multi-criteria decision method (MCDM) software. They showed that some fatty acids have major influences on the fuel properties whereas others have minimal influence. Based on the fuel properties and amounts of lipid content rank order is drawn by PROMETHEE-GAIA which helped to select the best algae species for biodiesel production in subtropical climates. Three species had fatty acid profiles that gave the best fuel properties although only one of these (Nannochloropsis oculata) is considered the best choice because of its higher lipid content.