Modelling of some biological materials using continuum mechanics

Continuum mechanics provides a mathematical framework for modelling the physical stresses experienced by a material. Recent studies show that physical stresses play an important role in a wide variety of biological processes, including dermal wound healing, soft tissue growth and morphogenesis. Thus, continuum mechanics is a useful mathematical tool for modelling a range of biological phenomena. Unfortunately, classical continuum mechanics is of limited use in biomechanical problems. As cells refashion the �bres that make up a soft tissue, they sometimes alter the tissue's fundamental mechanical structure. Advanced mathematical techniques are needed in order to accurately describe this sort of biological `plasticity'. A number of such techniques have been proposed by previous researchers. However, models that incorporate biological plasticity tend to be very complicated. Furthermore, these models are often di�cult to apply and/or interpret, making them of limited practical use. One alternative approach is to ignore biological plasticity and use classical continuum mechanics. For example, most mechanochemical models of dermal wound healing assume that the skin behaves as a linear viscoelastic solid. Our analysis indicates that this assumption leads to physically unrealistic results. In this thesis we present a novel and practical approach to modelling biological plasticity. Our principal aim is to combine the simplicity of classical linear models with the sophistication of plasticity theory. To achieve this, we perform a careful mathematical analysis of the concept of a `zero stress state'. This leads us to a formal de�nition of strain that is appropriate for materials that undergo internal remodelling. Next, we consider the evolution of the zero stress state over time. We develop a novel theory of `morphoelasticity' that can be used to describe how the zero stress state changes in response to growth and remodelling. Importantly, our work yields an intuitive and internally consistent way of modelling anisotropic growth. Furthermore, we are able to use our theory of morphoelasticity to develop evolution equations for elastic strain. We also present some applications of our theory. For example, we show that morphoelasticity can be used to obtain a constitutive law for a Maxwell viscoelastic uid that is valid at large deformation gradients. Similarly, we analyse a morphoelastic model of the stress-dependent growth of a tumour spheroid. This work leads to the prediction that a tumour spheroid will always be in a state of radial compression and circumferential tension. Finally, we conclude by presenting a novel mechanochemical model of dermal wound healing that takes into account the plasticity of the healing skin.

Molecular epidemiology of respiratory viruses in a paediatric cohort from Indonesia

Acute lower respiratory tract infections (ALRTIs) are a common cause of morbidity and mortality among children under 5 years of age and are found worldwide, with pneumonia as the most severe manifestation. Although the incidence of severe disease varies both between individuals and countries, there is still no clear understanding of what causes this variation. Studies of community-acquired pneumonia (CAP) have traditionally not focused on viral causes of disease due to a paucity of diagnostic tools. However, with the emergence of molecular techniques, it is now known that viruses outnumber bacteria as the etiological agents of childhood CAP, especially in children under 2 years of age. The main objective of this study was to investigate viruses contributing to disease severity in cases of childhood ALRTI, using a two year cohort study following 2014 infants and children enrolled in Bandung, Indonesia. A total of 352 nasopharyngeal washes collected from 256 paediatric ALRTI patients were used for analysis. A subset of samples was screened using a novel microarray pathogen detection method that identified respiratory syncytial virus (RSV), human metapneumovirus (hMPV) and human rhinovirus (HRV) in the samples. Real-time RT-PCR was used both for confirming and quantifying viruses found in the nasopharyngeal samples. Viral copy numbers were determined and normalised to the numbers of human cells collected with the use of 18S rRNA. Molecular epidemiology was performed for RSV A and hMPV using sequences to the glycoprotein gene and nucleoprotein gene respectively, to determine genotypes circulating in this Indonesian paediatric cohort. This study found that HRV (119/352; 33.8%) was the most common virus detected as the cause of respiratory tract infections in this cohort, followed by the viral pathogens RSV A (73/352; 20.7%), hMPV (30/352; 8.5%) and RSV B (12/352; 3.4%). Co-infections of more than two viruses were detected in 31 episodes (defined as an infection which occurred more than two weeks apart), accounting for 8.8% of the 352 samples tested or 15.4% of the 201 episodes with at least one virus detected. RSV A genotypes circulating in this population were predominantly GA2, GA5 and GA7, while hMPV genotypes circulating were mainly A2a (27/30; 90.0%), B2 (2/30; 6.7%) and A1 (1/30; 3.3%). This study found no evidence of disease severity associated either with a specific virus or viral strain, or with viral load. However, this study did find a significant association with co-infection of RSV A and HRV with severe disease (P = 0.006), suggesting that this may be a novel cause of severe disease.

Molecular characterisation of six badnavirus species associated with leaf streak disease of banana in East Africa

