Modeling ion channel dynamics through reflected stochastic differential equations

Ion channels are membrane proteins that open and close at random and play a vital role in the electrical dynamics of excitable cells. The stochastic nature of the conformational changes these proteins undergo can be significant, however current stochastic modeling methodologies limit the ability to study such systems. Discrete-state Markov chain models are seen as the "gold standard," but are computationally intensive, restricting investigation of stochastic effects to the single-cell level. Continuous stochastic methods that use stochastic differential equations (SDEs) to model the system are more efficient but can lead to simulations that have no biological meaning. In this paper we show that modeling the behavior of ion channel dynamics by a reflected SDE ensures biologically realistic simulations, and we argue that this model follows from the continuous approximation of the discrete-state Markov chain model. Open channel and action potential statistics from simulations of ion channel dynamics using the reflected SDE are compared with those of a discrete-state Markov chain method. Results show that the reflected SDE simulations are in good agreement with the discrete-state approach. The reflected SDE model therefore provides a computationally efficient method to simulate ion channel dynamics while preserving the distributional properties of the discrete-state Markov chain model and also ensuring biologically realistic solutions. This framework could easily be extended to other biochemical reaction networks. © 2012 American Physical Society.

Application of a human bone engineering platform to an in vitro and in vivo breast cancer metastasis model

Breast cancer in its advanced stage has a high predilection to the skeleton. Currently, treatment options of breast cancer-related bone metastasis are restricted to only palliative therapeutic modalities. This is due to the fact that mechanisms regarding the breast cancer celI-bone colonisation as well as the interactions of breast cancer cells with the bone microenvironment are not fully understood, yet. This might be explained through a lack of appropriate in vitro and in vivo models that are currently addressing the above mentioned issue. Hence the hypothesis that the translation of a bone tissue engineering platform could lead to improved and more physiological in vitro and in vivo model systems in order to investigate breast cancer related bone colonisation was embraced in this PhD thesis. Therefore the first objective was to develop an in vitro model system that mimics human mineralised bone matrix to the highest possible extent to examine the specific biological question, how the human bone matrix influences breast cancer cell behaviour. Thus, primary human osteoblasts were isolated from human bone and cultured under osteogenic conditions. Upon ammonium hydroxide treatment, a cell-free intact mineralised human bone matrix was left behind. Analyses revealed a similar protein and mineral composition of the decellularised osteoblast matrix to human bone. Seeding of a panel of breast cancer cells onto the bone mimicking matrix as well as reference substrates like standard tissue culture plastic and collagen coated tissue culture plastic revealed substrate specific differences of cellular behaviour. Analyses of attachment, alignment, migration, proliferation, invasion, as well as downstream signalling pathways showed that these cellular properties were influenced through the osteoblast matrix. The second objective of this PhD project was the development of a human ectopic bone model in NOD/SCID mice using medical grade polycaprolactone tricalcium phosphate (mPCL-TCP) scaffold. Human osteoblasts and mesenchymal stem cells were seeded onto an mPCL-TCP scaffold, fabricated using a fused deposition modelling technique. After subcutaneous implantation in conjunction with the bone morphogenetic protein 7, limited bone formation was observed due to the mechanical properties of the applied scaffold and restricted integration into the soft tissue of flank of NOD/SCID mice. Thus, a different scaffold fabrication technique was chosen using the same polymer. Electrospun tubular scaffolds were seeded with human osteoblasts, as they showed previously the highest amount of bone formation and implanted into the flanks of NOD/SCID mice. Ectopic bone formation with sufficient vascularisation could be observed. After implantation of breast cancer cells using a polyethylene glycol hydrogel in close proximity to the newly formed bone, macroscopic communication between the newly formed bone and the tumour could be observed. Taken together, this PhD project showed that bone tissue engineering platforms could be used to develop an in vitro and in vivo model system to study cancer cell colonisation in the bone microenvironment.

Talk and web searching in an early years classroom

