Serum KL-6 and the mortality of patients with connective tissue disease-associated interstitial lung disease: A meta-analysis.

Connective tissue disease-associated interstitial lung disease (CTD-ILD) is an important underlying cause of morbidity and mortality in patients with CTD. Serum Krebs von den Lungen-6 (KL-6) is an immune factor which has been related to the severity of ILD. This systematic review and meta-analysis aimed to evaluate the association between serum KL-6 and mortality of patients with CTD-ILD. Longitudinal studies relevant to the aim of the meta-analysis were retrieved by search of electronic databases including PubMed, Web of Science, and Embase. A random-effects model was used to combine the data by incorporating the influence of between-study heterogeneity. Fifteen cohorts involving 1737 patients with CTD-ILD were included. During a mean follow-up of 35.3 months, 430 (24.8%) patients died. Compared to those with a lower KL-6 at admission, patients with a higher KL-6 were associated with a higher mortality risk during follow-up (risk ratio: 2.18, 95% confidence interval: 1.66 to 2.87, P < 0.001; I2 = 20%). Subgroup analysis showed a significant association in studies from Asian countries, but not in those from non-Asian countries; in studies with cutoff of KL-6 derived in receiver operating characteristic (ROC) curve analysis, but not in those derived from other methods; in studies with multivariate analysis, but not in those with univariate analysis (P for subgroup difference all < 0.05). The association was not significantly affected by different CTDs or methods for measuring serum KL-6. In conclusion, a high serum KL-6 may be a risk factor of increased mortality in patients with CTD-ILD.

Assessing osteoporosis in postmenopausal women: preliminary results using a novel lumbar spine phantom-based MRI scoring method.

OBJECTIVE: To develop a novel magnetic resonance imaging (MRI) phantom for producing F-score (for fat) and W-score (for water) and to evaluate the performance of these scores in assessing osteoporosis and related vertebral fractures. MATERIALS AND METHODS: First, a real-time phantom consisting of oil and water tubes was manufactured. Then, 30 female volunteers (age: 62.3 ± 6.3 years) underwent lumbar spine examination with MRI (using a novel phantom) and dual-energy X-ray absorptiometry (DXA), following ethical approval. MRI phantom-based F-score and W-score were defined by normalizing the vertebral signal intensities (SIs) by the oil and water SIs of the phantom on T1- and T2-weighted images, respectively. The diagnostic performances of the new scores for assessing osteoporosis and vertebral fractures were examined using receiver operating characteristic analysis and compared with DXA-measured areal bone mineral density (DXA-aBMD). RESULTS: The F-score and W-score were greater in the osteoporotic patients (3.93 and 2.29) than the non-osteoporotic subjects (3.05 and 1.79) and achieved AUC values of 0.85 and 0.74 (p < 0.05), respectively, when detecting osteoporosis. Similarly, F-score and W-score had greater values for the fracture patients (3.94 and 2.53) than the non-fracture subjects (3.14 and 1.69) and produced better AUC values (0.90 for W-score and 0.79 for F-score) compared to DXA-aBMD (AUC: 0.27, p < 0.05). In addition, the F-score and W-score had a strong correlation (r = 0.77; p < 0.001). CONCLUSION: A novel real-time lumber spine MRI phantom was developed, based upon which newly defined F-score and W-score were able to detect osteoporosis and demonstrated an improved ability over DXA-aBMD in differentiating patients with vertebral fractures.

By-degree Health and Economic Impacts of Lyme Disease, Eastern and Midwestern United States.

Lyme disease (LD) is the most common vector-borne disease in the United States (U.S.). This paper assesses how climate change may influence LD incidence in the eastern and upper Midwestern U.S. and the associated economic burden. We estimated future Ixodes scapularis habitat suitability and LD incidence with a by-degree approach using variables from an ensemble of multiple climate models. We then applied estimates for present-day and projected habitat suitability for I. scapularis, present-day presence of Borrelia burgdorferi, and projected climatological variables to model reported LD incidence at the county level among adults, children, and the total population. Finally, we applied an estimate of healthcare expenses to project economic impacts. We show an overall increase in LD cases with regional variation. We estimate an increase in incidence in New England and the upper Midwestern U.S. and a concurrent decrease in incidence in Virginia and North Carolina. At 3°C of national warming from the 1986-2015 baseline climate, we project approximately 55,000 LD cases, a 38-percent increase from present-day estimates. At 6°C of warming, our most extreme scenario, we project approximately 92,000 LD cases in the region, an increase of 145 percent relative to current levels. Annual LD-related healthcare expenses at 3°C of warming are estimated to be $236 million (2021 dollars), approximately 38 percent greater than present-day. These results may inform decision-makers tasked with addressing climate risks, the public, and healthcare professionals preparing for treatment and prevention of LD.

