Design methods and geometric modelling of scaffolds for bone tissue engineering

Every year, thousand of surgical treatments are performed in order to fix up or completely substitute, where possible, organs or tissues affected by degenerative diseases. Patients with these kind of illnesses stay long times waiting for a donor that could replace, in a short time, the damaged organ or the tissue. The lack of biological alternates, related to conventional surgical treatments as autografts, allografts, e xenografts, led the researchers belonging to different areas to collaborate to find out innovative solutions. This research brought to a new discipline able to merge molecular biology, biomaterial, engineering, biomechanics and, recently, design and architecture knowledges. This discipline is named Tissue Engineering (TE) and it represents a step forward towards the substitutive or regenerative medicine. One of the major challenge of the TE is to design and develop, using a biomimetic approach, an artificial 3D anatomy scaffold, suitable for cells adhesion that are able to proliferate and differentiate themselves as consequence of the biological and biophysical stimulus offered by the specific tissue to be replaced. Nowadays, powerful instruments allow to perform analysis day by day more accurateand defined on patients that need more precise diagnosis and treatments.Starting from patient specific information provided by TC (Computed Tomography) microCT and MRI(Magnetic Resonance Imaging), an image-based approach can be performed in order to reconstruct the site to be replaced. With the aid of the recent Additive Manufacturing techniques that allow to print tridimensional objects with sub millimetric precision, it is now possible to practice an almost complete control of the parametrical characteristics of the scaffold: this is the way to achieve a correct cellular regeneration. In this work, we focalize the attention on a branch of TE known as Bone TE, whose the bone is main subject. Bone TE combines osteoconductive and morphological aspects of the scaffold, whose main properties are pore diameter, structure porosity and interconnectivity. The realization of the ideal values of these parameters represents the main goal of this work: here we'll a create simple and interactive biomimetic design process based on 3D CAD modeling and generative algorithmsthat provide a way to control the main properties and to create a structure morphologically similar to the cancellous bone. Two different typologies of scaffold will be compared: the first is based on Triply Periodic MinimalSurface (T.P.M.S.) whose basic crystalline geometries are nowadays used for Bone TE scaffolding; the second is based on using Voronoi's diagrams and they are more often used in the design of decorations and jewellery for their capacity to decompose and tasselate a volumetric space using an heterogeneous spatial distribution (often frequent in nature). In this work, we will show how to manipulate the main properties (pore diameter, structure porosity and interconnectivity) of the design TE oriented scaffolding using the implementation of generative algorithms: "bringing back the nature to the nature".

Hexose-6-phosphate dehydrogenase modulates 11beta-hydroxysteroid dehydrogenase type 1-dependent metabolism of 7-keto- and 7beta-hydroxy-neurosteroids

BACKGROUND: The role of 11beta-hydroxysteroid dehydrogenase type 1 (11beta-HSD1) in the regulation of energy metabolism and immune system by locally reactivating glucocorticoids has been extensively studied. Experiments determining initial rates of enzyme activity revealed that 11beta-HSD1 can catalyze both the reductase and the dehydrogenase reaction in cell lysates, whereas it predominantly catalyzes the reduction of cortisone to cortisol in intact cells that also express hexose-6-phosphate dehydrogenase (H6PDH), which provides cofactor NADPH. Besides its role in glucocorticoid metabolism, there is evidence that 11beta-HSD1 is involved in the metabolism of 7-keto- and 7-hydroxy-steroids; however the impact of H6PDH on this alternative function of 11beta-HSD1 has not been assessed. METHODOLOGY: We investigated the 11beta-HSD1-dependent metabolism of the neurosteroids 7-keto-, 7alpha-hydroxy- and 7beta-hydroxy-dehydroepiandrosterone (DHEA) and 7-keto- and 7beta-hydroxy-pregnenolone, respectively, in the absence or presence of H6PDH in intact cells. 3D-structural modeling was applied to study the binding of ligands in 11beta-HSD1. PRINCIPAL FINDINGS: We demonstrated that 11beta-HSD1 functions in a reversible way and efficiently catalyzed the interconversion of these 7-keto- and 7-hydroxy-neurosteroids in intact cells. In the presence of H6PDH, 11beta-HSD1 predominantly converted 7-keto-DHEA and 7-ketopregnenolone into their corresponding 7beta-hydroxy metabolites, indicating a role for H6PDH and 11beta-HSD1 in the local generation of 7beta-hydroxy-neurosteroids. 3D-structural modeling offered an explanation for the preferred formation of 7beta-hydroxy-neurosteroids. CONCLUSIONS: Our results from experiments determining the steady state concentrations of glucocorticoids or 7-oxygenated neurosteroids suggested that the equilibrium between cortisone and cortisol and between 7-keto- and 7-hydroxy-neurosteroids is regulated by 11beta-HSD1 and greatly depends on the coexpression with H6PDH. Thus, the impact of H6PDH on 11beta-HSD1 activity has to be considered for understanding both glucocorticoid and neurosteroid action in different tissues.

