The Behavior of Water in Collagen and Hydroxyapatite Sites of Cortical Bone: Fracture, Mechanical Wear, and Load Bearing Studies

Rapidly Screening the Correlation between the Rotational Mobility and the Hydrogen Bonding Strength of Confined Water

Automated Deuterium Relaxation-Ordered SpectroscopY in solution (ADROSYS), an automated two-dimensional deuterium NMR methodology, discriminates between D2O populations (as well as deuterium-labeled alcohol groups) whose properties differ as a result of being confined inside nanoscale volumes. In this contribution, a proof-of-principle demonstration on reverse micelles (RMs) yields the insight that as the length scale of the confinement decreases from several nanometers down to less than a nanometer, the position of the signal peak migrates through the two-dimensional (2D) spectrum, tracing out a distinctive path in the 2D space (of relaxation time vs chemical shift). The signals typically follow a relatively gentle linear path for water confined on the scale of several nanometers, before curving once the surfactants confine the water molecules to length scales smaller than 1-2 nm. The qualitative shape of this path, especially in the regime of strong confinement, can change with different choices of surfactants, i.e., a different choice of chemistry at the edges of the confining environment. An important facet of this research was to demonstrate the relatively wide applicability of these techniques by showing that both: (1) Standard modern NMR instrumentation is capable of deploying an automated measurement, even though the choice of a deuterium nucleus is nonstandard and frequently requires companion proton spectra in order to reference the chemical shifts; and (2) well-established (though underutilized) modern techniques can process the resulting signal even though it involves the somewhat unusual combination of chemical shifts along one dimension and a distribution of relaxation times along another dimension. In addition to demonstrating that this technique can be deployed across many samples of interest, detailed facts pertaining to the broadening or shifting of resulting signals upon inclusion of various guest molecules are also discussed.

Open-system force-elongation relationship of collagen in chemo-mechanical equilibrium with water

A significant deformation mechanism of collagen at low loads is molecular uncoiling and rearrangement. Although the effect of hydration and cross-linking has been investigated at larger loads when collagen undergoes molecular sliding, their effects on collagen molecular reorganization remain unclear. Here we develop two thermodynamic models that use the notion of open-system elasticity to elucidate the effect of swelling due to water uptake during deformation of collagen networks under low and high cross-linking conditions. With low crosslinking, entropic contributions dominate resulting in rejection of solvent from the polymer network leading to reduced collagen stiffness with increased loads. Contrarily, high cross-linking inhibits initial coiling and structural kinking and the mechanical behavior is dominated by elastic energy. In this configuration, the solvent content depends on the sign of the applied load resulting in a non-linear open-system stress-strain relationship. The models provide insight on the parameters that impact the stress-strain relationships of hydrated collagen and can inform the way collagenous matrices are treated both in medical and laboratory settings.

3D Tortuosity and Diffusion Characterization in the Human Mineralized Collagen Fibril Using a Random Walk Model

Bone tissue is mainly composed at the nanoscale of apatite minerals, collagen molecules and water that form the mineralized collagen fibril (MCF). In this work, we developed a 3D random walk model to investigate the influence of bone nanostructure on water diffusion. We computed 1000 random walk trajectories of water molecules within the MCF geometric model. An important parameter to analyse transport behaviour in porous media is tortuosity, computed as the ratio between the effective path length and the straight-line distance between initial and final points. The diffusion coefficient is determined from the linear fit of the mean squared displacement of water molecules as a function of time. To achieve more insight into the diffusion phenomenon within MCF, we estimated the tortuosity and diffusivity at different quotes in the longitudinal direction of the model. Tortuosity is characterized by increasing values in the longitudinal direction. As expected, the diffusion coefficient decreases as tortuosity increases. Diffusivity outcomes confirm the findings achieved by experimental investigations. The computational model provides insights into the relation between the MCF structure and mass transport behaviour that may contribute to the improvement of bone-mimicking scaffolds.

Exploring the microstructure of hydrated collagen hydrogels under scanning electron microscopy

