Development of musculoskeletal models for the design and the pre-clinical validation of hip resurfacing prosthesis

Background. The surgical treatment of dysfunctional hips is a severe condition for the patient and a costly therapy for the public health. Hip resurfacing techniques seem to hold the promise of various advantages over traditional THR, with particular attention to young and active patients. Although the lesson provided in the past by many branches of engineering is that success in designing competitive products can be achieved only by predicting the possible scenario of failure, to date the understanding of the implant quality is poorly pre-clinically addressed. Thus revision is the only delayed and reliable end point for assessment. The aim of the present work was to model the musculoskeletal system so as to develop a protocol for predicting failure of hip resurfacing prosthesis. Methods. Preliminary studies validated the technique for the generation of subject specific finite element (FE) models of long bones from Computed Thomography data. The proposed protocol consisted in the numerical analysis of the prosthesis biomechanics by deterministic and statistic studies so as to assess the risk of biomechanical failure on the different operative conditions the implant might face in a population of interest during various activities of daily living. Physiological conditions were defined including the variability of the anatomy, bone densitometry, surgery uncertainties and published boundary conditions at the hip. The protocol was tested by analysing a successful design on the market and a new prototype of a resurfacing prosthesis. Results. The intrinsic accuracy of models on bone stress predictions (RMSE < 10%) was aligned to the current state of the art in this field. The accuracy of prediction on the bone-prosthesis contact mechanics was also excellent (< 0.001 mm). The sensitivity of models prediction to uncertainties on modelling parameter was found below 8.4%. The analysis of the successful design resulted in a very good agreement with published retrospective studies. The geometry optimisation of the new prototype lead to a final design with a low risk of failure. The statistical analysis confirmed the minimal risk of the optimised design over the entire population of interest. The performances of the optimised design showed a significant improvement with respect to the first prototype (+35%). Limitations. On the authors opinion the major limitation of this study is on boundary conditions. The muscular forces and the hip joint reaction were derived from the few data available in the literature, which can be considered significant but hardly representative of the entire variability of boundary conditions the implant might face over the patients population. This moved the focus of the research on modelling the musculoskeletal system; the ongoing activity is to develop subject-specific musculoskeletal models of the lower limb from medical images. Conclusions. The developed protocol was able to accurately predict known clinical outcomes when applied to a well-established device and, to support the design optimisation phase providing important information on critical characteristics of the patients when applied to a new prosthesis. The presented approach does have a relevant generality that would allow the extension of the protocol to a large set of orthopaedic scenarios with minor changes. Hence, a failure mode analysis criterion can be considered a suitable tool in developing new orthopaedic devices.

Numerical simulations of magma chamber dynamics at Campi Flegrei, and associated seismicity, deformation and gravity changes

