N-site Phosphorylation Systems with 2N-1 Steady States

Multisite protein phosphorylation plays a prominent role in intracellular processes like signal transduction, cell-cycle control and nuclear signal integration. Many proteins are phosphorylated in a sequential and distributive way at more than one phosphorylation site. Mathematical models of \(n\) -site sequential distributive phosphorylation are therefore studied frequently. In particular, in Wang and Sontag (J Math Biol 57:29–52, 2008), it is shown that models of \(n\) -site sequential distributive phosphorylation admit at most \(2n-1\) steady states. Wang and Sontag furthermore conjecture that for odd \(n\) , there are at most \(n\) and that, for even \(n\) , there are at most \(n+1\) steady states. This, however, is not true: building on earlier work in Holstein et al. (Bull Math Biol 75(11):2028–2058, 2013), we present a scalar determining equation for multistationarity which will lead to parameter values where a \(3\) -site system has \(5\) steady states and parameter values where a \(4\) -site system has \(7\) steady states. Our results therefore are counterexamples to the conjecture of Wang and Sontag. We furthermore study the inherent geometric properties of multistationarity in \(n\) -site sequential distributive phosphorylation: the complete vector of steady state ratios is determined by the steady state ratios of free enzymes and unphosphorylated protein and there exists a linear relationship between steady state ratios of phosphorylated protein.     

A thermomechanical framework for reconciling the effects of ultraviolet radiation exposure time and wavelength on connective tissue elasticity

Augmentation of the mechanical properties of connective tissue using ultraviolet (UV) radiation—by targeting collagen cross-linking in the tissue at predetermined UV exposure time \((t)\) and wavelength \((\lambda )\) —has been proposed as a therapeutic method for supporting the treatment for structural-related injuries and pathologies. However, the effects of \(\lambda \) and \(t\) on the tissue elasticity, namely elastic modulus \((E)\) and modulus of resilience \((u_\mathrm{Y})\) , are not entirely clear. We present a thermomechanical framework to reconcile the \(t\) - and \(\lambda \) -related effects on \(E\) and \(u_\mathrm{Y}\) . The framework addresses (1) an energy transfer model to describe the dependence of the absorbed UV photon energy, \(\xi \) , per unit mass of the tissue on \(t\) and \(\lambda \) , (2) an intervening thermodynamic shear-related parameter, \(G\) , to quantify the extent of UV-induced cross-linking in the tissue, (3) a threshold model for the \(G\) versus \(\xi \) relationship, characterized by   \(t_\mathrm{C}\) —the critical \(t\) underpinning the association of \(\xi \) with \(G\) —and (4) the role of \(G\) in the tissue elasticity. We hypothesized that \(G\) regulates \(E\) (UV-stiffening hypothesis) and \(u_\mathrm{Y}\) (UV-resilience hypothesis). The framework was evaluated with the support from data derived from tensile testing on isolated ligament fascicles, treated with two levels of \(\lambda \) (365 and 254 nm) and three levels of \(t\) (15, 30 and 60 min). Predictions from the energy transfer model corroborated the findings from a two-factor analysis of variance of the effects of \(t\) and \(\lambda \) treatments. Student’s t test revealed positive change in \(E\) and \(u_\mathrm{Y}\) with increases in \(G\) —the findings lend support to the hypotheses, implicating the implicit dependence of UV-induced cross-links on \(t\) and \(\lambda \) for directing tissue stiffness and resilience. From a practical perspective, the study is a step in the direction to establish a UV irradiation treatment protocol for effective control of exogenous cross-linking in connective tissues.     

Approximation of sojourn-times via maximal couplings: motif frequency distributions