Collagen hydrogels are a rapidly expanding platform in bioengineering and soft materials engineering for novel applications focused on medical therapeutics, medical devices and biosensors. Observations linking microstructure to material properties and function enables rational design strategies to control this space. Visualisation of the microscale organisation of these soft hydrated materials presents unique technical challenges due to the relationship between hydration and the molecular organisation of a collagen gel. Scanning electron microscopy is a robust tool widely employed to visualise and explore materials on the microscale. However, investigation of collagen gel microstructure is difficult without imparting structural changes during preparation and/or observation. Electrons are poorly propagated within liquid-phase materials, limiting the ability of electron microscopy to interrogate hydrated gels. Sample preparation techniques to remove water induce artefactual changes in material microstructure particularly in complex materials such as collagen, highlighting a critical need to develop robust material handling protocols for the imaging of collagen hydrogels. Here a collagen hydrogel is fabricated, and the gel state explored under high-vacuum (10^-6  Pa) and low-vacuum (80-120 Pa) conditions, and in an environmental SEM chamber. Visualisation of collagen fibres is found to be dependent on the degree of sample hydration, with higher imaging chamber pressures and humidity resulting in decreased feature fidelity. Reduction of imaging chamber pressure is used to induce evaporation of gel water content, revealing collagen fibres of significantly larger diameter than observed in samples dehydrated prior to imaging. Rapid freezing and cryogenic handling of the gel material is found to retain a porous 3D structure following sublimation of the gel water content. Comparative analysis of collagen hydrogel materials demonstrates the care needed when preparing hydrogel samples for electron microscopy.

Low-field and variable-field NMR relaxation studies of H2O and D2O molecular dynamics in articular cartilage

Osteoarthritis (OA) as the main degenerative disease of articular cartilage in joints is accompanied by structural and compositional changes in the tissue. Degeneration is a consequence of a reduction of the amount of macromolecules, the so-called proteoglycans, and of a corresponding increase in water content, both leading to structural weakening of cartilage. NMR investigations of cartilage generally address only the relaxation properties of water. In this study, two-dimensional (T1-T2) measurements of bovine articular cartilage samples were carried out for different stages of hydration, complemented by molecular exchange with D2O and treatment by trypsin which simulates degeneration by OA. Two signal components were identified in all measurements, characterized by very different T2 which suggests liquid-like and solid-like dynamics. These measurements allow the quantification of separate hydrogen components and their assignment to defined physical pools which had been discussed repeatedly in the literature, i.e. bulk-like water and a combination of protein hydrogens and strongly bound water. The first determination of 2H relaxation dispersion in comparison to 1H dispersion suggests intramolecular interactions as the dominating source for the pronounced magnetic field dependence of the longitudinal relaxation time T1.

Environment-controlled water adsorption at hydroxyapatite/collagen interfaces

Water contributes to the structure of bone by coupling hydroxyapatite to collagen over the hierarchical levels of tissue organization. Bone water exists in two states, bound or mobile, each accomplishing different roles. Although many experimental studies show that the amount of bound water correlates with bone strength, a molecular understanding of the interactions between hydroxyapatite, collagen and water is missing. In this work, we unveil the water adsorption properties of bone tissues at the nanoscale using advanced density functional theory methods. We demonstrate that environmental factors such as collagen conformation or degree of confinement, rather than the surface itself, dictate the adsorption mode, strength and density of water on hydroxyapatite. While the results derived in this paper come from a simplified model of bone tissues, they are consistent with experimental observations and constitute a reasonable starting point for more realistic models of bone tissues. For example, we show that environmental changes expected in aging bone lead to reduced water adsorption capabilities, which is consistent with weaker bones at the macroscale. Our findings provide a new interpretation of molecular interactions in bone tissues with the potential to impact bone repair strategies.

3D random walk model of diffusion in human Hypo- and Hyper- mineralized collagen fibrils

Bone tissue is composed at the nanoscale of apatite minerals, collagen molecules and water that form the mineralized collagen fibril (MCF). Water has a crucial role in bone biomineralization. We developed a 3D random walk model to investigate the water diffusion process within the MCF for three different scenarios, namely low, intermediate and high mineral volume fraction. The MCF geometric model is obtained after applying 6·10⁶ translational and rotational perturbations to an ordered arrangement of mineral. Subsequently, we compute 300 random trajectories of water molecules within the MCF for each mineral volume fraction. Every trajectory is constituted of up to 500 k positions of the water particle. We determined the diffusion coefficient from the linear fit of the mean squared displacement of water molecules as a function of time. We investigate changes in the diffusivity values in relation to variation of bone mineral content. The analysis performed on the random walk data, for all mineralization conditions, leads to diffusion coefficients in good agreement with the diffusivity outcomes achieved from previous experimental studies. Thus, the 3D geometrical configuration adopted in this numerical study appears suitable for modelling the MCF with different volume fractions, from hypo- to hyper-mineralized conditions. We observed that low mineral content is associated with an increase of the water diffusion, while lower values of diffusivity are determined in hypermineralized conditions. In agreement with experimental data, our results highlight the influence of the structural alterations on the mass transport properties.