Study of particle movement in the nearshore region of Lake Superior with radioisotope tracers

Biogeochemical processes in the coastal region, including the coastal area of the Great Lakes, are of great importance due to the complex physical, chemical and biological characteristics that differ from those on either the adjoining land or open water systems. Particle-reactive radioisotopes, both naturally occurring (210Pb, 210Po and 7Be) and man-made (137Cs), have proven to be useful tracers for these processes in many systems. However, a systematic isotope study on the northwest coast of the Keweenaw Peninsula in Lake Superior has not yet been performed. In this dissertation research, field sampling, laboratory measurements and numerical modeling were conducted to understand the biogeochemistry of the radioisotope tracers and some particulate-related coastal processes. In the first part of the dissertation, radioisotope activities of 210Po and 210Pb in a variability of samples (dissolved, suspended particle, sediment trap materials, surficial sediment) were measured. A completed picture of the distribution and disequilibrium of this pair of isotopes was drawn. The application of a simple box model utilizing these field observations reveals short isotope residence times in the water column and a significant contribution of sediment resuspension (for both particles and isotopes). The results imply a highly dynamic coastal region. In the second part of this dissertation, this conclusion is examined further. Based on intensive sediment coring, the spatial distribution of isotope inventories (mainly 210Pb, 137Cs and 7Be) in the nearshore region was determined. Isotope-based focusing factors categorized most of the sampling sites as non- or temporary depositional zones. A twodimensional steady-state box-in-series model was developed and applied to individual transects with the 210Pb inventories as model input. The modeling framework included both water column and upper sediments down to the depth of unsupported 210Pb penetration. The model was used to predict isotope residence times and cross-margin fluxes of sediments and isotopes at different locations along each transect. The time scale for sediment focusing from the nearshore to offshore regions of the transect was on the order of 10 years. The possibility of sediment longshore movement was indicated by high inventory ratios of 137Cs: 210Pb. Local deposition of fine particles, including fresh organic carbon, may explain the observed distribution of benthic organisms such as Diporeia. In the last part of this dissertation, isotope tracers, 210Pb and 210Po, were coupled into a hydrodynamic model for Lake Superior. The model was modified from an existing 2-D finite difference physical-biological model which has previously been successfully applied on Lake Superior. Using the field results from part one of this dissertation as initial conditions, the model was used to predict the isotope distribution in the water column; reasonable results were achieved. The modeling experiments demonstrated the potential for using a hydrodynamic model to study radioisotope biogeochemistry in the lake, although further refinements are necessary.

Patient-specific finite-element simulation of the human cornea: A clinical validation study on cataract surgery

