Skeletal muscle accumulates intramuscular fat through age and obesity. Muscle quality, a measure of muscle strength per unit size, decreases in these conditions. It is not clear how fat influences this loss in performance. Changes to structural parameters (e.g. fibre pennation and connective tissue properties) affect the muscle quality. This study investigated the mechanisms that lead to deterioration in muscle performance due to changes in intramuscular fat, pennation and aponeurosis stiffness. A finite-element model of the human gastrocnemius was developed as a fibre-reinforced composite biomaterial containing contractile fibres within the base material. The base-material properties were modified to include intramuscular fat in five different ways. All these models with fat generated lower fibre stress and muscle quality than their lean counterparts. This effect is due to the higher stiffness of the tissue in the fatty models. The fibre deformations influence their interactions with the aponeuroses, and these change with fatty inclusions. Muscles with more compliant aponeuroses generated lower forces. The muscle quality was further reduced for muscles with lower pennation. This study shows that whole-muscle force is dependent on its base-material properties and changes to the base material due to fatty inclusions result in reductions to force and muscle quality.

During bone remodeling, bone-resorbing osteoclasts and bone-forming osteoblasts are organized in bone multicellular units (BMUs), which travel at a rate of 20-40 mm/d for 6-12 mo, maintaining a cylindrical structure. However, the interplay of local BMU geometry with biochemical regulation is poorly understood. We developed a mathematical model of BMU describing changes in time and space of the concentrations of proresorptive cytokine RANKL and its inhibitor osteoprotegerin (OPG), in osteoclast and osteoblast numbers, and in bone mass. We assumed that osteocytes surrounding a microfracture produce RANKL, which attracted osteoclasts. OPG and RANKL were produced by osteoblasts and diffused through bone, RANKL was eliminated by binding to OPG and RANK. Osteoblasts were coupled to osteoclasts through paracrine factors. The evolution of the BMU arising from this model was studied using numerical simulations. Our model recapitulated the spatio-temporal dynamics observed in vivo in a cross-section of bone. In response to a RANKL field, osteoclasts moved as a well-confined cutting cone. The coupling of osteoclasts to osteoblasts allowed for sufficient recruitment of osteoblasts to the resorbed surfaces. The RANKL field was the highest at the microfracture in front of the BMU, whereas the OPG field peaked at the back of the BMU, resulting in the formation of a RANKL/OPG gradient, which strongly affected the rate of BMU progression and its size. Thus, the spatial organization of a BMU provides important constraints on the roles of RANKL and OPG as well as possibly other regulators in determining the outcome of remodeling in the BMU.

During contraction the energy of muscle tissue increases due to energy from the hydrolysis of ATP. This energy is distributed across the tissue as strain-energy potentials in the contractile elements, strain-energy potential from the 3D deformation of the base-material tissue (containing cellular and extracellular matrix effects), energy related to changes in the muscle's nearly incompressible volume and external work done at the muscle surface. Thus, energy is redistributed through the muscle's tissue as it contracts, with only a component of this energy being used to do mechanical work and develop forces in the muscle's longitudinal direction. Understanding how the strain-energy potentials are redistributed through the muscle tissue will help enlighten why the mechanical performance of whole muscle in its longitudinal direction does not match the performance that would be expected from the contractile elements alone. Here we demonstrate these physical effects using a 3D muscle model based on the finite element method. The tissue deformations within contracting muscle are large, and so the mechanics of contraction were explained using the principles of continuum mechanics for large deformations. We present simulations of a contracting medial gastrocnemius muscle, showing tissue deformations that mirror observations from magnetic resonance imaging. This paper tracks the redistribution of strain-energy potentials through the muscle tissue during fixed-end contractions, and shows how fibre shortening, pennation angle, transverse bulging and anisotropy in the stress and strain of the muscle tissue are all related to the interaction between the material properties of the muscle and the action of the contractile elements.

Geometric integration theory can be employed when numerically solving ODEs or PDEs with constraints. In this paper, we present several one-step algorithms of various orders for ODEs on a collection of spheres. To demonstrate the versatility of these algorithms, we present representative calculations for reduced free rigid body motion (a conservative ODE) and a discretization of micromagnetics (a dissipative PDE). We emphasize the role of isotropy in geometric integration and link numerical integration schemes to modern differential geometry through the use of partial connection forms; this theoretical framework generalizes moving frames and connections on principal bundles to manifolds with nonfree actions.Comment: This paper appeared in prin

Skeletal muscle can contain neuromuscular compartments that are spatially distinct regions that can receive relatively independent levels of activation. This study tested how the magnitude and direction of the force developed by a whole muscle would change when the muscle activity was regionalized within the muscle. A 3D finite element model of a muscle with its bounding aponeurosis was developed for the lateral gastrocnemius, and isometric contractions were simulated for a series of conditions with either a uniform activation pattern, or regionally distinct activation patterns: in all cases the mean activation from all fibers within the muscle reached 10%. The models showed emergent features of the fiber geometry that matched physiological characteristics: with fibers shortening, rotating to greater pennation, adopting curved trajectories in 3D and changes in the thickness and width of the muscle belly. Simulations were repeated for muscle with compliant, normal and stiff aponeurosis and the aponeurosis stiffness affected the changes to the fiber geometry and the resultant muscle force. Changing the regionalization of the activity resulted to changes in the magnitude, direction and center of the force vector from the whole muscle. Regionalizing the muscle activity resulted in greater muscle force than the simulation with uniform activity across the muscle belly. The study shows how the force from a muscle depends on the complex interactions between the muscle fibers and connective tissues and the region of muscle that is active.

Abstract. The mechanical properties of vertebrate bone are largely determined by a process which involves the complex interplay of three different cell types. This process is called bone remodeling, and occurs asynchronously at multiple sites in the mature skeleton. The cells involved are bone resorbing osteoclasts, bone matrix producing osteoblasts and mechanosensing osteocytes. These cells communicate with each other by means of autocrine and paracrine signaling factors and operate in complex entities, the so-called bone multicellular units (BMU). To investigate the BMU dynamics in silico, we develop a novel mathematical model resulting in a system of nonlinear partial differential equations with time delays. The model describes the osteoblast and osteoclast populations together with the dynamics of the key messenger molecule RANKL and its decoy receptor OPG. Scaling theory is used to address parameter sensitivity and predict the emergence of pathological remodeling regimes. The model is studied numerically in one and two space dimensions using finite difference schemes in space and explicit delay equation solvers in time. The computational results are in agreement with in vivo observations and provide new insights into the role of the RANKL/OPG pathway in the spatial regulation of bone remodeling.

No abstract

In this article we consider regularizations of the Dirac delta distribution with applications to prototypical elliptic and hyperbolic partial differential equations (PDEs). We study the convergence of a sequence of distributions S H to a singular term S as a parameter H (associated with the support size of S H ) shrinks to zero. We characterize this convergence in both the weak- * topology of distributions, as well as in a weighted Sobolev norm. These notions motivate a framework for constructing regularizations of the delta distribution that includes a large class of existing methods in the literature. This framework allows different regularizations to be compared. The convergence of solutions of PDEs with these regularized source terms is then studied in various topologies such as pointwise convergence on a deleted neighborhood and weighted Sobolev norms. We also examine the lack of symmetry in tensor product regularizations and effects of dissipative error in hyperbolic problems.