Growth of confined cancer spheroids : a combined experimental marker, critical modelling approach

A critical step in the dissemination of ovarian cancer is the formation of multicellular spheroids from cells shed from the primary tumour. The objectives of this study were to apply bioengineered three-dimensional (3D) microenvironments for culturing ovarian cancer spheroids in vitro and simultaneously to build on a mathematical model describing the growth of multicellular spheroids in these biomimetic matrices. Cancer cells derived from human epithelial ovarian carcinoma were embedded within biomimetic hydrogels of varying stiffness and grown for up to 4 weeks. Immunohistochemistry, imaging and growth analyses were used to quantify the dependence of cell proliferation and apoptosis on matrix stiffness, long-term culture and treatment with the anti-cancer drug paclitaxel. The mathematical model was formulated as a free boundary problem in which each spheroid was treated as an incompressible porous medium. The functional forms used to describe the rates of cell proliferation and apoptosis were motivated by the experimental work and predictions of the mathematical model compared with the experimental output. This work aimed to establish whether it is possible to simulate solid tumour growth on the basis of data on spheroid size, cell proliferation and cell death within these spheroids. The mathematical model predictions were in agreement with the experimental data set and simulated how the growth of cancer spheroids was influenced by mechanical and biochemical stimuli including matrix stiffness, culture duration and administration of a chemotherapeutic drug. Our computational model provides new perspectives on experimental results and has informed the design of new 3D studies of chemoresistance of multicellular cancer spheroids.

Porohyperelastic analysis to explore mechanical properties of chondrocytes using numerical modeling and experiments : a finite element study

There are many continuum mechanical models have been developed such as liquid drop models, solid models, and so on for single living cell biomechanics studies. However, these models do not give a fully approach to exhibit a clear understanding of the behaviour of single living cells such as swelling behaviour, drag effect, etc. Hence, the porohyperelastic (PHE) model which can capture those aspects would be a good candidature to study cells behaviour (e.g. chondrocytes in this study). In this research, an FEM model of single chondrocyte cell will be developed by using this PHE model to simulate Atomic Force Microscopy (AFM) experimental results with the variation of strain rate. This material model will be compared with viscoelastic model to demonstrate the advantages of PHE model. The results have shown that the maximum value of force applied of PHE model is lower at lower strain rates. This is because the mobile fluid does not have enough time to exude in case of very high strain rate and also due to the lower permeability of the membrane than that of the protoplasm of chondrocyte. This behavior is barely observed in viscoelastic model. Thus, PHE model is the better model for cell biomechanics studies.

Feasibility of agent-based modelling of articular cartilage including a conceptual representation of its structure

Articular cartilage is a complex structure with an architecture in which fluid-swollen proteoglycans constrained within a 3D network of collagen fibrils. Because of the complexity of the cartilage structure, the relationship between its mechanical behaviours at the macroscale level and its components at the micro-scale level are not completely understood. The research objective in this thesis is to create a new model of articular cartilage that can be used to simulate and obtain insight into the micro-macro-interaction and mechanisms underlying its mechanical responses during physiological function. The new model of articular cartilage has two characteristics, namely: i) not use fibre-reinforced composite material idealization ii) Provide a framework for that it does probing the micro mechanism of the fluid-solid interaction underlying the deformation of articular cartilage using simple rules of repartition instead of constitutive / physical laws and intuitive curve-fitting. Even though there are various microstructural and mechanical behaviours that can be studied, the scope of this thesis is limited to osmotic pressure formation and distribution and their influence on cartilage fluid diffusion and percolation, which in turn governs the deformation of the compression-loaded tissue. The study can be divided into two stages. In the first stage, the distributions and concentrations of proteoglycans, collagen and water were investigated using histological protocols. Based on this, the structure of cartilage was conceptualised as microscopic osmotic units that consist of these constituents that were distributed according to histological results. These units were repeated three-dimensionally to form the structural model of articular cartilage. In the second stage, cellular automata were incorporated into the resulting matrix (lattice) to simulate the osmotic pressure of the fluid and the movement of water within and out of the matrix; following the osmotic pressure gradient in accordance with the chosen rule of repartition of the pressure. The outcome of this study is the new model of articular cartilage that can be used to simulate and study the micromechanical behaviours of cartilage under different conditions of health and loading. These behaviours are illuminated at the microscale level using the socalled neighbourhood rules developed in the thesis in accordance with the typical requirements of cellular automata modelling. Using these rules and relevant Boundary Conditions to simulate pressure distribution and related fluid motion produced significant results that provided the following insight into the relationships between osmotic pressure gradient and associated fluid micromovement, and the deformation of the matrix. For example, it could be concluded that: 1. It is possible to model articular cartilage with the agent-based model of cellular automata and the Margolus neighbourhood rule. 2. The concept of 3D inter connected osmotic units is a viable structural model for the extracellular matrix of articular cartilage. 3. Different rules of osmotic pressure advection lead to different patterns of deformation in the cartilage matrix, enabling an insight into how this micromechanism influences macromechanical deformation. 4. When features such as transition coefficient were changed, permeability (representing change) is altered due to the change in concentrations of collagen, proteoglycans (i.e. degenerative conditions), the deformation process is impacted. 5. The boundary conditions also influence the relationship between osmotic pressure gradient and fluid movement at the micro-scale level. The outcomes are important to cartilage research since we can use these to study the microscale damage in the cartilage matrix. From this, we are able to monitor related diseases and their progression leading to potential insight into drug-cartilage interaction for treatment. This innovative model is an incremental progress on attempts at creating further computational modelling approaches to cartilage research and other fluid-saturated tissues and material systems.