Recent Australian early childhood policy and curriculum guidelines promoting the use of technologies invite investigations of young children’s practices in classrooms. This study examined the practices of one preparatory year classroom, to show teacher and child interactions as they engaged in Web searching. The study investigated the in situ practices of the teacher and children to show how they accomplished the Web search. The data corpus consists of eight hours of videorecorded interactions over three days where children and teachers engaged in Web searching. One episode was selected that showed a teacher and two children undertaking a Web search. The episode is shown to consist of four phases: deciding on a new search subject, inputting the search query, considering the result options, and exploring the selected result. The sociological perspectives of ethnomethodology and conversation analysis were employed as the conceptual and methodological frameworks of the study, to analyse the video-recorded teacher and child interactions as they co-constructed a Web search. Ethnomethodology is concerned with how people make ‘sense’ in everyday interactions, and conversation analysis focuses on the sequential features of interaction to show how the interaction unfolds moment by moment. This extended single case analysis showed how the Web search was accomplished over multiple turns, and how the children and teacher collaboratively engaged in talk. There are four main findings. The first was that Web searching featured sustained teacher-child interaction, requiring a particular sort of classroom organisation to enable the teacher to work in this sustained way. The second finding was that the teacher’s actions recognised the children’s interactional competence in situ, orchestrating an interactional climate where everyone was heard. The third finding was that the teacher drew upon a range of interactional resources designed to progress the activity at hand, that of accomplishing the Web search. The teacher drew upon the interactional resources of interrogatives, discourse markers, and multi-unit turns during the Web search, and these assisted the teacher and children to co-construct their discussion, decide upon and co-ordinate their future actions, and accomplish the Web search in a timely way. The fourth finding explicates how particular social and pedagogic orders are accomplished through talk, where children collaborated with each other and with the teacher to complete the Web search. The study makes three key recommendations for the field of early childhood education. The study’s first recommendation is that fine-grained transcription and analysis of interaction aids in understanding interactional practices of Web searching. This study offers material for use in professional development, such as using transcribed and videorecorded interactions to highlight how teachers strategically engage with children, that is, how talk works in classroom settings. Another strategy is to focus on the social interactions of members engaging in Web searches, which is likely to be of interest to teachers as they work to engage with children in an increasingly online environment. The second recommendation involves classroom organisation; how teachers consider and plan for extended periods of time for Web searching, and how teachers accommodate children’s prior knowledge of Web searching in their classrooms. The third recommendation is in relation to future empirical research, with suggested possible topics focusing on the social interactions of children as they engage with peers as they Web search, as well as investigations of techno-literacy skills as children use the Internet in the early years.

A coupled SPH-DEM model for fluid and solid mechanics of apple parenchyma cells during drying

A coupled SPH-DEM based two-dimensional (2-D) micro-scale single cell model is developed to predict basic cell-level shrinkage effects of apple parenchyma cells during air drying. In this newly developed drying model, Smoothed Particle Hydrodynamics (SPH) is used to model the low Reynolds Number fluid motions of the cell protoplasm, and a Discrete Element Method (DEM) is employed to simulate the polymer-like cell wall. Simulations results reasonably agree with published experimental drying results on cellular shrinkage properties such as cellular area, diameter and perimeter. These preliminary results indicate that the model is effective for the modelling and simulation of apple parenchyma cells during air drying.

Dormant but migratory tumour cells in desmoplastic stroma of invasive ductal carcinomas

Mortality in breast cancer is linked to metastasis and recurrence yet there is no acceptable biological model for cancer relapse. We hypothesise that there might exist primary tumour cells capable of escaping surgery by migration and resisting radiotherapy and chemotherapy to cause cancer recurrence. We investigated this possibility in invasive ductal carcinoma (IDC) tissue and observed the presence of solitary primary tumour cells (SPCs) in the dense collagen stroma that encapsulates intratumoural cells (ICs). In IDC tissue sections, collagen was detected with either Masson's Trichrome or by second harmonics imaging. Cytokeratin-19 (CK-19) and vimentin (VIM) antibodies were, respectively, used to identify epithelial-derived tumour cells and to indicate epithelial to mesenchymal transition (EMT). Confocal/multiphoton microscopy showed that ICs from acini were mainly CK-19 +ve and were encapsulated by dense stromal collagen. Within the stroma, SPCs were detected by their staining for both CK-19 and VIM (confirming EMT). ICs and SPCs were subsequently isolated by laser capture microdissection followed by multiplex tandem-PCR studies. SPCs were found to be enriched for pro-migratory and anti-proliferative genes relative to ICs. In vitro experiments using collagen matrices at 20 mg/cm 3, similar in density to tumour matrices, demonstrated that SPC-like cells were highly migratory but dormant, phenotypes that recapitulated the genotypes of SPCs in clinical tissue. These data suggest that SPCs located at the breast cancer perimeter are invasive and dormant such that they may exceed surgical margins and resist local and adjuvant therapies. This study has important connotations for a role of SPCs in local recurrence.