Sojourn-times provide a versatile framework to assess the statistical significance of motifs in genome-wide searches even under non-Markovian background models. However, the large state spaces encountered in genomic sequence analyses make the exact calculation of sojourn-time distributions computationally intractable in long sequences. Here, we use coupling and analytic combinatoric techniques to approximate these distributions in the general setting of Polish state spaces, which encompass discrete state spaces. Our approximations are accompanied with explicit, easy to compute, error bounds for total variation distance. Broadly speaking, if \({\mathsf{T}}_n\) is the random number of times a Markov chain visits a certain subset \({\mathsf{T}}\) of states in its first \(n\) transitions, then we can usually approximate the distribution of \({\mathsf{T}}_n\) for \(n\) of order \((1-\alpha )^{-m}\) , where \(m\) is the largest integer for which the exact distribution of \({\mathsf{T}}_m\) is accessible and \(0\le \alpha \le 1\) is an ergodicity coefficient associated with the probability transition kernel of the chain. This gives access to approximations of sojourn-times in the intermediate regime where \(n\) is perhaps too large for exact calculations, but too small to rely on Normal approximations or stationarity assumptions underlying Poisson and compound Poisson approximations. As proof of concept, we approximate the distribution of the number of matches with a motif in promoter regions of C. elegans. Mathematical properties of the proposed ergodicity coefficients and connections with additive functionals of homogeneous Markov chains as well as ergodicity of non-homogeneous Markov chains are also explored.     

Pareto genealogies arising from a Poisson branching evolution model with selection

We study a class of coalescents derived from a sampling procedure out of $N$ i.i.d. Pareto $\left( \alpha \right) $ random variables, normalized by their sum, including $\beta $ –size-biasing on total length effects ( $\beta <\alpha $ ). Depending on the range of $\alpha ,$ we derive the large $N$ limit coalescents structure, leading either to a discrete-time Poisson-Dirichlet $ \left(\alpha ,-\beta \right) \Xi -$ coalescent ( $\alpha \in \left[ 0,1\right) $ ), or to a family of continuous-time Beta $\left( 2-\alpha ,\alpha -\beta \right) \Lambda -$ coalescents ( $\alpha \in \left[ 1,2\right) $ ), or to the Kingman coalescent ( $\alpha \ge 2$ ). We indicate that this class of coalescent processes (and their scaling limits) may be viewed as the genealogical processes of some forward in time evolving branching population models including selection effects. In such constant-size population models, the reproduction step, which is based on a fitness-dependent Poisson Point Process with scaling power-law $\left( \alpha \right) $ intensity, is coupled to a selection step consisting of sorting out the $N$ fittest individuals issued from the reproduction step.     

Approximating the dynamics of communicating cells in a diffusive medium by ODEs—homogenization with localization

Bacteria may change their behavior depending on the population density. Here we study a dynamical model in which cells of radius $R$ within a diffusive medium communicate with each other via diffusion of a signalling substance produced by the cells. The model consists of an initial boundary value problem for a parabolic PDE describing the exterior concentration $u$ of the signalling substance, coupled with $N$ ODEs for the masses $a_i$ of the substance within each cell. We show that for small $R$ the model can be approximated by a hierarchy of models, namely first a system of $N$ coupled delay ODEs, and in a second step by $N$ coupled ODEs. We give some illustrations of the dynamics of the approximate model.     

Sensory uncertainty and stick balancing at the fingertip

The effects of sensory input uncertainty, $\varepsilon $ , on the stability of time-delayed human motor control are investigated by calculating the minimum stick length, $\ell _\mathrm{crit}$ , that can be stabilized in the inverted position for a given time delay, $\tau $ . Five control strategies often discussed in the context of human motor control are examined: three time-invariant controllers [proportional–derivative, proportional–derivative–acceleration (PDA), model predictive (MP) controllers] and two time-varying controllers [act-and-wait (AAW) and intermittent predictive controllers]. The uncertainties of the sensory input are modeled as a multiplicative term in the system output. Estimates based on the variability of neural spike trains and neural population responses suggest that $\varepsilon \approx 7$ –13 %. It is found that for this range of uncertainty, a tapped delay-line type of MP controller is the most robust controller. In particular, this controller can stabilize inverted sticks of the length balanced by expert stick balancers (0.25–0.5 m when $\tau \approx 0.08$  s). However, a PDA controller becomes more effective when $\varepsilon > 15\,\%$ . A comparison between $\ell _\mathrm{crit}$ for human stick balancing at the fingertip and balancing on the rubberized surface of a table tennis racket suggest that friction likely plays a role in balance control. Measurements of $\ell _\mathrm{crit},\,\tau $ , and a variability of the fluctuations in the vertical displacement angle, an estimate of $\varepsilon $ , may make it possible to study the changes in control strategy as motor skill develops.     