Non-invasive fractional flow reserve derived from reduced-order coronary model and machine learning prediction of stenosis flow resistance.

BACKGROUND AND OBJECTIVE: Recently, computational fluid dynamics enables the non-invasive calculation of fractional flow reserve (FFR) based on 3D coronary model, but it is time-consuming. Currently, machine learning technique has emerged as an efficient and reliable approach for prediction, which allows saving a lot of analysis time. This study aimed at developing a simplified FFR prediction model for rapid and accurate assessment of functional significance of stenosis. METHODS: A reduced-order lumped parameter model (LPM) of coronary system and cardiovascular system was constructed for rapidly simulating coronary flow, in which a machine learning model was embedded for accurately predicting stenosis flow resistance at a given flow from anatomical features of stenosis. Importantly, the LPM was personalized in both structures and parameters according to coronary geometries from computed tomography angiography and physiological measurements such as blood pressure and cardiac output for personalized simulations of coronary pressure and flow. Coronary lesions with invasive FFR ≤ 0.80 were defined as hemodynamically significant. RESULTS: A total of 91 patients (93 lesions) who underwent invasive FFR were involved in FFR derived from machine learning (FFRML) calculation. Of the 93 lesions, 27 lesions (29.0%) showed lesion-specific ischemia. The average time of FFRML simulation was about 10 min. On a per-vessel basis, the FFRML and FFR were significantly correlated (r = 0.86, p < 0.001). The diagnostic accuracy, sensitivity, specificity, positive predictive value and negative predictive value were 91.4%, 92.6%, 90.9%, 80.6% and 96.8%, respectively. The area under the receiver-operating characteristic curve of FFRML was 0.984. CONCLUSION: In this selected cohort of patients, the FFRML improves the computational efficiency and ensures the accuracy. The favorable performance of FFRML approach greatly facilitates its potential application in detecting hemodynamically significant coronary stenosis in future routine clinical practice.

A simplified coronary model for diagnosis of ischemia-causing coronary stenosis.

BACKGROUND AND OBJECTIVE: The functional assessment of the severity of coronary stenosis from coronary computed tomography angiography (CCTA)-derived fractional flow reserve (FFR) has recently attracted interest. However, existing algorithms run at high computational cost. Therefore, this study proposes a fast calculation method of FFR for the diagnosis of ischemia-causing coronary stenosis. METHODS: We combined CCTA and machine learning to develop a simplified single-vessel coronary model for rapid calculation of FFR. First, a zero-dimensional model of single-vessel coronary was established based on CCTA, and microcirculation resistance was determined through the relationship between coronary pressure and flow. In addition, a coronary stenosis model based on machine learning was introduced to determine stenosis resistance. Computational FFR (cFFR) was then obtained by combining the zero-dimensional model and the stenosis model with inlet boundary conditions for resting (cFFRr) and hyperemic (cFFRh) aortic pressure, respectively. We retrospectively analyzed 75 patients who underwent clinically invasive FFR (iFFR), and verified the model accuracy by comparison of cFFR with iFFR. RESULTS: The average computing time of cFFR was less than 2 s. The correlations between cFFRr and cFFRh with iFFR were r = 0.89 (p < 0.001) and r = 0.90 (p < 0.001), respectively. Diagnostic accuracy, sensitivity, specificity, positive predictive value, negative predictive value, positive likelihood ratio, negative likelihood ratio for cFFRr and cFFRh were 90.7%, 95.0%, 89.1%, 76.0%, 98.0%, 8.7, 0.1 and 92.0%, 95.0%, 90.9%, 79.2%, 98.0%, 10.5, 0.1, respectively. CONCLUSIONS: The proposed model enables rapid prediction of cFFR and exhibits high diagnostic performance in selected patient cohorts. The model thus provides an accurate and time-efficient computational tool to detect ischemia-causing stenosis and assist with clinical decision-making.