The planning of refractive surgical interventions is a challenging task. Numerical modeling has been proposed as a solution to support surgical intervention and predict the visual acuity, but validation on patient specific intervention is missing. The purpose of this study was to validate the numerical predictions of the post-operative corneal topography induced by the incisions required for cataract surgery. The corneal topography of 13 patients was assessed preoperatively and postoperatively (1-day and 30-day follow-up) with a Pentacam tomography device. The preoperatively acquired geometric corneal topography – anterior, posterior and pachymetry data – was used to build patient-specific finite element models. For each patient, the effects of the cataract incisions were simulated numerically and the resulting corneal surfaces were compared to the clinical postoperative measurements at one day and at 30-days follow up. Results showed that the model was able to reproduce experimental measurements with an error on the surgically induced sphere of 0.38D one day postoperatively and 0.19D 30 days postoperatively. The standard deviation of the surgically induced cylinder was 0.54D at the first postoperative day and 0.38D 30 days postoperatively. The prediction errors in surface elevation and curvature were below the topography measurement device accuracy of ±5μm and ±0.25D after the 30-day follow-up. The results showed that finite element simulations of corneal biomechanics are able to predict post cataract surgery within topography measurement device accuracy. We can conclude that the numerical simulation can become a valuable tool to plan corneal incisions in cataract surgery and other ophthalmosurgical procedures in order to optimize patients' refractive outcome and visual function.

Identifiability of diffusion and sorption parameters from in situ diffusion experiments by using simultaneously tracer dilution and claystone data

In situ diffusion experiments are performed in geological formations at underground research laboratories to overcome the limitations of laboratory diffusion experiments and investigate scale effects. Tracer concentrations are monitored at the injection interval during the experiment (dilution data) and measured from host rock samples around the injection interval at the end of the experiment (overcoring data). Diffusion and sorption parameters are derived from the inverse numerical modeling of the measured tracer data. The identifiability and the uncertainties of tritium and Na-22(+) diffusion and sorption parameters are studied here by synthetic experiments having the same characteristics as the in situ diffusion and retention (DR) experiment performed on Opalinus Clay. Contrary to previous identifiability analyses of in situ diffusion experiments, which used either dilution or overcoring data at approximate locations, our analysis of the parameter identifiability relies simultaneously on dilution and overcoring data, accounts for the actual position of the overcoring samples in the claystone, uses realistic values of the standard deviation of the measurement errors, relies on model identification criteria to select the most appropriate hypothesis about the existence of a borehole disturbed zone and addresses the effect of errors in the location of the sampling profiles. The simultaneous use of dilution and overcoring data provides accurate parameter estimates in the presence of measurement errors, allows the identification of the right hypothesis about the borehole disturbed zone and diminishes other model uncertainties such as those caused by errors in the volume of the circulation system and the effective diffusion coefficient of the filter. The proper interpretation of the experiment requires the right hypothesis about the borehole disturbed zone. A wrong assumption leads to large estimation errors. The use of model identification criteria helps in the selection of the best model. Small errors in the depth of the overcoring samples lead to large parameter estimation errors. Therefore, attention should be paid to minimize the errors in positioning the depth of the samples. The results of the identifiability analysis do not depend on the particular realization of random numbers. (C) 2012 Elsevier B.V. All rights reserved.

Upscaling of anisotropy in unsaturated Miller-similar porous media

Geological and pedological processes rarely form isotropic media as is usually assumed in transport studies. Anisotropy at the Darcy or field scale may be detected directly by measuring flow parameters or may become indirectly evident from movement and shape of solute plumes. Anisotropic behavior of a soil at one scale may, in many cases, be related to the presence of lower-scale directional structures. Miller similitude with different pore-scale geometries of the basic element is used to model macroscopic flow and transport behavior. Analytical expressions for the anisotropic conductivity tensor are derived based on the dynamic law that governs the flow problem at the pore scale. The effects of anisotropy on transport parameters are estimated by numerical modeling.

Comet 1P/Halley multifluid MHD model for the Giotto fly-by

