Vascular deposition patterns for nanoparticles in an inflamed patient-specific arterial tree

Inflammation, a precursor to many diseases including cancer and atherosclerosis, induces differential surface expression of specific vascular molecules. Blood-borne nanoparticles (NPs), loaded with therapeutic and imaging agents, can recognize and use these molecules as vascular docking sites. Here, a computational model is developed within the isogeometric analysis framework to understand and predict the vascular deposition of NPs within an inflamed arterial tree. The NPs have a diameter ranging from 0.1 to $2.0\,\upmu $ m and are decorated with antibodies directed toward three endothelial adhesion molecules, namely intravascular cell adhesion molecule-1 (ICAM-1), vascular cell adhesion molecule-1 (VCAM-1), and E-selectin, whose surface density depends on the local wall shear stress. Results indicate VCAM-1 targeted NPs adhere more, with ICAM-1 directed NPs adhering least efficiently, resulting in approximately an order-of-magnitude lower average particle surface density. ICAM-1 and E-selectin directed $0.5\,\upmu $ m NPs are distributed more uniformly (heterogeneity index $\approx $  0.9 and 1.0, respectively) over the bifurcating vascular branches compared to their VCAM-1 counterparts (heterogeneity index $\approx $  1.4). When the NPs are coated with antibodies for VCAM-1 and E-selectin in equal proportions, a more uniform vascular distribution is achieved compared with VCAM-1-only targeted particles, thus demonstrating the advantage of NP multivalency in vascular targeting. Furthermore, the larger NPs ( $2\,\upmu $ m) adhere more ( $\approx $  200 %) in the lower branches compared to the upper branch. This computational framework provides insights into how size, ligand type, density, and multivalency can be manipulated to enhance NP vascular adhesion in an individual patient.     

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.     

Basic Features of a Cell Electroporation Model: Illustrative Behavior for Two Very Different Pulses

Science increasingly involves complex modeling. Here we describe a model for cell electroporation in which membrane properties are dynamically modified by poration. Spatial scales range from cell membrane thickness (5 nm) to a typical mammalian cell radius (10  \(\upmu\) m), and can be used with idealized and experimental pulse waveforms. The model consists of traditional passive components and additional active components representing nonequilibrium processes. Model responses include measurable quantities: transmembrane voltage, membrane electrical conductance, and solute transport rates and amounts for the representative “long” and “short” pulses. The long pulse—1.5 kV/cm, 100  \(\upmu\) s—evolves two pore subpopulations with a valley at \({\sim}\) 5 nm, which separates the subpopulations that have peaks at \({\sim}\) 1.5 and \({\sim}\) 12 nm radius. Such pulses are widely used in biological research, biotechnology, and medicine, including cancer therapy by drug delivery and nonthermal physical tumor ablation by causing necrosis. The short pulse—40 kV/cm, 10 ns—creates 80-fold more pores, all small ( \(<\) 3 nm; \(\sim\) 1 nm peak). These nanosecond pulses ablate tumors by apoptosis. We demonstrate the model’s responses by illustrative electrical and poration behavior, and transport of calcein and propidium. We then identify extensions for expanding modeling capability. Structure-function results from MD can allow extrapolations that bring response specificity to cell membranes based on their lipid composition. After a pulse, changes in pore energy landscape can be included over seconds to minutes, by mechanisms such as cell swelling and pulse-induced chemical reactions that slowly alter pore behavior.     

Principal trabecular structural orientation predicted by quantitative ultrasound is strongly correlated with \upmu FEA determined anisotropic apparent stiffness