Quantifying interspecific spatial overlap in aquatic macrophyte communities

Levins’s asymmetrical α index quantifies between species overlap over resources more realistically than similar-purpose single-value indices. The associated community-wide \(\bar \alpha\) index expresses the degree of “species packing”. Both indices were formulated upon competing animal (i.e., mobile) organisms and are independent of population densities. However, overlap over resources for nonmobile organisms such as plants may have an impact even below carrying capacity. The proposed \(\hat \alpha\) index, based on Levins’s α index, quantifies spatial overlap for plants integrating information on species spatial distribution and crowding conditions. The \(\hat \alpha\) index is specifically designed for plant distribution data collected in discrete plots with density expressed as percent coverage (%cover) of substratum. We also propose a community-wide \({\hat \alpha_{\text{c}}}\) index, conceptually analogous to \(\bar \alpha\) , but furnished with a measure of dispersion (se \({\hat \alpha_{\text{c}}}\) ). Species importance within the community is inferred from comparisons of pairwise \(\hat \alpha\) ’s with \({\hat \alpha_{\text{c}}}\) . The \(\hat \alpha\) and \({\hat \alpha_{\text{c}}}\) indices correlate closely and exponentially with plant density, and correct apparent over- and underestimations of interaction intensity at low and very high crowding by Levins’s α and \(\bar \alpha\) , respectively. Index application to aquatic plant communities gave results consistent with within-community and general ecological patterns, suggesting a high potential of the proposed \(\hat \alpha\) and \({\hat \alpha_{\text{c}}}\) indices in basic and applied macrophyte ecological studies and management.     

Modulation of the cAMP Response by G\alpha _i and G\beta \gamma : A Computational Study of G Protein Signaling in Immune Cells

Cyclic AMP is important for the resolution of inflammation, as it promotes anti-inflammatory signaling in several immune cell lines. In this paper, we present an immune cell specific model of the cAMP signaling cascade, paying close attention to the specific isoforms of adenylyl cyclase (AC) and phosphodiesterase that control cAMP production and degradation, respectively, in these cells. The model describes the role that G protein subunits, including G \(\alpha _s\) , G \(\alpha _i\) , and G \(\beta \gamma \) , have in regulating cAMP production. Previously, G \(\alpha _i\) activation has been shown to increase the level of cAMP in certain immune cell types. This increase in cAMP is thought to be mediated by \(\beta \gamma \) subunits which are released upon G \(\alpha \) activation and can directly stimulate specific isoforms of AC. We conduct numerical experiments in order to explore the mechanisms through which G \(\alpha _i\) activation can increase cAMP production. An important conclusion of our analysis is that the relative abundance of different G protein subunits is an essential determinant of the cAMP profile in immune cells. In particular, our model predicts that limited availability of \(\beta \gamma \) subunits may both \((i)\) enable immune cells to link inflammatory G \(\alpha _i\) signaling to anti-inflammatory cAMP production thereby creating a balanced immune response to stimulation with low concentrations of PGE2, and \((ii)\) prohibit robust anti-inflammatory cAMP signaling in response to stimulation with high concentrations of PGE2.     

Individual-Based Model for Quorum Sensing with Background Flow

