Incorporation of exogenous RGD peptide and inter-species blending as strategies for enhancing human corneal limbal epithelial cell growth on Bombyx mori silk fibroin membranes

While fibroin isolated from the cocoons of domesticated silkworm Bombyx mori supports growth of human corneal limbal epithelial (HLE) cells, the mechanism of cell attachment remains unclear. In the present study we sought to enhance the attachment of HLE cells to membranes of Bombyx mori silk fibroin (BMSF) through surface functionalization with an arginine-glycine-aspartic acid (RGD)-containing peptide. Moreover, we have examined the response of HLE cells to BMSF when blended with the fibroin produced by a wild silkworm, Antheraea pernyi, which is known to contain RGD sequences within its primary structure. A procedure to isolate A. pernyi silk fibroin (APSF) from the cocoons was established, and blends of the two fibroins were prepared at five different BMSF/APSF ratios. In another experiment, BMSF surface was modified by binding chemically the GRGDSPC peptide using a water-soluble carbodiimide. Primary HLE were grown in the absence of serum on membranes made of BMSF, APSF, and their blends, as well as on RGD-modified BMSF. There was no statistically significant enhancing effect on the cell attachment due to the RGD presence. This suggests that the adhesion through RGD ligands may have a complex mechanism, and the investigated strategies are of limited value unless the factors contributing to this mechanism become better known.

Numerical solution of an optimal control model of dendritic cell treatment of a growing tumour

A new optimal control model of the interactions between a growing tumour and the host immune system along with an immunotherapy treatment strategy is presented. The model is based on an ordinary differential equation model of interactions between the growing tu- mour and the natural killer, cytotoxic T lymphocyte and dendritic cells of the host immune system, extended through the addition of a control function representing the application of a dendritic cell treat- ment to the system. The numerical solution of this model, obtained from a multi species Runge–Kutta forward-backward sweep scheme, is described. We investigate the effects of varying the maximum al- lowed amount of dendritic cell vaccine administered to the system and find that control of the tumour cell population is best effected via a high initial vaccine level, followed by reduced treatment and finally cessation of treatment. We also found that increasing the strength of the dendritic cell vaccine causes an increase in the number of natural killer cells and lymphocytes, which in turn reduces the growth of the tumour.

Humanised xenograft models of bone metastasis revisited : novel insights into species-specific mechanisms of cancer cell osteotropism

The determinants and key mechanisms of cancer cell osteotropism have not been identified, mainly due to the lack of reproducible animal models representing the biological, genetic and clinical features seen in humans. An ideal model should be capable of recapitulating as many steps of the metastatic cascade as possible, thus facilitating the development of prognostic markers and novel therapeutic strategies. Most animal models of bone metastasis still have to be derived experimentally as most syngeneic and transgeneic approaches do not provide a robust skeletal phenotype and do not recapitulate the biological processes seen in humans. The xenotransplantation of human cancer cells or tumour tissue into immunocompromised murine hosts provides the possibility to simulate early and late stages of the human disease. Human bone or tissue-engineered human bone constructs can be implanted into the animal to recapitulate more subtle, species-specific aspects of the mutual interaction between human cancer cells and the human bone microenvironment. Moreover, the replication of the entire "organ" bone makes it possible to analyse the interaction between cancer cells and the haematopoietic niche and to confer at least a partial human immunity to the murine host. This process of humanisation is facilitated by novel immunocompromised mouse strains that allow a high engraftment rate of human cells or tissue. These humanised xenograft models provide an important research tool to study human biological processes of bone metastasis.

Dosimetric effects of a high-density spinal implant

In this study, a treatment plan for a spinal lesion, with all beams transmitted though a titanium vertebral reconstruction implant, was used to investigate the potential effect of a high-density implant on a three-dimensional dose distribution for a radiotherapy treatment. The BEAMnrc/DOSXYZnrc and MCDTK Monte Carlo codes were used to simulate the treatment using both a simplified, recltilinear model and a detailed model incorporating the full complexity of the patient anatomy and treatment plan. The resulting Monte Carlo dose distributions showed that the commercial treatment planning system failed to accurately predict both the depletion of dose downstream of the implant and the increase in scattered dose adjacent to the implant. Overall, the dosimetric effect of the implant was underestimated by the commercial treatment planning system and overestimated by the simplified Monte Carlo model. The value of performing detailed Monte Carlo calculations, using the full patient and treatment geometry, was demonstrated.