Biomechanical Study and Analysis for Cardiovascular/Skeletal Materials and Devices.

The Special Issue entitled "Biomechanical Study and Analysis for Cardiovascular/Skeletal Materials and Devices" addresses biological functional materials and devices relevant to cardiovascular diseases and orthopedic conditions [...].

Craniocaudal toggling increases the risk of screw loosening in osteoporotic vertebrae.

BACKGROUND AND OBJECTIVE: Screw loosening remains a prominent problem for osteoporotic patients undergoing pedicle screw fixation surgeries but its underlying mechanisms are not fully understood. This study sought to examine the interactive effect of craniocaudal or axial cyclic loading (toggling) and osteoporosis on screw fixation. METHODS: QCT-based finite element models of normal (n = 7; vBMD = 156 ± 13 mg/cm3) and osteoporotic vertebrae (n = 7; vBMD = 72 ± 6 mg/cm3) were inserted with pedicle screws and loaded with or without craniocaudal toggling. Among them, a representative normal vertebra (age: 55; BMD: 140 mg/cm3) and an osteoporotic vertebra (age: 64; BMD: 79 mg/cm3) were also loaded with or without axial toggling. The individual and interactive effects of craniocaudal toggling and osteoporosis on screw fixation strength (the force when the pull-up displacement of the screw head reached 1 mm) and bone tissue failure (characterized by equivalent plastic strain) were examined by repeated measure ANOVA. RESULTS: A significant interactive effect between craniocaudal toggling and osteoporosis on screw fixation strength was detected (p = 0.008). Specifically, craniocaudal toggling led to a marked decrease in the fixation strength (68%, p < 0.05) and stiffness (83%, p < 0.05) only in the osteoporotic vertebrae but had no effect on screw fixation strength and stiffness of the normal vertebrae (p > 0.05). Likewise, most of the bone tissues around the screw in the osteoporotic vertebrae yielded following craniocaudal toggling whereas this result was not seen in the normal vertebrae. The axial toggling had no significant effect on bone tissue failure as well as pedicle screw fixation in normal or osteoporotic vertebrae. CONCLUSIONS: Craniocaudal toggling substantially reduces the screw fixation strength of the osteoporotic vertebrae by progressively increasing tissue failure around the screw, and therefore may contribute to the higher rates of screw loosening in osteoporotic compared to normal patients, whereas axial toggling is not a risk factor for pedicle screw loosening in normal or osteoporotic patients.

Effects of dynamization timing and degree on bone healing of different fracture types.

Dynamization, that is, increasing interfragmentary movement (IFM) by reducing fixation stiffness from a rigid to a more flexible state, has been successfully used in clinical practice to promote fracture healing. However, it remains unclear how dynamization timing and degree affect bone healing of different fracture types. Finite element models of tibial fractures based on the OTA/AO classification (Simple: A1-Spiral, A2-Oblique, A3-Transverse; Wedge: B2-Spiral, B3-Fragmented; Complex: C2-Segment, C3-Irregular), in combination with fuzzy logic-based mechano-regulatory tissue differentiation algorithms, were used to simulate the healing process when dynamization of varied degrees (dynamization coefficient or DC = 0-0.9; 0.9 represents 90% reduction in the fixation stiffness relative to a rigid fixation) were applied at different time points after fracture. The fuzzy logic-based algorithms have been validated with a preclinical animal model. The results showed that the healing responses of type A fractures were more sensitive to the changes in dynamization degree and timing comparing with type B or C fractures. Additionally, the optimal dynamization regime for each fracture type was different. For type A fractures, a moderate dynamization degree (e.g., DC = 0.5) applied after Week 1 promoted the recovery of biomechanical integrity. For type B and C fractures, the effective dynamization included a greater dynamization degree (DC = 0.7) applied after Week 2. Our results further demonstrated that the fracture morphology affected interfragmentary strain environments within the callus, leading to varied healing results for different fracture types. These results suggest that the effects of dynamization are highly dependent of the fracture types. Therefore, specific dynamization strategies should be chosen for different fracture types to achieve optimal healing outcomes.

