The effect of osteocyte apoptosis on signalling in the osteocyte and bone lining cell network: a computer simulation

Osteocytes play a critical role in the regulation of bone remodelling by translating strain due to mechanical loading into biochemical signals transmitted through the interconnecting lacuno-canalicular network to bone lining cells (BLCs) on the bone surface. This work aims to examine the effects of disruption of that intercellular communication by simulation of osteocyte apoptosis in the bone matrix. A model of a uniformly distributed osteocyte network has been developed that simulates the signalling through the network to the BLCs based on strain level. Bi-directional and asymmetric communication between neighbouring osteocytes and BLCs is included. The effect of osteocyte apoptosis is examined by preventing signalling at and through the affected cells. The simulation shows that apoptosis of only 3% of the osteocyte cells leads to a significant reduction in the peak signal at the BLCs. Furthermore, experiments with the model confirm how important the location and density of the apoptotic osteocytes are to the signalling received at the bone surface. With 5% and 9% osteocyte apoptosis, the mean peak BLC levels were reduced by 25% and 37% respectively. Such a significant reduction in the signal at the BLCs may explain a possible mechanism that leads to the increased remodelling and eventual bone loss observed with osteoporosis. More generally, it provides a unique framework for a broader exploration of the role of osteocyte and bi-directional and asymmetric cell-cell communication in mechanotransduction, and the effects of disruption to that communication.

Masticatory loadings and cranial deformation in Macaca fascicularis: a finite element analysis sensitivity study

Biomechanical analyses are commonly conducted to investigate how craniofacial form relates to function, particularly in relation to dietary adaptations. However, in the absence of corresponding muscle activation patterns, incomplete muscle data recorded experimentally for different individuals during different feeding tasks are frequently substituted. This study uses finite element analysis (FEA) to examine the sensitivity of the mechanical response of a Macaca fascicularis cranium to varying muscle activation patterns predicted via multibody dynamic analysis. Relative to the effects of varying bite location, the consequences of simulated variations in muscle activation patterns and of the inclusion/exclusion of whole muscle groups were investigated. The resulting cranial deformations were compared using two approaches; strain maps and geometric morphometric analyses. The results indicate that, with bite force magnitude controlled, the variations among the mechanical responses of the cranium to bite location far outweigh those observed as a consequence of varying muscle activations. However, zygomatic deformation was an exception, with the activation levels of superficial masseter being most influential in this regard. The anterior portion of temporalis deforms the cranial vault, but the remaining muscles have less profound effects. This study for the first time systematically quantifies the sensitivity of an FEA model of a primate skull to widely varying masticatory muscle activations and finds that, with the exception of the zygomatic arch, reasonable variants of muscle loading for a second molar bite have considerably less effect on cranial deformation and the resulting strain map than does varying molar bite point. The implication is that FEA models of biting crania will generally produce acceptable estimates of deformation under load as long as muscle activations and forces are reasonably approximated. In any one FEA study, the biological significance of the error in applied muscle forces is best judged against the magnitude of the effect that is being investigated.

Shearing mechanics and the influence of a flexible symphysis during oral food processing in Sphenodon (Lepidosauria: Rhynchocephalia)

The New Zealand tuatara, Sphenodon, has a specialized feeding system in which the teeth of the lower jaw close between two upper tooth rows before sliding forward to slice food apart like a draw cut saw. This shearing action is unique amongst living amniotes but has been compared with the chewing power stroke of mammals. We investigated details of the jaw movement using multibody dynamics analysis of an anatomically accurate three-dimensional computer model constructed from computed tomography scans. The model predicts that a flexible symphysis is necessary for changes in the intermandibular angle that permits prooral movement. Models with the greatest symphysial flexibility allow the articulation surface of the articular to follow the quadrate cotyle with the least restriction, and suggest that shearing is accompanied by a long axis rotation of the lower jaws. This promotes precise point loading between the cutting edges of particular teeth, enhancing the effectiveness of the shearing action. Given that Sphenodon is a relatively inactive reptile, we suggest that the link between oral food processing and endothermy has been overstated. Food processing improves feeding efficiency, a consideration of particular importance when food availability is unpredictable. Although this feeding mechanism is today limited to Sphenodon, a survey of fossil rhynchocephalians suggests that it was once more widespread.

Functional evolution of the feeding system in rodents