The microarchitecture and alignment of trabecular bone adapts to the particular mechanical milieu applied to it. Due to this anisotropic mechanical property, measurement orientation has to be taken into consideration when assessing trabecular bone quality and fracture risk prediction. Quantitative ultrasound (QUS) has demonstrated the ability in predicting the principal structural orientation (PSO) of trabecular bone. Although the QUS prediction for PSO is very close to that of \(\upmu \) CT, certain angle differences still exist. It remains unknown whether this angle difference can induce significant differences in mechanical properties or not. The objective of this study was to evaluate the mechanical properties in different PSOs predicted using different methods, QUS and \(\upmu \) CT, thus to investigate the ability of QUS as a means to predict the PSO of trabecular bone noninvasively. By validating the ability of QUS to predict the PSO of trabecular bone, it is beneficial for future QUS applications because QUS measurements in the PSO can provide information more correlated with the mechanical properties than with other orientations. In this study, seven trabecular bone balls from distal bovine femurs were used to generate finite element models based on the 3-dimensional \(\upmu \) CT images. Uniaxial compressive loading was performed on the bone ball models in the finite element analysis (FEA) in six different orientations (three anatomical orientations, two PSOs predicted by QUS and the longest vector of mean intercept length (MIL) tensor calculated by \(\upmu \) CT). The stiffness was calculated based on the reaction force of the bone balls under loading, and the von Mises stress results showed that both the mechanical properties in the PSOs predicted by QUS are significantly higher than the anatomical orientations and comparatively close to the longest vector of MIL tensor. The stiffness in the PSOs predicted by QUS is also highly correlated with the stiffness in the MIL tensor orientation (ATTmax vs. MIL, \(R^{2}\)  = 0.98, \(p<001\) ; UVmax vs. MIL, \(R^{2}\)  = 0.92, \(p<001\) ). These results were validated by in vitro mechanical testing on the bone ball samples. This study demonstrates that the PSO of trabecular bone predicted by QUS has an equally strong apparent stiffness with the orientation predicted by \(\upmu \) CT.     

Increased Damping of Plasmon Resonances in Gold Nanoparticles Due to Broadening of the Band Structure

For the first time, systematic investigations of the damping parameter A of gold nanoparticles as a function of photon energy are presented. A is an essential parameter that quantifies the size-dependent optical properties of metal nanoparticles in the dielectric function. To determine the damping parameter, the dephasing time T ₂ of gold nanoparticles has been systematically determined under ultrahigh vacuum conditions as a function of photon energy. Dephasing times ranging from  $T_2 = 5$ fs to $T_2 = 17$ fs were measured, and subsequently, the damping parameter has been extracted. We found a strong resonance-like damping of the plasmon resonance in the vicinity of the onset of the interband transition. While the damping parameter scatters statistically around a value of  $A = 0.19$ nm/fs for photon energies below  $h\nu = 1.70$ eV, it increases rapidly to 0.32 nm/fs for  $h\nu = 1.85$ eV. For higher photon energies, A decreases steadily to $A = 0.24$ nm/fs at $h\nu = 2.15$ eV. A comparison to former measurements as well as to theoretical predictions reveals surface scattering and a discretizing and broadening of the band structure that influences the interband transition as the most dominant size-dependent damping mechanisms. The latter, i.e., a damping via increased interband transitions, assumes a coherent damping process of the oscillating electrons and, as a consequence, the plasmon is treated as a two-level system. Thus, the results deliver new physical insight to the fundamental understanding of surface plasmons.     

Release of non-exchangeable {}^{{{\text{15}}}}{\text{NH}}^{{\text{ + }}}_{{\text{4}}} from subgrade, decomposed granite substrates and uptake by non-mycorrhizal and mycorrhizal California native annual grass, Vulpia microstachys