Quorum sensing is a wide-spread mode of cell–cell communication among bacteria in which cells release a signalling substance at a low rate. The concentration of this substance allows the bacteria to gain information about population size or spatial confinement. We consider a model for \(N\) cells which communicate with each other via a signalling substance in a diffusive medium with a background flow. The model consists of an initial boundary value problem for a parabolic PDE describing the exterior concentration \(u\) of the signalling substance, coupled with \(N\) ODEs for the masses \(a_i\) of the substance within each cell. The cells are balls of radius \(R\) in \(\mathbb {R} ^3\) , and under some scaling assumptions we formally derive an effective system of \(N\) ODEs describing the behaviour of the cells. The reduced system is then used to study the effect of flow on communication in general, and in particular for a number of geometric configurations.     

Bone refilling in cortical basic multicellular units: insights into tetracycline double labelling from a computational model

Bone remodelling is carried out by ‘bone multicellular units’ ( $\text{ BMU }$ s) in which active osteoclasts and active osteoblasts are spatially and temporally coupled. The refilling of new bone by osteoblasts towards the back of the $\text{ BMU }$ occurs at a rate that depends both on the number of osteoblasts and on their secretory activity. In cortical bone, a linear phenomenological relationship between matrix apposition rate and $\text{ BMU }$ cavity radius is found experimentally. How this relationship emerges from the combination of complex, nonlinear regulations of osteoblast number and secretory activity is unknown. Here, we extend our previous mathematical model of cell development within a single cortical $\text{ BMU }$ to investigate how osteoblast number and osteoblast secretory activity vary along the $\text{ BMU }$ ’s closing cone. The mathematical model is based on biochemical coupling between osteoclasts and osteoblasts of various maturity and includes the differentiation of osteoblasts into osteocytes and bone lining cells, as well as the influence of $\text{ BMU }$ cavity shrinkage on osteoblast development and activity. Matrix apposition rates predicted by the model are compared with data from tetracycline double labelling experiments. We find that the linear phenomenological relationship observed in these experiments between matrix apposition rate and $\text{ BMU }$ cavity radius holds for most of the refilling phase simulated by our model, but not near the start and end of refilling. This suggests that at a particular bone site undergoing remodelling, bone formation starts and ends rapidly, supporting the hypothesis that osteoblasts behave synchronously. Our model also suggests that part of the observed cross-sectional variability in tetracycline data may be due to different bone sites being refilled by $\text{ BMU }$ s at different stages of their lifetime. The different stages of a $\text{ BMU }$ ’s lifetime (such as initiation stage, progression stage, and termination stage) depend on whether the cell populations within the $\text{ BMU }$ are still developing or have reached a quasi-steady state whilst travelling through bone. We find that due to their longer lifespan, active osteoblasts reach a quasi-steady distribution more slowly than active osteoclasts. We suggest that this fact may locally enlarge the Haversian canal diameter (due to a local lack of osteoblasts compared to osteoclasts) near the $\text{ BMU }$ ’s point of origin.     

Do calcium buffers always slow down the propagation of calcium waves?

Calcium buffers are large proteins that act as binding sites for free cytosolic calcium. Since a large fraction of cytosolic calcium is bound to calcium buffers, calcium waves are widely observed under the condition that free cytosolic calcium is heavily buffered. In addition, all physiological buffered excitable systems contain multiple buffers with different affinities. It is thus important to understand the properties of waves in excitable systems with the inclusion of buffers. There is an ongoing controversy about whether or not the addition of calcium buffers into the system always slows down the propagation of calcium waves. To solve this controversy, we incorporate the buffering effect into the generic excitable system, the FitzHugh–Nagumo model, to get the buffered FitzHugh–Nagumo model, and then to study the effect of the added buffer with large diffusivity on traveling waves of such a model in one spatial dimension. We can find a critical dissociation constant ( $K=K(a)$ ) characterized by system excitability parameter $a$ such that calcium buffers can be classified into two types: weak buffers ( $K\in (K(a),\infty )$ ) and strong buffers ( $K\in (0,K(a))$ ). We analytically show that the addition of weak buffers or strong buffers but with its total concentration $b_0^1$ below some critical total concentration $b_{0,c}^1$ into the system can generate a traveling wave of the resulting system which propagates faster than that of the origin system, provided that the diffusivity $D_1$ of the added buffers is sufficiently large. Further, the magnitude of the wave speed of traveling waves of the resulting system is proportional to $\sqrt{D_1}$ as $D_1\rightarrow \infty $ . In contrast, the addition of strong buffers with the total concentration $b_0^1>b_{0,c}^1$ into the system may not be able to support the formation of a biologically acceptable wave provided that the diffusivity $D_1$ of the added buffers is sufficiently large.     