Role of intratumoural heterogeneity in cancer drug resistance : molecular and clinical perspectives

Drug resistance continues to be a major barrier to the delivery of curative therapies in cancer. Historically, drug resistance has been associated with over-expression of drug transporters, changes in drug kinetics or amplification of drug targets. However, the emergence of resistance in patients treated with new-targeted therapies has provided new insight into the complexities underlying cancer drug resistance. Recent data now implicate intratumoural heterogeneity as a major driver of drug resistance. Single cell sequencing studies that identified multiple genetically distinct variants within human tumours clearly demonstrate the heterogeneous nature of human tumours. The major contributors to intratumoural heterogeneity are (i) genetic variation, (ii) stochastic processes, (iii) the microenvironment and (iv) cell and tissue plasticity. Each of these factors impacts on drug sensitivity. To deliver curative therapies to patients, modification of current therapeutic strategies to include methods that estimate intratumoural heterogeneity and plasticity will be essential.

DNA damage in response to an Ironman triathlon

The major aims of this study were to investigate the effect of an Ironman triathlon on DNA migration in the single cell gel electrophoresis assay, apoptosis and necrosis in the cytokinesis-block micronucleus cytome assay with lymphocytes and on changes of total antioxidant capacity in plasma. Blood samples were taken 2 days (d) before, within 20 min, 1 d, 5 d and 19 d post-race. The level of strand breaks decreased (p<0.05) immediately after the race, then increased (p<0.01) 1 d post-race and declined (p<0.01) until 19 d post-race. Apoptotic and necrotic cells decreased (p<0.01) and the total antioxidant status increased (p<0.01) immediately after the race. The results indicate that ultra-endurance exercise does not cause prolonged DNA damage in well-trained male athletes.

Estimate of Lagrangian integral scales in shallow tidal water using high resolution GPS-tracked drifters

In estuaries and natural water channels, the estimate of velocity and dispersion coefficients is critical to the knowledge of scalar transport and mixing. This estimate is rarely available experimentally at sub-tidal time scale in shallow water channels where high frequency is required to capture its spatio-temporal variation. This study estimates Lagrangian integral scales and autocorrelation curves, which are key parameters for obtaining velocity fluctuations and dispersion coefficients, and their spatio-temporal variability from deployments of Lagrangian drifters sampled at 10 Hz for a 4-hour period. The power spectral densities of the velocities between 0.0001 and 0.8 Hz were well fitted with a slope of 5/3 predicted by Kolmogorov’s similarity hypothesis within the inertial subrange, and were similar to the Eulerian power spectral previously observed within the estuary. The result showed that large velocity fluctuations determine the magnitude of the integral time scale, TL. Overlapping of short segments improved the stability of the estimate of TL by taking advantage of the redundant data included in the autocorrelation function. The integral time scales were about 20 s and varied by up to a factor of 8. These results are essential inputs for spatial binning of velocities, Lagrangian stochastic modelling and single particle analysis of the tidal estuary.

Exclusion of angiotensinogen gene in molecular basis of human hypertension: Sibpair linkage and association analyses in Australian anglo-caucasians