Release rates of recently fixed $ {\text{NH}}^{{\text{ + }}}_{{\text{4}}} $ from non-exchangeable interlayer sites in 2:1 silicate minerals were determined for decomposed granite (DG) saprolites from three locations in California, USA. Recently-fixed $ {\text{NH}}^{{\text{ + }}}_{{\text{4}}} $ release from the DG substrate was quantified by extracting diffused $ {\text{NH}}^{{\text{ + }}}_{{\text{4}}} $ with H-resin, as well as a native, annual grass Vulpia microstachys. The $ {\text{NH}}^{{\text{ + }}}_{{\text{4}}} $ release data varied with via the method of extraction, which included H-resin pre-treatments (Na⁺ or H⁺) and V. microstachys uptake (mycorrhizal inoculated or uninoculated). After 6 weeks (1008 h), more $ {\text{NH}}^{{\text{ + }}}_{{\text{4}}} $ was recovered from fixed interlayer positions by the H-resins as compared to uptake by V. microstachys. The H⁺ treated H-resins recovered more released $ {\text{NH}}^{{\text{ + }}}_{{\text{4}}} $ (≈94 mg ${\text{NH}}^{{\text{ + }}}_{{\text{4}}} - {\text{N}}\;{\text{kg}}^{1} $ or (12%) of total fixed $ {\text{NH}}^{{\text{ + }}}_{{\text{4}}} $ ) in two of the three DG samples as compared to the Na⁺ treated resins, (which recovered ≈70–78 mg ${\text{NH}}^{{\text{ + }}}_{{\text{4}}} - {\text{N}}\;{\text{kg}}^{{{\text{ - 1}}}} $ (or 9–10%) of the total fixed $ {\text{NH}}^{{\text{ + }}}_{{\text{4}}} $ ). The V. microstachys assimilated 8–9% of the total fixed $ {\text{NH}}^{{\text{ + }}}_{{\text{4}}} $ with mycorrhizal inoculum as compared to only 2% without a mycorrhizal inoculum, over the same time period. The fixed $ {\text{NH}}^{{\text{ + }}}_{{\text{4}}} $ release kinetics from the H-resin experiments were most accurately described by first order and power function models, and can be characterized as biphasic using a heterogeneous diffusion model. Uptake of both the ¹⁵N and ambient, unlabelled N from the soils was closely related to plant biomass. There was no significant difference in percent of N per unit of biomass between the control and mycorrhizal treatments. The findings presented here indicate that observed, long-term $ {\text{NH}}^{{\text{ + }}}_{{\text{4}}} $ release rates from DG in studies utilizing resins, may overestimate the levels of fixed $ {\text{NH}}^{{\text{ + }}}_{{\text{4}}} $ made available to plants and microorganisms. Additionally, the study suggested that mycorrhizae facilitate the acquisition and plant uptake of fixed $ {\text{NH}}^{{\text{ + }}}_{{\text{4}}} $ , resulting in markedly increased plant biomass production.     

Fluid flow in the osteocyte mechanical environment: a fluid–structure interaction approach

Osteocytes are believed to be the primary sensor of mechanical stimuli in bone, which orchestrate osteoblasts and osteoclasts to adapt bone structure and composition to meet physiological loading demands. Experimental studies to quantify the mechanical environment surrounding bone cells are challenging, and as such, computational and theoretical approaches have modelled either the solid or fluid environment of osteocytes to predict how these cells are stimulated in vivo. Osteocytes are an elastic cellular structure that deforms in response to the external fluid flow imposed by mechanical loading. This represents a most challenging multi-physics problem in which fluid and solid domains interact, and as such, no previous study has accounted for this complex behaviour. The objective of this study is to employ fluid–structure interaction (FSI) modelling to investigate the complex mechanical environment of osteocytes in vivo. Fluorescent staining of osteocytes was performed in order to visualise their native environment and develop geometrically accurate models of the osteocyte in vivo. By simulating loading levels representative of vigorous physiological activity ( $3,000\,\upmu \upvarepsilon $ compression and 300 Pa pressure gradient), we predict average interstitial fluid velocities $(\sim 60.5\,\upmu \text{ m/s })$ and average maximum shear stresses $(\sim 11\, \text{ Pa })$ surrounding osteocytes in vivo. Interestingly, these values occur in the canaliculi around the osteocyte cell processes and are within the range of stimuli known to stimulate osteogenic responses by osteoblastic cells in vitro. Significantly our results suggest that the greatest mechanical stimulation of the osteocyte occurs in the cell processes, which, cell culture studies have indicated, is the most mechanosensitive area of the cell. These are the first computational FSI models to simulate the complex multi-physics mechanical environment of osteocyte in vivo and provide a deeper understanding of bone mechanobiology.     

Integrating qPLM and biomechanical test data with an anisotropic fiber distribution model and predictions of TGF-\upbeta 1 and IGF-1 regulation of articular cartilage fiber modulus

A continuum mixture model with distinct collagen (COL) and glycosaminoglycan elastic constituents was developed for the solid matrix of immature bovine articular cartilage. A continuous COL fiber volume fraction distribution function and a true COL fiber elastic modulus ( $E^\mathrm{f})$ were used. Quantitative polarized light microscopy (qPLM) methods were developed to account for the relatively high cell density of immature articular cartilage and used with a novel algorithm that constructs a 3D distribution function from 2D qPLM data. For specimens untreated and cultured in vitro, most model parameters were specified from qPLM analysis and biochemical assay results; consequently, $E^\mathrm{f}$ was predicted using an optimization to measured mechanical properties in uniaxial tension and unconfined compression. Analysis of qPLM data revealed a highly anisotropic fiber distribution, with principal fiber orientation parallel to the surface layer. For untreated samples, predicted $E^\mathrm{f}$ values were 175 and 422 MPa for superficial (S) and middle (M) zone layers, respectively. TGF- $\upbeta $ 1 treatment was predicted to increase and decrease $E^\mathrm{f}$ values for the S and M layers to 281 and 309 MPa, respectively. IGF-1 treatment was predicted to decrease $E^\mathrm{f}$ values for the S and M layers to 22 and 26 MPa, respectively. A novel finding was that distinct native depth-dependent fiber modulus properties were modulated to nearly homogeneous values by TGF- $\upbeta $ 1 and IGF-1 treatments, with modulated values strongly dependent on treatment.     

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.     