Biomechanical MRI detects reduced bone strength in subjects with vertebral fractures.

Vertebral fracture is one of the most serious consequences of osteoporosis. Estimation of vertebral strength from magnetic resonance imaging (MRI) scans may provide a new approach for the prediction of vertebral fractures. To that end, we sought to establish a biomechanical MRI (BMRI) method to compute vertebral strength and test its ability to distinguish fracture from non-fracture subjects. This case-control study included 30 subjects without vertebral fractures and 15 subjects with vertebral fractures. All subjects underwent MRI with a mDIXON-Quant sequence and quantitative computed tomography (QCT), from which proton fat fraction-based bone marrow adipose tissue (BMAT) content and volumetric bone mineral density (vBMD) were measured, respectively. Nonlinear finite element analysis was applied to MRI and QCT scans of L2 vertebrae to compute vertebral strength (BMRI- and BCT-strength). The differences in BMAT content, vBMD, BMRI-strength and BCT-strength between the two groups were examined by t-tests. Receiver operating characteristic (ROC) analysis was performed to assess the ability of each measured parameter to distinguish fracture from non-fracture subjects. Results showed that the fracture group had 23 % lower BMRI-strength (P < .001) and 19 % higher BMAT content (P < .001) than the non-fracture group, whereas no significant difference in vBMD was detected between the two groups. A poor correlation was found between vBMD and BMRI-strength (R2 = 0.33). Compared to vBMD and BMAT content, BMRI- and BCT-strength had the larger area under the curve (0.82 and 0.84, respectively) and provided better sensitivity and specificity in separating fracture from non-fracture subjects. In conclusion, BMRI is capable of detecting reduced bone strength in patients with vertebral fracture, and may serve as a new approach for risk assessment of vertebral fracture.

Osteocyte Estrogen Receptor ß (Ot-ERß) Regulates Bone Turnover and Skeletal Adaptive Response to Mechanical Loading Differently in Male and Female Growing and Adult Mice.

Age-related bone loss is a failure of balanced bone turnover and diminished skeletal mechanoadaptation. Estrogen receptors, ERα and ERß, play critical roles in osteoprotective regulation activated by estrogen and mechanical signals. Previous studies mainly focused on ERα and showed that osteocyte-ERα (Ot-ERα) regulated trabecular, but not cortical bone, and played a minor role in load-induced cortical adaptation. However, the role of Ot-ERß in bone mass regulation remains unrevealed. To address this issue, we characterized bone (re)modeling and gene expression in male and female mice with Ot-ERß deletion (ERß-dOT) and littermate control (LC) at 10 weeks (young) or 28 weeks (adult) of age, as well as their responses to in vivo tibial compressive loading. Increased cancellous bone mass appeared in the L4 vertebral body of young male ERß-dOT mice. At the same time, femoral cortical bone gene expression showed signs consistent with elevated osteoblast and osteoclast activities (type-I collagen, Cat K, RANKL). Upregulated androgen receptor (AR) expression was observed in young male ERß-dOT mice relative to LC, suggesting a compensatory effect of testosterone on male bone protection. In contrast, bone mass in L4 decreased in adult male ERß-dOT mice, attributed to potentially increased bone resorption activity (Cat K) with no change in bone formation. There was no effect of ERß-dOT on bone mass or gene expression in female mice. Sex-dependent regulation of Ot-ERß also appeared in load-induced cortical responsiveness. Young female ERß-dOT mice showed an enhanced tibial cortical anabolic adaptation compared with LC. In contrast, an attenuated cortical anabolic response presented at the proximal tibia in male ERß-dOT mice at both ages. For the first time, our findings suggest that Ot-ERß regulates bone (re)modeling and the response to mechanical signals through different mechanisms in males and females. © 2022 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

A Novel Screw Modeling Approach to Study the Effects of Screw Parameters on Pullout Strength.

