Search

Biological tissues acquire various characteristic shapes through morphogenesis. Tissue shapes result from the spatiotemporally heterogeneous cellular activities influenced by mechanical and biochemical environments. To investigate multicellular tissue morphogenesis, this study aimed to develop a novel multiscale method that can connect each cellular activity to the mechanical behaviors of the whole tissue by constructing continuum-based particle models of cellular activities. This study proposed mechanical models of cell growth and proliferation that are expressed as volume expansion and cell division by extending the material point method. By simulating cell hypertrophy and proliferation under both free and constraint conditions, the proposed models demonstrated potential for evaluating the mechanical state and tracing cells throughout tissue morphogenesis. Moreover, the effect of a cell size checkpoint was incorporated into the cell proliferation model to investigate the mechanical behaviors of the whole tissue depending on the condition of cellular activities. Consequently, the accumulation of strain energy density was suppressed because of the influence of the checkpoint. In addition, the whole tissues acquired different shapes depending on the influence of the checkpoint. Thus, the models constructed herein enabled us to investigate the change in the mechanical behaviors of the whole tissue according to each cellular activity depending on the mechanical state of the cells during morphogenesis., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Size control in biological tissues involves multicellular communication via mechanical forces during development. Although fundamental cellular behaviours in response to mechanical stimuli underlie size maintenance during morphogenetic processes, the mechanisms underpinning the cellular mechano-response system that maintains size along an axis of a polarized tissue remain elusive. Here, we show how the diameter of an epithelial tube is maintained during murine epididymal development by combining quantitative imaging, mechanical perturbation and mathematical modelling. We found that epithelial cells counteract compressive forces caused by cell division exclusively along the circumferential axis of the tube to produce polarized contractile forces, eventually leading to an oriented cell rearrangement. Moreover, a mathematical model that includes the polarized mechano-responsive regime explains how the diameter of proliferating tubes is maintained. Our findings pave the way for an improved understanding of the cellular response to mechanical forces that involves collective multicellular behaviours for organizing diverse tissue morphologies., Competing Interests: Competing interestsThe authors declare no competing or financial interests., (© 2019. Published by The Company of Biologists Ltd.)