Influence of grain boundaries on the vibrational properties of silver nanowires

Metal and semiconductor nanowires (NWs) have been widely employed as the building blocks of the nanoelectromechanical systems, which usually acted a resonant beam. Recent researches reported that nanowires are often polycrystalline, which contains grain boundaries (GBs) that transect the whole nanowire into a bamboo like structure. Based on the larger-scale molecular dynamics (MD) simulations, a comprehensive investigation of the influence from grain boundaries on the vibrational properties of doubly clamped Ag NWs is conducted. It is found that, the presence of grain boundary will result in significant energy dissipation during the resonance of polycrystalline NWs, which leads a great deterioration to the quality factor. Further investigation reveals that the energy dissipation is originated from the plastic deformation of polycrystalline NWs in the form of the nucleation of partial dislocations or the generation of micro stacking faults around the GBs and the micro stacking faults is found to keep almost intact during the whole vibration process. Moreover, it is observed that the closer of the grain boundary getting to the regions with the highest strain state, the more energy dissipation will be resulted from the plastic deformation. In addition, either the increase of the number of grain boundaries or the decrease of the distance between the grain boundary and the highest strain state region is observed to induce a lower first resonance frequency. This work sheds lights on the better understanding of the mechanical properties of polycrystalline NWs, which benefits the increasing utilities of NWs in diverse nano-electronic devices.

Assessment of the molecular structure of the borate mineral sakhaite Ca12Mg4(BO3)7(CO3)4Cl(OH)2·H2O using vibrational spectroscopy

The structure of the borate mineral sakhaite Ca12Mg4(BO3)7(CO3)4Cl(OH)2·H2O, a borate–carbonate of calcium and magnesium has been assessed using vibrational spectroscopy. Assignment of bands is undertaken by comparison with the data from other published results. Intense Raman band at 1134 cm−1 with a shoulder at 1123 cm−1 is assigned to the symmetric stretching mode. The Raman spectrum displays bands at 1479, 1524 and 1560 cm−1 which are assigned to the antisymmetric stretching vibrations. The observation of multiple carbonate stretching bands supports the concept that the carbonate units are non-equivalent. The Raman band at 968 cm−1 with a shoulder at 950 cm−1 is assigned to the symmetric stretching mode of trigonal boron. Raman bands at 627 and 651 cm−1 are assigned to the out-of-plane bending modes of trigonal and tetrahedral boron. Raman spectroscopy coupled with infrared spectroscopy enables the molecular structure of the mineral sakhaite to be assessed.

Assessment of the molecular structure of the borate mineral boracite Mg3B7O13Cl using vibrational spectroscopy

Boracite is a magnesium borate mineral with formula: Mg3B7O13Cl and occurs as blue green, colorless, gray, yellow to white crystals in the orthorhombic – pyramidal crystal system. An intense Raman band at 1009 cm−1 was assigned to the BO stretching vibration of the B7O13 units. Raman bands at 1121, 1136, 1143 cm−1 are attributed to the in-plane bending vibrations of trigonal boron. Four sharp Raman bands observed at 415, 494, 621 and 671 cm−1 are simply defined as trigonal and tetrahedral borate bending modes. The Raman spectrum clearly shows intense Raman bands at 3405 and 3494 cm−1, thus indicating that some Cl anions have been replaced with OH units. The molecular structure of a natural boracite has been assessed by using vibrational spectroscopy.

Modelling the flow of juice through a mill

A juice flow model has been developed to estimate the juice expression at the four nips of a sixroller mill. An extended volumetric theory was applied to determine the juice expressed at each nip. The model was applied to a first and final mill, using typical mill settings and an empirical equation to estimate reabsorption. Results of using the model for typical heavy-duty pressure feeder settings show that most of the juice is expressed at the pressure feeder nip. Since the pressure feeders are remote from the mill, a significant portion of the juice is expressed before the bagasse enters the mill.