Banana leaf streak disease, caused by several species of Banana streak virus (BSV), is widespread in East Africa. We surveyed for this disease in Uganda and Kenya, and used rolling-circle amplification (RCA) to detect the presence of BSV in banana. Six distinct badnavirus sequences, three from Uganda and three from Kenya, were amplified for which only partial sequences were previously available. The complete genomes were sequenced and characterised. The size and organisation of all six sequences was characteristic of other badnaviruses, including conserved functional domains present in the putative polyprotein encoded by open reading frame (ORF) 3. Based on nucleotide sequence analysis within the reverse transcriptase/ribonuclease H-coding region of open reading frame 3, we propose that these sequences be recognised as six new species and be designated as Banana streak UA virus, Banana streak UI virus, Banana streak UL virus, Banana streak UM virus, Banana streak CA virus and Banana streak IM virus. Using PCR and species-specific primers to test for the presence of integrated sequences, we demonstrated that sequences with high similarity to BSIMV only were present in several banana cultivars which had tested negative for episomal BSV sequences.

Molecular analysis of the myosin gene family in Arabidopsis thaliana

Myosin is believed to act as the molecular motor for many actin-based motility processes in eukaryotes. It is becoming apparent that a single species may possess multiple myosin isoforms, and at least seven distinct classes of myosin have been identified from studies of animals, fungi, and protozoans. The complexity of the myosin heavy-chain gene family in higher plants was investigated by isolating and characterizing myosin genomic and cDNA clones from Arabidopsis thaliana. Six myosin-like genes were identified from three polymerase chain reaction (PCR) products (PCR1, PCR11, PCR43) and three cDNA clones (ATM2, MYA2, MYA3). Sequence comparisons of the deduced head domains suggest that these myosins are members of two major classes. Analysis of the overall structure of the ATM2 and MYA2 myosins shows that they are similar to the previously-identified ATM1 and MYA1 myosins, respectively. The MYA3 appears to possess a novel tail domain, with five IQ repeats, a six-member imperfect repeat, and a segment of unique sequence. Northern blot analyses indicate that some of the Arabidopsis myosin genes are preferentially expressed in different plant organs. Combined with previous studies, these results show that the Arabidopsis genome contains at least eight myosin-like genes representing two distinct classes.

An advanced implicit meshless approach for fractional partial differential equation in computational mechanics

Recently, the numerical modelling and simulation for fractional partial differential equations (FPDE), which have been found with widely applications in modern engineering and sciences, are attracting increased attentions. The current dominant numerical method for modelling of FPDE is the explicit Finite Difference Method (FDM), which is based on a pre-defined grid leading to inherited issues or shortcomings. This paper aims to develop an implicit meshless approach based on the radial basis functions (RBF) for numerical simulation of time fractional diffusion equations. The discrete system of equations is obtained by using the RBF meshless shape functions and the strong-forms. The stability and convergence of this meshless approach are then discussed and theoretically proven. Several numerical examples with different problem domains are used to validate and investigate accuracy and efficiency of the newly developed meshless formulation. The results obtained by the meshless formations are also compared with those obtained by FDM in terms of their accuracy and efficiency. It is concluded that the present meshless formulation is very effective for the modelling and simulation for FPDE.

Atomistic numerical investigation of single-crystal copper nanowire with surface defects

Based on the embedded atom method (EAM), a molecular dynamics (MD) simulation is performed to study the single-crystal copper nanowire with surface defects through tension. The tension simulations for nanowire without defect are first carried out under different temperatures, strain rates and time steps and then surface defect effects for nanowire are investigated. The stress-strain curves obtained by the MD simulations of various strain rates show a rate below 1 x 10(9) s-1 will exert less effect on the yield strength and yield point, and the Young's modulus is independent of strain rate. a time step below 5 fs is recommend for the atomic model during the MD simulation. It is observed that high temperature leads to low Young's modulus, as well as the yield strength. The surface defects on nanowires are systematically studied in considering different defect orientations. It is found that the surface defect serves as a dislocation source, and the yield strength shows 34.20% decresse with 45 degree surface defect. Both yield strength and yield point are significantly influenced by the surface defects, except the Young's modulus.

Molecular evolutionary dynamics of Ross River virus and implications for vaccine efficacy

Ross River virus (RRV) is a mosquito-borne member of the genus Alphavirus that causes epidemic polyarthritis in humans, costing the Australian health system at least US$10 million annually. Recent progress in RRV vaccine development requires accurate assessment of RRV genetic diversity and evolution, particularly as they may affect the utility of future vaccination. In this study, we provide novel RRV genome sequences and investigate the evolutionary dynamics of RRV from time-structured E2 gene datasets. Our analysis indicates that, although RRV evolves at a similar rate to other alphaviruses (mean evolutionary rate of approx. 8x10(-4) nucleotide substitutions per site year(-1)), the relative genetic diversity of RRV has been continuously low through time, possibly as a result of purifying selection imposed by replication in a wide range of natural host and vector species. Together, these findings suggest that vaccination against RRV is unlikely to result in the rapid antigenic evolution that could compromise the future efficacy of current RRV vaccines.

Molecular structure of Mg-Al, Mn-Al and Zn-Al halotrichites : an infrared spectroscopic study