Screw loosening remains a prominent problem for osteoporotic patients undergoing pedicle screw fixation surgeries and is affected by screw parameters (e.g., diameter, pitch, and thread angle). However, the individual and interactive effects of these parameters on screw fixation are not fully understood. Furthermore, the current finite element modeling of a threaded screw is less computationally efficient. To address these issues, we (1) explored a novel "simulated threaded screw" approach (virtual threads assigned to the contact elements of a simplified screw) and compared its performance with threaded and simplified screws, and (2) examined this approach the individual and interactive effects of altering screw diameter (5.5-6.5 mm), pitch (1-2 mm) and half-thread angle (20-30 deg) on pullout strength of normal vertebrae. Results demonstrated that the "simulated threaded screw" approach equivalently predicted pullout strength compared to the "threaded screw" approach (R2 = 0.99, slope = 1). We further found that the pullout strength was most sensitive to the change in screw diameter, followed by thread angle, pitch, and interactions of diameter*pitch or diameter*angle. In conclusion, the "simulated threaded screw" approach can achieve the same predictive capability compared to threaded modeling of the screw. The current findings may serve as useful references for planning of screw parameters, so as to improve the complication of screw loosening.

Biomechanical CT-computed bone strength predicts the risk of subsequent vertebral fracture.

Following primary fractures and percutaneous kyphoplasty (PKP), patients have a high risk of incurring a subsequent vertebral fracture (SVF). Given that SVF is a consequence of mechanical deterioration of the vertebra, we sought to examine whether vertebral strength derived from QCT-based finite element analysis (i.e., BCT) can predict the risk of SVF. Sixty-six patients who underwent PKP were categorized into two groups: control or non-SVF group (age: 70 ± 7 years; n = 40) and SVF group (age: 69 ± 8 years; n = 26). BCT was performed on L4 or L3 vertebrae to noninvasively measure vertebral strength. Vertebral strength was also estimated based upon the geometry and material properties of the vertebra. Additionally, trabecular volumetric bone mineral density (vBMD) and L1 Hounsfield unit (HU) were measured. t-Test, χ2 test or Mann Whitney U test were used to compare differences in these parameters between the two groups. The predictive abilities of BCT strength and other measured parameters were evaluated using the receiver operating characteristic (ROC) analysis. Results showed no significant difference in either vBMD or L1 HU between the control and SVF groups (p > 0.05), whereas BCT-computed and estimated vertebral strength values were significantly reduced by 33 % and 24 % for the SVF group relative to the non-SVF group, respectively. ROC curve indicated that BCT strength had the largest area under the curve, compared to other parameters. These results suggest that BCT-computed vertebral strength may serve as a surrogate for assessing risk of SVF.

Effect of the accordion technique on bone regeneration during distraction osteogenesis: A computational study.

BACKGROUND AND OBJECTIVE: Distraction osteogenesis (DO), a bone lengthening technique, is widely employed to treat congenital and acquired limb length discrepancies and large segmental bone defects. However, a major issue of DO is the prolonged consolidation phase (10-36 months) during which patients must wear a cumbersome external fixator. Attempts have been made to accelerate the healing process of DO by an alternating distraction and compression mode (so-called "accordion" technique or AT). However, it remains unclear how varied AT parameters affect DO outcomes and what the most effective AT mode is. METHODS: Based on an experimentally-verified mechanobiological model, we performed a parametric analysis via in silico simulation of the bone regeneration process of DO under different AT modes, including combinations of varied application times (AT began at week 1-8 of the consolidation phase), durations (AT was used continuously for 1 week, 2 weeks or 4 weeks) and rates (distraction or compression at 0.25, 0.5, 0.75, and 1 mm/12 h). The control group had no AT applied during the consolidation phase. RESULTS: Compared with the control group (no AT), AT applied at an early consolidation stage (e.g. week 1 of the consolidation phase) significantly enhanced bone formation and reduced the overall healing time. However, the effect of AT on bone healing was dependent on its duration and rate. Specifically, a moderate rate of AT (e.g. 0.5 mm/12 h) lasting for two weeks promoted blood perfusion recovery and bone regeneration, ultimately shortening the healing time. Conversely, over-high rates (e.g. 1 mm/12 h) and longer durations (e.g. 4 weeks) of AT adversely affected bone regeneration and blood perfusion recovery, thereby delaying bone bridging. CONCLUSIONS: These results suggest that the therapeutic effects of AT on DO are highly dependent of the AT parameters of choice. Under appropriate durations and rates, the AT applied at an early consolidation phase is beneficial for blood recovery and bone regeneration. These results may provide a basis for selecting effective AT modes to accelerate consolidation and reduce the overall treatment period of DO.