Understanding the complex relationships between quantities measured by volcanic monitoring network and shallow magma processes is a crucial headway for the comprehension of volcanic processes and a more realistic evaluation of the associated hazard. This question is very relevant at Campi Flegrei, a volcanic quiescent caldera immediately north-west of Napoli (Italy). The system activity shows a high fumarole release and periodic ground slow movement (bradyseism) with high seismicity. This activity, with the high people density and the presence of military and industrial buildings, makes Campi Flegrei one of the areas with higher volcanic hazard in the world. In such a context my thesis has been focused on magma dynamics due to the refilling of shallow magma chambers, and on the geophysical signals detectable by seismic, deformative and gravimetric monitoring networks that are associated with this phenomenologies. Indeed, the refilling of magma chambers is a process frequently occurring just before a volcanic eruption; therefore, the faculty of identifying this dynamics by means of recorded signal analysis is important to evaluate the short term volcanic hazard. The space-time evolution of dynamics due to injection of new magma in the magma chamber has been studied performing numerical simulations with, and implementing additional features in, the code GALES (Longo et al., 2006), recently developed and still on the upgrade at the Istituto Nazionale di Geofisica e Vulcanologia in Pisa (Italy). GALES is a finite element code based on a physico-mathematical two dimensional, transient model able to treat fluids as multiphase homogeneous mixtures, compressible to incompressible. The fundamental equations of mass, momentum and energy balance are discretised both in time and space using the Galerkin Least-Squares and discontinuity-capturing stabilisation technique. The physical properties of the mixture are computed as a function of local conditions of magma composition, pressure and temperature.The model features enable to study a broad range of phenomenologies characterizing pre and sin-eruptive magma dynamics in a wide domain from the volcanic crater to deep magma feeding zones. The study of displacement field associated with the simulated fluid dynamics has been carried out with a numerical code developed by the Geophysical group at the University College Dublin (O’Brien and Bean, 2004b), with whom we started a very profitable collaboration. In this code, the seismic wave propagation in heterogeneous media with free surface (e.g. the Earth’s surface) is simulated using a discrete elastic lattice where particle interactions are controlled by the Hooke’s law. This method allows to consider medium heterogeneities and complex topography. The initial and boundary conditions for the simulations have been defined within a coordinate project (INGV-DPC 2004-06 V3_2 “Research on active volcanoes, precursors, scenarios, hazard and risk - Campi Flegrei”), to which this thesis contributes, and many researchers experienced on Campi Flegrei in volcanological, seismic, petrological, geochemical fields, etc. collaborate. Numerical simulations of magma and rock dynamis have been coupled as described in the thesis. The first part of the thesis consists of a parametric study aimed at understanding the eect of the presence in magma of carbon dioxide in magma in the convection dynamics. Indeed, the presence of this volatile was relevant in many Campi Flegrei eruptions, including some eruptions commonly considered as reference for a future activity of this volcano. A set of simulations considering an elliptical magma chamber, compositionally uniform, refilled from below by a magma with volatile content equal or dierent from that of the resident magma has been performed. To do this, a multicomponent non-ideal magma saturation model (Papale et al., 2006) that considers the simultaneous presence of CO2 and H2O, has been implemented in GALES. Results show that the presence of CO2 in the incoming magma increases its buoyancy force promoting convection ad mixing. The simulated dynamics produce pressure transients with frequency and amplitude in the sensitivity range of modern geophysical monitoring networks such as the one installed at Campi Flegrei . In the second part, simulations more related with the Campi Flegrei volcanic system have been performed. The simulated system has been defined on the basis of conditions consistent with the bulk of knowledge of Campi Flegrei and in particular of the Agnano-Monte Spina eruption (4100 B.P.), commonly considered as reference for a future high intensity eruption in this area. The magmatic system has been modelled as a long dyke refilling a small shallow magma chamber; magmas with trachytic and phonolitic composition and variable volatile content of H2O and CO2 have been considered. The simulations have been carried out changing the condition of magma injection, the system configuration (magma chamber geometry, dyke size) and the resident and refilling magma composition and volatile content, in order to study the influence of these factors on the simulated dynamics. Simulation results allow to follow each step of the gas-rich magma ascent in the denser magma, highlighting the details of magma convection and mixing. In particular, the presence of more CO2 in the deep magma results in more ecient and faster dynamics. Through this simulations the variation of the gravimetric field has been determined. Afterward, the space-time distribution of stress resulting from numerical simulations have been used as boundary conditions for the simulations of the displacement field imposed by the magmatic dynamics on rocks. The properties of the simulated domain (rock density, P and S wave velocities) have been based on data from literature on active and passive tomographic experiments, obtained through a collaboration with A. Zollo at the Dept. of Physics of the Federici II Univeristy in Napoli. The elasto-dynamics simulations allow to determine the variations of the space-time distribution of deformation and the seismic signal associated with the studied magmatic dynamics. In particular, results show that these dynamics induce deformations similar to those measured at Campi Flegrei and seismic signals with energies concentrated on the typical frequency bands observed in volcanic areas. The present work shows that an approach based on the solution of equations describing the physics of processes within a magmatic fluid and the surrounding rock system is able to recognise and describe the relationships between geophysical signals detectable on the surface and deep magma dynamics. Therefore, the results suggest that the combined study of geophysical data and informations from numerical simulations can allow in a near future a more ecient evaluation of the short term volcanic hazard.