A vibrational spectroscopic study of philipsbornite PbAl3(AsO4)2(OH)5⋅H2O-molecular structural implications and relationship to the crandallite subgroup arsenates

The presence of arsenic in the environment is a hazard. The accumulation of arsenate by a range of cations in the formation of minerals provides a mechanism for the remediation of arsenate contamination. The formation of the crandallite group of minerals provides a mechanism for arsenate accumulation. Among the crandallite minerals are philipsbornite, arsenocrandallite and arsenogoyazite. Raman spectroscopy complimented with infrared spectroscopy has enabled aspects of the structure of philipsbornite to be studied. The Raman spectrum of philipsbornite displays an intense band at around 840 cm−1 attributed to the overlap of the symmetric and antisymmetric stretching modes. Raman bands observed at 325, 336, 347, 357, 376 and 399 cm−1 are assigned to the ν2 (AsO4)3− symmetric bending vibration (E) and to the ν4 bending vibration (F2). The observation of multiple bending modes supports the concept of a reduction in symmetry of the arsenate anion in philipsbornite. Evidence for phosphate in the mineral is provided. By using an empirical formula, hydrogen bond distances for the OH units in philipsbornite of 2.8648 Å, 2.7864 Å, 2.6896 Å cm−1 and 2.6220 were calculated.

The molecular structure of matioliite – NaMgAl5(PO4)4(OH)6·2(H2O) : a pegmatite mineral from Minas Gerais, Brazil

Detailed spectroscopic and chemical investigation of matioliite, including infrared and Raman spectroscopy, scanning electron microscopy and electron probe microanalysis has been carried out on homogeneous samples from the Gentil pegmatite, Mendes Pimentel, Minas Gerais, Brazil. The chemical composition is (wt.%): FeO 2.20, CaO 0.05, Na2O 1.28, MnO 0.06, Al2O3 39.82, P2O5 42.7, MgO 4.68, F 0.02 and H2O 9.19; total 100.00. The mineral crystallize in the monoclinic crystal system, C2/c space group, with a = 25.075(1) Å, b = 5.0470(3) Å, c = 13.4370(7) Å, β = 110.97(3)°, V = 1587.9(4) Å3, Z = 4. Raman spectroscopy coupled with infrared spectroscopy supports the concept of phosphate, hydrogen phosphate and dihydrogen phosphate units in the structure of matioliite. Infrared and Raman bands attributed to water and hydroxyl stretching modes are identified. Vibrational spectroscopy adds useful information to the molecular structure of matioliite.

The molecular structure of the phosphate mineral chalcosiderite – a vibrational spectroscopic study

The mineral chalcosiderite with formula CuFe6(PO4)4(OH)8⋅4H2O has been studied by Raman spectroscopy and by infrared spectroscopy. A comparison of the chalcosiderite spectra is made with the spectra of turquoise. The spectra of the mineral samples are very similar in the 1200–900 cm−1 region but strong differences are observed in the 900–100 cm−1 region. The effect of substitution of Fe for Al in chalcosiderite shifts the bands to lower wave numbers. Factor group analysis (FGA) implies four OH stretching vibrations for both the water and hydroxyl units. Two bands ascribed to water are observed at 3276 and 3072 cm−1. Three hydroxyl stretching vibrations are observed. Calculations using a Libowitzky type formula show that the hydrogen bond distances of the water molecules are 2.745 and 2.812 Å which are considerably shorter than the values for the hydroxyl units 2.896, 2.917 and 2.978 Å. Two phosphate stretching vibrations at 1042 and 1062 cm−1 in line with the two independent phosphate units in the structure of chalcosiderite. Three bands are observed at 1102, 1159 and 1194 cm−1 assigned to the phosphate antisymmetric stretching vibrations. FGA predicts six bands but only three are observed due to accidental degeneracy. Both the ν2 and ν4 bending regions are complex. Four Raman bands observed at 536, 580, 598 and 636 cm−1 are assigned to the ν4 bending modes. Raman bands at 415, 420, 475 and 484 cm−1are assigned to the phosphate ν2 bending modes. Vibrational spectroscopy enables aspects of the molecular structure of chalcosiderite to be assessed.

Numerical analysis of shape transition in graphene nanoribbons

Graphene nanoribbon (GNR) with free edges can exhibit non-flat morphologies due to pre-existing edge stress. Using molecular dynamics (MD) simulations, we investigate the free-edge effect on the shape transition in GNRs with different edge types, including regular (armchair and zigzag), armchair terminated with hydrogen and reconstructed armchair. The results show that initial edge stress and energy are dependent on the edge configurations. It is confirmed that pre-strain on the free edges is a possible way to limit the random shape transition of GNRs. In addition, the influence of surface attachment on the shape transition is also investigated in this work. It is found that surface attachment can lead to periodic ripples in GNRs, dependent on the initial edge configurations.

