1
|
Zhai Y, Han G, Gao Q, Bai W. Analysis of Transient Cutting Forces in Cortical Bone During Ultrasonically Assisted Cutting. J Biomech Eng 2025; 147:061003. [PMID: 40192649 DOI: 10.1115/1.4068371] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2024] [Indexed: 04/20/2025]
Abstract
Ultrasonically assisted cutting (UAC), a process characterized by high-performance material removal and enhanced surface finish, is widely employed in orthopedic surgery. However, variability in the mechanical properties of cortical bone may lead to unstable fractures and fluctuating cutting force during material removal, particularly under high-frequency vibration cutting. This study introduces a transient shear strength model that utilizes strain rate fluctuations to estimate cutting forces in the UAC process. The impact of varying osteon orientations and strain rate ranges on the yield strength of cortical bone is analyzed to elucidate changes in its mechanical properties under UAC conditions. Additionally, strain rates from conventional cutting (CC) and UAC, measured through digital image correlation (DIC), are compared with model predictions. The results demonstrate that the proposed model accurately predicts cutting forces and associated changes in thrust. This research offers a fresh insight into the dynamics of fluctuating forces during UAC, potentially inspiring advancements in orthopedic surgical instruments.
Collapse
Affiliation(s)
- Yuhao Zhai
- Mechanical Engineering and Electronic Information, China University of Geosciences, Wuhan 430074, China
| | - Guangchao Han
- Mechanical Engineering and Electronic Information, China University of Geosciences, Wuhan 430074, China; Shenzhen Research Institute, China University of Geosciences, Shenzhen 518057, China
- Department of Medicine, The University of Hong Kong-Shenzhen Hospital, Shenzhen 518053, China
| | - Qingpeng Gao
- Department of Medicine, The University of Hong Kong-Shenzhen Hospital, Shenzhen 518053, China
- University of Hong Kong - Shenzhen Hospital
| | - Wei Bai
- Mechanical Engineering and Electronic Information, China University of Geosciences, Wuhan 430074, China; Shenzhen Research Institute, China University of Geosciences, Shenzhen 518057, China
- Department of Medicine, The University of Hong Kong-Shenzhen Hospital, Shenzhen 518053, China
| |
Collapse
|
2
|
Skinner J, Langley N, Fahrenholtz S, Shanavas Y, Waletzki B, Brown R, Herrick J, Shyamsunder L, Goguen P, Rajan S. Microscopic characteristics of peri- and postmortem fracture surfaces. Forensic Sci Int 2024; 365:112223. [PMID: 39437525 DOI: 10.1016/j.forsciint.2024.112223] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2024] [Revised: 08/07/2024] [Accepted: 09/09/2024] [Indexed: 10/25/2024]
Abstract
This study investigated if microscopic surface features captured with a scanning electron microscope (SEM) effectively discriminate fracture timing. We hypothesized that microscopic fracture characteristics, including delamination, osteon pullout, and microcracks, may vary as bone elasticity decreases, elucidating perimortem and postmortem events more reliably than macroscopic analyses. Thirty-seven unembalmed, defleshed human femoral shafts from males (n=18) and females (n=2) aged 33-81 years were fractured at experimentally simulated postmortem intervals (PMIs) ranging from 1 to 60 warm weather days (250-40,600 ADH). A gravity convection oven was used to approximate tissue decomposition at 37 C and 27 C, and the resulting heat-time unit (accumulated degree hours, or ADH) was used to examine fractures in elastic/wet versus brittle/dry bone. The bones were fractured with a drop test frame using a three-point bending setup, sensors were used to calculate fracture energy, and high-speed photography documented fracture events. The following data were collected to relate fracture appearance to the biomechanical properties of bone: PMI (postmortem interval) length in ADH, temperature, humidity, collagen percentage, water loss, bone mineral density, cortical bone thickness, fracture energy, age, sex, cause of death, and microscopic fracture feature scores. SEM micrographs were collected from the primary tension zones of each fracture surface, and three microscopic fracture characteristics were scored from a region of interest in the center of the tension zone: percentage of delaminated osteons, percent osteon pullout, and number of microcracks. Multiple linear regression showed that microscopic fracture surface features are strong predictors of ADH (adjusted R-squared=0.67 for the 0 - 40,000 ADH samples; adjusted R-squared=0.92 for the 0-16,000 ADH samples). Osteon pullout is the single best predictor of ADH. Additionally, water loss is the primary driver of bone elasticity changes in low ADH samples, while collagen fibers appear to remain intact until later in the postmortem interval (approximately 40,000 ADH in this study). The results of this study indicate microscopic fracture surface analysis detects the biomechanical effects of decreased elasticity more reliably and with greater sensitivity than macroscopic analysis.
Collapse
Affiliation(s)
- Jessica Skinner
- Mayo Clinic Arizona, 13400 E Shea Blvd, Scottsdale, AZ 85259, USA.
| | - Natalie Langley
- Mayo Clinic Arizona, 13400 E Shea Blvd, Scottsdale, AZ 85259, USA
| | | | - Yuktha Shanavas
- State University of New York Upstate Medical University, 750 E Adams St, Syracuse, NY 13210, USA
| | - Brian Waletzki
- MilliporeSigma, 900 Middlesex Turnpike, Billerica, MA 01821, USA
| | - Robert Brown
- Mayo Clinic Rochester, 200 1st St SW, Rochester, MN 55902, USA
| | - James Herrick
- Mayo Clinic Rochester, 200 1st St SW, Rochester, MN 55902, USA
| | | | - Peter Goguen
- Arizona State University, 850 S McAllister Ave, Tempe, AZ 85281, USA
| | - Subramaniam Rajan
- Arizona State University, 850 S McAllister Ave, Tempe, AZ 85281, USA
| |
Collapse
|
3
|
Waidi YO, Debnath S, Datta S, Chatterjee K. 3D-Printed Silk Proteins for Bone Tissue Regeneration and Associated Immunomodulation. Biomacromolecules 2024; 25:5512-5540. [PMID: 39133748 DOI: 10.1021/acs.biomac.4c00540] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/10/2024]
Abstract
Current bone repair methods have limitations, prompting the exploration of innovative approaches. Tissue engineering emerges as a promising solution, leveraging biomaterials to craft scaffolds replicating the natural bone environment, facilitating cell growth and differentiation. Among fabrication techniques, three-dimensional (3D) printing stands out for its ability to tailor intricate scaffolds. Silk proteins (SPs), known for their mechanical strength and biocompatibility, are an excellent choice for engineering 3D-printed bone tissue engineering (BTE) scaffolds. This article comprehensively reviews bone biology, 3D printing, and the unique attributes of SPs, specifically detailing criteria for scaffold fabrication such as composition, structure, mechanics, and cellular responses. It examines the structural, mechanical, and biological attributes of SPs, emphasizing their suitability for BTE. Recent studies on diverse 3D printing approaches using SPs-based for BTE are highlighted, alongside advancements in their 3D and four-dimensional (4D) printing and their role in osteo-immunomodulation. Future directions in the use of SPs for 3D printing in BTE are outlined.
Collapse
Affiliation(s)
- Yusuf Olatunji Waidi
- Department of Materials Engineering, Indian Institute of Science, C. V. Raman Avenue, Bangalore 560012, India
| | - Souvik Debnath
- Department of Materials Engineering, Indian Institute of Science, C. V. Raman Avenue, Bangalore 560012, India
| | - Sudipto Datta
- Department of Materials Engineering, Indian Institute of Science, C. V. Raman Avenue, Bangalore 560012, India
| | - Kaushik Chatterjee
- Department of Bioengineering, Indian Institute of Science, C. V. Raman Avenue, Bangalore 560012, India
- Department of Materials Engineering, Indian Institute of Science, C. V. Raman Avenue, Bangalore 560012, India
| |
Collapse
|
4
|
Alijani H, Vaughan TJ. Exploring the hierarchical structure of lamellar bone and its impact on fracture behaviour: A computational study using a phase field damage model. J Mech Behav Biomed Mater 2024; 153:106471. [PMID: 38458079 DOI: 10.1016/j.jmbbm.2024.106471] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2023] [Revised: 02/02/2024] [Accepted: 02/16/2024] [Indexed: 03/10/2024]
Abstract
Bone is a naturally occurring composite material composed of a stiff mineral phase and a compliant organic matrix of collagen and non-collagenous proteins (NCP). While diverse mineral morphologies such as platelets and grains have been documented, the precise role of individual constituents, and their morphology, remains poorly understood. To understand the role of constituent morphology on the fracture behaviour of lamellar bone, a damage based representative volume element (RVE) was developed, which considered various mineral morphologies and mineralised collagen fibril (MCF) configurations. This model framework incorporated a novel phase-field damage model to predict the onset and evolution of damage at mineral-mineral and mineral-MCF interfaces. It was found that platelet-based mineral morphologies had superior mechanical performance over their granular counterparts, owing to their higher load-bearing capacity, resulting from a higher aspect ratio. It was also found that MCFs had a remarkable capacity for energy dissipation under axial loading, with these fibrillar structures acting as barriers to crack propagation, thereby enhancing overall elongation and toughness. Interestingly, the presence of extrafibrillar platelet-based minerals also provided an additional toughening through a similar mechanism, whereby these structures also inhibited crack propagation. These findings demonstrate that the two primary constituent materials of lamellar bone play a key role in its toughening behaviour, with combined effect by both mineral and MCFs to inhibit crack propagation at this scale. These results have provided novel insight into the fracture behaviour of lamellar bone, enhancing our understanding of microstructure-property relationships at the sub-tissue level.
Collapse
Affiliation(s)
- Hamid Alijani
- Biomechanics Research Centre (BioMEC), Biomedical Engineering, School of Engineering, College of Science and Engineering, University of Galway, Ireland
| | - Ted J Vaughan
- Biomechanics Research Centre (BioMEC), Biomedical Engineering, School of Engineering, College of Science and Engineering, University of Galway, Ireland.
| |
Collapse
|
5
|
Vercher-Martínez A, Megías R, Belda R, Vargas P, Giner E. Estimation of the in-plane ultimate stress of lamellar tissue as a function of bone mineral density and osteocyte lacunae porosity. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2024; 248:108120. [PMID: 38492277 DOI: 10.1016/j.cmpb.2024.108120] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/09/2024] [Revised: 03/05/2024] [Accepted: 03/06/2024] [Indexed: 03/18/2024]
Abstract
BACKGROUND AND OBJECTIVE Detailed finite element models based on medical images (μ-CT) are commonly used to analyze the mechanical behavior of bone at microscale. In order to simulate the tissue failure onset, isotropic failure criteria of lamellar tissue are often used, despite its non-isotropic and heterogeneous nature. The main goal of the present work is to estimate the in-plane ultimate stress of lamellar bone, considering the influence of mineral content and the porosity due to the osteocyte lacunae concentration. METHODS To this aim, a representative volume cell of lamellar tissue is modeled numerically, including: (1) non-isotropic elastic properties of tissue as a function of the bone mineral density and (2) explicit modeling of the osteocyte lacunae, considering the range of porosity content, size and orientation of ellipsoid-shaped lacunae. Firstly, the element size for the finite element models have been defined from a preliminary convergence analysis. Bounds on the ultimate stress of non-porous lamellar tissue are estimated for two values of bone mineral density, considering the results of tensile and compressive tests of wet osteons from the literature. Subsequently, the ultimate stress of lamellar tissue considering several values of micro-porosity is addressed. RESULTS Results obtained in this work show that the strength of lamellar bone decreases exponentially with the increase of lacunae porosity concentration. Ultimate stress of non-porous tissue (p=0%) increases with high mineral content, reaching a value of S¯transc=355.40±39.80 MPa for compression in the transversal direction of the fiber bundles, being BMD=1.246g/cm3. The mean value for the longitudinal to transverse strength ratio evaluated for porosity p=0%,1% and 5% and a mineral content BMD=1.2g/cm3, is 2.47:1 for tension and 1.55:1 for compression. These values are in agreement with literature. CONCLUSIONS Osteocyte lacunae act as stress concentrators, acting as potential stimulus for the bone regeneration process. A novel micromechanical model for the in-plane ultimate stress of lamellar tissue as a function of mineral content and lacunae concentration is presented. Additional considerations about the intralamellar shear stress evolution are also given.
Collapse
Affiliation(s)
- Ana Vercher-Martínez
- Dept. de Ingeniería Mecánica y de Materiales, Instituto de Ingeniería Mecánica y Biomecánica de Valencia - I2MB, Universitat Politècnica de València, Camino de Vera, Building 5E, 46022 Valencia, Spain.
| | - Raquel Megías
- Dept. de Ingeniería Mecánica y de Materiales, Instituto de Ingeniería Mecánica y Biomecánica de Valencia - I2MB, Universitat Politècnica de València, Camino de Vera, Building 5E, 46022 Valencia, Spain
| | - Ricardo Belda
- Dept. de Ingeniería Mecánica y de Materiales, Instituto de Ingeniería Mecánica y Biomecánica de Valencia - I2MB, Universitat Politècnica de València, Camino de Vera, Building 5E, 46022 Valencia, Spain
| | - Pablo Vargas
- Dept. de Ingeniería Mecánica y de Materiales, Instituto de Ingeniería Mecánica y Biomecánica de Valencia - I2MB, Universitat Politècnica de València, Camino de Vera, Building 5E, 46022 Valencia, Spain
| | - Eugenio Giner
- Dept. de Ingeniería Mecánica y de Materiales, Instituto de Ingeniería Mecánica y Biomecánica de Valencia - I2MB, Universitat Politècnica de València, Camino de Vera, Building 5E, 46022 Valencia, Spain
| |
Collapse
|
6
|
Alijani H, Vaughan TJ. Micromechanical modelling of transverse fracture behaviour of lamellar bone using a phase-field damage model: The role of non-collagenous proteins and mineralised collagen fibrils. J Mech Behav Biomed Mater 2024; 153:106472. [PMID: 38432183 DOI: 10.1016/j.jmbbm.2024.106472] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2023] [Revised: 12/20/2023] [Accepted: 02/18/2024] [Indexed: 03/05/2024]
Abstract
At the tissue-scale and above, there are now well-established structure-property relationships that provide good approximations of the biomechanical performance of bone through, for example, power-law relationships that relate tissue mineral density to elastic properties. However, below the tissue-level, the individual role of the constituents becomes prominent and these simple relationships tend to break down, with more detailed theoretical and computational models are required to describe the mechanical response. In this study, a two-dimensional micromechanics damage-based representative volume element (RVE) of lamellar bone was developed, which included a novel implementation of a phase-field damage model to describe the behaviour of non-collagenous proteins at mineral-mineral and mineral-fibril interface regions. It was found that, while the stiffness of the tissue was governed by the relative proportion of extra-fibrillar mineral and mineralised collagen fibrils, the strength and toughness of the tissue in transverse direction relied on the interactions occurring at mineral-mineral and mineral-fibril interfaces, highlighting the prominence of non-collagenous proteins in determine fracture-based processes at this scale. While fractures tended to initiate in mineral rich areas of the extra-fibrillar mineral matrix, it was found that the presence of mineralised collagen fibrils at low density did not provide a substantial contribution to crack propagation behaviour under transverse loading. However, at physiological volume fraction (VfMCF=50%), different scenarios could arise depending on the relative strength value of the interphase around the MCFs ( [Formula: see text] ) to the interphase between individual minerals ( [Formula: see text] ): (i) When [Formula: see text] , MCFs appear to facilitate crack propagation with MCF-mineral debonding being the dominant failure mode; (ii) once γ>1, the MCFs hinder the microcracks, leading to inhibition of crack propagation, which can be regarded as an energy dissipation mechanism. The effective fracture properties of the tissue also experience a sudden increase in fracture work density (J-integral) once the crack is arrested by MCFs or severely deflected. Collectively, the predicted behaviour of the model compared well to those reported through experimental and computational methods, highlighting its potential to provide further understanding into the mechanistic response of bone ultrastructure alterations related to the structural and compositional changes resulting from disease and aging.
Collapse
Affiliation(s)
- Hamid Alijani
- Biomechanics Research Centre (BioMEC), Biomedical Engineering, School of Engineering, College of Science and Engineering, University of Galway, Ireland
| | - Ted J Vaughan
- Biomechanics Research Centre (BioMEC), Biomedical Engineering, School of Engineering, College of Science and Engineering, University of Galway, Ireland.
| |
Collapse
|
7
|
Hu H, Yan M, Jiang J, Huang A, Cai S, Lan L, Ye K, Chen D, Tang K, Zuo Q, Zeng Y, Tang W, Fu J, Jiang C, Wang Y, Yan Z, He X, Qiao L, Zhao Y. A state-of-the-art review on biomass-derived carbon materials for supercapacitor applications: From precursor selection to design optimization. THE SCIENCE OF THE TOTAL ENVIRONMENT 2024; 912:169141. [PMID: 38072258 DOI: 10.1016/j.scitotenv.2023.169141] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/17/2023] [Revised: 11/22/2023] [Accepted: 12/04/2023] [Indexed: 12/27/2023]
Abstract
Biomass-derived carbon materials have the characteristics of a wide range of precursor sources, controllable carbon nano-dimension, large specific surface area and abundant heteroatoms doping. At present, biomass-derived carbon materials have been widely used in electrochemical energy storage devices, especially the research and development of biomass-derived carbon materials for supercapacitors has become mature and in-depth. Therefore, it is of importance to summarize the advanced technologies and strategies for optimizing biomass-derived carbon materials for supercapacitors, which will effectively promote the further development of high-performance supercapacitors. In this review, the recent research progress of biomass-derived carbon materials is provided in detail, including the selection of biomass precursors, the design of carbon nano-dimension and the theory of heteroatom doping. Besides, the preparation methods of biomass-derived carbon materials and the related processes of optimizing the electrochemical performance are also summarized. This review ends with the perspectives for future research directions and challenges in the field of biomass-derived carbon materials for electrochemical applications. This review aims to provide helpful reference information for the nano-dimensional design and electrochemical performance optimization of biomass-derived carbon materials for the practical application of supercapacitors.
Collapse
Affiliation(s)
- Hengyuan Hu
- College of Water Conservancy and Hydropower Engineering, Sichuan Agricultural University, Ya'an 625014, China
| | - Minglei Yan
- College of Water Conservancy and Hydropower Engineering, Sichuan Agricultural University, Ya'an 625014, China.
| | - Jietao Jiang
- Key Laboratory of Mesoscopic Chemistry of MOE and Jiangsu Provincial Lab for Nanotechnology, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China.
| | - Ankui Huang
- College of Water Conservancy and Hydropower Engineering, Sichuan Agricultural University, Ya'an 625014, China
| | - Sicheng Cai
- College of Water Conservancy and Hydropower Engineering, Sichuan Agricultural University, Ya'an 625014, China
| | - Linxuan Lan
- College of Water Conservancy and Hydropower Engineering, Sichuan Agricultural University, Ya'an 625014, China
| | - Kewei Ye
- College of Water Conservancy and Hydropower Engineering, Sichuan Agricultural University, Ya'an 625014, China
| | - Danlei Chen
- College of Water Conservancy and Hydropower Engineering, Sichuan Agricultural University, Ya'an 625014, China
| | - Kewen Tang
- College of Water Conservancy and Hydropower Engineering, Sichuan Agricultural University, Ya'an 625014, China
| | - Qin Zuo
- College of Water Conservancy and Hydropower Engineering, Sichuan Agricultural University, Ya'an 625014, China
| | - Yun Zeng
- College of Water Conservancy and Hydropower Engineering, Sichuan Agricultural University, Ya'an 625014, China
| | - Wei Tang
- College of Water Conservancy and Hydropower Engineering, Sichuan Agricultural University, Ya'an 625014, China
| | - Junheng Fu
- College of Water Conservancy and Hydropower Engineering, Sichuan Agricultural University, Ya'an 625014, China
| | - Chenglu Jiang
- College of Water Conservancy and Hydropower Engineering, Sichuan Agricultural University, Ya'an 625014, China
| | - Yong Wang
- School of Physics, University of Electronic Science and Technology of China, Chengdu 610054, China.
| | - Zhenhua Yan
- Key Laboratory of Advanced Energy Materials Chemistry (Ministry of Education), Nankai University, Tianjin 300071, China
| | - Xi He
- Shaanxi Joint Laboratory of Artificial Intelligence, Shaanxi University of Science and Technology, Xi'an 710021, China
| | - Liang Qiao
- School of Physics, University of Electronic Science and Technology of China, Chengdu 610054, China
| | - Yang Zhao
- College of Water Conservancy and Hydropower Engineering, Sichuan Agricultural University, Ya'an 625014, China.
| |
Collapse
|
8
|
Vellwock AE, Libonati F. XFEM for Composites, Biological, and Bioinspired Materials: A Review. MATERIALS (BASEL, SWITZERLAND) 2024; 17:745. [PMID: 38591618 PMCID: PMC10856485 DOI: 10.3390/ma17030745] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/27/2023] [Revised: 01/09/2024] [Accepted: 01/29/2024] [Indexed: 04/10/2024]
Abstract
The eXtended finite element method (XFEM) is a powerful tool for structural mechanics, assisting engineers and designers in understanding how a material architecture responds to stresses and consequently assisting the creation of mechanically improved structures. The XFEM method has unraveled the extraordinary relationships between material topology and fracture behavior in biological and engineered materials, enhancing peculiar fracture toughening mechanisms, such as crack deflection and arrest. Despite its extensive use, a detailed revision of case studies involving XFEM with a focus on the applications rather than the method of numerical modeling is in great need. In this review, XFEM is introduced and briefly compared to other computational fracture models such as the contour integral method, virtual crack closing technique, cohesive zone model, and phase-field model, highlighting the pros and cons of the methods (e.g., numerical convergence, commercial software implementation, pre-set of crack parameters, and calculation speed). The use of XFEM in material design is demonstrated and discussed, focusing on presenting the current research on composites and biological and bioinspired materials, but also briefly introducing its application to other fields. This review concludes with a discussion of the XFEM drawbacks and provides an overview of the future perspectives of this method in applied material science research, such as the merging of XFEM and artificial intelligence techniques.
Collapse
Affiliation(s)
- Andre E. Vellwock
- B CUBE—Center for Molecular Bioengineering, Technische Universität Dresden, 01307 Dresden, Germany;
| | - Flavia Libonati
- Department of Mechanical, Energy, Management and Transportation Engineering, University of Genoa, 16145 Genoa, Italy
| |
Collapse
|
9
|
Indermaur M, Casari D, Kochetkova T, Willie BM, Michler J, Schwiedrzik J, Zysset P. Tensile Mechanical Properties of Dry Cortical Bone Extracellular Matrix: A Comparison Among Two Osteogenesis Imperfecta and One Healthy Control Iliac Crest Biopsies. JBMR Plus 2023; 7:e10826. [PMID: 38130764 PMCID: PMC10731133 DOI: 10.1002/jbm4.10826] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 02/27/2023] [Revised: 08/22/2023] [Accepted: 09/14/2023] [Indexed: 12/23/2023] Open
Abstract
Osteogenesis imperfecta (OI) is a genetic, collagen-related bone disease that increases the incidence of bone fractures. Still, the origin of this brittle mechanical behavior remains unclear. The extracellular matrix (ECM) of OI bone exhibits a higher degree of bone mineralization (DBM), whereas compressive mechanical properties at the ECM level do not appear to be inferior to healthy bone. However, it is unknown if collagen defects alter ECM tensile properties. This study aims to quantify the tensile properties of healthy and OI bone ECM. In three transiliac biopsies (healthy n = 1, OI type I n = 1, OI type III n = 1), 23 microtensile specimens (gauge dimensions 10 × 5 × 2 μm3) were manufactured and loaded quasi-statically under tension in vacuum condition. The resulting loading modulus and ultimate strength were extracted. Interestingly, tensile properties in OI bone ECM were not inferior compared to controls. All specimens revealed a brittle failure behavior. Fracture surfaces were graded according to their mineralized collagen fibers (MCF) orientation into axial, mixed, and transversal fracture surface types (FST). Furthermore, tissue mineral density (TMD) of the biopsy cortices was extracted from micro-computed tomogra[hy (μCT) images. Both FST and TMD are significant factors to predict loading modulus and ultimate strength with an adjusted R 2 of 0.556 (p = 2.65e-05) and 0.46 (p = 2.2e-04), respectively. The influence of MCF orientation and DBM on the mechanical properties of the neighboring ECM was further verified with quantitative polarized Raman spectroscopy (qPRS) and site-matched nanoindentation. MCF orientation and DBM were extracted from the qPRS spectrum, and a second mechanical model was developed to predict the indentation modulus with MCF orientation and DBM (R 2 = 67.4%, p = 7.73e-07). The tensile mechanical properties of the cortical bone ECM of two OI iliac crest biopsies are not lower than the one from a healthy and are primarily dependent on MCF orientation and DBM. © 2023 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.
Collapse
Affiliation(s)
- Michael Indermaur
- ARTORG Center for Biomedical EngineeringUniversity of BernBernSwitzerland
| | - Daniele Casari
- ARTORG Center for Biomedical EngineeringUniversity of BernBernSwitzerland
- Swiss Federal Laboratories for Material Science and TechnologyThunSwitzerland
| | - Tatiana Kochetkova
- Swiss Federal Laboratories for Material Science and TechnologyThunSwitzerland
| | - Bettina M. Willie
- Research Centre, Shriners Hospital for Children‐Canada, Department of Pediatric SurgeryMcGill UniversityMontrealQCCanada
| | - Johann Michler
- Swiss Federal Laboratories for Material Science and TechnologyThunSwitzerland
| | - Jakob Schwiedrzik
- Swiss Federal Laboratories for Material Science and TechnologyThunSwitzerland
| | - Philippe Zysset
- ARTORG Center for Biomedical EngineeringUniversity of BernBernSwitzerland
| |
Collapse
|
10
|
Kochetkova T, Hanke MS, Indermaur M, Groetsch A, Remund S, Neuenschwander B, Michler J, Siebenrock KA, Zysset P, Schwiedrzik J. Composition and micromechanical properties of the femoral neck compact bone in relation to patient age, sex and hip fracture occurrence. Bone 2023; 177:116920. [PMID: 37769956 DOI: 10.1016/j.bone.2023.116920] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 07/04/2023] [Revised: 09/22/2023] [Accepted: 09/25/2023] [Indexed: 10/03/2023]
Abstract
Current clinical methods of bone health assessment depend to a great extent on bone mineral density (BMD) measurements. However, these methods only act as a proxy for bone strength and are often only carried out after the fracture occurs. Besides BMD, composition and tissue-level mechanical properties are expected to affect the whole bone's strength and toughness. While the elastic properties of the bone extracellular matrix (ECM) have been extensively investigated over the past two decades, there is still limited knowledge of the yield properties and their relationship to composition and architecture. In the present study, morphological, compositional and micropillar compression bone data was collected from patients who underwent hip arthroplasty. Femoral neck samples from 42 patients were collected together with anonymous clinical information about age, sex and primary diagnosis (coxarthrosis or hip fracture). The femoral neck cortex from the inferomedial region was analyzed in a site-matched manner using a combination of micromechanical testing (nanoindentation, micropillar compression) together with micro-CT and quantitative polarized Raman spectroscopy for both morphological and compositional characterization. Mechanical properties, as well as the sample-level mineral density, were constant over age. Only compositional properties demonstrate weak dependence on patient age: decreasing mineral to matrix ratio (p = 0.02, R2 = 0.13, 2.6 % per decade) and increasing amide I sub-peak ratio I∼1660/I∼1683 (p = 0.04, R2 = 0.11, 1.5 % per decade). The patient's sex and diagnosis did not seem to influence investigated bone properties. A clear zonal dependence between interstitial and osteonal cortical zones was observed for compositional and elastic bone properties (p < 0.0001). Site-matched microscale analysis confirmed that all investigated mechanical properties except yield strain demonstrate a positive correlation with the mineral fraction of bone. The output database is the first to integrate the experimentally assessed microscale yield properties, local tissue composition and morphology with the available patient clinical information. The final dataset was used for bone fracture risk prediction in-silico through the principal component analysis and the Naïve Bayes classification algorithm. The analysis showed that the mineral to matrix ratio, indentation hardness and micropillar yield stress are the most relevant parameters for bone fracture risk prediction at 70 % model accuracy (0.71 AUC). Due to the low number of samples, further studies to build a universal fracture prediction algorithm are anticipated with the higher number of patients (N > 200). The proposed classification algorithm together with the output dataset of bone tissue properties can be used for the future comparison of existing methods to evaluate bone quality as well as to form a better understanding of the mechanisms through which bone tissue is affected by aging or disease.
Collapse
Affiliation(s)
- Tatiana Kochetkova
- Empa, Swiss Federal Laboratories for Materials Science and Technology, Thun, Switzerland.
| | - Markus S Hanke
- Department of Orthopedic Surgery, Inselspital, University of Bern, Switzerland
| | - Michael Indermaur
- ARTORG Center for Biomedical Engineering Research, University of Bern, Switzerland
| | - Alexander Groetsch
- Empa, Swiss Federal Laboratories for Materials Science and Technology, Thun, Switzerland
| | - Stefan Remund
- Institute for Applied Laser, Photonics and Surface Technologies (ALPS), Bern University of Applied Sciences, Burgdorf, Switzerland
| | - Beat Neuenschwander
- Institute for Applied Laser, Photonics and Surface Technologies (ALPS), Bern University of Applied Sciences, Burgdorf, Switzerland
| | - Johann Michler
- Empa, Swiss Federal Laboratories for Materials Science and Technology, Thun, Switzerland
| | - Klaus A Siebenrock
- Department of Orthopedic Surgery, Inselspital, University of Bern, Switzerland
| | - Philippe Zysset
- ARTORG Center for Biomedical Engineering Research, University of Bern, Switzerland
| | - Jakob Schwiedrzik
- Empa, Swiss Federal Laboratories for Materials Science and Technology, Thun, Switzerland.
| |
Collapse
|
11
|
Zhang Y, Cheng Z, Liu Z, Shen X, Cai C, Li M, Luo Z. Functionally Tailored Metal-Organic Framework Coatings for Mediating Ti Implant Osseointegration. ADVANCED SCIENCE (WEINHEIM, BADEN-WURTTEMBERG, GERMANY) 2023; 10:e2303958. [PMID: 37705110 PMCID: PMC10582459 DOI: 10.1002/advs.202303958] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/15/2023] [Revised: 07/24/2023] [Indexed: 09/15/2023]
Abstract
Owing to their mechanical resilience and non-toxicity, titanium implants are widely applied as the major treatment modality for the clinical intervention against bone fractures. However, the intrinsic bioinertness of Ti and its alloys often impedes the effective osseointegration of the implants, leading to severe adverse complications including implant loosening, detachment, and secondary bone damage. Consequently, new Ti implant engineering strategies are urgently needed to improve their osseointegration after implantation. Remarkably, metalorganic frameworks (MOFs) are a class of novel synthetic material consisting of coordinated metal species and organic ligands, which have demonstrated a plethora of favorable properties for modulating the interfacial properties of Ti implants. This review comprehensively summarizes the recent progress in the development of MOF-coated Ti implants and highlights their potential utility for modulating the bio-implant interface to improve implant osseointegration, of which the discussions are outlined according to their physical traits, chemical composition, and drug delivery capacity. A perspective is also provided in this review regarding the current limitations and future opportunities of MOF-coated Ti implants for orthopedic applications. The insights in this review may facilitate the rational design of more advanced Ti implants with enhanced therapeutic performance and safety.
Collapse
Affiliation(s)
- Yuan Zhang
- Joint Disease & Sport Medicine CentreDepartment of OrthopaedicsXinqiao HospitalArmy Medical UniversityChongqing400038China
| | - Zhuo Cheng
- School of Life ScienceChongqing UniversityChongqing400044China
| | - Zaiyang Liu
- Joint Disease & Sport Medicine CentreDepartment of OrthopaedicsXinqiao HospitalArmy Medical UniversityChongqing400038China
| | - Xinkun Shen
- Department of OrthopaedicsRuian People's HospitalThe Third Affiliated Hospital of Wenzhou Medical UniversityWenzhou325016China
| | - Chunyuan Cai
- Department of OrthopaedicsRuian People's HospitalThe Third Affiliated Hospital of Wenzhou Medical UniversityWenzhou325016China
| | - Menghuan Li
- School of Life ScienceChongqing UniversityChongqing400044China
| | - Zhong Luo
- School of Life ScienceChongqing UniversityChongqing400044China
| |
Collapse
|
12
|
Marty AG, Barbone PE, Morgan EF. Multiscale theoretical model shows that aging-related mechanical degradation of cortical bone is driven by microstructural changes in addition to porosity. J Mech Behav Biomed Mater 2023; 145:106029. [PMID: 37499524 PMCID: PMC10528045 DOI: 10.1016/j.jmbbm.2023.106029] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2023] [Revised: 07/12/2023] [Accepted: 07/14/2023] [Indexed: 07/29/2023]
Abstract
This study aims to gain mechanistic understanding of how aging-related changes in the microstructure of cortical bone drive mechanical consequences at the macroscale. To that end, cortical bone was modeled as a bundle of elastic-plastic, parallel fibers, which represented osteons and interstitial tissue, loaded in uniaxial tension. Distinct material properties were assigned to each fiber in either the osteon or interstitial fiber "families." Models representative of mature (20-60 yrs.) bone, and elderly (60+) bone were created by modeling aging via the following changes to the input parameters: (i) increasing porosity from 5% to 15%, (ii) increasing the ratio of the number of osteon fibers relative to interstitial fibers from 40% to 50%, and (iii) changing the fiber material properties from representing mature bone samples to representing elderly bone samples (i.e., increased strength and decreased toughness of interstitial fibers together with decreased toughness of osteon fibers). To understand the respective contributions of these changes, additional models isolating one or two of each of these were also created. From the computed stress-strain curve for the fiber bundle, the yield point (ϵy, σy), ultimate point (ϵu, σu), and toughness (UT) for the bundle as a whole were measured. We found that changes to all three input parameters were required for the model to capture the aging-related decline in cortical bone mechanical properties consistent with those previously reported in the literature. In both mature and elderly bundles, rupture of the interstitial fibers drove the initial loss of strength following the ultimate point. Plasticity and more gradual rupture of the osteons drove the remainder of the response. Both the onset and completion of interstitial fiber rupture occurred at lower strains in the elderly vs. mature case. These findings point to the importance of studying microstructural changes beyond porosity, such as the area fraction of osteons and the material properties of osteon and interstitial tissue, in order to further understanding of aging-related changes in bone.
Collapse
Affiliation(s)
- André Gutiérrez Marty
- Mechanical Engineering, Boston University, 110 Cummington Mall, Boston, 02115, MA, USA; Center for Multiscale and Translational Mechanobiology, Boston University, 110 Cummington Mall, Boston, 02115, MA, USA.
| | - Paul E Barbone
- Mechanical Engineering, Boston University, 110 Cummington Mall, Boston, 02115, MA, USA; Center for Multiscale and Translational Mechanobiology, Boston University, 110 Cummington Mall, Boston, 02115, MA, USA.
| | - Elise F Morgan
- Mechanical Engineering, Boston University, 110 Cummington Mall, Boston, 02115, MA, USA; Center for Multiscale and Translational Mechanobiology, Boston University, 110 Cummington Mall, Boston, 02115, MA, USA; Biomedical Engineering, Boston University, 44 Cummington Mall, Boston, 02115, MA, USA.
| |
Collapse
|
13
|
Fernández MP, Schwiedrzik J, Bürki A, Peyrin F, Michler J, Zysset PK, Wolfram U. In situ synchrotron radiation μCT indentation of cortical bone: Anisotropic crack propagation, local deformation, and fracture. Acta Biomater 2023:S1742-7061(23)00237-4. [PMID: 37127075 DOI: 10.1016/j.actbio.2023.04.038] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2023] [Revised: 04/21/2023] [Accepted: 04/25/2023] [Indexed: 05/03/2023]
Abstract
The development of treatment strategies for skeletal diseases relies on the understanding of bone mechanical properties in relation to its structure at different length scales. At the microscale, indention techniques can be used to evaluate the elastic, plastic, and fracture behaviour of bone tissue. Here, we combined in situ high-resolution SRμCT indentation testing and digital volume correlation to elucidate the anisotropic crack propagation, deformation, and fracture of ovine cortical bone under Berkovich and spherical tips. Independently of the indenter type we observed significant dependence of the crack development due to the anisotropy ahead of the tip, with lower strains and smaller crack systems developing in samples indented in the transverse material direction, where the fibrillar bone ultrastructure is largely aligned perpendicular to the indentation direction. Such alignment allows to accommodate the strain energy, inhibiting crack propagation. Higher tensile hoop strains generally correlated with regions that display significant cracking radial to the indenter, indicating a predominant Mode I fracture. This was confirmed by the three-dimensional analysis of crack opening displacements and stress intensity factors along the crack front obtained for the first time from full displacement fields in bone tissue. The X-ray beam significantly influenced the relaxation behaviour independent of the tip. Raman analyses did not show significant changes in specimen composition after irradiation compared to non-irradiated tissue, suggesting an embrittlement process that may be linked to damage of the non-fibrillar organic matrix. This study highlights the importance of three-dimensional investigation of bone deformation and fracture behaviour to explore the mechanisms of bone failure in relation to structural changes due to aging or disease. STATEMENT OF SIGNIFICANCE: : Characterising the three-dimensional deformation and fracture behaviour of bone remains essential to decipher the interplay between structure, function, and composition with the aim to improve fracture prevention strategies. The experimental methodology presented here, combining high-resolution imaging, indentation testing and digital volume correlation, allows us to quantify the local deformation, crack propagation, and fracture modes of cortical bone tissue. Our results highlight the anisotropic behaviour of osteonal bone and the complex crack propagation patterns and fracture modes initiating by the intricate stress states beneath the indenter tip. This is of wide interest not only for the understanding of bone fracture but also to understand other architectured (bio)structures providing an effective way to quantify their toughening mechanisms in relation to their main mechanical function.
Collapse
Affiliation(s)
- Marta Peña Fernández
- School of Engineering and Physical Science, Institute for Mechanical, Process and Energy Engineering, Heriot-Watt University, Edinburgh, United Kingdom
| | - Jakob Schwiedrzik
- EMPA, Swiss Federal Laboratories for Materials Science and Technology, Laboratory for Mechanics of Materials and Nanostructures, Thun, Switzerland
| | - Alexander Bürki
- ARTORG Centre for Biomedical Engineering Research, University of Bern, Bern, Switzerland
| | - Françoise Peyrin
- Université de Lyon, INSA Lyon, Université Claude Bernard Lyon 1, UJM-Saint Etienne, CNRS UMR 5220, Inserm U1294, CREATIS, Lyon, France
| | - Johann Michler
- EMPA, Swiss Federal Laboratories for Materials Science and Technology, Laboratory for Mechanics of Materials and Nanostructures, Thun, Switzerland
| | - Philippe K Zysset
- ARTORG Centre for Biomedical Engineering Research, University of Bern, Bern, Switzerland
| | - Uwe Wolfram
- School of Engineering and Physical Science, Institute for Mechanical, Process and Energy Engineering, Heriot-Watt University, Edinburgh, United Kingdom.
| |
Collapse
|
14
|
Xu M, An B, Zhang D. Competing mechanisms in fracture of staggered mineralized collagen fibril arrays. J Mech Behav Biomed Mater 2023; 141:105761. [PMID: 36905708 DOI: 10.1016/j.jmbbm.2023.105761] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2022] [Revised: 01/17/2023] [Accepted: 03/05/2023] [Indexed: 03/08/2023]
Abstract
Mineralized collagen fibril (MCF) arrays are important structural elements involved in inelastic deformation and fracture process of bone. Recent experiments have shown that MCF breakage has an influence on toughening of bone. Motivated by the experiments, we carry out the analyses of fracture in staggered MCF arrays. The plastic deformation of extrafibrillar matrix (EFM), debonding of the MCF-EFM interface, plastic deformation of MCFs and MCF fracture are accounted for in the calculations. It is found that the fracture of MCF arrays is controlled by competition between MCF breakage and debonding of the MCF-EFM interface. The MCF-EFM interface with high shear strength and large shear fracture energy is capable of activating MCF breakage, which promotes plastic energy dissipation of MCF arrays. In the absence of MCF breakage, damage energy dissipation is higher than plastic energy dissipation and debonding of the MCF-EFM interface provides the major contribution to toughening of bone. We have further revealed that the relative contributions of interfacial debonding mechanism and plastic deformation of MCF arrays are dependent on fracture properties of the MCF-EFM interface in the normal direction. The high normal strength gives rise to enhanced damage energy dissipation and amplified plastic deformation of MCF arrays; while high normal fracture energy of the interface suppresses plastic deformation of MCFs.
Collapse
Affiliation(s)
- Min Xu
- Shanghai Key Laboratory of Mechanics in Energy Engineering, Shanghai Institute of Applied Mathematics and Mechanics, School of Mechanics and Engineering Science, Shanghai University, Shanghai, 200444, China
| | - Bingbing An
- Shanghai Key Laboratory of Mechanics in Energy Engineering, Shanghai Institute of Applied Mathematics and Mechanics, School of Mechanics and Engineering Science, Shanghai University, Shanghai, 200444, China; Shaoxing Institute of Technology, Shanghai University, Shaoxing, 312074, China.
| | - Dongsheng Zhang
- Shanghai Key Laboratory of Mechanics in Energy Engineering, Shanghai Institute of Applied Mathematics and Mechanics, School of Mechanics and Engineering Science, Shanghai University, Shanghai, 200444, China; Shaoxing Institute of Technology, Shanghai University, Shaoxing, 312074, China
| |
Collapse
|
15
|
Alizadeh E, Omairey S, Zysset P. Investigating the post-yield behavior of mineralized bone fibril arrays using a 3D non-linear finite element unit-cell model. J Mech Behav Biomed Mater 2023; 139:105660. [PMID: 36638635 DOI: 10.1016/j.jmbbm.2023.105660] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2022] [Revised: 12/24/2022] [Accepted: 01/04/2023] [Indexed: 01/09/2023]
Abstract
In this study, we propose a 3D non-linear finite element (FE) unit-cell model to investigate the post-yield behavior of mineralized collagen fibril arrays (FAY). We then compare the predictions of the model with recent micro-tensile and micropillar compression tests in both axial and transverse directions. The unit cell consists of mineralized collagen fibrils (MCFs) embedded in an extrafibrillar matrix (EFM), and the FE mesh is equipped with cohesive interactions and a custom plasticity model. The simulation results confirm that MCF plays a dominant role in load bearing prior to yielding under axial tensile loading. Damage was initiated via debonding in shear and progressive sliding at the MCF/EFM interface, and resulted in MCF pull-out until brittle failure. In transverse tensile loading, EFM carried most of the load in pre-yield deformation, and then mixed normal/shear debonding between MCF and EFM began to form, which eventually produced brittle delamination of the two phases. The loading/unloading FE analysis in compression along both axial and transverse directions demonstrated perfect plasticity without any reduction in elastic modulus, i.e., damage due to the interfaces as seen in micropillar compression. Beyond the brittle and ductile nature of the stress-strain curves, in tensile and compressive loading, the simulated post-yield behavior and failure mechanism are in good quantitative agreement with the experimental observations. Our rather simple but efficient unit-cell FE model can reproduce qualitatively and quantitatively the mechanical behavior of bone ECM under tensile and compressive loading along the two main orientations. The model's integration into higher length scales may be useful in describing the macroscopic post-yield and failure behavior of trabecular and cortical bone in greater detail.
Collapse
Affiliation(s)
- Elham Alizadeh
- ARTORG Centre for Biomedical Engineering Research, University of Bern, Bern, Switzerland.
| | - Sadik Omairey
- Brunel Composites Centre, College of Engineering, Design and Physical Sciences, Brunel University London, London, UK
| | - Philippe Zysset
- ARTORG Centre for Biomedical Engineering Research, University of Bern, Bern, Switzerland
| |
Collapse
|
16
|
Lamellar thickness measurements in control and osteogenesis imperfecta human bone, with development of a method of automated thickness averaging to simplify quantitation. Matrix Biol 2023; 116:85-101. [PMID: 36592737 DOI: 10.1016/j.matbio.2022.12.006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2022] [Revised: 12/18/2022] [Accepted: 12/29/2022] [Indexed: 01/01/2023]
Abstract
Lamellar bone that forms in moderate and severe osteogenesis imperfecta (OI) is composed of structurally irregular lamellae compared to those in control bone. OI and control cortical bone fragments were prepared for light microscopy in standardized fashion: decalcified, embedded in plastic, sectioned and stained with toluidine blue. Polarization light microscopy (PLM) was used to demonstrate and quantify bright and dark lamellar thicknesses in cortical bone fragments from 5 patients with moderate to severe OI in whom type I collagen structural/molecular defects were detected and in control bone from 5 patients. Rigid selection criteria identified lamellar regions for quantification. Thicknesses of bright and dark lamellae were measured manually at 20X magnification using a histomorphometric image analysis system. A method of automated thickness averaging was developed to determine lamellar thicknesses from PLM images to make measurement faster. Our study demonstrates, for the first time, that in OI bone from patients with type I collagen structural/molecular defects mean lamellar thickness measurements (along with the bright and dark lamellar thicknesses) were less than those in control bone by statistically highly significant differences. The mean value for bright lamellae was less than that for dark lamellae in both control and OI bone. The ratio of mean values for bright/dark lamellar thicknesses was the same in control and OI bone. The automated method obtained similar results to the manual method. Lamellar bone in moderate and severe OI with type I collagen defects is composed of thinner and less structurally regular lamellae than those in control bone. This finding indicates that lamellar thickness measurements can be helpful in assessing the effect of specific collagen and collagen-related mutations on OI bone synthesis and warrant inclusion in research and clinical histomorphometric assessments.
Collapse
|
17
|
Groetsch A, Stelzl S, Nagel Y, Kochetkova T, Scherrer NC, Ovsianikov A, Michler J, Pethö L, Siqueira G, Nyström G, Schwiedrzik J. Microscale 3D Printing and Tuning of Cellulose Nanocrystals Reinforced Polymer Nanocomposites. SMALL (WEINHEIM AN DER BERGSTRASSE, GERMANY) 2023; 19:e2202470. [PMID: 36449596 DOI: 10.1002/smll.202202470] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/08/2022] [Revised: 11/04/2022] [Indexed: 06/17/2023]
Abstract
The increasing demand for functional materials and an efficient use of sustainable resources makes the search for new material systems an ever growing endeavor. With this respect, architected (meta-)materials attract considerable interest. Their fabrication at the micro- and nanoscale, however, remains a challenge, especially for composites with highly different phases and unmodified reinforcement fillers. This study demonstrates that it is possible to create a non-cytotoxic nanocomposite ink reinforced by a sustainable phase, cellulose nanocrystals (CNCs), to print and tune complex 3D architectures using two-photon polymerization, thus, advancing the state of knowledge toward the microscale. Micro-compression, high-res scanning electron microscopy, (polarised) Raman spectroscopy, and composite modeling are used to study the structure-property relationships. A 100% stiffness increase is observed already at 4.5 wt% CNC while reaching a high photo-polymerization degree of ≈80% for both neat polymers and CNC-composites. Polarized Raman and the Halpin-Tsai composite-model suggest a random CNC orientation within the polymer matrix. The microscale approach can be used to tune arbitrary small scale CNC-reinforced polymer-composites with comparable feature sizes. The new insights pave the way for future applications where the 3D printing of small structures is essential to improve performances of tissue-scaffolds, extend bio-electronics applications or tailor microscale energy-absorption devices.
Collapse
Affiliation(s)
- Alexander Groetsch
- Laboratory for Mechanics of Materials and Nanostructures Department of Advanced Materials and Surfaces, Empa - Swiss Federal Laboratories for Materials Science and Technology, Feuerwerkerstrasse 39, Thun, 3602, Switzerland
| | - Samuel Stelzl
- Research Group 3D Printing and Biofabrication Institute of Materials Science and Technology, TU Wien, Getreidemarkt 9, Vienna, 1060, Austria
| | - Yannick Nagel
- Laboratory for Cellulose & Wood Materials, Empa - Swiss Federal Laboratories for Materials Science and Technology, Ueberlandstrasse 129, Dübendorf, 8600, Switzerland
| | - Tatiana Kochetkova
- Laboratory for Mechanics of Materials and Nanostructures Department of Advanced Materials and Surfaces, Empa - Swiss Federal Laboratories for Materials Science and Technology, Feuerwerkerstrasse 39, Thun, 3602, Switzerland
| | - Nadim C Scherrer
- Bern University of Applied Sciences, HKB, Bern, 3027, Switzerland
| | - Aleksandr Ovsianikov
- Research Group 3D Printing and Biofabrication Institute of Materials Science and Technology, TU Wien, Getreidemarkt 9, Vienna, 1060, Austria
| | - Johann Michler
- Laboratory for Mechanics of Materials and Nanostructures Department of Advanced Materials and Surfaces, Empa - Swiss Federal Laboratories for Materials Science and Technology, Feuerwerkerstrasse 39, Thun, 3602, Switzerland
| | - Laszlo Pethö
- Laboratory for Mechanics of Materials and Nanostructures Department of Advanced Materials and Surfaces, Empa - Swiss Federal Laboratories for Materials Science and Technology, Feuerwerkerstrasse 39, Thun, 3602, Switzerland
| | - Gilberto Siqueira
- Laboratory for Cellulose & Wood Materials, Empa - Swiss Federal Laboratories for Materials Science and Technology, Ueberlandstrasse 129, Dübendorf, 8600, Switzerland
| | - Gustav Nyström
- Laboratory for Cellulose & Wood Materials, Empa - Swiss Federal Laboratories for Materials Science and Technology, Ueberlandstrasse 129, Dübendorf, 8600, Switzerland
- Department of Health Sciences and Technology, ETH Zurich, Universitätstrasse 2, Zurich, 8092, Switzerland
| | - Jakob Schwiedrzik
- Laboratory for Mechanics of Materials and Nanostructures Department of Advanced Materials and Surfaces, Empa - Swiss Federal Laboratories for Materials Science and Technology, Feuerwerkerstrasse 39, Thun, 3602, Switzerland
| |
Collapse
|
18
|
Alunni Cardinali M, Di Michele A, Mattarelli M, Caponi S, Govoni M, Dallari D, Brogini S, Masia F, Borri P, Langbein W, Palombo F, Morresi A, Fioretto D. Brillouin-Raman microspectroscopy for the morpho-mechanical imaging of human lamellar bone. J R Soc Interface 2022; 19:20210642. [PMID: 35104431 PMCID: PMC8807060 DOI: 10.1098/rsif.2021.0642] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2021] [Accepted: 12/21/2021] [Indexed: 12/12/2022] Open
Abstract
Bone has a sophisticated architecture characterized by a hierarchical organization, starting at the sub-micrometre level. Thus, the analysis of the mechanical and structural properties of bone at this scale is essential to understand the relationship between its physiology, physical properties and chemical composition. Here, we unveil the potential of Brillouin-Raman microspectroscopy (BRaMS), an emerging correlative optical approach that can simultaneously assess bone mechanics and chemistry with micrometric resolution. Correlative hyperspectral imaging, performed on a human diaphyseal ring, reveals a complex microarchitecture that is reflected in extremely rich and informative spectra. An innovative method for mechanical properties analysis is proposed, mapping the intermixing of soft and hard tissue areas and revealing the coexistence of regions involved in remodelling processes, nutrient transportation and structural support. The mineralized regions appear elastically inhomogeneous, resembling the pattern of the osteons' lamellae, while Raman and energy-dispersive X-ray images through scanning electron microscopy show an overall uniform distribution of the mineral content, suggesting that other structural factors are responsible for lamellar micromechanical heterogeneity. These results, besides giving an important insight into cortical bone tissue properties, highlight the potential of BRaMS to access the origin of anisotropic mechanical properties, which are almost ubiquitous in other biological tissues.
Collapse
Affiliation(s)
- M. Alunni Cardinali
- Department of Physics and Geology, University of Perugia, Via A. Pascoli, Perugia 06123, Italy
| | - A. Di Michele
- Department of Physics and Geology, University of Perugia, Via A. Pascoli, Perugia 06123, Italy
| | - M. Mattarelli
- Department of Physics and Geology, University of Perugia, Via A. Pascoli, Perugia 06123, Italy
| | - S. Caponi
- Istituto Officina Dei Materiali, National Research Council (IOM-CNR), Unit of Perugia, c/o Department of Physics and Geology, University of Perugia, Via A. Pascoli, Perugia 06123, Italy
| | - M. Govoni
- Reconstructive Orthopaedic Surgery and Innovative Techniques – Musculoskeletal Tissue Bank, IRCCS Istituto Ortopedico Rizzoli, Via G.C. Pupilli 1, Bologna 40136, Italy
| | - D. Dallari
- Reconstructive Orthopaedic Surgery and Innovative Techniques – Musculoskeletal Tissue Bank, IRCCS Istituto Ortopedico Rizzoli, Via G.C. Pupilli 1, Bologna 40136, Italy
| | - S. Brogini
- Complex Structure of Surgical Sciences and Technologies, IRCCS Istituto Ortopedico Rizzoli, Via Di Barbiano 1/10, Bologna 40136, Italy
| | - F. Masia
- School of Biosciences, Cardiff University, Museum Avenue, Cardiff CF10 3AX, UK
| | - P. Borri
- School of Biosciences, Cardiff University, Museum Avenue, Cardiff CF10 3AX, UK
| | - W. Langbein
- School of Physics and Astronomy, Cardiff University, The Parade, Cardiff CF24 3AA, UK
| | - F. Palombo
- School of Physics and Astronomy, University of Exeter, Exeter EX4 4QL, UK
| | - A. Morresi
- Department of Chemistry, Biology and Biotechnology, University of Perugia, Via Elce di Sotto 8, Perugia 06123, Italy
| | - D. Fioretto
- Department of Physics and Geology, University of Perugia, Via A. Pascoli, Perugia 06123, Italy
- CEMIN - Center of Excellence for Innovative Nanostructured Material, Via Elce di Sotto 8, Perugia 06123, Italy
| |
Collapse
|
19
|
A multiscale finite element investigation on the role of intra- and extra-fibrillar mineralisation on the elastic properties of bone tissue. J Mech Behav Biomed Mater 2022; 129:105139. [DOI: 10.1016/j.jmbbm.2022.105139] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2021] [Revised: 02/01/2022] [Accepted: 02/19/2022] [Indexed: 11/24/2022]
|
20
|
Abstract
Understanding the properties of bone is of both fundamental and clinical relevance. The basis of bone’s quality and mechanical resilience lies in its nanoscale building blocks (i.e., mineral, collagen, non-collagenous proteins, and water) and their complex interactions across length scales. Although the structure–mechanical property relationship in healthy bone tissue is relatively well characterized, not much is known about the molecular-level origin of impaired mechanics and higher fracture risks in skeletal disorders such as osteoporosis or Paget’s disease. Alterations in the ultrastructure, chemistry, and nano-/micromechanics of bone tissue in such a diverse group of diseased states have only been briefly explored. Recent research is uncovering the effects of several non-collagenous bone matrix proteins, whose deficiencies or mutations are, to some extent, implicated in bone diseases, on bone matrix quality and mechanics. Herein, we review existing studies on ultrastructural imaging—with a focus on electron microscopy—and chemical, mechanical analysis of pathological bone tissues. The nanometric details offered by these reports, from studying knockout mice models to characterizing exact disease phenotypes, can provide key insights into various bone pathologies and facilitate the development of new treatments.
Collapse
|
21
|
Microscale compressive behavior of hydrated lamellar bone at high strain rates. Acta Biomater 2021; 131:403-414. [PMID: 34245895 DOI: 10.1016/j.actbio.2021.07.005] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2021] [Revised: 06/30/2021] [Accepted: 07/02/2021] [Indexed: 11/20/2022]
Abstract
The increased risk of fracture in the elderly associated with metabolic conditions like osteoporosis poses a significant strain on health care systems worldwide. Due to bone's hierarchical nature, it is necessary to study its mechanical properties and failure mechanisms at several length scales. We conducted micropillar compression experiments on ovine cortical bone to assess the anisotropic mechanical response at the lamellar scale over a wide range of strain rates (10-4 to 8·102 s-1). At the microscale, lamellar bone exhibits a strain rate sensitivity similar to what is reported at the macroscale suggesting that it is an intrinsic property of the extracellular matrix. Significant shear band thickening was observed at high strain rates by HRSEM and STEM imaging. This is likely caused by the material's inability to accommodate the imposed deformation by propagation of thin kink bands and shear cracks at high strain rates, leading to shear band thickening and nucleation. The post-yield behavior is strain rate and direction dependent: hardening was observed for transverse oriented micropillars and hardening modulus increases with strain rate by a factor of almost 2, while axially oriented micropillars showed strain softening and an increase of the softening peak width and work to ultimate stress as a function of strain rate. This suggests that for compression at the micrometer scale, energy absorption in bone increases with strain rate. This study highlights the importance of investigating bone strength and post-yield behavior at lower length scales, under hydrated conditions and at clinically relevant strain rates. STATEMENT OF SIGNIFICANCE: We performed micropillar compression experiments of ovine cortical bone at two different orientations and over seven orders of magnitude of strain rate. Experiments were performed under humid condition to mimic the natural conditions of bone in a human body using a newly developed micro-indenter setup. The strain rate sensitivity was found to be of a similar magnitude to what has been reported for higher length scales, suggesting that the strain rate sensitivity is an intrinsic property of the bone extracellular matrix. In addition, localized shear deformation in thick bands was observed for the first time at high strain rates, highlighting the importance of investigating bone under conditions representative of an accident or fall at several length scales.
Collapse
|
22
|
Casari D, Kochetkova T, Michler J, Zysset P, Schwiedrzik J. Microtensile failure mechanisms in lamellar bone: Influence of fibrillar orientation, specimen size and hydration. Acta Biomater 2021; 131:391-402. [PMID: 34175475 DOI: 10.1016/j.actbio.2021.06.032] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2021] [Revised: 06/17/2021] [Accepted: 06/18/2021] [Indexed: 11/24/2022]
Abstract
A mechanistic understanding of bone fracture is indispensable for developing improved fracture risk assessment in clinics. Since bone is a hierarchically structured material, gaining such knowledge requires analysis at multiple length scales. Here, the tensile response of cortical bone is characterized at the lamellar length scale under dry and hydrated conditions with the aim of investigating the influence of bone's microstructure and hydration on its microscale strength and toughness. For individual lamellae, bone strength strongly correlates with the underlying mineralized collagen fibrils orientation and shows a 2.3-fold increase compared to the macroscale. When specimen size is increased to a few lamellae, the influence of fibril orientation and the size effect on strength are significantly reduced. These findings highlight the critical influence of defects, such as canaliculi and interlamellar interfaces, when assessing larger volumes. Hydration leads up to a 3-fold strength decrease but activates several toughening mechanisms enabling inelastic deformation. In axial specimens, toughening is seen through fibril bridging and crack kinking. In transverse specimens, water presence leads to a progressive but stable crack growth parallel to the fibril orientation, suggesting crack-tip plasticity at the fibrillar interfaces. This work offers a better understanding of the role of interfaces, porosity, and hydration in crack initiation under tensile loading, which is a crucial step towards improved clinical management of disease-related bone fractures through multiscale modeling approaches. STATEMENT OF SIGNIFICANCE: Bone features a complex hierarchical structure which gives rise to several toughening mechanisms across several length scales. To better understand bone fracture, particularly the changes associated with age and disease, it is essential to investigate bone mechanical response at different levels of its hierarchical structure. For the first time, we were able to observe the nucleation of a single crack in hydrated bone lamellae under well-controlled uniaxial tensile loading conditions. These experiments highlight the role of water, interfaces, defects, and the ratio of defect to specimen size on bone's apparent strength and toughness. Such knowledge can be used in the future to develop multiscale models enabling improved clinical management of disease-related bone fractures.
Collapse
|
23
|
Indermaur M, Casari D, Kochetkova T, Peruzzi C, Zimmermann E, Rauch F, Willie B, Michler J, Schwiedrzik J, Zysset P. Compressive Strength of Iliac Bone ECM Is Not Reduced in Osteogenesis Imperfecta and Increases With Mineralization. J Bone Miner Res 2021; 36:1364-1375. [PMID: 33740286 PMCID: PMC8359849 DOI: 10.1002/jbmr.4286] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 07/17/2020] [Revised: 03/03/2021] [Accepted: 03/07/2021] [Indexed: 01/13/2023]
Abstract
Osteogenesis imperfecta (OI) is an inheritable, genetic, and collagen-related disorder leading to an increase in bone fragility, but the origin of its "brittle behavior" is unclear. Because of its complex hierarchical structure, bone behaves differently at various length scales. This study aims to compare mechanical properties of human OI bone with healthy control bone at the extracellular matrix (ECM) level and to quantify the influence of the degree of mineralization. Degree of mineralization and mechanical properties were analyzed under dry conditions in 12 fixed and embedded transiliac crest biopsies (control n = 6, OI type I n = 3, OI type IV n = 2, and OI type III n = 1). Mean degree of mineralization was measured by microcomputed tomography at the biopsy level and the mineral-to-matrix ratio was assessed by Raman spectroscopy at the ECM level. Both methods revealed that the degree of mineralization is higher for OI bone compared with healthy control. Micropillar compression is a novel technique for quantifying post-yield properties of bone at the ECM level. Micropillars (d = 5 μm, h = 10 μm) were fabricated using focused ion beam milling and quasi-statically compressed to capture key post-yield properties such as ultimate strength. The qualitative inspection of the stress-strain curves showed that both OI and healthy control bone have a ductile response at the ECM level. The quantitative results showed that compressive strength is not reduced in OI bone and is increasing with OI severity. Nanoindentation measurements revealed that OI bone tends to have a higher Young's modulus, hardness, and dissipated energy compared with healthy bone. Micropillar strength and indentation modulus increased linearly and significantly (p < .0001) with mineral-to-matrix ratio. In conclusion, this study indicates that compressive mechanical properties of dry OI bone at the iliac crest are not inferior to healthy control at the ECM level and increase with mineralization. © 2021 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).
Collapse
Affiliation(s)
- Michael Indermaur
- ARTORG Center for Biomedical Engineering ResearchUniversity of BernBernSwitzerland
| | - Daniele Casari
- Swiss Federal Laboratories for Material Science and Technology, EmpaThunSwitzerland
| | - Tatiana Kochetkova
- Swiss Federal Laboratories for Material Science and Technology, EmpaThunSwitzerland
| | - Cinzia Peruzzi
- Swiss Federal Laboratories for Material Science and Technology, EmpaThunSwitzerland
| | | | - Frank Rauch
- Shriners Hospital for ChildrenMontrealCanada
- McGill UniversityMontrealCanada
| | - Bettina Willie
- Shriners Hospital for ChildrenMontrealCanada
- McGill UniversityMontrealCanada
| | - Johann Michler
- Swiss Federal Laboratories for Material Science and Technology, EmpaThunSwitzerland
| | - Jakob Schwiedrzik
- Swiss Federal Laboratories for Material Science and Technology, EmpaThunSwitzerland
| | - Philippe Zysset
- ARTORG Center for Biomedical Engineering ResearchUniversity of BernBernSwitzerland
| |
Collapse
|
24
|
Sayin E, Baran ET, Elsheikh A, Mudera V, Cheema U, Hasirci V. Evaluating Oxygen Tensions Related to Bone Marrow and Matrix for MSC Differentiation in 2D and 3D Biomimetic Lamellar Scaffolds. Int J Mol Sci 2021; 22:4010. [PMID: 33924614 PMCID: PMC8068918 DOI: 10.3390/ijms22084010] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2021] [Revised: 04/02/2021] [Accepted: 04/06/2021] [Indexed: 02/06/2023] Open
Abstract
The physiological O2 microenvironment of mesenchymal stem cells (MSCs) and osteoblasts and the dimensionality of a substrate are known to be important in regulating cell phenotype and function. By providing the physiologically normoxic environments of bone marrow (5%) and matrix (12%), we assessed their potential to maintain stemness, induce osteogenic differentiation, and enhance the material properties in the micropatterned collagen/silk fibroin scaffolds that were produced in 2D or 3D. Expression of osterix (OSX) and vascular endothelial growth factor A (VEGFA) was significantly enhanced in the 3D scaffold in all oxygen environments. At 21% O2, OSX and VEGFA expressions in the 3D scaffold were respectively 13,200 and 270 times higher than those of the 2D scaffold. Markers for assessing stemness were significantly more pronounced on tissue culture polystyrene and 2D scaffold incubated at 5% O2. At 21% O2, we measured significant increases in ultimate tensile strength (p < 0.0001) and Young's modulus (p = 0.003) of the 3D scaffold compared to the 2D scaffold, whilst 5% O2 hindered the positive effect of cell seeding on tensile strength. In conclusion, we demonstrated that the 3D culture of MSCs in collagen/silk fibroin scaffolds provided biomimetic cues for bone progenitor cells toward differentiation and enhanced the tensile mechanical properties.
Collapse
Affiliation(s)
- Esen Sayin
- Department of Biotechnology, Middle East Technical University, 06800 Ankara, Turkey;
| | - Erkan Türker Baran
- Department of Tissue Engineering, University of Health Sciences, 34668 Istanbul, Turkey;
| | - Ahmed Elsheikh
- School of Engineering, The University of Liverpool, Liverpool L69 3GH, UK;
| | - Vivek Mudera
- UCL Centre for 3D Models of Health and Disease, Division of Surgery and Interventional Science, University College London, 43-45 Foley Street, Fitzrovia, London W1W 7TY, UK; (V.M.); (U.C.)
| | - Umber Cheema
- UCL Centre for 3D Models of Health and Disease, Division of Surgery and Interventional Science, University College London, 43-45 Foley Street, Fitzrovia, London W1W 7TY, UK; (V.M.); (U.C.)
| | - Vasif Hasirci
- Department of Biotechnology, Middle East Technical University, 06800 Ankara, Turkey;
- Department of Medical Engineering, Acibadem Mehmet Ali Aydinlar University, 34752 Istanbul, Turkey
| |
Collapse
|
25
|
|
26
|
Kochetkova T, Peruzzi C, Braun O, Overbeck J, Maurya AK, Neels A, Calame M, Michler J, Zysset P, Schwiedrzik J. Combining polarized Raman spectroscopy and micropillar compression to study microscale structure-property relationships in mineralized tissues. Acta Biomater 2021; 119:390-404. [PMID: 33122147 DOI: 10.1016/j.actbio.2020.10.034] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2020] [Revised: 10/07/2020] [Accepted: 10/22/2020] [Indexed: 02/07/2023]
Abstract
Bone is a natural composite possessing outstanding mechanical properties combined with a lightweight design. The key feature contributing to this unusual combination of properties is the bone hierarchical organization ranging from the nano- to the macro-scale. Bone anisotropic mechanical properties from two orthogonal planes (along and perpendicular to the main bone axis) have already been widely studied. In this work, we demonstrate the dependence of the microscale compressive mechanical properties on the angle between loading direction and the mineralized collagen fibril orientation in the range between 0° and 82°. For this, we calibrated polarized Raman spectroscopy for quantitative collagen fibril orientation determination and validated the method using widely used techniques (small angle X-ray scattering, micro-computed tomography). We then performed compression tests on bovine cortical bone micropillars with known mineralized collagen fibril angles. A strong dependence of the compressive micromechanical properties of bone on the fibril orientation was found with a high degree of anisotropy for both the elastic modulus (Ea/Et=3.80) and the yield stress (σay/σty=2.54). Moreover, the post-yield behavior was found to depend on the MCF orientation with a transition between softening to hardening behavior at approximately 50°. The combination of methods described in this work allows to reliably determine structure-property relationships of bone at the microscale, which may be used as a measure of bone quality.
Collapse
|
27
|
Wang Y, Ural A. A three-dimensional multiscale finite element model of bone coupling mineralized collagen fibril networks and lamellae. J Biomech 2020; 112:110041. [DOI: 10.1016/j.jbiomech.2020.110041] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2020] [Revised: 08/23/2020] [Accepted: 09/01/2020] [Indexed: 02/06/2023]
|