Cells generate mechanical forces primarily from interactions between F-actin, cross-linking proteins, myosin motors, and other actin-binding proteins in the cytoskeleton. To understand how molecular interactions between the cytoskeletal elements generate forces, a number of in vitro experiments have been performed but are limited in their ability to accurately reproduce the diversity of motor mobility. In myosin motility assays, myosin heads are fixed on a surface and glide F-actin. By contrast, in reconstituted gels, the motion of both myosin and F-actin is unrestricted. Because only these two extreme conditions have been used, the importance of mobility of motors for network behaviors has remained unclear. In this study, to illuminate the impacts of motor mobility on the contractile behaviors of the actin cytoskeleton, we employed an agent-based computational model based on Brownian dynamics. We find that if motors can bind to only one F-actin like myosin I, networks are most contractile at intermediate mobility. In this case, less motor mobility helps motors stably pull F-actins to generate tensile forces, whereas higher motor mobility allows F-actins to aggregate into larger clustering structures. The optimal intermediate motor mobility depends on the stall force and affinity of motors that are regulated by mechanochemical rates. In addition, we find that the role of motor mobility can vary drastically if motors can bind to a pair of F-actins. A network can exhibit large contraction with high motor mobility because motors bound to antiparallel pairs of F-actins can exert similar forces regardless of their mobility. Results from this study imply that the mobility of molecular motors may critically regulate contractile behaviors of actin networks in cells., (Copyright © 2019 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

This study demonstrates computational simulations of multicellular deformation coupled with chemical patterning in the three-dimensional (3D) space. To address these aspects, we proposes a novel mathematical model, where a reaction-diffusion system is discretely expressed at a single cell level and combined with a 3D vertex model. To investigate complex phenomena emerging from the coupling of patterning and deformation, as an example, we employed an activator-inhibitor system and converted the activator concentration of individual cells into their growth rate. Despite the simplicity of the model, by growing a monolayer cell vesicle, the coupling system provided rich morphological dynamics such as undulation, tubulation, and branching. Interestingly, the morphological variety depends on the difference in time scales between patterning and deformation, and can be partially understood by the intrinsic hysteresis in the activator-inhibitor system with domain growth. Importantly, the model can be applied to 3D multicellular dynamics that couple the reaction-diffusion patterning with various cell behaviors, such as deformation, rearrangement, division, apoptosis, differentiation, and proliferation. Thus, the results demonstrate the significant advantage of the proposed model as well as the biophysical importance of exploring spatiotemporal dynamics of the coupling phenomena of patterning and deformation in 3D space.

To understand Wolff's law, bone adaptation by remodeling at the cellular and tissue levels has been discussed extensively through experimental and simulation studies. For the clinical application of a bone remodeling simulation, it is significant to establish a macroscopic model that incorporates clarified microscopic mechanisms. In this study, we proposed novel macroscopic models based on the microscopic mechanism of osteocytic mechanosensing, in which the flow of fluid in the lacuno-canalicular porosity generated by fluid pressure gradients plays an important role, and theoretically evaluated the proposed models, taking biological rationales of bone adaptation into account. The proposed models were categorized into two groups according to whether the remodeling equilibrium state was defined globally or locally, i.e., the global or local uniformity models. Each remodeling stimulus in the proposed models was quantitatively evaluated through image-based finite element analyses of a swine cancellous bone, according to two introduced criteria associated with the trabecular volume and orientation at remodeling equilibrium based on biological rationales. The evaluation suggested that nonuniformity of the mean stress gradient in the local uniformity model, one of the proposed stimuli, has high validity. Furthermore, the adaptive potential of each stimulus was discussed based on spatial distribution of a remodeling stimulus on the trabecular surface. The theoretical consideration of a remodeling stimulus based on biological rationales of bone adaptation would contribute to the establishment of a clinically applicable and reliable simulation model of bone remodeling.

Epithelial invagination is one of the fundamental deformation modes during morphogenesis, and is essential for deriving the three-dimensional shapes of organs from a flat epithelial sheet. Invagination occurs in an orderly manner according to the spatial pattern of the contractile cells; however, it remains elusive how tissue deformation can be caused by cellular activity in the patterned region. In this study, we investigated the mechanical role of the spatial patterns of the contractile cells in invagination of growing tissue using multicellular dynamics simulations. We found that cell proliferation and apical constriction were responsible for expanding the degree of tissue deformation and determining the location of the deformation, respectively. The direction of invagination depended on the spatial pattern of the contractile cells. Further, comparing the simulation results of surface and line contractions as possible modes of apical constriction, we found that the direction of invagination differed between these two modes even if the spatial pattern was the same. These results indicate that the buckling of the epithelial cell sheet caused by cell proliferation causes the invagination, with the direction and location determined by the configuration of the wedge-shaped cells given by the spatial pattern of the contractile cells., (© 2017 The Authors Development, Growth & Differentiation published by John Wiley & Sons Australia, Ltd on behalf of Japanese Society of Developmental Biologists.)

Morphogenesis in multicellular organisms is accompanied by apoptotic cell behaviors: cell shrinkage and cell disappearance. The mechanical effects of these behaviors are spatiotemporally regulated within multicellular dynamics to achieve proper tissue sizes and shapes in three-dimensional (3D) space. To analyze 3D multicellular dynamics, 3D vertex models have been suggested, in which a reversible network reconnection (RNR) model has successfully expressed 3D cell rearrangements during large deformations. To analyze the effects of apoptotic cell behaviors on 3D multicellular morphogenesis, we modeled cell apoptosis based on the RNR model framework. Cell shrinkage was modeled by the potential energy as a function of individual cell times during the apoptotic phase. Cell disappearance was modeled by merging neighboring polyhedrons at their boundary surface according to the topological rules of the RNR model. To establish that the apoptotic cell behaviors could be expressed as modeled, we simulated morphogenesis driven by cell apoptosis in two types of tissue topology: 3D monolayer cell sheet and 3D compacted cell aggregate. In both types of tissue topology, the numerical simulations successfully illustrated that cell aggregates gradually shrank because of successive cell apoptosis. During tissue shrinkage, the number of cells in aggregates decreased while maintaining individual cell size and shape. Moreover, in case of localizing apoptotic cells within a part of the 3D monolayer cell aggregate, the cell apoptosis caused the global tissue bending by pulling on surrounding cells. In case of localizing apoptotic cells on the surface of the 3D compacted cell aggregate, the cell apoptosis caused successive, directional cell rearrangements from the inside to the surface. Thus, the proposed model successfully provided a basis for expressing apoptotic cell behaviors during 3D multicellular morphogenesis based on an RNR model framework.

Tissue morphogenesis in multicellular organisms is accompanied by proliferative cell behaviors: cell division (increase in cell number after each cell cycle) and cell growth (increase in cell volume during each cell cycle). These proliferative cell behaviors can be regulated by multicellular dynamics to achieve proper tissue sizes and shapes in three-dimensional (3D) space. To analyze multicellular dynamics, a reversible network reconnection (RNR) model has been suggested, in which each cell shape is expressed by a single polyhedron. In this study, to apply the RNR model to simulate tissue morphogenesis involving proliferative cell behaviors, we model cell proliferation based on a RNR model framework. In this model, cell division was expressed by dividing a polyhedron at a planar surface for which cell division behaviors were characterized by three quantities: timing, intracellular position, and normal direction of the dividing plane. In addition, cell growth was expressed by volume growth as a function of individual cell times within their respective cell cycles. Numerical simulations using the proposed model showed that tissues grew during successive cell divisions with several cell cycle times. During these processes, the cell number in tissues increased while maintaining individual cell size and shape. Furthermore, tissue morphology dramatically changed based on different regulations of cell division directions. Thus, the proposed model successfully provided a basis for expressing proliferative cell behaviors during morphogenesis based on a RNR model framework.

Morphogenesis of tissues in organ development is accompanied by large three-dimensional (3D) deformations, in which mechanical interactions among multiple cells are spatiotemporally regulated. To reveal the deformation mechanisms, in this study, we developed the reversible network reconnection (RNR) model. The model is developed on the basis of 3D vertex model, which expresses a multicellular aggregate as a network composed of vertices. 3D vertex models have successfully simulated morphogenetic dynamics by expressing cellular rearrangements as network reconnections. However, the network reconnections in 3D vertex models can cause geometrical irreversibility, energetic inconsistency, and topological irreversibility, therefore inducing unphysical results and failures in simulating large deformations. To resolve these problems, we introduced (1) a new definition of the shapes of polygonal faces between cellular polyhedrons, (2) an improved condition for network reconnections, (3) a new condition for potential energy functions, and (4) a new constraint condition for the shapes of polygonal faces that represent cell-cell boundaries. Mathematical and computational analyses demonstrated that geometrical irreversibility, energetic inconsistency, and topological irreversibility were resolved by suppressing the geometrical gaps in the network and avoiding the generation of irreversible network patterns in reconnections. Lastly, to demonstrate the applicability of the RNR model, we simulated tissue deformation of growing cell sheets and showed that our model can simulate large tissue deformations, in which large changes occur in the local curvatures and layer formations of tissues. Thus, the RNR model enables in silico recapitulation of complex tissue morphogenesis.

In general, a rapidly growing epithelial sheet during tissue morphogenesis shows a smooth and continuous curvature on both inner cavity (apical) and basement membrane (basal) sides. For instance, epithelia of the neural tube and optic vesicle in the early embryo maintain continuous curvatures in their local domains, even during their rapid growth. However, given that cell divisions, which substantially perturb the local force balance, frequently and successively occur in an uncoordinated manner, it is not self-evident to explain how the tissue keeps a continuous curvature at large. In the majority of developing embryonic epithelia with smooth surfaces, their curvatures are apically concave, because of the presence of strong tangential contractile force on the apical side. In this numerical study, we demonstrate that tangential contractile forces on the apical surface play a critical role in the maintenance of smooth curvatures in the epithelium and reduce irregular undulations caused by uncoordinated generation of local pushing force. Using a reversible network reconnection (RNR) model, which we previously developed to make numerical analyses highly reproducible even under rapid tissue-growth conditions, we performed simulations for morphodynamics to examine the effect of apical contractile forces on the continuity of curvatures. Interestingly, the presence of apical contractile forces suppressed irregular undulations not only on the apical side but also on the basal surface. These results indicate that cellular contractile forces on the apical surface control not only the shape at a single cell level but also at a tissue level as a result of emergent mechanical coordination., (Copyright © 2013 Elsevier Ltd. All rights reserved.)

The process of bone remodeling is regulated by metabolic activities of many bone cells. While osteoclasts and osteoblasts are responsible for bone resorption and formation, respectively, activities of these cells are believed to be controlled by a mechanosensory system of osteocytes embedded in the extracellular bone matrix. Several experimental and theoretical studies have suggested that the strain-derived interstitial fluid flow in lacuno-canalicular porosity serves as the prime mover for bone remodeling. Previously, we constructed a mathematical model for trabecular bone remodeling that interconnects the microscopic cellular activities with the macroscopic morphological changes in trabeculae through the mechanical hierarchy. This model assumes that fluid-induced shear stress acting on osteocyte processes is a driving force for bone remodeling. The validity of this model has been demonstrated with a remodeling simulation using a two-dimensional trabecular model. In this study, to investigate the effects of loading frequency, which is thought to be a significant mechanical factor in bone remodeling, we simulated morphological changes of a three-dimensional single trabecula under cyclic uniaxial loading with various frequencies. The results of the simulation show the trabecula reoriented to the loading direction with the progress of bone remodeling. Furthermore, as the imposed loading frequency increased, the diameter of the trabecula in the equilibrium state was enlarged by remodeling. These results indicate that our simulation model can successfully evaluate the relationship between loading frequency and trabecular bone remodeling., (Copyright © 2011 Elsevier Ltd. All rights reserved.)

Actomyosin contractility is a major force-generating mechanism that drives rearrangement of actomyosin networks; it is fundamental to cellular functions such as cellular reshaping and movement. Thus, to clarify the mechanochemical foundation of the emergence of cellular functions, understanding the relationship between actomyosin contractility and rearrangement of actomyosin networks is crucial. For this purpose, in this study, we present a new particulate-based model for simulating the motions of actin, non-muscle myosin II, and α-actinin. To confirm the model's validity, we successfully simulated sliding and bending motions of actomyosin filaments, which are observed as fundamental behaviors in dynamic rearrangement of actomyosin networks in migrating keratocytes. Next, we simulated the dynamic rearrangement of actomyosin networks. Our simulation results indicate that an increase in the density fraction of myosin induces a higher-order structural transition of actomyosin filaments from networks to bundles, in addition to increasing the force generated by actomyosin filaments in the network. We compare our simulation results with experimental results and confirm that actomyosin bundles bridging focal adhesions and the characteristics of myosin-dependent rearrangement of actomyosin networks agree qualitatively with those observed experimentally., (Copyright © 2011 Elsevier Ltd. All rights reserved.)

Membrane protrusion is a mechanochemical process of active membrane deformation driven by actin polymerization. Previously, Brownian ratchet (BR) was modeled on the basis of the underlying molecular mechanism. However, because the BR requires a priori load that cannot be determined without information of the cell shape, it cannot be effective in studies in which resultant shapes are to be solved. Other cellular-scale models describing the protrusion have also been suggested for modeling a whole cell; however, these models were not developed on the basis of coarse-grained physics representing the underlying molecular mechanism. Therefore, to express the membrane protrusion on the cellular scale, we propose a novel mathematical model, the coarse-grained BR (CBR), which is derived on the basis of nonequilibrium thermodynamics theory. The CBR can reproduce the BR within the limit of the quasistatic process of membrane protrusion and can estimate the protrusion velocity consistently with an effective elastic constant that represents the state of the energy of the membrane. Finally, to demonstrate the applicability of the CBR, we attempt to perform a cellular-scale simulation of migrating keratocyte in which the proposed CBR is used for the membrane protrusion model on the cellular scale. The results show that the experimentally observed shapes of the leading edge are well reproduced by the simulation. In addition, The trend of dependences of the protrusion velocity on the curvature of the leading edge, the temperature, and the substrate stiffness also agreed with the other experimental results. Thus, the CBR can be considered an appropriate cellular-scale model to express the membrane protrusion on the basis of its underlying molecular mechanism.

In bone functional adaptation by remodelling, osteocytes in the lacuno-canalicular system are believed to play important roles in the mechanosensory system. Under dynamic loading, bone matrix deformation generates an interstitial fluid flow in the lacuno-canalicular system; this flow induces shear stress on the osteocytic process membrane that is known to stimulate the osteocytes. In this sense, the osteocytes behave as mechanosensors and deliver mechanical information to neighbouring cells through the intercellular communication network. In this study, bone remodelling is assumed to be regulated by the mechanical signals collected by the osteocytes. From the viewpoint of multi-scale biomechanics, we propose a mathematical model of trabecular bone remodelling that takes into account the osteocytic mechanosensory network system. Based on this model, a computational simulation of trabecular bone remodelling was conducted for a single trabecula under cyclic uniaxial loading, demonstrating functional adaptation to the applied mechanical loading as a load-bearing construct.

Load-induced interstitial fluid flow in lacuno-canalicular porosity is believed to play an important role in cellular activities regulating adaptive bone remodeling. To investigate interstitial fluid behavior based on poroelasticity, it is important to determine the anisotropic permeability tensor reflecting the morphological features of the lacuno-canalicular porosity as fluid channels. In this study, we presented an estimation method of trabecular permeability by describing the analytical relationship between the volume orientation (VO) fabric tensor, which represents the canalicular orientation, and the permeability tensor. The relationship showed that the trabecular permeability tensor is proportional to the product of the volume fraction of the interstitial fluid and the VO fabric tensor of the canaliculi. We applied the proposed method to a two-dimensional fluorescent image of a trabecular cross section to quantify the canalicular anisotropy and the trabecular permeability tensor. The results indicated that the canaliculi are predominantly oriented in the radial direction of the trabecula, and the permeability depends strongly on the canalicular morphology., ((c) 2009 Elsevier Ltd. All rights reserved.)

The actin filament, which is the most abundant component of the cytoskeleton, plays important roles in fundamental cellular activities such as shape determination, cell motility, and mechanosensing. In each activity, the actin filament dynamically changes its structure by polymerization, depolymerization, and severing. These phenomena occur on the scales ranging from the dynamics of actin molecules to filament structural changes with its deformation due to the various forces, for example, by the membrane and solvent. To better understand the actin filament dynamics, it is important to focus on these scales and develop its mathematical model. Thus, the objectives of this study were to model and simulate actin filament polymerization, depolymerization, and severing based on the Brownian dynamics method. In the model, the actin monomers and the solvent were considered as globular particles and a continuum, respectively. The motion of the actin molecules was assumed to follow the Langevin equation. The polymerization, which increases the filament length, was determined by the distance between the center of the actin particle at the barbed end and actin particles in the solvent. The depolymerization, which decreases the filament length, was modeled such that the number of dissociation particles from the filament end per unit time was constant. In addition, the filament severing, in which one filament divides into two, was modeled to occur at an equal rate along the filament. Then, we simulated the actin filament dynamics using the developed model, and analyzed the filament elongation rate, its turnover, and the effects of filament severing on the polymerization and depolymerization. Results indicated that the model reproduced the linear dependence of the filament elongation on time, filament turnover process by polymerization and depolymerization, and acceleration of the polymerization and depolymerization by severing, which qualitatively agreed with those observed in experiments.

In adaptive bone remodeling, it is believed that bone cells such as osteoblasts, osteocytes and osteoclasts can sense mechanical stimuli and modulate their remodeling activities. However, the mechanosensing mechanism by which these cells sense mechanical stimuli and transduce mechanical signals into intracellular biochemical signals is still not clearly understood. From the viewpoint of cell biomechanics, it is important to clarify the mechanical conditions under which the cellular mechanosensing mechanism is activated. The aims of this study were to evaluate a mechanical condition, that is, the local strain on the cell membrane, at the initiation point of the intracellular calcium signaling response to the applied mechanical stimulus in osteoblast-like MC3T3-E1 cells, and to investigate the effect of deformation velocity on the characteristics of the cellular response. To apply a local deformation to a single cell, a glass microneedle was directly indented to the cell and moved horizontally on the cell membrane. To observe the cellular response and the deformation of the cell membrane, intracellular calcium ions and the cell membrane were labeled using fluorescent dyes and simultaneously observed by confocal laser scanning microscopy. The strain distribution on the cell membrane attributable to the applied local deformation and the strain magnitude at the initiation point of the calcium signaling responses were analyzed using obtained fluorescence images. From two-dimensionally projected images, it was found that there is a local compressive strain at the initiation point of calcium signaling. Moreover, the cellular response revealed velocity dependence, that is, the cells seemed to respond with a higher sensitivity to a higher deformation velocity. From the viewpoint of cell biomechanics, these results provide us a fundamental understanding of the mechanosensing mechanism of osteoblast-like cells.

In bone tissue engineering using a biodegradable scaffold, geometry of the porous scaffold microstructure is a key factor for controlling mechanical function of the bone-scaffold system in the regeneration process as well as after the regeneration. In this study, we propose a framework for the optimal design of the porous scaffold microstructure by three-dimensional computational simulation of bone tissue regeneration that consists of scaffold degradation and new bone formation. The rate of scaffold degradation due to hydrolysis, that leads to decrease in mechanical properties, was simply assumed to relate to the water content diffused from the surface to the bulk material. For the new bone formation on both bone and scaffold surfaces, the rate equation of trabecular surface remodeling driven by mechanical stimulation was applied. Solving these two phenomena in the same time frame, the bone regeneration process in the bone-scaffold system was predicted by computational simulation using a voxel finite element method. The change in the mechanical function of the bone-scaffold system during the regeneration process was quantitatively evaluated by measuring the change in total strain energy, and this was used for the evaluation function to optimize the scaffold microstructure that provides the desired mechanical function during and after the bone regeneration process. A case study conducted for the scaffold with a simple microstructure demonstrated that the proposed simulation method could be applied to the design of a porous scaffold microstructure. In addition, the regeneration process was found to be very complex even though the simple rate equations for scaffold regeneration and new bone formation were used because of the coupling effects of these phenomena.

Fabric and compliance tensors of a cube of cancellous bone with a complicated three-dimensional trabecular structure were obtained for trabecular surface remodeling by using a digital image-based model combined with a large-scale finite element method. Using mean intercept length and a homogenization method, the fabric and compliance tensors were determined for the trabecular structure obtained in the computer remodeling simulation. The tensorial quantities obtained indicated that anisotropic structural changes occur in cancellous bone adapting to the compressive loading condition. There were good correlations between the fabric tensor, bone volume fraction, and compliance tensor in the remodeling process. The result demonstrates that changes in the structural and mechanical properties of cancellous bone are essentially anisotropic and should be expressed by tensorial quantities.

Computational simulation of trabecular surface remodeling was conducted to investigate the effects of a spinal fixation screw on trabecular structural changes in a vertebral body. By using voxel-based finite elements, computational models of the bone and screw were constructed in two structural scales of a vertebral body with an implanted screw and a bone-screw interface. In the vertebral body, the implantation of the fixation screw caused changes in the mechanical environment in cancellous bone, leading to trabecular structural changes at the cancellous level. The effects of the screw on trabecular orientation were greater in the regions above and below the screw than in those in front of the screw. In the case of the bone-screw interface, trabecular structural changes depended on the direction of load applied to the screw. It was suggested that the bone resorption predicted in the pull-out loading case is a candidate cause of the loosening of the screw. The results indicate that the effects of the implanted screw on trabecular structural changes are more important for longer-term fixation.

Two-dimensional simulation of trabecular surface remodeling was conducted for a human proximal femur to investigate the structural change of cancellous bone toward a uniform stress state. Considering that a local mechanical stimulus plays an important role in cellular activities in bone remodeling, local stress nonuniformity was assumed to drive trabecular structural change to seek a uniform stress state. A large-scale pixel-based finite element model was used to simulate structural changes of individual trabeculae over the entire bone. As a result, the initial structure of trabeculae changed from isotropic to anisotropic due to trabecular microstructural changes caused by surface remodeling according to the mechanical environment in the proximal femur. Under a single-loading condition, it was shown that the apparent structural property evaluated by fabric ellipses corresponded to the apparent stress state in cancellous bone. As is observed in the actual bone, a distributed trabecular structure was obtained under a multiple-loading condition. Through these studies, it was concluded that trabecular surface remodeling toward a local uniform stress state at the trabecular level could naturally bring about functional adaptation phenomenon at the apparent tissue level. The proposed simulation model would be capable of providing insight into the hierarchical mechanism of trabecular surface remodeling at the microstructural level up to the apparent tissue level.