Linkage with essential hypertension has been claimed for a microsatellite marker near the angiotensinogen gene (AGT; chromosome 1q42), as has association for the AGT variants M235T, G(-6)A and A(-20)C. To more rigorously evaluate AGT as a candidate gene for hypertension we performed sibpair analysis with multiple microsatellite markers surrounding this locus and using more sophisticated analysis programs. We also performed an association study of the AGT variants in unrelated subjects with a strong family history (two affected parents). For the linkage study, single and multiplex polymerase chain reaction (PCRs) and automated genescan analysis were conducted on DNA from 175 Australian Anglo-Celtic Caucasian hypertensives for the following markers: D1S2880-(2.1 cM)-D1S213-(2.8 cM)-D1S251-(6.5 cM)-AGT-(2.0 cM) -D1S235. Statistical evaluation of genotype data by nonparametric methods resulted in the following scores: Single-point analysis - SPLINK, P > 0.18; APM method, P > 0.25; ASPEX, MLOD < 0.28; SIB-PAIR, P > 0. 24; Multipoint analysis - MAPMAKER/SIBS, MLOD < 0.24; GENEHUNTER, P > 0.35. Exclusion scores of Lod -4.1 to -5.1 were obtained for these markers using MAPMAKER/SIBS for a lambda(s) of 1.6. The association study of G(-6)A, A(-20)C and M235T variants in 111 hypertensives with strong family history and 190 normotensives with no family history showed significant linkage disequilibrium between particular haplotypes, but we could find no association with hypertension. The present study therefore excludes AGT in the etiology of hypertension, at least in the population of Australian Anglo-Celtic Caucasians studied.

Computation of confined swirling flow using RNG and k-epsilon turbulence models

In this work we numerically model isothermal turbulent swirling flow in a cylindrical burner. Three versions of the RNG k-epsilon model are assessed against performance of the standard k-epsilon model. Sensitivity of numerical predictions to grid refinement, differing convective differencing schemes and choice of (unknown) inlet dissipation rate, were closely scrutinised to ensure accuracy. Particular attention is paid to modelling the inlet conditions to within the range of uncertainty of the experimental data, as model predictions proved to be significantly sensitive to relatively small changes in upstream flow conditions. We also examine the characteristics of the swirl--induced recirculation zone predicted by the models over an extended range of inlet conditions. Our main findings are: - (i) the standard k-epsilon model performed best compared with experiment; - (ii) no one inlet specification can simultaneously optimize the performance of the models considered; - (iii) the RNG models predict both single-cell and double-cell IRZ characteristics, the latter both with and without additional internal stagnation points. The first finding indicates that the examined RNG modifications to the standard k-e model do not result in an improved eddy viscosity based model for the prediction of swirl flows. The second finding suggests that tuning established models for optimal performance in swirl flows a priori is not straightforward. The third finding indicates that the RNG based models exhibit a greater variety of structural behaviour, despite being of the same level of complexity as the standard k-e model. The plausibility of the predicted IRZ features are discussed in terms of known vortex breakdown phenomena.

Multiple human single-stranded DNA binding proteins function in genome maintenance : structural, biochemical and functional analysis

DNA exists predominantly in a duplex form that is preserved via specific base pairing. This base pairing affords a considerable degree of protection against chemical or physical damage and preserves coding potential. However, there are many situations, e.g. during DNA damage and programmed cellular processes such as DNA replication and transcription, in which the DNA duplex is separated into two singlestranded DNA (ssDNA) strands. This ssDNA is vulnerable to attack by nucleases, binding by inappropriate proteins and chemical attack. It is very important to control the generation of ssDNA and protect it when it forms, and for this reason all cellular organisms and many viruses encode a ssDNA binding protein (SSB). All known SSBs use an oligosaccharide/oligonucleotide binding (OB)-fold domain for DNA binding. SSBs have multiple roles in binding and sequestering ssDNA, detecting DNA damage, stimulating strand-exchange proteins and helicases, and mediation of protein–protein interactions. Recently two additional human SSBs have been identified that are more closely related to bacterial and archaeal SSBs. Prior to this it was believed that replication protein A, RPA, was the only human equivalent of bacterial SSB. RPA is thought to be required for most aspects of DNA metabolism including DNA replication, recombination and repair. This review will discuss in further detail the biological pathways in which human SSBs function.

Highly discriminatory single-nucleotide polymorphism interrogation of Escherichia coli by use of allele-specific real-time PCR and eBURST analysis