The masticatory musculature of rodents has evolved to enable both gnawing at the incisors and chewing at the molars. In particular, the masseter muscle is highly specialised, having extended anteriorly to originate from the rostrum. All living rodents have achieved this masseteric expansion in one of three ways, known as the sciuromorph, hystricomorph and myomorph conditions. Here, we used finite element analysis (FEA) to investigate the biomechanical implications of these three morphologies, in a squirrel, guinea pig and rat. In particular, we wished to determine whether each of the three morphologies is better adapted for either gnawing or chewing. Results show that squirrels are more efficient at muscle-bite force transmission during incisor gnawing than guinea pigs, and that guinea pigs are more efficient at molar chewing than squirrels. This matches the known diet of nuts and seeds that squirrels gnaw, and of grasses that guinea pigs grind down with their molars. Surprisingly, results also indicate that rats are more efficient as well as more versatile feeders than both the squirrel and guinea pig. There seems to be no compromise in biting efficiency to accommodate the wider range of foodstuffs and the more general feeding behaviour adopted by rats. Our results show that the morphology of the skull and masticatory muscles have allowed squirrels to specialise as gnawers and guinea pigs as chewers, but that rats are high-performance generalists, which helps explain their overwhelming success as a group.

Intraluminal thrombus has a selective influence on matrix metalloproteinases and their inhibitors (tissue inhibitors of matrix metalloproteinases) in the wall of abdominal aortic aneurysms

BACKGROUND

The influence of intraluminal thrombus (ILT) on the proteolytic environment within the wall of an abdominal aortic aneurysm (AAA) is unknown. This is the first study to examine the correlation between ILT thickness and the levels of matrix metalloproteinases (MMPs) and their natural inhibitors (tissue inhibitors of matrix metalloproteinases [TIMPs]) within the adjacent AAA wall.

METHODS

Thirty-five patients undergoing elective repair of AAAs were studied. A single full-thickness infrarenal aortic sample was obtained uniformly from the arteriotomy site from each patient. All samples were snap frozen and analyzed for total and active MMP 2, 8, and 9 and TIMP 1 and 2. Thrombus thickness at the specimen site was measured on the preoperative contrast computed tomographic angiograms.

RESULTS

There was a statistically significant correlation between ILT thickness, concentration of TIMP 1, and active concentration of MMP 9. MMP 2 (active and total) and TIMP 2 demonstrated a positive correlation with ILT thickness, although not statistically significant.

CONCLUSION

In this novel study, we found a significant positive correlation of ILT thickness with active MMP 9 and TIMP 1 concentration in the adjacent AAA wall, and this may have implications for AAA expansion and eventual rupture.

Functional relationship between skull form and feeding mechanics in Sphenodon, and implications for diapsid skull development

The vertebrate skull evolved to protect the brain and sense organs, but with the appearance of jaws and associated forces there was a remarkable structural diversification. This suggests that the evolution of skull form may be linked to these forces, but an important area of debate is whether bone in the skull is minimised with respect to these forces, or whether skulls are mechanically "over-designed" and constrained by phylogeny and development. Mechanical analysis of diapsid reptile skulls could shed light on this longstanding debate. Compared to those of mammals, the skulls of many extant and extinct diapsids comprise an open framework of fenestrae (window-like openings) separated by bony struts (e.g., lizards, tuatara, dinosaurs and crocodiles), a cranial form thought to be strongly linked to feeding forces. We investigated this link by utilising the powerful engineering approach of multibody dynamics analysis to predict the physiological forces acting on the skull of the diapsid reptile Sphenodon. We then ran a series of structural finite element analyses to assess the correlation between bone strain and skull form. With comprehensive loading we found that the distribution of peak von Mises strains was particularly uniform throughout the skull, although specific regions were dominated by tensile strains while others were dominated by compressive strains. Our analyses suggest that the frame-like skulls of diapsid reptiles are probably optimally formed (mechanically ideal: sufficient strength with the minimal amount of bone) with respect to functional forces; they are efficient in terms of having minimal bone volume, minimal weight, and also minimal energy demands in maintenance.

Finite element modelling of squirrel, guinea pig and rat skulls: using geometric morphometrics to assess sensitivity