Mechanical testing and biomechanical CT analysis to assess vertebral flexion strength of Chinese cadavers.

Biomechanical CT (BCT), i.e., quantitative computed tomography-based finite element analysis (QCT-FEA), promises an improved technique over bone mineral density (BMD) in predicting bone strength and the risk of osteoporotic vertebral fractures. However, most of the BCT models only consider a uniform compressive loading condition and they have not been validated for Chinese subjects. This study examined the ability of BCT to predict wedge fracture-related vertebral flexion strength in a cohort of Chinese cadaveric vertebrae. Twelve human vertebrae were scanned with dual energy X-ray absorptiometry (DXA) and QCT to measure areal and volumetric BMD, respectively. To produce wedge fractures, the cadaveric vertebrae were experimentally loaded until failure under a 15° flexion. Vertebral flexion stiffness and strength were measured from the force-displacement curve. Voxel-based heterogeneous FE models of the vertebrae were created and virtually tested in uniform compression and 15° flexion to compute compressive and flexion strength (and stiffness), respectively. The predictions of vertebral flexion strength with BMD or BCT measures were evaluated with linear regression analyses. Results showed weak correlations between experimentally-measured flexion strength vs. DXA-aBMD (R2 = 0.26) or QCT-vBMD (R2 = 0.39). However, there were strong correlations between experimentally-measured flexion strength vs. BCT-computed vertebral strength under either flexion (R2 = 0.71) or compression (R2 = 0.70) loading conditions, although flexion reduced the BCT-computed vertebral strength by 9.2%. These results suggest that, regardless of whether a uniform compression or a flexion loading is simulated, BCT can predict in vitro vertebral flexion strength better than BMD.

Zinc-Based Biodegradable Materials for Orthopaedic Internal Fixation.

Traditional inert materials used in internal fixation have caused many complications and generally require removal with secondary surgeries. Biodegradable materials, such as magnesium (Mg)-, iron (Fe)- and zinc (Zn)-based alloys, open up a new pathway to address those issues. During the last decades, Mg-based alloys have attracted much attention by researchers. However, the issues with an over-fast degradation rate and release of hydrogen still need to be overcome. Zn alloys have comparable mechanical properties with traditional metal materials, e.g., titanium (Ti), and have a moderate degradation rate, potentially serving as a good candidate for internal fixation materials, especially at load-bearing sites of the skeleton. Emerging Zn-based alloys and composites have been developed in recent years and in vitro and in vivo studies have been performed to explore their biodegradability, mechanical property, and biocompatibility in order to move towards the ultimate goal of clinical application in fracture fixation. This article seeks to offer a review of related research progress on Zn-based biodegradable materials, which may provide a useful reference for future studies on Zn-based biodegradable materials targeting applications in orthopedic internal fixation.

Physical and Chemical Characterization of Biomineralized Collagen with Different Microstructures.

Mineralized collagen is the basic unit in hierarchically organized natural bone with different structures. Polyacrylic acid (PAA) and periodic fluid shear stress (FSS) are the most common chemical and physical means to induce intrafibrillar mineralization. In the present study, non-mineralized collagen, extrafibrillar mineralized (EM) collagen, intrafibrillar mineralized (IM) collagen, and hierarchical intrafibrillar mineralized (HIM) collagen induced by PAA and FSS were prepared, respectively. The physical and chemical properties of these mineralized collagens with different microstructures were systematically investigated afterwards. Transmission electron microscopy (TEM) and scanning electron microscopy (SEM) showed that mineralized collagen with different microstructures was prepared successfully. The pore density of the mineralized collagen scaffold is higher under the action of periodic FSS. Fourier transform infrared spectroscopy (FTIR) analysis showed the formation of the hydroxyapatite (HA) crystal. A significant improvement in the pore density, hydrophilicity, enzymatic stability, and thermal stability of the mineralized collagen indicated that the IM collagen under the action of periodic FSS was beneficial for maintaining collagen activity. HIM collagen fibers, which are prepared under the co-action of periodic FSS and sodium tripolyphosphate (TPP), may pave the way for new bone substitute material applications.

Interactive effects of various loading parameters on the fluid dynamics within the lacunar-canalicular system for a single osteocyte.