The interaction of comets with the solar wind has been the focus of many studies including numerical modeling. We compare the results of our multifluid MHD simulation of comet 1P/Halley to data obtained during the flyby of the European Space Agency's Giotto spacecraft in 1986. The model solves the full set of MHD equations for the individual fluids representing the solar wind protons, the cometary light and heavy ions, and the electrons. The mass loading, charge-exchange, dissociative ion-electron recombination, and collisional interactions between the fluids are taken into account. The computational domain spans over several million kilometers, and the close vicinity of the comet is resolved to the details of the magnetic cavity. The model is validated by comparison to the corresponding Giotto observations obtained by the Ion Mass Spectrometer, the Neutral Mass Spectrometer, the Giotto magnetometer experiment, and the Johnstone Plasma Analyzer instrument. The model shows the formation of the bow shock, the ion pile-up, and the diamagnetic cavity and is able to reproduce the observed temperature differences between the pick-up ion populations and the solar wind protons. We give an overview of the global interaction of the comet with the solar wind and then show the effects of the Lorentz force interaction between the different plasma populations.

Redistribution of particles across the nucleus of comet 67P/Churyumov-Gerasimenko

Context. We present an investigation of the surface properties of areas on the nucleus of comet 67P/Churyumov-Gerasimenko. Aims. We aim to show that transport of material from one part of the cometary nucleus to another is a significant mechanism that influences the appearance of the nucleus and the surface thermal properties. Methods. We used data from the OSIRIS imaging system onboard the Rosetta spacecraft to identify surface features on the nucleus that can be produced by various transport mechanisms. We used simple calculations based on previous works to establish the plausibility of dust transport from one part of the nucleus to another. Results. We show by observation and modeling that "airfall" as a consequence of non-escaping large particles emitted from the neck region of the nucleus is a plausible explanation for the smooth thin deposits in the northern hemisphere of the nucleus. The consequences are also discussed. We also present observations of aeolian ripples and ventifacts. We show by numerical modeling that a type of saltation is plausible even under the rarified gas densities seen at the surface of the nucleus. However, interparticle cohesive forces present difficulties for this model, and an alternative mechanism for the initiation of reptation and creep may result from the airfall mechanism. The requirements on gas density and other parameters of this alternative make it a more attractive explanation for the observations. The uncertainties and implications are discussed.

Age determination and stable isotope record of foraminifera of sediment core SO82_5-2

Surface and deepwater paleoclimate records in Irminger Sea core SO82-5 (59°N, 31°W) and Icelandic Sea core PS2644 (68°N, 22°W) exhibit large fluctuations in thermohaline circulation (THC) from 60 to 18 calendar kyr B.P., with a dominant periodicity of 1460 years from 46 to 22 calendar kyr B.P., matching the Dansgaard-Oeschger (D-O) cycles in the Greenland Ice Sheet Project 2 (GISP2) temperature record [Grootes and Stuiver, 1997, doi:10.1029/97JC00880]. During interstadials, summer sea surface temperatures (SSTsu) in the Irminger Sea averaged to 8°C, and sea surface salinities (SSS) averaged to ~36.5, recording a strong Irminger Current and Atlantic THC. During stadials, SSTsu dropped to 2°-4°C, in phase with SSS drops by ~1-2. They reveal major meltwater injections along with the East Greenland Current, which turned off the North Atlantic deepwater convection and hence the heat advection to the north, in harmony with various ocean circulation and ice models. On the basis of the IRD composition, icebergs came from Iceland, east Greenland, and perhaps Svalbard and other northern ice sheets. However, the southward drifting icebergs were initially jammed in the Denmark Strait, reaching the Irminger Sea only with a lag of 155-195 years. We also conclude that the abrupt stadial terminations, the D-O warming events, were tied to iceberg melt via abundant seasonal sea ice and brine water formation in the meltwater-covered northwestern North Atlantic. In the 1/1460-year frequency band, benthic ?18O brine water spikes led the temperature maxima above Greenland and in the Irminger Sea by as little as 95 years. Thus abundant brine formation, which was induced by seasonal freezing of large parts of the northwestern Atlantic, may have finally entrained a current of warm surface water from the subtropics and thereby triggered the sudden reactivation of the THC. In summary, the internal dynamics of the east Greenland ice sheet may have formed the ultimate pacemaker of D-O cycles.