Rodents are defined by a uniquely specialized dentition and a highly complex arrangement of jaw-closing muscles. Finite element analysis (FEA) is an ideal technique to investigate the biomechanical implications of these specializations, but it is essential to understand fully the degree of influence of the different input parameters of the FE model to have confidence in the model's predictions. This study evaluates the sensitivity of FE models of rodent crania to elastic properties of the materials, loading direction, and the location and orientation of the models' constraints. Three FE models were constructed of squirrel, guinea pig and rat skulls. Each was loaded to simulate biting on the incisors, and the first and the third molars, with the angle of the incisal bite varied over a range of 45°. The Young's moduli of the bone and teeth components were varied between limits defined by findings from our own and previously published tests of material properties. Geometric morphometrics (GMM) was used to analyse the resulting skull deformations. Bone stiffness was found to have the strongest influence on the results in all three rodents, followed by bite position, and then bite angle and muscle orientation. Tooth material properties were shown to have little effect on the deformation of the skull. The effect of bite position varied between species, with the mesiodistal position of the biting tooth being most important in squirrels and guinea pigs, whereas bilateral vs. unilateral biting had the greatest influence in rats. A GMM analysis of isolated incisor deformations showed that, for all rodents, bite angle is the most important parameter, followed by elastic properties of the tooth. The results here elucidate which input parameters are most important when defining the FE models, but also provide interesting glimpses of the biomechanical differences between the three skulls, which will be fully explored in future publications.

Why do humans have chins? Testing the mechanical significance of modern human symphyseal morphology with finite element analysis

The modern human mandibular symphysis differs from those of all other primates in being vertically orientated and possessing a chin, but the functional significance of this unique morphology is not well understood. Some hypotheses propose that it is an adaptation to specific loads occurring during masticatory function. This study uses finite element analysis to examine these symphyseal loads in a model of a modern human mandible. By modifying the symphyseal cross-sectional form, the mechanical significance of the presence of the chin and symphyseal orientation is tested, and modern human and Neanderthal symphyseal cross-sections are compared with regard to their ability to withstand different loads. The results show that changes in symphyseal form have profound effects on the strains. The presence of a chin leads to lower symphyseal strains overall, whereas a vertical orientation of the symphysis results in higher strains under wishboning, but not under vertical bending in the coronal plane and dorsoventral shear. Compared to Neanderthals, the modern human symphysis shows higher strains during dorsoventral shear and wishboning, but is as effective as the Neanderthal symphysis in resisting vertical bending in the coronal plane and the loads resulting from simulated incision and unilateral molar biting. In general, the results of this study corroborate prior hypotheses about the mechanical effects of the human chin and vertical symphyseal orientation and support the idea that the relative importance of wishboning and vertical bending in the coronal plane might have played a role in the evolution of modern human symphyseal morphology.

Combining geometric morphometrics and functional simulation: an emerging toolkit for virtual functional analyses

The development of virtual methods for anatomical reconstruction and functional simulation of skeletal structures offers great promise in evolutionary and ontogenetic investigations of form-function relationships. Key developments reviewed here include geometric morphometric methods for the analysis and visualization of variations in form (size and shape), finite element methods for the prediction of mechanical performance of skeletal structures under load and multibody dynamics methods for the simulation and prediction of musculoskeletal function. These techniques are all used in studies of form and function in biology, but only recently have they been combined in novel ways to facilitate biomechanical modelling that takes account of variations in form, can statistically compare performance, and relate performance to form and its covariates. Here we provide several examples that illustrate how these approaches can be combined and we highlight areas that require further investigation and development before we can claim a mature theory and toolkit for a statistical biomechanical framework that unites these methods.

Isolated word recognition of silent speech using magnetic implants and sensors

There are a number of situations where individuals wish to communicate verbally but are unable to use conventional means-so called 'silent speech'. These include speakers in noisy and covert situations as well as patients who have lost their voice as a result of a laryngectomy or similar procedure. This paper focuses on those who are unable to speak following a laryngectomy and assesses the possibility of speech recognition based on a magnetic implant/sensors system. Permanent magnets are placed on the tongue and lips and the changes in magnetic field resulting from movement during speech are monitored using a set of magnetic sensors. The sensor signals are compared to sets of pre-recorded templates using the dynamic time warping (DTW) method, and the best match is identified. Experimental trials are reported for subjects with intact larynx, typically using 500-1000 utterances used for speaker dependant training and testing. It is shown that recognition rates of over 90% are achievable for vocabularies of at least 57 isolated words: sufficient to drive command-and-control applications.

Validation of a Morphometric Reconstruction Technique Applied to a Juvenile Pelvis