A computational study of mechanical properties of collagen-based bio-composites

Studying changes in collagen deformation behavior at the nanoscale due to variations in mineralization and hydration is important for characterizing and developing collagen-based bio-composites. Recent studies also find that carbon nanotubes (CNTs) show promise as a reinforcing material for collagenous bio-composites. Currently, the effects of variation in mineral, water, and CNT content on collagen gap and overlap region mechanics during compression is unexplored. We use molecular dynamics simulations to investigate how variations in mineral, water, and CNT contents of collagen bio-composites in compression change their deformation behavior and thermal properties. Results indicate that variations in mineral and water content affect the collagen structure due to expansion or contraction of the gap and overlap regions. The deformation mechanisms of the gap and overlap regions also change. The presence of CNTs in non-mineralized collagen reduces the deformation of the gap region and increases the bio-composite elastic modulus to ranges comparable to mineralized collagen. The collagen/CNT bio-composites are also determined to have a higher specific heat than the studied mineralized collagen bio-composites, making them more likely to be resistant to thermal damage that could occur during implantation or functional use of a collagen collagen/CNT bio-composite biomaterial.

Technical Note: Post-burial alteration of bones: Quantitative characterization with solid-state 1H MAS NMR

The identification of markers of the modifications occurring in human bones after death and of the sedimentary and post-sedimentary processes affecting their state of preservation, is of interest for several scientific disciplines. A new index, obtained from spectral deconvolution of the 1H MAS NMR spectra of bones, relating the number of organic protons to the amount of hydrogen nuclei in the OH- groups of bioapatite, is proposed as indicator of the state of preservation of the organic fraction. In the osteological material from three different archaeological sites, this index resulted positively correlated with the extent of collagen loss derived from infrared spectroscopy. Its sensitivity to changes in the physical and chemical characteristics of bone allows to identify distinct diagenetic pathways specific to each site and to distinguish different trajectories within the same site.

Characterization of bony anatomic regions in pediatric and adult healthy volunteers using diffuse optical spectroscopic imaging

SIGNIFICANCE

Diffuse optical spectroscopic imaging (DOSI) measures quantitative functional and molecular information in thick tissue in a noninvasive manner using near-infrared light. DOSI may be useful for diagnosis and prognosis of bone pathologies including osteosarcoma and Ewing's sarcoma, but little is currently known about DOSI-derived parameters in bony anatomic locations where this disease occurs.

AIM

Our goal is to quantify the optical characteristics and chromophore content of bony anatomic locations of healthy volunteers and assess differences due to anatomic region, age, sex, ethnicity, race, and body fat.

APPROACH

Fifty-five healthy volunteers aged 4 to 72 were enrolled in the study. The optical properties and quantitative tissue concentrations of oxyhemoglobin, deoxyhemoglobin, water, and lipids were assessed at the distal humerus, distal femur, and proximal tibia. Body fat was assessed using skinfold calipers. One volunteer underwent a more comprehensive body scan from neck to foot to explore chromophore distributions within an individual. Regression analysis was used to identify the most important sources of variation in the measured data set.

RESULTS

Statistical differences between bony locations were found for scattering, water, and lipids, but not for hemoglobin. All chromophores had statistical differences with sex, but there were no significant age-related correlations. Regression analysis revealed that body fat measured with skinfold calipers was the most important predictor of oxy-, deoxy-, total hemoglobin, and lipids. Hemoglobin and lipid levels were highly correlated (ρ  ≥  0.7) over the subject population and within the single-subject body scan.

CONCLUSIONS

DOSI can successfully measure bony anatomic sites where osteosarcomas and Ewing's sarcomas commonly occur. Future studies of bone pathology using DOSI should account for the variation caused by anatomic region, sex, race and ethnicity, and body fat as these cause substantial variations in DOSI-derived metrics.

Water-lipid interactions in native bone by high-resolution solid-state NMR spectroscopy