Cell migration and proliferation on homogeneous and non-homogeneous domains : modelling on the scale of individuals and populations

Cell migration is a behaviour critical to many key biological effects, including wound healing, cancerous cell invasion and morphogenesis, the development of an organism from an embryo. However, given that each of these situations is distinctly different and cells are extremely complicated biological objects, interest lies in more basic experiments which seek to remove conflating factors and present a less complex environment within which cell migration can be experimentally examined. These include in vitro studies like the scratch assay or circle migration assay, and ex vivo studies like the colonisation of the hindgut by neural crest cells. The reduced complexity of these experiments also makes them much more enticing as problems to mathematically model, like done here. The primary goal of the mathematical models used in this thesis is to shed light on which cellular behaviours work to generate the travelling waves of invasion observed in these experiments, and to explore how variations in these behaviours can potentially predict differences in this invasive pattern which are experimentally observed when cell types or chemical environment are changed. Relevant literature has already identified the difficulty of distinguishing between these behaviours when using traditional mathematical biology techniques operating on a macroscopic scale, and so here a sophisticated individual-cell-level model, an extension of the Cellular Potts Model (CPM), is been constructed and used to model a scratch assay experiment. This model includes a novel mechanism for dealing with cell proliferations that allowed for the differing properties of quiescent and proliferative cells to be implemented into their behaviour. This model is considered both for its predictive power and used to make comparisons with the travelling waves which result in more traditional macroscopic simulations. These comparisons demonstrate a surprising amount of agreement between the two modelling frameworks, and suggest further novel modifications to the CPM that would allow it to better model cell migration. Considerations of the model’s behaviour are used to argue that the dominant effect governing cell migration (random motility or signal-driven taxis) likely depends on the sort of invasion demonstrated by cells, as easily seen by microscopic photography. Additionally, a scratch assay simulated on a non-homogeneous domain consisting of a ’fast’ and ’slow’ region is also used to further differentiate between these different potential cell motility behaviours. A heterogeneous domain is a novel situation which has not been considered mathematically in this context, nor has it been constructed experimentally to the best of the candidate’s knowledge. Thus this problem serves as a thought experiment used to test the conclusions arising from the simulations on homogeneous domains, and to suggest what might be observed should this non-homogeneous assay situation be experimentally realised. Non-intuitive cell invasion patterns are predicted for diffusely-invading cells which respond to a cell-consumed signal or nutrient, contrasted with rather expected behaviour in the case of random-motility-driven invasion. The potential experimental observation of these behaviours is demonstrated by the individual-cell-level model used in this thesis, which does agree with the PDE model in predicting these unexpected invasion patterns. In the interest of examining such a case of a non-homogeneous domain experimentally, some brief suggestion is made as to how this could be achieved.

Adsorption and dissociation of ammonia borane outside and inside single-walled carbon nanotubes : A density functional theory study

Amonia borane (AB) has been identified as a potential candidate highcapacity hydrogen storage material. This work probes the adsorption and dissociation of AB inside and outside single-walled carbon nanotubes (SWCNTs) within the framework of density functional theory. The dissociation barriers of AB have been calculated and compared with that of the isolated AB molecule. On the basis of the present calculations, no notable improvement results from SWCNT confinement; on the contrary, the dissociation barrier slightly increases with respect to isolated AB.

A density functional theory study on CO2 capture and activation by graphene-like boron nitride with boron vacancy

First principle calculations for a hexagonal (graphene-like) boron nitride (g-BN) monolayer sheet in the presence of a boron-atom vacancy show promising properties for capture and activation of carbon dioxide. CO2 is found to decompose to produce an oxygen molecule via an intermediate chemisorption state on the defect g-BN sheet. The three stationary states and two transition states in the reaction pathway are confirmed by minimum energy pathway search and frequency analysis. The values computed for the two energy barriers involved in this catalytic reaction after enthalpy correction indicate that the catalytic reaction should proceed readily at room temperature.

A density functional theory study of CO2 and N2 adsorption on aluminium nitride single walled nanotubes

Strong binding of isolated carbon dioxide (CO2) on aluminium nitride (AlN) single walled nanotubes is verified using two different functionals. Two optimized configurations corresponding to physisorption and chemisorption are linked by a low energy barrier, such that the chemisorbed state is accessible and thermodynamically favored at low temperatures. In contrast, N2 is found only to form a physisorbed complex with the AlN nanotube, suggesting the potential application of aluminium nitride based materials for CO2 fixation. The effect of nanotube diameter on gas adsorption properties is also discussed. The diameter is found to have an important effect on the chemisorption of CO2, but has little effect on the physisorption of either CO2 or N2.