Effects of collagen deposition on passive and active mechanical properties of large pulmonary arteries in hypoxic pulmonary hypertension

Proximal pulmonary artery (PA) stiffening is a strong predictor of mortality in pulmonary hypertension. Collagen accumulation is mainly responsible for PA stiffening in hypoxia-induced pulmonary hypertension (HPH) in mouse models. We hypothesized that collagen cross-linking and the type I isoform are the main determinants of large PA mechanical changes during HPH, which we tested by exposing mice that resist type I collagen degradation (Col1a1 $^\mathrm{R/R})$ and littermate controls (Col1a1 $^{+/+})$ to hypoxia for 10 days with or without $\beta $ -aminopropionitrile (BAPN) treatment to prevent cross-link formation. Static and dynamic mechanical tests were performed on isolated PAs with smooth muscle cells (SMC) in passive and active states. Percentages of type I and III collagen were quantified by histology; total collagen content and cross-linking were measured biochemically. In the SMC passive state, for both genotypes, hypoxia tended to increase PA stiffness and damping capacity, and BAPN treatment limited these increases. These changes were correlated with collagen cross-linking ( $p<0.05$ ). In the SMC active state, hypoxia increased PA dynamic stiffness and BAPN had no effect in Col1a1 $^{+/+}$ mice ( $p<0.05$ ). PA stiffness did not change in Col1a1 $^\mathrm{R/R}$ mice. Similarly, damping capacity did not change for either genotype. Type I collagen accumulated more in Col1a1 $^{+/+}$ mice, whereas type III collagen increased more in Col1a1 $^\mathrm{R/R}$ mice during HPH. In summary, PA passive mechanical properties (both static and dynamic) are related to collagen cross-linking. Type I collagen turnover is critical to large PA remodeling during HPH when collagen metabolism is not mutated and type III collagen may serve as a reserve.     

A locally convergent rotationally invariant particle swarm optimization algorithm

Several well-studied issues in the particle swarm optimization algorithm are outlined and some earlier methods that address these issues are investigated from the theoretical and experimental points of view. These issues are the: stagnation of particles in some points in the search space, inability to change the value of one or more decision variables, poor performance when the swarm size is small, lack of guarantee to converge even to a local optimum (local optimizer), poor performance when the number of dimensions grows, and sensitivity of the algorithm to the rotation of the search space. The significance of each of these issues is discussed and it is argued that none of the particle swarm optimizers we are aware of can address all of these issues at the same time. To address all of these issues at the same time, a new general form of velocity update rule for the particle swarm optimization algorithm that contains a user-definable function \(f\) is proposed. It is proven that the proposed velocity update rule guarantees to address all of these issues if the function \(f\) satisfies the following two conditions: (i) the function \(f\) is designed in such a way that for any input vector \(\vec {y}\) in the search space, there exists a region \(A\) which contains \(\vec {y}\) and \( f\!\left( {\vec {y}} \right) \) can be located anywhere in \(A\) , and (ii) \(f\) is invariant under any affine transformation. An example of function \(f\) is provided that satisfies these conditions and its performance is examined through some experiments. The experiments confirm that the proposed algorithm (with an appropriate function \(f)\) can effectively address all of these issues at the same time. Also, comparisons with earlier methods show that the overall ability of the proposed method for solving benchmark functions is significantly better.     

Topological classification and enumeration of RNA structures by genus