Geochemistry of lavas from the Lau-Tonga arc and back-arc system

Detailed comparison of mineralogy, and major and trace geochemistry are presented for the modern Lau Basin spreading centers, the Sites 834-839 lavas, the modern Tonga-Kermadec arc volcanics, the northern Tongan boninites, and the Lau Ridge volcanics. The data clearly confirm the variations from near normal mid-ocean-ridge basalt (N-MORB) chemistries (e.g., Site 834, Central Lau Spreading Center) to strongly arc-like (e.g., Site 839, Valu Fa), the latter closely comparable to the modern arc volcanoes. Sites 835 and 836 and the East Lau Spreading Center represent transitional chemistries. Bulk compositions range from andesitic to basaltic, but lavas from Sites 834 and 836 and the Central Lau Spreading Center extend toward more silica-undersaturated compositions. The Valu Fa and modern Tonga-Kermadec arc lavas, in contrast, are dominated by basaltic andesites. The phenocryst and groundmass mineralogies show the strong arc-like affinities of the Site 839 lavas, which are also characterized by the existence of very magnesian olivines (up to Fo90-92) and Cr-rich spinels in Units 3 and 6, and highly anorthitic plagioclases in Units 2 and 9. The regional patterns of mineralogical and geochemical variations are interpreted in terms of two competing processes affecting the inferred magma sources: (1) mantle depletion processes, caused by previous melt extractions linked to backarc magmatism, and (2) enrichment in large-ion-lithophile elements, caused by a subduction contribution. A general trend of increasing depletion is inferred both eastward across the Lau Basin toward the modern arc, and northward along the Tongan (and Kermadec) Arc. Numerical modeling suggests that multistage magma extraction can explain the low abundances (relative to N-MORB) of elements such as Nb, Ta, and Ti, known to be characteristic of island arc magmas. It is further suggested that a subduction jump following prolonged slab rollback could account for the initiation of the Lau Basin opening, plausibly allowing a later influx of new mantle, as required by the recognition of a two-stage opening of the Lau Basin.

Heat flow values in the Sea of Okhotsk region

Tectonic structure and anomalous distributions of geophysical fields of the Sea of Okhotsk region are considered; the lack of reliable data on age of the lithosphere beneath basins of various origin in the Sea of Okhotsk is noted. Model calculations based on geological and geophysical data yielded 65 Ma (Cretaceous-Paleocene boundary) age for the Central Okhotsk rise underlain by the continental lithosphere. This estimate agrees with the age (the end of Cretaceous) derived from seismostratigraphic data. A comparative analysis of theoretical and measured heat flows in the Akademii Nauk Rise, underlain by the thinned continental crust, is performed. The analysis points to a higher (by 20%) value of the measured thermal background of the rise, which is consistent with high negative gradient of gravity anomalies in this area. Calculations yielded 36 Ma (Early Oligocene) age and lithosphere thickness of 50 km for the South Okhotsk depression, whose seafloor was formed by processes of back-arc spreading. The estimated age of the depression is supported by kinematic data on the region; the calculated thickness of the lithosphere coincides with the value estimated from data of magnetotelluric sounding here. This indicates that formation time (36 Ma) of the South Okhotsk depression was estimated correctly. Numerical modeling performed for determination of the basement age of rifting basins in the Sea of Okhotsk gave the following estimates: 18 Ma (Early Miocene) for the Deryugin Basin, 12 Ma (Middle Miocene) for the TINRO Basin, and 23 Ma (Late Oligocene) for the West Kamchatka Trough. These estimates agree with formation time (Oligocene-Quaternary) of the sedimentary cover in rifting basins of the Sea of Okhotsk derived from geological and geophysical data. Model temperature estimates are obtained for lithologic and stratigraphic boundaries of the sedimentary cover in the Deryugin and TINRO Basins and the West Kamchatka Trough; the temperature analysis indicates that the latter two structures are promising for oil and hydrocarbon gas generation; the West Kamchatka Trough possesses better reservoir properties compared to the TINRO and Deryugin Basins. The latter is promising for generation of hydrocarbon gas. Paleogeodynamic reconstructions of the Sea of Okhotsk region evolution are obtained for times of 90, 66, and 36 Ma on the base of kinematic, geomagnetic, structural, tectonic, geothermal, and other geological and geophysical data.