The assessment of change in rural landscape and the analysis of the driving forces. Proposal of a spatial model for the rural built environment

The presented study carried out an analysis on rural landscape changes. In particular the study focuses on the understanding of driving forces acting on the rural built environment using a statistical spatial model implemented through GIS techniques. It is well known that the study of landscape changes is essential for a conscious decision making in land planning. From a bibliography review results a general lack of studies dealing with the modeling of rural built environment and hence a theoretical modelling approach for such purpose is needed. The advancement in technology and modernity in building construction and agriculture have gradually changed the rural built environment. In addition, the phenomenon of urbanization of a determined the construction of new volumes that occurred beside abandoned or derelict rural buildings. Consequently there are two types of transformation dynamics affecting mainly the rural built environment that can be observed: the conversion of rural buildings and the increasing of building numbers. It is the specific aim of the presented study to propose a methodology for the development of a spatial model that allows the identification of driving forces that acted on the behaviours of the building allocation. In fact one of the most concerning dynamic nowadays is related to an irrational expansion of buildings sprawl across landscape. The proposed methodology is composed by some conceptual steps that cover different aspects related to the development of a spatial model: the selection of a response variable that better describe the phenomenon under study, the identification of possible driving forces, the sampling methodology concerning the collection of data, the most suitable algorithm to be adopted in relation to statistical theory and method used, the calibration process and evaluation of the model. A different combination of factors in various parts of the territory generated favourable or less favourable conditions for the building allocation and the existence of buildings represents the evidence of such optimum. Conversely the absence of buildings expresses a combination of agents which is not suitable for building allocation. Presence or absence of buildings can be adopted as indicators of such driving conditions, since they represent the expression of the action of driving forces in the land suitability sorting process. The existence of correlation between site selection and hypothetical driving forces, evaluated by means of modeling techniques, provides an evidence of which driving forces are involved in the allocation dynamic and an insight on their level of influence into the process. GIS software by means of spatial analysis tools allows to associate the concept of presence and absence with point futures generating a point process. Presence or absence of buildings at some site locations represent the expression of these driving factors interaction. In case of presences, points represent locations of real existing buildings, conversely absences represent locations were buildings are not existent and so they are generated by a stochastic mechanism. Possible driving forces are selected and the existence of a causal relationship with building allocations is assessed through a spatial model. The adoption of empirical statistical models provides a mechanism for the explanatory variable analysis and for the identification of key driving variables behind the site selection process for new building allocation. The model developed by following the methodology is applied to a case study to test the validity of the methodology. In particular the study area for the testing of the methodology is represented by the New District of Imola characterized by a prevailing agricultural production vocation and were transformation dynamic intensively occurred. The development of the model involved the identification of predictive variables (related to geomorphologic, socio-economic, structural and infrastructural systems of landscape) capable of representing the driving forces responsible for landscape changes.. The calibration of the model is carried out referring to spatial data regarding the periurban and rural area of the study area within the 1975-2005 time period by means of Generalised linear model. The resulting output from the model fit is continuous grid surface where cells assume values ranged from 0 to 1 of probability of building occurrences along the rural and periurban area of the study area. Hence the response variable assesses the changes in the rural built environment occurred in such time interval and is correlated to the selected explanatory variables by means of a generalized linear model using logistic regression. Comparing the probability map obtained from the model to the actual rural building distribution in 2005, the interpretation capability of the model can be evaluated. The proposed model can be also applied to the interpretation of trends which occurred in other study areas, and also referring to different time intervals, depending on the availability of data. The use of suitable data in terms of time, information, and spatial resolution and the costs related to data acquisition, pre-processing, and survey are among the most critical aspects of model implementation. Future in-depth studies can focus on using the proposed model to predict short/medium-range future scenarios for the rural built environment distribution in the study area. In order to predict future scenarios it is necessary to assume that the driving forces do not change and that their levels of influence within the model are not far from those assessed for the time interval used for the calibration.

Cross-Talk tra Bifidobacterium e intestino umano: impatto sulle attività health promoting