Adsorption of carbon dioxide and nitrogen on single-layer aluminum nitride nanostructures studied by density functional theory

The adsorption of carbon dioxide and nitrogen molecules on aluminum nitride (AlN) nanostructures has been explored using first-principle computational methods. Optimized configurations corresponding to physisorption and, subsequentially, chemisorption of CO2 are identified, in contrast to N2, for which only a physisorption structure is found. Transition-state searches imply a low energy barrier between the physisorption and chemisorption states for CO2 such that the latter is accessible and thermodynamically favored at room temperature. The effective binding energy of the optimized chemisorption structure is apparently larger than those for other CO2 adsorptive materials, suggesting the potential for application of aluminum nitride nanostructures for carbon dioxide capture and storage.

The effect of Fe doping on adsorption of CO2/N2within carbon nanotubes : a density functional theory study with dispersion corrections

An ab initio density functional theory (DFT) study with correction for dispersive interactions was performed to study the adsorption of N2 and CO2 inside an (8, 8) single-walled carbon nanotube. We find that the approach of combining DFT and van der Waals correction is very effective for describing the long-range interaction between N2/CO2 and the carbon nanotube (CNT). Surprisingly, exohedral doping of an Fe atom onto the CNT surface will only affect the adsorption energy of the quadrupolar CO2 molecule inside the CNT (20–30%), and not that of molecular N2. Our results suggest the feasibility of enhancement of CO2/N2 separation in CNT-based membranes by using exohedral doping of metal atoms.

The catalytic role of an isolated-Ti atom in the hydrogenation of Ti-doped Al(001) surface : An ab initio density functional theory calculation

Recent work [S. Chaudhuri, J.T. Muckerman, J. Phys. Chem. B 109 (2005) 6952] reported that two Ti-substituted atoms on an Al(0 0 1) surface can form a catalytically active site for the dissociation of H2, but the diffusion barrier of atomic H away from Ti site is as high as 1.57 eV. By using ab initio density functional calculations, we found that two hydrogen molecules can dissociate on isolated-Ti atom doped Al(0 0 1) surface with small activation barriers (0.21 and 0.235 eV for first and second H2, respectively). Additionally, the diffusion barrier of atomic H away from Ti site is also moderate (0.47 eV). These results contribute further towards understanding the improved kinetics observed in recycling of hydrogen with Ti-doped NaAlH4.

Role of charge in destabilizing AlH4 and BH4 complex anions for hydrogen storage applications : Ab initio density functional calculations

NaAlH4 and LiBH4 are potential candidate materials for mobile hydrogen storage applications, yet they have the drawback of being highly stable and desorbing hydrogen only at elevated temperatures. In this letter, ab initio density functional theory calculations reveal how the stabilities of the AlH4 and BH4 complex anions will be affected by reducing net anionic charge. Tetrahedral AlH4 and BH4 complexes are found to be distorted with the decrease of negative charge. One H-H distance becomes smaller and the charge density will overlap between them at a small anion charge. The activation energies to release of H2 from AlH4 and BH4 complexes are thus greatly decreased. We demonstrate that point defects such as neutral Na vacancies or substitution of a Na atom with Ti on the NaAlH4(001) surface can potentially cause strong distortion of neighboring AlH4 complexes and even induce spontaneous dehydrogenation. Our results help to rationalize the conjecture that the suppression of charge transfer to AlH4 and BH4 anion as a consequence of surface defects should be very effective for improving the recycling performance of H2 in NaAlH4 and LiBH4. The understanding gained here will aid in the rational design and development of hydrogen storage materials based on these two systems.

Van der Waals-corrected density functional theory : benchmarking for hydrogen–nanotube and nanotube–nanotube interactions

Density functional theory (DFT) is a powerful approach to electronic structure calculations in extended systems, but suffers currently from inadequate incorporation of long-range dispersion, or Van der Waals (VdW) interactions. VdW-corrected DFT is tested for interactions involving molecular hydrogen, graphite, single-walled carbon nanotubes (SWCNTs), and SWCNT bundles. The energy correction, based on an empirical London dispersion term with a damping function at short range, allows a reasonable physisorption energy and equilibrium distance to be obtained for H2 on a model graphite surface. The VdW-corrected DFT calculation for an (8, 8) nanotube bundle reproduces accurately the experimental lattice constant. For H2 inside or outside an (8, 8) SWCNT, we find the binding energies are respectively higher and lower than that on a graphite surface, correctly predicting the well known curvature effect. We conclude that the VdW correction is a very effective method for implementing DFT calculations, allowing a reliable description of both short-range chemical bonding and long-range dispersive interactions. The method will find powerful applications in areas of SWCNT research where empirical potential functions either have not been developed, or do not capture the necessary range of both dispersion and bonding interactions.