The study of structural and dynamical properties of lipid and its associated interaction with different components of bone is essential to understand its role at a different level of bone homeostasis such as bone mineralization and bone metabolism. In this article, we present water-dependent dynamical changes observed in lipids (triglycerides) in its absolute native environment inside bone by high-resolution ¹H solid-state nuclear magnetic resonance spectroscopy (ssNMR). Relaxation measurement (T₂ measurement) ssNMR experiments were performed at different levels of water network induced by dehydration and H/D exchange in native bone. Our measurements reflect the changes in the local environment and dynamical properties of triglyceride due to different hydration levels. The present study explains the role of water in stabilizing the structural properties of triglycerides in bone hence will help understand its pathological role associated with bone physiology and bone disorders.

Static Solid Relaxation Ordered Spectroscopy: SS-ROSY

A two-dimensional pulse sequence is introduced for correlating nuclear magnetic resonance anisotropic chemical shifts to a relaxation time (e.g., T₁) in solids under static conditions. The sequence begins with a preparatory stage for measuring relaxation times, and is followed by a multiple pulse sequence for homonuclear dipolar decoupling. Data analysis involves the use of Fourier transform, followed by a one-dimensional inverse Laplace transform for each frequency index. Experimental results acquired on solid samples demonstrate the general approach, and additional variations involving heteronuclear decoupling and magic angle spinning are discussed.

The effect of NaOCl and heat treatment on static and dynamic mechanical properties and chemical changes of dentine

OBJECTIVES

To determine the effect of heat on flexural strength (FS), maximum strain (MS), storage modulus (SM), tan delta (TD) and chemical changes through micro-Raman spectroscopy of dentine exposed to 2.5% NaOCl or saline.

METHOD

ology: Dentine bars were randomly allocated to 8 test groups. Half (groups 2,4,6,8) were treated with NaOCl for 20 min; the rest (groups 1,3,5,7) remained in saline. FS/MS were measured in groups 1-4 (n = 15) (3/4 were also heated to 200 °C & re-hydrated in saline). Micro-Raman spectroscopy was performed on bars from groups 1-4. SM/TD were measured in 5-8: in 5/6 (n = 10), repeated after heating (200 °C), then following re-hydration; in 7/8 (n = 3) after heating to 25-185 °C.

RESULTS

Increase in MS on heat and FS/MS on heat + NaOCl was not significant (P > 0.05). SM increased (P = 0.06) after heat treatment but reduced to initial state after rehydration (P = 0.03). TD did not change (P = 0.4) after heat (200 °C) treatment but rehydration increased it compared with pre-treatment state (P = 0.001). For dentine bars pre-treated with NaOCl, SM did not change (P = 0.6) after heat (200 °C) treatment or rehydration but TD significantly increased (P = 0.02) upon re-hydration compared with pre- (P=0.007), or post- (P = 0.03) heat-treatment states. SM and TD varied between 25-185 °C with no consistent trend amongst the NaOCl pre-treated bars. Micro-Raman only detected chemical changes following NaOCl treatment in the mineral phase.

CONCLUSIONS

Exposure of dentine bars to heat and NaOCl produced only moderate changes to quasi-static but marked changes to viscoelastic properties, which may be explained by chemical alterations.

Alterations of elastin in female reproductive tissues arising from advancing parity

Female reproductive tissues undergo significant alterations during pregnancy, which may compromise the structural integrity of extracellular matrix proteins. Here, we report on modifications of elastic fibers, which are primarily composed of elastin and believed to provide a scaffold to the reproductive tissues, due to parity and parturition. Elastic fibers from the upper vaginal wall of virgin Sprague Dawley rats were investigated and compared to rats having undergone one, three, or more than five pregnancies. Optical microscopy was used to study fiber level changes. Mass spectrometry, ¹³C and ²H NMR, was applied to study alterations of elastin from the uterine horns. Spectrophotometry was used to measure matrix metalloproteinases-2,9 and tissue inhibitor of metalloproteinase-1 concentration changes in the uterine horns. Elastic fibers were found to exhibit increase in tortuosity and fragmentation with increased pregnancies. Surprisingly, secondary structure, dynamics, and crosslinking of elastin from multiparous cohorts appear similar to healthy mammalian tissues, despite fragmentation observed at the fiber level. In contrast, elastic fibers from virgin and single pregnancy cohorts are less fragmented and comprised of elastin exhibiting structure and dynamics distinguishable from multiparous groups, with reduced crosslinking. These alterations were correlated to matrix metalloproteinases-2,9 and tissue inhibitor of metalloproteinase-1 concentrations. This work indicates that fiber level alterations resulting from pregnancy and/or parturition, such as fragmentation, rather than secondary structure (e.g. elastin crosslinking density), appear to govern scaffolding characteristics in the female reproductive tissues.