The osteocyte lacunar-canalicular system (LCS) serves as a mechanotransductive core where external loading applied to the skeleton is transduced into mechanical signals (e.g., fluid shear) that can be sensed by mechanosensors (osteocytes). The fluid velocity and shear stress within the LCS are affected by various loading parameters. However, the interactive effect of distinct loading parameters on the velocity and shear stress in the LCS remains unclear. To address this issue, we developed a multiscale modeling approach, combining a poroelastic finite element (FE) model with a single osteocytic LCS unit model to calculate the flow velocity and shear stress within the LCS. Next, a sensitivity analysis was performed to investigate individual and interactive effects of strain magnitude, strain rate, number of cycles, and intervening short rests between loading cycles on the velocity and shear stress around the osteocyte. Lastly, we developed a relatively simple regression model to predict those outcomes. Our results demonstrated that the strain magnitude or rate alone were the main factors affecting the velocity and shear stress; however, the combination of these two was not directly additive, and addition of a short rest between cycles could enhance the combination of these two related factors. These results show highly interactive effects of distinct loading parameters on fluid velocity and shear stress in the LCS. Specifically, our results suggest that an enhanced fluid dynamics environment in the LCS can be achieved with a brief number of load cycles combined with short rest insertion and high strain magnitude and rate.

In vivo and in silico monitoring bone regeneration during distraction osteogenesis of the mouse femur.

BACKGROUND AND OBJECTIVE: Distraction osteogenesis (DO) is a mechanobiological process of producing new bone by gradual and controlled distraction of the surgically separated bone segments. Mice have been increasingly used to study the role of relevant biological factors in regulating bone regeneration during DO. However, there remains a lack of in silico DO models and related mechano-regulatory tissue differentiation algorithms for mouse bone. This study sought to establish an in silico model based on in vivo experimental data to simulate the bone regeneration process during DO of the mouse femur. METHODS: In vivo micro-CT, including time-lapse morphometry was performed to monitor the bone regeneration in the distraction gap. A 2D axisymmetric finite element model, with a geometry originating from the experimental data, was created. Bone regeneration was simulated with a fuzzy logic-based two-stage (distraction and consolidation) mechano-regulatory tissue differentiation algorithm, which was adjusted from that used for fracture healing based on our in vivo experimental data. The predictive potential of the model was further tested with varied distraction frequencies and distraction rates. RESULTS: The computational simulations showed similar bone regeneration patterns to those observed in the micro-CT data from the experiment throughout the DO process. This consisted of rapid bone formation during the first 10 days of the consolidation phase, followed by callus reshaping via bone remodeling. In addition, the computational model predicted a faster and more robust bone healing response as the model's distraction frequency was increased, whereas higher or lower distraction rates were not conducive to healing. CONCLUSIONS: This in silico model could be used to investigate basic mechanobiological mechanisms involved in bone regeneration or to optimize DO strategies for potential clinical applications.

Orthophosphate and alkaline phosphatase induced the formation of apatite with different multilayered structures and mineralization balance.

Mineralized collagen is a natural organic-inorganic composite. The combination of organic collagen and inorganic apatite to form different nanostructures is the key to producing bone substitutes with biomechanical properties that are as identical to normal bone as possible. However, the formation of apatite with different nanostructures during collagen mineralization is unexplored. Here, pyrophosphate (Pyro-P), as an important hydrolysate of adenosine triphosphate in the body, was introduced to prepare mineralized collagen under the regulation of alkaline phosphatase (ALP) and orthophosphate (Ortho-P). Scanning electron microscopy (SEM) and transmission electron microscopy (TEM) results showed that mineralized collagen, which combined with different crystallinities and multilayered structured apatite, was successfully prepared. A combination of ion chromatography (IC), Fourier transform infrared (FTIR) spectroscopy, circular dichroism (CD), and thermogravimetry (TG) analyses revealed the crucial role of Ortho-P in the formation of multilayered flower-shaped apatite with different crystallinities and in the maintenance of mineralization balance. Mineralization balance is of great significance for maintaining normal bone morphology during bone regeneration. Overall, our results provide a promising method to produce new bone substitute materials for the repair of large bone defects and a deeper insight into the mechanisms of biomineralization.