Adatom CCV Auger rates via the local density of states

In this letter the core-core-valence Auger transitions of an atomic impurity, both in bulk or adsorbed on a jellium-like surface, are computed within a DFT framework. The Auger rates calculated by the Fermi golden rule are compared with those determined by an approximate and simpler expression. This is based on the local density of states (LDOS) with a core hole present, in a region around the impurity nucleus. Different atoms, Na and Mg, solids, Al and Ag, and several impurity locations are considered. We obtain an excellent agreement between KL1V and KL23V rates worked out with the two approaches. The radius of the sphere in which we calculate the LDOS is the relevant parameter of the simpler approach. Its value only depends on the atomic species regardless of the location of the impurity and the type of substrate. (C) 2003 Elsevier B.V. All rights reserved.

Effect of CT electron density identification on Monte Carlo simulations of radiotherapeutic treatments

Introduction: The accurate identification of tissue electron densities is of great importance for Monte Carlo (MC) dose calculations. When converting patient CT data into a voxelised format suitable for MC simulations, however, it is common to simplify the assignment of electron densities so that the complex tissues existing in the human body are categorized into a few basic types. This study examines the effects that the assignment of tissue types and the calculation of densities can have on the results of MC simulations, for the particular case of a Siemen’s Sensation 4 CT scanner located in a radiotherapy centre where QA measurements are routinely made using 11 tissue types (plus air). Methods: DOSXYZnrc phantoms are generated from CT data, using the CTCREATE user code, with the relationship between Hounsfield units (HU) and density determined via linear interpolation between a series of specified points on the ‘CT-density ramp’ (see Figure 1(a)). Tissue types are assigned according to HU ranges. Each voxel in the DOSXYZnrc phantom therefore has an electron density (electrons/cm3) defined by the product of the mass density (from the HU conversion) and the intrinsic electron density (electrons /gram) (from the material assignment), in that voxel. In this study, we consider the problems of density conversion and material identification separately: the CT-density ramp is simplified by decreasing the number of points which define it from 12 down to 8, 3 and 2; and the material-type-assignment is varied by defining the materials which comprise our test phantom (a Supertech head) as two tissues and bone, two plastics and bone, water only and (as an extreme case) lead only. The effect of these parameters on radiological thickness maps derived from simulated portal images is investigated. Results & Discussion: Increasing the degree of simplification of the CT-density ramp results in an increasing effect on the resulting radiological thickness calculated for the Supertech head phantom. For instance, defining the CT-density ramp using 8 points, instead of 12, results in a maximum radiological thickness change of 0.2 cm, whereas defining the CT-density ramp using only 2 points results in a maximum radiological thickness change of 11.2 cm. Changing the definition of the materials comprising the phantom between water and plastic and tissue results in millimetre-scale changes to the resulting radiological thickness. When the entire phantom is defined as lead, this alteration changes the calculated radiological thickness by a maximum of 9.7 cm. Evidently, the simplification of the CT-density ramp has a greater effect on the resulting radiological thickness map than does the alteration of the assignment of tissue types. Conclusions: It is possible to alter the definitions of the tissue types comprising the phantom (or patient) without substantially altering the results of simulated portal images. However, these images are very sensitive to the accurate identification of the HU-density relationship. When converting data from a patient’s CT into a MC simulation phantom, therefore, all possible care should be taken to accurately reproduce the conversion between HU and mass density, for the specific CT scanner used. Acknowledgements: This work is funded by the NHMRC, through a project grant, and supported by the Queensland University of Technology (QUT) and the Royal Brisbane and Women's Hospital (RBWH), Brisbane, Australia. The authors are grateful to the staff of the RBWH, especially Darren Cassidy, for assistance in obtaining the phantom CT data used in this study. The authors also wish to thank Cathy Hargrave, of QUT, for assistance in formatting the CT data, using the Pinnacle TPS. Computational resources and services used in this work were provided by the HPC and Research Support Group, QUT, Brisbane, Australia.