Sr-Nd isotopic analyses of sand-sized sediment from the San Francisco Bay coastal system

A diverse suite of geochemical tracers, including 87Sr/86Sr and 143Nd/144Nd isotope ratios, the rare earth elements (REEs), and select trace elements were used to determine sand-sized sediment provenance and transport pathways within the San Francisco Bay coastal system. This study complements a large interdisciplinary effort (Barnard et al., 2012) that seeks to better understand recent geomorphic change in a highly urbanized and dynamic estuarine-coastal setting. Sand-sized sediment provenance in this geologically complex system is important to estuarine resource managers and was assessed by examining the geographic distribution of this suite of geochemical tracers from the primary sources (fluvial and rock) throughout the bay, adjacent coast, and beaches. Due to their intrinsic geochemical nature, 143Nd/144Nd isotopic ratios provide the most resolved picture of where sediment in this system is likely sourced and how it moves through this estuarine system into the Pacific Ocean. For example, Nd isotopes confirm that the predominant source of sand-sized sediment to Suisun Bay, San Pablo Bay, and Central Bay is the Sierra Nevada Batholith via the Sacramento River, with lesser contributions from the Napa and San Joaquin Rivers. Isotopic ratios also reveal hot-spots of local sediment accumulation, such as the basalt and chert deposits around the Golden Gate Bridge and the high magnetite deposits of Ocean Beach. Sand-sized sediment that exits San Francisco Bay accumulates on the ebb-tidal delta and is in part conveyed southward by long-shore currents. Broadly, the geochemical tracers reveal a complex story of multiple sediment sources, dynamic intra-bay sediment mixing and reworking, and eventual dilution and transport by energetic marine processes. Combined geochemical results provide information on sediment movement into and through San Francisco Bay and further our understanding of how sustained anthropogenic activities which limit sediment inputs to the system (e.g., dike and dam construction) as well as those which directly remove sediments from within the Bay, such as aggregate mining and dredging, can have long-lasting effects.

Methane concentration profiles and isotopic composition of sediment cores from the Black Sea

We present high resolution profiles for the methane concentration and the carbon isotope composition of methane from surface sediments and from the sediment-water transition in the Black Sea. At shallow water sites methane migrates from the sediment into the water column, and the magnitude of this upward migrating flux depends on the depth of the sulfate-methane transition (SMT) in the sediment. The isotope data reveal that the sediments at shallow water sites are a source for methane depleted in 13C relative to the isotope composition of methane in the water column. At deep water sites the methane concentration first decreases with depth in the sediment to reach lowest values at the Unit I to Unit II transition. Below this transition the concentration increases again. Numerical modeling of methane concentration and isotope data shows that high methane oxidation rates occur in the surface sediment layer, indicating that the removal of methane in the surface sediments is not related to the anaerobic oxidation of methane coupled to sulfate reduction that occurs a few meters deep in the sediment, at the SMT. Instead, near-surface methane consumption in the euxinic Black Sea sediments appears to be related to lithological stratification. Furthermore, a map of the diffusive methane fluxes in the Black Sea surface sediments indicates that approximately half of the Black Sea seafloor acts as a sink for methane and thus limits the flux of methane to the atmosphere.