Bifidobacteria constitute up to 3% of the total microbiota and represent one of the most important healthpromoting bacterial groups of the human intestinal microflora. The presence of Bifidobacterium in the human gastrointestinal tract has been directly related to several health-promoting activities; however, to date, no information about the specific mechanisms of interaction with the host is available. The first health-promoting activities studied in these job was the oxalate-degrading activity. Oxalic acid occurs extensively in nature and plays diverse roles, especially in pathological processes. Due to its highly oxidizing effects, hyper absorption or abnormal synthesis of oxalate can cause serious acute disorders in mammals and be lethal in extreme cases. Intestinal oxalate-degrading bacteria could therefore be pivotal in maintaining oxalate homeostasis, reducing the risk of kidney stone development. In this study, the oxalate-degrading activity of 14 bifidobacterial strains was measured by a capillary electrophoresis technique. The oxc gene, encoding oxalyl-CoA decarboxylase, a key enzyme in oxalate catabolism, was isolated by probing a genomic library of B. animalis subsp. lactis BI07, which was one of the most active strains in the preliminary screening. The genetic and transcriptional organization of oxc flanking regions was determined, unravelling the presence of other two independently transcribed open reading frames, potentially responsible for B. animalis subsp. lactis ability to degrade oxalate. Transcriptional analysis, using real-time quantitative reverse transcription PCR, revealed that these genes were highly induced in cells first adapted to subinhibitory concentrations of oxalate and then exposed to pH 4.5. Acidic conditions were also a prerequisite for a significant oxalate degradation rate, which dramatically increased in oxalate pre-adapted cells, as demonstrated in fermentation experiments with different pH-controlled batch cultures. These findings provide new insights in the characterization of oxalate-degrading probiotic bacteria and may support the use of B. animalis subsp. lactis as a promising adjunct for the prophylaxis and management of oxalate-related kidney disease. In order to provide some insight into the molecular mechanisms involved in the interaction with the host, in the second part of the job, we investigated whether Bifidobacterium was able to capture human plasminogen on the cell surface. The binding of human plasminogen to Bifidobacterium was dependent on lysine residues of surface protein receptors. By using a proteomic approach, we identified six putative plasminogen-binding proteins in the cell wall fraction of three strain of Bifidobacterium. The data suggest that plasminogen binding to Bifidobactrium is due to the concerted action of a number of proteins located on the bacterial cell surface, some of which are highly conserved cytoplasmic proteins which have other essential cellular functions. Our findings represent a step forward in understanding the mechanisms involved in the Bifidobacterium-host interaction. In these job w studied a new approach based on to MALDI-TOF MS to measure the interaction between entire bacterial cells and host molecular target. MALDI-TOF (Matrix Assisted Laser Desorption Ionization-Time of Flight)—mass spectrometry has been applied, for the first time, in the investigation of whole Bifidobacterium cells-host target proteins interaction. In particular, by means of this technique, a dose dependent human plasminogen-binding activity has been shown for Bifidobacterium. The involvement of lysine binding sites on the bacterial cell surface has been proved. The obtained result was found to be consistent with that from well-established standard methodologies, thus the proposed MALDI-TOF approach has the potential to enter as a fast alternative method in the field of biorecognition studies involving in bacterial cells and proteins of human origin.

Jasmonates and abscisic acid influence fruit ripening and plant water use: practical, physiological and morphological aspects