Effects of hydration and mineralization on the deformation mechanisms of collagen fibrils in bone at the nanoscale

Bone is a biomaterial with a structural load-bearing function. Investigating the biomechanics of bone at the nanoscale is important in application to tissue engineering, the development of bioinspired materials, and for characterizing factors such as age, trauma, or disease. At the nanoscale, bone is composed of fibrils that are primarily a composite of collagen, apatite crystals (mineral), and water. Though several studies have been done characterizing the mechanics of fibrils, the effects of variation and distribution of water and mineral content in fibril gap and overlap regions are unexplored. We investigate how the deformation mechanisms of collagen fibrils change as a function of mineral and water content. We use molecular dynamics to study the mechanics of collagen fibrils of 0 wt%, 20 wt%, and 40 wt% mineralization and 0 wt%, 2 wt%, and 4 wt% hydration under applied tensile stresses. We observe that the stress-strain behavior becomes more nonlinear with an increase in hydration, and an increase in mineral content for hydrated fibrils under tensile stress reduces the nonlinear stress versus strain behavior caused by hydration. The Young's modulus of both non-mineralized and mineralized fibrils decreases as the water content increases. As the water content increases, the gap/overlap ratio increases by approximately 40% for the non-mineralized cases and 16% for the highly mineralized cases. Our results indicate that variations in mineral and water content change the distribution of water in collagen fibrils and that the water distribution changes the deformation of gap and overlap regions under tensile loading.

Non-destructive NIR spectral imaging assessment of bone water: Comparison to MRI measurements

Bone fracture risk increases with age, disease states, and with use of certain therapeutics, such as acid-suppressive drugs, steroids and high-dose bisphosphonates. Historically, investigations into factors that underlie bone fracture risk have focused on evaluation of bone mineral density (BMD). However, numerous studies have pointed to factors other than BMD that contribute to fragility, including changes in bone collagen and water. The goal of this study is to investigate the feasibility of using near infrared spectral imaging (NIRSI) to determine the spatial distribution and relative amount of water and organic components in whole cross-sections of bone, and to compare those results to those obtained using magnetic resonance imaging (MRI) methods. Cadaver human whole-section tibiae samples harvested from 18 donors of ages 27-97years underwent NIRSI and ultrashort echo time (UTE) MRI. As NIRSI data is comprised of broad absorbances, second derivative processing was evaluated as a means to narrow peaks and obtain compositional information. The (inverted) second derivative peak heights of the NIRSI absorbances correlated significantly with the mean peak integration of the water, collagen and fat NIR absorbances, respectively, indicating that either processing method could be used for compositional assessment. The 5797cm^-1 absorbance was validated as arising from the fat present in bone marrow, as it completely disappeared after ultrasonication. The MRI UTE-determined bound water content in tibial cortical bone samples ranged from 62 to 91%. The NIRSI water peaks at 5152cm^-1 and at 7008cm^-1 correlated significantly with the UTE data, with r=0.735, p=0.016, and r=0.71, p=0.0096, respectively. There was also a strong correlation between the intensity of the NIRSI water peak at 7008cm^-1 and the intensity of the collagen peak at 4608cm^-1 (r=0.69, p=0.004). Since NIRSI requires minimal to no sample preparation, this approach has great potential to become a gold standard modality for the investigation of changes in water content, distribution, and environment in pre-clinical studies of bone pathology and therapeutics.

Myocardial T1 mapping and extracellular volume quantification: an overview of technical and biological confounders

Novel tissue biomarkers based on the spin–lattice relaxation time T1, a fundamental property in the theory of magnetic resonance physics, have emerged as a new approach for myocardial tissue characterization with many validated clinical applications. This article is intended as an overview of the physical and physiological mechanisms underlying the interpretation and the accuracy of any practical measurement of T1, or derived biomarkers such as extravascular volume fraction, and also includes a discussion of potential pitfalls. Numerous caveats und knowledge gaps related to the precise interpretation of T1-based biomarkers remain, which are being addressed incrementally through ongoing research. Equally important, further careful standardization will pave the way for a wider clinical translation of these novel T1-based biomarkers of tissue remodeling, which have been well validated for their sensitivity to pathophysiological changes, though for the most part in single-center studies.

Biosynthetic PCL-graft-Collagen Bulk Material for Tissue Engineering Applications