Near infrared (NIR), X-ray diffraction (XRD) and infrared (IR) spectroscopy have been applied to halotrichites of the formula MgAl2(SO4)4∙22H2O, MnAl2(SO4)4∙22H2O and ZnAl2(SO4)4∙22H2O. Comparison of the halotrichites in different spectral regions has shown that the incorporation of a divalent transition metal into the halotrichite structure causes a shift in OH stretching band positions to lower wavenumbers. Therefore, an increase in hydrogen bonded water is observed for divalent cations with a larger molecular mass. XRD has confirmed the formation of halotrichite for all three samples and characteristic peaks of halotrichite have been identified at 18.5 and 24.5° 2θ, along with a group of six peaks between 5 and 15° 2θ. It has been observed that Mg-Al and Mn-Al halotrichite are very similar in structure, while Zn-Al showed several differences particularly in the NIR spectra. This work has shown that halotrichite structures can be synthesised and characterised by infrared and NIR spectroscopy.

The molecular mechanisms utilised by mesenchymal stem cells in migration to acute myocardial infarction

Heart damage caused by acute myocardial infarction (AMI) is a leading cause of death and disability in Australia. Novel therapies are still required for the treatment of this condition due to the poor reparative ability of the heart. As such, cellular therapies that assist in the recovery of heart muscle are of great current interest. Culture expanded mesenchymal stem cells (MSC) represent a stem and progenitor cell population that has been shown to promote tissue recovery in pre-clinical studies of AMI. For MSC-based therapies in the clinic, an intravenous route of administration would ideally be used due to the low cost, ease of delivery and relative safety. The study of MSC migration is therefore clinically relevant for a minimally invasive cell therapy to promote regeneration of damaged tissue. C57BL/6, UBI-GFP-BL/6 and CD44-/-/GFP+/+ mice were utilised to investigate mMSC migration. To assist in murine models of MSC migration, a novel method was used for the isolation of murine MSC (mMSC). These mMSC were then expanded in culture and putative mMSC were positive for Sca-1, CD90.2, and CD44 and were negative for CD45 and CD11b. Furthermore, mMSC from C57BL/6 and UBI-GFP-BL/6 mice were shown to differentiate into cells of the mesodermal lineage. Cells from CD44-/-/GFP+/+ mice were positive for Sca-1 and CD90.2, and negative for CD44, CD45 and CD11b however, these cells were unable to differentiate into adipocytes and chondrocytes and express lineage specific genes, PLIN and ACAN. Analysis of mMSC chemokine receptor (CR) expression showed that although mMSC do express chemokine receptors, (including those specific for chemokines released after AMI), these were low or undetectable by mRNA. However, protein expression could be detected, which was predominantly cytoplasmic. It was further shown that in both healthy (unperturbed) and inflamed tissues, mMSC had very little specific migration and engraftment after intravenous injection. To determine if poor mMSC migration was due to the inability of mMSC to respond to chemotactic stimuli, chemokine expression in bone marrow, skin injury and hearts (healthy and after AMI) was analysed at various time points by quantitative real-time PCR (qRT PCR). Many chemokines were up-regulated after skin biopsy and AMI, but the highest acute levels were found for CXCL12 and CCL7. Due to their high expression in infarcted hearts, the chemokines CXCL12 and CCL7 were tested for their effect on mMSC migration. Despite CR expression at both protein and mRNA levels, migration in response to CXCL12 and CCL7 was low in mMSC cultured on Nunclon plastic. A novel tissue culture plastic technology (UpCellTM) was then used that allowed gentle non-enzymatic dissociation of mMSC, thus preserving surface expression of the CRs. Despite this the in vitro data indicated that CXCL12 fails to induce significant migration ability of mMSC, while CCL7 induces significant, but low-level migration. We speculated this may be because of low levels of surface expression of chemokine receptors. In a strategy to increase cell surface expression of mMSC chemokine receptors and enhance their in vitro and in vivo migration capacity, mMSC were pre-treated with pro-inflammatory cytokines. Increased levels of both mRNA and surface protein expression were found for CRs by pre-treating mMSC with pro-inflammatory cytokines including TNF-á, IFN-ã, IL-1á and IL-6. Furthermore, the chemotactic response of mMSC to CXCL12 and CCL7 was significantly higher with these pretreated cells. Finally, the effectiveness of this type of cell manipulation was demonstrated in vivo, where mMSC pre-treated with TNF-á and IFN-ã showed significantly increased migration in skin injury and AMI models. Therefore this thesis has demonstrated, using in vitro and in vivo models, the potential for prior manipulation of MSC as a possible means for increasing the utility of intravenously delivery for MSC-based cellular therapies.

Molecular dynamics study of dynamic buckling properties of nanowires with defects

Based on the embedded atom method (EAM) and molecular dynamics (MD) method, the deformation properties of Cu nanowires with different single defects under dynamic compression have been studied. The mechanical behaviours of the perfect nanowire are first studied, and the critical stress decreases with the increase of the nanowire’s length, which is well agreed with the modified Euler theory. We then consider the effects to the buckling phenomenon resulted from different defects. It is found that obvious decrease of the critical stress is resulted from different defects, and the largest decrease is found in nanowire with the surface vertical defect. Surface defects are found exerting larger influence than internal defects. The buckling duration is found shortened due to different defects except the nanowire with surface horizon defect, which is also found possessing the largest deflection. Different deflections are also observed for different defected nanowires. It is find that due to surface defects, only deflection in one direction is happened, but for internal defects, more complex deflection circumstances are observed.