Obtain Quadruple Intense Plasmonic Resonances from Multilayered Gold Nanoshells by Silver Coating: Application in Multiplex Sensing

Four intense and separate localized surface plasmon resonance (LSPR) absorption peaks have been obtained in the gold-dielectric–gold–silver multilayer nanoshells. The silver coating on the gold shell results in a new LSPR peak at about 400 nm corresponding to the $ {{\left| {\omega_{+}^{-}} \right\rangle}_{Ag }} $ mode. The intense local electric field concentrated in the silver shell at the wavelength of 400 nm indicates that this new plasmonic band is coming from the symmetric coupling between the antibonding silver shell plasmon mode and the inner sphere plasmon. Increasing the silver shell thickness also leads to the intensity increasing of the $ {{\left| {\omega_{+}^{-}} \right\rangle}_{Au }} $ mode and blue shift of $ \left| {\omega_{-}^{+}} \right\rangle $ and $ \left| {\omega_{-}^{-}} \right\rangle $ modes. Therefore, quadruple intense plasmonic resonances in the visible region could be achieved in gold-dielectric–gold–silver multilayer nanoshells by tuning the geometrical parameters. And the quadruple intense plasmonic resonances in the visible region provide well potential for multiplex biosensing based on LSPR.     

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.     

Deriving fluorometer-specific values of relative PSI fluorescence intensity from quenching of F 0 fluorescence in leaves of Arabidopsis thaliana and Zea mays

The effect of stepwise increments of red light intensities on pulse-amplitude modulated (PAM) chlorophyll (Chl) fluorescence from leaves of A. thaliana and Z. mays was investigated. Minimum and maximum fluorescence were measured before illumination (F ₀ and F _M, respectively) and at the end of each light step ( $ F^{\prime}_{0} $ and $ F^{\prime}_{\text{M}} $ , respectively). Calculated $ F^{\prime}_{0} $ values derived from F ₀, F _M and $ F^{\prime}_{\text{M}} $ fluorescence according to Oxborough and Baker (1997) were lower than the corresponding measured $ F^{\prime}_{0} $ values. Based on the concept that calculated $ F^{\prime}_{0} $ values are under-estimated because the underlying theory ignores PSI fluorescence, a method was devised to gain relative PSI fluorescence intensities from differences between calculated and measured $ F^{\prime}_{0} $ . This method yields fluorometer-specific PSI data as its input data (F ₀, F _M, $ F^{\prime}_{0} $ and $ F^{\prime}_{\text{M}} $ ) depend solely on the spectral properties of the fluorometer used. Under the present conditions, the PSI contribution to F ₀ fluorescence was 0.24 in A. thaliana and it was independent on the light acclimation status; the corresponding value was 0.50 in Z. mays. Correction for PSI fluorescence affected Z. mays most: the linear relationship between PSI and PSII photochemical yields was clearly shifted toward the one-to-one proportionality line and maximum electron transport was increased by 50 %. Further, correction for PSI fluorescence increased the PSII reaction center-specific parameter, 1/F ₀ ? 1/F _M, up to 50 % in A. thaliana and up to 400 % in Z. mays.     

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.     

Numerical Modeling of Bi-polar (AC) Pulse Electroporation of Single Cell in Microchannel to Create Nanopores on its Membrane