Cohesive zone modelling of mineralized collagen fibril arrays in bending

Hard biological materials such as bone possess superior material properties of high stiffness and toughness. Two unique characteristics of bone microstructure are a large aspect ratio of mineralized collagen fibrils (MCF), and an extremely thin and large area of extrafibrillar protein matrix located between the MCF. The objective of this study is to investigate the effects of: (1) MCF aspect ratio, and (2) energy dissipation in extrafibrillar protein matrix on the mechanical behaviour of MCF arrays. In this study, notched specimens of MCF arrays in extrafibrillar protein matrix are subjected to bending. Cohesive zone model was implemented to simulate the failure of extrafibrillar protein matrix. The study reveals that the MCF array with a higher MCF aspect ratio and the MCF array with a higher protein energy dissipation in the interface direction are able to sustain a higher bending force and dissipate higher energy.

Molecular dynamics simulation of fracture strength and morphology of defective graphene

Different types of defects can be introduced into graphene during material synthesis, and significantly influence the properties of graphene. In this work, we investigated the effects of structural defects, edge functionalisation and reconstruction on the fracture strength and morphology of graphene by molecular dynamics simulations. The minimum energy path analysis was conducted to investigate the formation of Stone-Wales defects. We also employed out-of-plane perturbation and energy minimization principle to study the possi-ble morphology of graphene nanoribbons with edge-termination. Our numerical results show that the fracture strength of graphene is dependent on defects and environmental temperature. However, pre-existing defects may be healed, resulting in strength recovery. Edge functionalization can induce compressive stress and ripples in the edge areas of gra-phene nanoribbons. On the other hand, edge reconstruction contributed to the tensile stress and curved shape in the graphene nanoribbons.

Risk modelling in quality clinical registries : Monitoring lesion treatment failure rate in percutaneous coronary interventions

Aims: This paper describes the development of a risk adjustment (RA) model predictive of individual lesion treatment failure in percutaneous coronary interventions (PCI) for use in a quality monitoring and improvement program. Methods and results: Prospectively collected data for 3972 consecutive revascularisation procedures (5601 lesions) performed between January 2003 and September 2011 were studied. Data on procedures to September 2009 (n = 3100) were used to identify factors predictive of lesion treatment failure. Factors identified included lesion risk class (p < 0.001), occlusion type (p < 0.001), patient age (p = 0.001), vessel system (p < 0.04), vessel diameter (p < 0.001), unstable angina (p = 0.003) and presence of major cardiac risk factors (p = 0.01). A Bayesian RA model was built using these factors with predictive performance of the model tested on the remaining procedures (area under the receiver operating curve: 0.765, Hosmer–Lemeshow p value: 0.11). Cumulative sum, exponentially weighted moving average and funnel plots were constructed using the RA model and subjectively evaluated. Conclusion: A RA model was developed and applied to SPC monitoring for lesion failure in a PCI database. If linked to appropriate quality improvement governance response protocols, SPC using this RA tool might improve quality control and risk management by identifying variation in performance based on a comparison of observed and expected outcomes.

Numerical Characterization of Nanowires

In this chapter, we will present a contemporary review of the hitherto numerical characterization of nanowires (NWs). The bulk of the research reported in the literatures concern metallic NWs including Al, Cu, Au, Ag, Ni, and their alloys NWs. Research has also been reported for the investigation of some nonmetallic NWs, such as ZnO, GaN, SiC, SiO2. A plenty of researches have been conducted regarding the numerical investigation of NWs. Issues analyzed include structural changes under different loading situations, the formation and propagation of dislocations, and the effect of the magnitude of applied loading on deformation mechanics. Efforts have also been made to correlate simulation results with experimental measurements. However, direct comparisons are difficult since most simulations are carried out under conditions of extremely high strain/loading rates and small simulation samples due to computational limitations. Despite of the immense numerical studies of NWs, a significant work still lies ahead in terms of problem formulation, interpretation of results, identification and delineation of deformation mechanisms, and constitutive characterization of behavior. In this chapter, we present an introduction of the commonly adopted experimental and numerical approaches in studies of the deformation of NWs in Section 1. An overview of findings concerning perfect NWs under different loading situations, such as tension, compression, torsion, and bending are presented in Section 2. In Section 3, we will detail some recent results from the authors’ own work with an emphasis on the study of influences from different pre-existing defect on NWs. Some thoughts on future directions of the computational mechanics of NWs together with Conclusions will be given in the last section.