Biosynthetic materials have emerged as one of the most exciting and productive fields in polymer chemistry due to their widespread adoption and potential applications in tissue engineering (TE) research. In this work, we report the synthesis of a poly(ε-caprolactone)-graft-collagen (PCL-g-Coll) copolymer. We combine its good mechanical and biodegradable PCL properties with the great biological properties of type I collagen as a functional material for TE. PCL, previously dissolved in dimethylformamide/dichloromethane mixture, and reacted with collagen using carbodiimide coupling chemistry. The synthesised material was characterised physically, chemically and biologically, using pure PCL and PCL/Coll blend samples as control. Infrared spectroscopy evidenced the presence of amide I and II peaks for the conjugated material. Similarly, XPS evidenced the presence of C-N and N-C=O bonds (8.96 ± 2.02% and 8.52 ± 0.63%; respectively) for PCL-g-Coll. Static contact angles showed a slight decrease in the conjugated sample. However, good biocompatibility and metabolic activity was obtained on PCL-g-Coll films compared to PCL and blend controls. After 3 days of culture, fibroblasts exhibited a spindle-like morphology, spreading homogeneously along the PCL-g-Coll film surface. We have engineered a functional biosynthetic polymer that can be processed by electrospinning.

Pregnancy-Induced Dynamical and Structural Changes of Reproductive Tract Collagen

The tissues and organs of the female reproductive tract and pelvic floor undergo significant remodeling and alterations to allow for fetal growth and birth. In this work, we report on a study of the alterations of murine reproductive tract collagen resulting from pregnancy and parturition by spectrophotometry, histology, and (13)C, (2)H nuclear magnetic resonance (NMR). Four different cohorts of rats were investigated that included virgin, multiparous, two- and fourteen-day postpartum primiparous rats. (13)C CPMAS NMR revealed small chemical shift differences across the cohorts. The measured H-C internuclear correlation times indicated differences in dynamics of some motifs. However, the dynamics of the major amino acids, e.g., Gly, remained unaltered with respect to parity. (2)H NMR relaxation measurements revealed an additional water reservoir in the postpartum and multiparous cohorts pointing to redistribution of water due to pregnancy and/or parturition. Spectrophotometric measurements indicated that the collagen content in virgin rats was highest. Histological analysis of the upper vaginal wall indicated a signature of collagen fiber dissociation with smooth muscle and a change in the density of collagen fibers in multiparous rats.

Separation of collagen-bound and porous bone water transverse relaxation in mice: proposal of a multi-step approach

The separation and quantification of collagen-bound water (CBW) and pore water (PW) components of the cortical bone signal are important because of their different contribution to bone mechanical properties. Ultrashort TE (UTE) imaging can be used to exploit the transverse relaxation from CBW and PW, allowing their quantification. We tested, for the first time, the feasibility of UTE measurements in mice for the separation and quantification of the transverse relaxation of CBW and PW in vivo using three different approaches for T2 * determination. UTE sequences were acquired at 4.7 T in six mice with 10 different TEs (50-5000 μs). The transverse relaxation time T2 * of CBW (T2 *cbw ) and PW (T2 *pw ) and the CBW fraction (bwf) were computed using a mono-exponential (i), a standard bi-exponential (ii) and a new multi-step bi-exponential (iii) approach. Regions of interest were drawn at multiple levels of the femur and vertebral body cortical bone for each mouse. The sum of the normalized squared residuals (Res) and the homogeneity of variance were tested to compare the different methods. In the femur, approach (i) yielded mean T2 * ± standard deviation (SD) of 657 ± 234 μs. With approach (ii), T2 *cbw , T2 *pw and bwf were 464 ± 153 μs, 15 777 ± 10 864 μs and 57.6 ± 9.9%, respectively. For approach (iii), T2 *cbw , T2 *pw and bwf were 387 ± 108 μs, 7534 ± 2765 μs and 42.5 ± 6.2%, respectively. Similar values were obtained from vertebral bodies. Res with approach (ii) was lower than with the two other approaches (p < 0.007), but T2 *pw and bwf variance was lower with approach (iii) than with approach (ii) (p < 0.048). We demonstrated that the separation and quantification of cortical bone water components with UTE sequences is feasible in vivo in mouse models. The direct bi-exponential approach exhibited the best approximation to the measured signal curve with the lowest residuals; however, the newly proposed multi-step algorithm resulted in substantially lower variability of the computed parameters. Copyright © 2016 John Wiley & Sons, Ltd.