AC electroporation of a single cell in a microchannel was numerically studied. A \(15\,\upmu\) m diameter cell was considered in a microchannel \(25\,\upmu\) m in height and the influences of AC electric pulse on its membrane were numerically investigated. The cell was assumed to be suspended between two electroporative electrodes embedded on the walls of a microchannel. An amplitude and a time span of applied electric pulse were chosen to be 80 kV/m and \(10\,\upmu\) s, respectively. For different frequency values (50, 100, 200, and 500 kHz), simulations were performed to show how the cell membrane was electroporated and the creation of nanopores. Obtained numerical results show that the most and the largest nanopores are created around poles of cell (nearest points of cell membrane to the electrodes). The numerical simulations also demonstrate that increased frequency will slightly decrease electroporated area of the cell membrane; additionally, growth of the created nanopores will be stabilized. It has also been proven that size and number of the created nanopores will be decreased by moving from the poles to the equator of the cell. There is almost no nanopore created in the vicinity of the equator. Frequency affects the rate of generation of nanopores. In case of AC electroporation, creation of nanopores has two phases that periodically repeat over time. In each period, the pore density sharply increases and then becomes constant. Enhancement of the frequency will result in decrease in time span of the periods. In each period, size of the created nanopores sharply increases and then slightly decreases. However, until the AC electric pulse is present, overall trends of creation and development of nanopores will be ascending. Variation of the size and number of created nanopores can be explained by considering time variation of transmembrane potential (difference of electric potential on two sides of cell membrane) which is clear in the results presented in this study.     

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.     

Cherry Picking: A Characterization of the Temporal Hybridization Number for a Set of Phylogenies

Recently, we have shown that calculating the minimum–temporal-hybridization number for a set ${\mathcal{P}}$ of rooted binary phylogenetic trees is NP-hard and have characterized this minimum number when ${\mathcal{P}}$ consists of exactly two trees. In this paper, we give the first characterization of the problem for ${\mathcal{P}}$ being arbitrarily large. The characterization is in terms of cherries and the existence of a particular type of sequence. Furthermore, in an online appendix to the paper, we show that this new characterization can be used to show that computing the minimum–temporal hybridization number for two trees is fixed-parameter tractable.     

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.     

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.     

A network with tunable clustering, degree correlation and degree distribution, and an epidemic thereon

A random network model which allows for tunable, quite general forms of clustering, degree correlation and degree distribution is defined. The model is an extension of the configuration model, in which stubs (half-edges) are paired to form a network. Clustering is obtained by forming small completely connected subgroups, and positive (negative) degree correlation is obtained by connecting a fraction of the stubs with stubs of similar (dissimilar) degree. An SIR (Susceptible $\rightarrow $ Infective $\rightarrow $ Recovered) epidemic model is defined on this network. Asymptotic properties of both the network and the epidemic, as the population size tends to infinity, are derived: the degree distribution, degree correlation and clustering coefficient, as well as a reproduction number $R_*$ , the probability of a major outbreak and the relative size of such an outbreak. The theory is illustrated by Monte Carlo simulations and numerical examples. The main findings are that (1) clustering tends to decrease the spread of disease, (2) the effect of degree correlation is appreciably greater when the disease is close to threshold than when it is well above threshold and (3) disease spread broadly increases with degree correlation $\rho $ when $R_*$ is just above its threshold value of one and decreases with $\rho $ when $R_*$ is well above one.     

Towards a conceptualization of ETL and physical storage of semantic data warehouses as a service

The data warehouse technology has become the incontestable tool for businesses and organizations to make strategic decisions to ensure their competitively. The construction of a data warehouse ( $\mathcal{D}\mathcal{W}$ ) passes by selecting relevant information sources, extracting relevant data and loading them into the $\mathcal{D}\mathcal{W}$ . These processes require a precise expertise from designers related to logical and physical implementations of information sources, which is not usually an easy task. The diversity and heterogeneity of information sources makes the construction process of the $\mathcal{D}\mathcal{W}$ complex and time consuming. Domain ontologies have been proposed to reduce heterogeneity between sources, platforms, services, etc. They resolve syntax and semantic conflicts. The phenomenon of adopting domain ontologies by organizations creates a new type of databases, called semantic databases ( $\mathcal{S}\mathcal{D}\mathcal{B}$ ). As a consequence, they become a candidate for building the semantic $\mathcal{D}\mathcal{W}$ ( $\mathcal{S}\mathcal{D}\mathcal{W}$ ). To handle the diversity of information sources and hide the implementations aspects of sources, proposing a generic framework for constructing $\mathcal{D}\mathcal{W}$ becomes a necessity. In this paper, we first proposed an ontology-based approach for designing $\mathcal{S}\mathcal {D}\mathcal{B}$ . Secondly, ETL phases are defined at ontological level to hide the implementation details. Thirdly, a storage service for ontologies and its associated data is given. Finally, our proposal is validated through a case study and a tool.