Three-dimensional reconstructions of bone geometry from microCT (computed tomography) data are frequently used in biomechanical and finite element analyses. Digitization of bone models is usually a simple process for specimens with a complete geometry, but in instances of damage or disarticulation it can be very challenging. Subsequent to digitization, further imaging techniques are often required to estimate the geometry of missing bone or connecting cartilage. This paper presents an innovative approach to the reconstruction of incomplete scan data, to reproduce proper anatomical arrangements of bones, including absent connecting cartilaginous elements. Utilizing geometric morphometric tools, the reconstruction technique is validated through comparison of a reconstructed 9 year old pelvis, to the original CT data. A principal component analysis and an overlay of the two pelves provide a measure of the accuracy of the reconstructed model. Future work aims to investigate the biomechanical effects of any minor positional error on the bone's predicted structural properties through the use of finite element analysis.

Predicting muscle activation patterns from motion and anatomy: modelling the skull of Sphenodon (Diapsida: Rhynchocephalia)

The relationship between skull shape and the forces generated during feeding is currently under widespread scrutiny and increasingly involves the use of computer simulations such as finite element analysis. The computer models used to represent skulls are often based on computed tomography data and thus are structurally accurate; however, correctly representing muscular loading during food reduction remains a major problem. Here, we present a novel approach for predicting the forces and activation patterns of muscles and muscle groups based on their known anatomical orientation (line of action). The work was carried out for the lizard-like reptile Sphenodon (Rhynchocephalia) using a sophisticated computer-based model and multi-body dynamics analysis. The model suggests that specific muscle groups control specific motions, and that during certain times in the bite cycle some muscles are highly active whereas others are inactive. The predictions of muscle activity closely correspond to data previously recorded from live Sphenodon using electromyography. Apparent exceptions can be explained by variations in food resistance, food size, food position and lower jaw motions. This approach shows considerable promise in advancing detailed functional models of food acquisition and reduction, and for use in other musculoskeletal systems where no experimental determination of muscle activity is possible, such as in rare, endangered or extinct species.

Biomechanical assessment of evolutionary changes in the lepidosaurian skull

The lepidosaurian skull has long been of interest to functional morphologists and evolutionary biologists. Patterns of bone loss and gain, particularly in relation to bars and fenestrae, have led to a variety of hypotheses concerning skull use and kinesis. Of these, one of the most enduring relates to the absence of the lower temporal bar in squamates and the acquisition of streptostyly. We performed a series of computer modeling studies on the skull of Uromastyx hardwickii, an akinetic herbivorous lizard. Multibody dynamic analysis (MDA) was conducted to predict the forces acting on the skull, and the results were transferred to a finite element analysis (FEA) to estimate the pattern of stress distribution. In the FEA, we applied the MDA result to a series of models based on the Uromastyx skull to represent different skull configurations within past and present members of the Lepidosauria. In this comparative study, we found that streptostyly can reduce the joint forces acting on the skull, but loss of the bony attachment between the quadrate and pterygoid decreases skull robusticity. Development of a lower temporal bar apparently provided additional support for an immobile quadrate that could become highly stressed during forceful biting.

Assessment of the role of sutures in a lizard skull: a computer modelling study

Sutures form an integral part of the functioning skull, but their role has long been debated among vertebrate morphologists and palaeontologists. Furthermore, the relationship between typical skull sutures, and those involved in cranial kinesis, is poorly understood. In a series of computational modelling studies, complex loading conditions obtained through multibody dynamics analysis were imposed on a finite element model of the skull of Uromastyx hardwickii, an akinetic herbivorous lizard. A finite element analysis (FEA) of a skull with no sutures revealed higher patterns of strain in regions where cranial sutures are located in the skull. From these findings, FEAs were performed on skulls with sutures (individual and groups of sutures) to investigate their role and function more thoroughly. Our results showed that individual sutures relieved strain locally, but only at the expense of elevated strain in other regions of the skull. These findings provide an insight into the behaviour of sutures and show how they are adapted to work together to distribute strain around the skull. Premature fusion of one suture could therefore lead to increased abnormal loading on other regions of the skull causing irregular bone growth and deformities. This detailed investigation also revealed that the frontal-parietal suture of the Uromastyx skull played a substantial role in relieving strain compared with the other sutures. This raises questions about the original role of mesokinesis in squamate evolution.

Combined finite element and multibody dynamics analysis of biting in a Uromastyx hardwickii lizard skull