To an RNA pseudoknot structure is naturally associated a topological surface, which has its associated genus, and structures can thus be classified by the genus. Based on earlier work of Harer–Zagier, we compute the generating function $\mathbf{D}_{g,\sigma }(z)=\sum _{n}\mathbf{d}_{g,\sigma }(n)z^n$ for the number $\mathbf{d}_{g,\sigma }(n)$ of those structures of fixed genus $g$ and minimum stack size $\sigma $ with $n$ nucleotides so that no two consecutive nucleotides are basepaired and show that $\mathbf{D}_{g,\sigma }(z)$ is algebraic. In particular, we prove that $\mathbf{d}_{g,2}(n)\sim k_g\,n^{3(g-\frac{1}{2})} \gamma _2^n$ , where $\gamma _2\approx 1.9685$ . Thus, for stack size at least two, the genus only enters through the sub-exponential factor, and the slow growth rate compared to the number of RNA molecules implies the existence of neutral networks of distinct molecules with the same structure of any genus. Certain RNA structures called shapes are shown to be in natural one-to-one correspondence with the cells in the Penner–Strebel decomposition of Riemann’s moduli space of a surface of genus $g$ with one boundary component, thus providing a link between RNA enumerative problems and the geometry of Riemann’s moduli space.     

Studies on mass attenuation coefficient, effective atomic number and electron density of some vitamins

Mass attenuation coefficient, $ \mu_{m} $ , atomic cross-section, $ \sigma_{i} $ , electronic cross-section, $ \sigma_{e} $ , effective atomic number, $ Z_{\text{eff}} $ and effective electron density, $ N_{\text{el}} $ , were determined experimentally and theoretically for some vitamins (retinol, beta-carotene, thiamine, riboflavin, niacinamide, pantothenic acid, pyridoxine, biotin, folic acid, cyanocobalamin, ascorbic acid, cholecalciferol, alpha-tocopherol, ketamine, hesperidin) at 30.82, 59.54, 80.99, 356.61, 661.66 and 1,408.01?keV photon energies using a NaI(Tl) scintillation detector. The theoretical mass attenuation coefficients were estimated using mixture rules. The calculated values were compared with the experimental values for all vitamins.     

Elastic behavior of a red blood cell with the membrane’s nonuniform natural state: equilibrium shape,motion transition under shear flow,and elongation during tank-treading motion

Direct numerical simulations of the mechanics of a single red blood cell (RBC) were performed by considering the nonuniform natural state of the elastic membrane. A RBC was modeled as an incompressible viscous fluid encapsulated by an elastic membrane. The in-plane shear and area dilatation deformations of the membrane were modeled by Skalak constitutive equation, while out-of-plane bending deformation was formulated by the spring model. The natural state of the membrane with respect to in-plane shear deformation was modeled as a sphere ( \(\alpha =0\) ), biconcave disk shape ( \(\alpha =1\) ) and their intermediate shapes ( \(0<\alpha <1\) ) with the nonuniformity parameter \(\alpha \) , while the natural state with respect to out-of-plane bending deformation was modeled as a flat plane. According to the numerical simulations, at an experimentally measured in-plane shear modulus of \(2.5\times 10^{-6}\,\hbox {N}/\hbox {m}\) and an out-of-plane bending rigidity of \(2.0\times 10^{-19}\,\hbox {N}\cdot \hbox {m}\) of the cell membrane, the following results were obtained. (i) The RBC shape at equilibrium was biconcave discoid for \(\alpha >0.22\) and cupped otherwise; (ii) the experimentally measured fluid shear stress at the transition between tumbling and tank-treading motions under shear flow was reproduced for \(0.05<\alpha <0.34\) ; (iii) the elongation deformation of the RBC during tank-treading motion from the simulation was consistent with that from in vitro experiments, irrespective of the \(\alpha \) value. Based on our RBC modeling, the three phenomena (i), (ii), and (iii) were mechanically consistent for \(0.22<\alpha <0.34\) . The condition \(0.05<\alpha <0.22\) precludes a biconcave discoid shape at equilibrium (i); however, it gives appropriate fluid shear stress at the motion transition under shear flow (ii), suggesting that a combined effect of \(\alpha \) and the natural state with respect to out-of-plane bending deformation is necessary for understanding details of the RBC mechanics at equilibrium. Our numerical results demonstrate that moderate nonuniformity in a membrane’s natural state with respect to in-plane shear deformation plays a key role in RBC mechanics.     