The aim of the present thesis was to better understand the physiological role of the phytohormones jasmonates (JAs) and abscisic acid (ABA) during fruit ripening in prospect of a possible field application of JAs and ABA to improve fruit yield and quality. In particular, the effects of exogenous application of these substances at different fruit developmental stages and under different experimental conditions were evaluated. Some aspects of the water relations upon ABA treatment were also analysed. Three fruit species, peach (Prunus persica L. Batsch), golden (Actinidia chinensis) and green kiwifruit (Actinidia deliciosa), and several of their cvs, were used for the trials. Different experimental models were adopted: fruits in planta, detached fruit, detached branches with fruit, girdled branches and micropropagated plants. The work was structured into four sets of experiments as follows: (i) Pre-harvest methyl jasmonate (MJ) application was performed at S3/S4 transition under field conditions in Redhaven peach; ethylene production, ripening index, fruit quality and shelf-life were assessed showing that MJ-treated fruit were firmer and thus less ripe than controls as confirmed by the Index of Absorbance Difference (IAD), but exhibited a shorter shelf-life due to an increase in ethylene production. Moreover, the time course of the expression of ethylene-, auxin- and other ripening-related genes was determined. Ripening-related ACO1 and ACS1 transcript accumulation was inhibited though transiently by MJ, and gene expression of the ethylene receptor ETR2 and of the ethylene-related transcription factor ERF2 was also altered. The time course of the expression of several auxin-related genes was strongly affected by MJ suggesting an increase in auxin biosynthesis, altered auxin conjugation and release as well as perception and transport; the need for a correct ethylene/auxin balance during ripening was confirmed. (ii) Pre- and post-harvest ABA applications were carried out under field conditions in Flaminia and O’Henry peach and Stark Red Gold nectarine fruit; ethylene production, ripening index, fruit quality and shelf-life were assessed. Results show that pre-harvest ABA applications increase fruit size and skin color intensity. Also post-harvest ABA treatments alter ripening-related parameters; in particular, while ethylene production is impaired in ABA-treated fruit soluble solids concentration (SSC) is enhanced. Following field ABA applications stem water potential was modified since ABA-treated peach trees retain more water. (iii) Pre- and post-harvest ABA and PDJ treatments were carried out in both kiwifruit species under field conditions at different fruit developmental stages and in post-harvest. Ripening index, fruit quality, plant transpiration, photosynthesis and stomatal conductance were assessed. Pre-harvest treatments enhance SSC in the two cvs and flesh color development in golden kiwifruit. Post-harvest applications of either ABA or ABA plus PDJ lead to increased SSC. In addition, ABA reduces gas exchanges in A. deliciosa. (iv) Spray, drench and dipping ABA treatments were performed in micropropagated peach plants and in peach and nectarine detached branches; plant water use and transpiration, biomass production and fruit dehydration were determined. In both plants and branches ABA significantly reduces water use and fruit dehydration. No negative effects on biomass production were detected. The present information, mainly arising from plant growth regulator application in a field environment, where plants have to cope with multiple biotic and abiotic stresses, may implement the perspectives for the use of these substances in the control of fruit ripening.

Involvement of the opioid system and BDNF in neuroplastic alterations occurring in neuropathic pain conditions

Chronic pain affects one in five adults, reducing quality of life and increasing risk of developing co-morbidities such as depression. Neuropathic pain results by lesions to the nervous system that alter its structure and function leading to spontaneous pain and amplified responses to noxious and innocuous stimuli. The Opioid System is probably the most important system involved in control of nociceptive transmission. Dynorphin and nociceptin systems have been suggested key mediators of some neuropathic pain aspects. An important role also for BDNF has been recently suggested since its involvement in the peripheral and central sensitization phenomena is known. We studied neuroplastic alterations occurring in chronic pain in mice subjected to the chronic constriction injury (CCI). We investigated gene expression alterations of both BDNF and Opioid System at spinal level at different intervals of time. A transient upregulation of pBDNF and pDYN was observed in spinal cord, while increasing upregulation of ppN/OFQ was found in the DRGs of injured mice. Development of neuropathic behavioral signs has been observed in ICR/CD-1 and BDNF+/+ mice, subjected to CCI. A different development of these signs was observed in BDNF+/-. We also studied gene expression changes of investigated systems in different brain areas fourteen days after surgery. We found pBDNF, pDYN, pKOP, ppN/OFQ and pNOP gene expression alterations in several areas of CCI mice. In the same brain regions we also determined bioactive nociceptin peptide levels, and elevated N/OFQ levels were observed in the amygdala area. Histone modifications studies have been performed in BDNF and DYN gene promoters of CCI animal spinal cord showing selected alterations in pDYN gene promoter. In addition, a preliminary characterization of the innovative NOP-EGFP mice was performed. Overall, our results could be useful to understand which and how neuropeptidergic systems are involved in neuroplastic mechanism occurring in neuropathic pain.