Lizard skulls vary greatly in shape and construction, and radical changes in skull form during evolution have made this an intriguing subject of research. The mechanics of feeding have surely been affected by this change in skull form, but whether this is the driving force behind the change is the underlying question that we are aiming to address in a programme of research. Here we have implemented a combined finite element analysis (FEA) and multibody dynamics analysis (MDA) to assess skull biomechanics during biting. A skull of Uromastyx hardwickii was assessed in the present study, where loading data (such as muscle force, bite force and joint reaction) for a biting cycle were obtained from an MDA and applied to load a finite element model. Fifty load steps corresponding to bilateral biting towards the front, middle and back of the dentition were implemented. Our results show the importance of performing MDA as a preliminary step to FEA, and provide an insight into the variation of stress during biting. Our findings show that higher stress occurs in regions where cranial sutures are located in functioning skulls, and as such support the hypothesis that sutures may play a pivotal role in relieving stress and producing a more uniform pattern of stress distribution across the skull. Additionally, we demonstrate how varying bite point affects stress distributions and relate stress distributions to the evolution of metakinesis in the amniote skull.

Rigid-body analysis of a lizard skull: modelling the skull of Uromastyx hardwickii

Lizard skulls vary greatly in their detailed morphology. Theoretical models and practical studies have posited a definite relationship between skull morphology and bite performance, but this can be difficult to demonstrate in vivo. Computer modelling provides an alternative approach, as long as hard and soft tissue components can be integrated and the model can be validated. An anatomically accurate three-dimensional computer model of an Uromastyx hardwickii skull was developed for rigid-body dynamic analysis. The Uromastyx jaw was first opened under motion control, and then muscle forces were applied to produce biting simulations where bite forces and joint forces were calculated. Bite forces comparable to those reported in the literature were predicted, and detailed muscular force information was produced along with additional information on the stabilizing role of temporal ligaments in late jaw closing.

Assessing mechanical function of the zygomatic region in macaques: validation and sensitivity testing of finite element models

Crucial to the interpretation of the results of any finite element analysis of a skeletal system is a test of the validity of the results and an assessment of the sensitivity of the model parameters. We have therefore developed finite element models of two crania of Macaca fascicularis and investigated their sensitivity to variations in bone material properties, the zygomatico-temporal suture and the loading regimen applied to the zygomatic arch. Maximum principal strains were validated against data derived from ex vivo strain gauge experiments using non-physiological loads applied to the macaque zygomatic arch. Elastic properties of the zygomatic arch bone and the zygomatico-temporal suture obtained by nanoindentation resulted in a high degree of congruence between experimental and simulated strains. The findings also indicated that the presence of a zygomatico-temporal suture in the model produced strains more similar to experimental values than a completely separated or fused arch. Strains were distinctly higher when the load was applied through the modelled superficial masseter compared with loading an array of nodes on the arch. This study demonstrates the importance of the accurate selection of the material properties involved in predicting strains in a finite element model. Furthermore, our findings strongly highlight the influence of the presence of craniofacial sutures on strains experienced in the face. This has implications when investigating craniofacial growth and masticatory function but should generally be taken into account in functional analyses of the craniofacial system of both extant and extinct species.

Three dimensional stereolithography models of cancellous bone structures from muCT data: testing and validation of finite element results

Stereolithography (STL) models of complex cancellous bone structures have been produced from three-dimensional micro-computed tomography data sets of human cancellous bone histological samples from four skeletal sites. The STL models have been mechanically tested and the derived stiffness compared with that predicted by finite element analysis. The results show a strong correlation (R2 = 0.941) between the predicted and calculated stiffnesses of the structures and show promise for the use of STL as an additional technique to complement the use of finite element models, for the assessment of the mechanical properties of complex cancellous bone structures.

Measurements of tortuosity in stereolithographical bone replicas using audiofrequency pulses

The tortuosity of five air-filled stereolithographical cancellous bone replicas has been obtained from measurements using audiofrequency pulses in a rectangular waveguide. The data obtained from the replicas yields information about anisotropy with respect to orthogonal axes of the passages that would be marrow filled in vivo. A strong relationship has been found between the acoustically measured tortuosity and the independently measured porosity. Use of stereolithographical bone replicas has the potential to simulate perforation and thinning of cancellous bone and hence evaluate the dependence of acoustic properties on cancellous bone microstructure. As an "extreme" illustration of such use, "inverses" of the original replicas have been manufactured and acoustic measurements have been made on them. The data reveal significantly greater tortuosity of the passages that are geometrically equivalent to the original solid bone structures.