The Basic Reproduction Number for Cellular SIR Networks

The basic reproduction number \(R_0\) is the average number of new infections produced by a typical infective individual in the early stage of an infectious disease, following the introduction of few infective individuals in a completely susceptible population. If \(R_0<1\) , then the disease dies, whereas for \(R_0>1\) the infection can invade the host population and persist. This threshold quantity is well studied for SIR compartmental or mean field models based on ordinary differential equations, and a general method for its computation has been proposed by van den Driessche and Watmough. We concentrate here on SIR epidemiological models that take into account the contact network N underlying the transmission of the disease. In this context, it is generally admitted that \(R_{0}\) can be approximated by the average number \(R_{2,3}\) of infective individuals of generation three produced by an infective of generation two. We give here a simple analytic formula of \(R_{2,3}\) for SIR cellular networks. Simulations on two-dimensional cellular networks with von Neumann and Moore neighborhoods show that \(R_{2,3}\) can be used to capture a threshold phenomenon related the dynamics of SIR cellular network and confirm the good quality of the simple approach proposed recently by Aparicio and Pascual for the particular case of Moore neighborhood.     

Modeling analysis of the lipid bilayer–cytoskeleton coupling in erythrocyte membrane

Studies of thermal fluctuations in discocytes, echinocytes, and spherocytes suggest that the coupling between lipid bilayer and cytoskeleton can affect viscoelastic behavior of single erythrocyte membranes. To test this hypothesis, we developed a 3D constitutive model describing viscoelastic behavior of erythrocyte membranes, at long relaxation times \(t \in [0.20\,\mathrm {s}, 1.05\,\mathrm {s}]\) and short relaxation times \(t \in [0.03\,\mathrm {s}, 0.20\,\mathrm {s}]\) . The model was constructed using combination of spring and spring pot rheological elements arranged in parallel. The rearrangement of cytoskeleton induced by changing the bending state of lipid bilayer was described by a modified Eyring model. The model predictions point to an anomalous nature of energy dissipation and an ordered harmonic nature of the coupling mechanism described by a series of fractional derivatives of the order n \(\alpha \) (where \( n \in [- 1, 2]\) ). As a result, the stress generated within the lipid bilayer is related to the rate of change of the irreversible stress within the cytoskeleton.     

Dynamic permeability of the lacunar–canalicular system in human cortical bone

A new method for the experimental determination of the permeability of a small sample of a fluid-saturated hierarchically structured porous material is described and applied to the determination of the lacunar–canalicular permeability \((K_\mathrm{LC})\) in bone. The interest in the permeability of the lacunar–canalicular pore system (LCS) is due to the fact that the LCS is considered to be the site of bone mechanotransduction due to the loading-driven fluid flow over cellular structures. The permeability of this space has been estimated to be anywhere from \(10^{-17}\;\) to \(10^{-25}\; \hbox {m}^{2}\) . However, the vascular pore system and LCS are intertwined, rendering the permeability of the much smaller-dimensioned LCS challenging to measure. In this study, we report a combined experimental and analytical approach that allowed the accurate determination of the \(K_\mathrm{LC}\) to be on the order of \(10^{-22}\; \hbox {m}^{2}\) for human osteonal bone. It was found that the \(K_\mathrm{LC}\) has a linear dependence on loading frequency, decreasing at a rate of \(2 \times 10^{-24}\; \hbox {m}^{2}\) /Hz from 1 to 100 Hz, and using the proposed model, the porosity alone was able to explain 86 % of the \(K_\mathrm{LC}\) variability.