In total, 782 Escherichia coli strains originating from various host sources have been analyzed in this study by using a highly discriminatory single-nucleotide polymorphism (SNP) approach. A set of eight SNPs, with a discrimination value (Simpson's index of diversity [D]) of 0.96, was determined using the Minimum SNPs software, based on sequences of housekeeping genes from the E. coli multilocus sequence typing (MLST) database. Allele-specific real-time PCR was used to screen 114 E. coli isolates from various fecal sources in Southeast Queensland (SEQ). The combined analysis of both the MLST database and SEQ E. coli isolates using eight high-D SNPs resolved the isolates into 74 SNP profiles. The data obtained suggest that SNP typing is a promising approach for the discrimination of host-specific groups and allows for the identification of human-specific E. coli in environmental samples. However, a more diverse E. coli collection is required to determine animal- and environment-specific E. coli SNP profiles due to the abundance of human E. coli strains (56%) in the MLST database.

Single-agent versus combination chemotherapy in patients with advanced non-small cell lung cancer and a performance status of 2 : prognostic factors and treatment selection based on two large randomized clinical trials

Purpose: Data from two randomized phase III trials were analyzed to evaluate prognostic factors and treatment selection in the first-line management of advanced non-small cell lung cancer patients with performance status (PS) 2. Patients and Methods: Patients randomized to combination chemotherapy (carboplatin and paclitaxel) in one trial and single-agent therapy (gemcitabine or vinorelbine) in the second were included in these analyses. Both studies had identical eligibility criteria and were conducted simultaneously. Comparison of efficacy and safety was performed between the two cohorts. A regression analysis identified prognostic factors and subgroups of patients that may benefit from combination or single-agent therapy. Results: Two hundred one patients were treated with combination and 190 with single-agent therapy. Objective responses were 37 and 15%, respectively. Median time to progression was 4.6 months in the combination arm and 3.5 months in the single-agent arm (p < 0.001). Median survival imes were 8.0 and 6.6 months, and 1-year survival rates were 31 and 26%, respectively. Albumin <3.5 g, extrathoracic metastases, lactate dehydrogenase ≥200 IU, and 2 comorbid conditions predicted outcome. Patients with 0-2 risk factors had similar outcomes independent of treatment, whereas patients with 3-4 factors had a nonsignificant improvement in median survival with combination chemotherapy. Conclusion: Our results show that PS2 non-small cell lung cancer patients are a heterogeneous group who have significantly different outcomes. Patients treated with first-line combination chemotherapy had a higher response and longer time to progression, whereas overall survival did not appear significantly different. A prognostic model may be helpful in selecting PS 2 patients for either treatment strategy. © 2009 by the International Association for the Study of Lung Cancer.

Porohyperelastic analysis of single osteocyte using AFM and inverse FEA

Osteocytes are the mature cells and perform as mechanosensors within the bone. The mechanical property of osteocytes plays an important role to fulfill these functions. However, little researches have been done to investigate the mechanical deformation properties of single osteocytes. Atomic Force Microscopy (AFM) is a state-of-art experimental facility for high resolution imaging of tissues, cells and any surfaces as well as for probing mechanical properties of the samples both qualitatively and quantitatively. In this paper, the experimental study based on AFM is firstly used to obtain forceindentation curves of single round osteocytes. The porohyperelastic (PHE) model of a single osteocyte is then developed by using the inverse finite element analysis (FEA) to identify and extract mechanical properties from the experiment results. It has been found that the PHE model is a good candidature for biomechanics studies of osteocytes.

Porohyperelastic analysis of single chondrocyte using AFM and inverse FEA

The aim of this paper is to determine the creep and relaxation responses of single chondrocytes in vitro. Firstly, Atomic Force Microscopy (AFM) was used to obtain the force-indentation curves of single chondrocytes at the strain-rate of 7.05 s-1. This result was then employed in inverse finite element analysis (FEA) using porohyperelastic (PHE) idealization of the cells to determine their mechanical properties. The PHE model results agreed well with AFM experimental data. This PHE model was then utilized to study chondrocyte’s creep and relaxation behaviors. The results revealed that the effect of fluid was predominant for cell’s mechanical behaviors and that the PHE is a good model for biomechanics studies of chondrocytes.