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.     

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.     

Effect of cyclic loading on the nanoscale deformation of hydroxyapatite and collagen fibrils in bovine bone

Cyclic compressive loading tests were carried out on bovine femoral bones at body temperature $(37\,^{\circ }\hbox {C})$ , with varying mean stresses ( $-55$ to $-80$  MPa) and loading frequencies (0.5–5 Hz). At various times, the cyclic loading was interrupted to carry out high-energy X-ray scattering measurements of the internal strains developing in the hydroxyapatite (HAP) platelets and the collagen fibrils. The residual strains upon unloading were always tensile in the HAP and compressive in the fibrils, and each increases in magnitude with loading cycles, which can be explained from damage at the HAP–collagen interface and accumulation of plastic deformation within the collagen phase. The samples tested at a higher mean stress and stress amplitude, and at lower loading frequencies exhibit greater plastic deformation and damage accumulation, which is attributed to greater contribution of creep. Synchrotron microcomputed tomography of some of the specimens showed that cracks are produced during cyclic loading and that they mostly occur concentric with Haversian canals.     

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.     

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.     

Theoretical study of the pH-dependent antioxidant properties of vitamin C

Molecules acting as antioxidants capable of scavenging reactive oxygen species (ROS) are of utmost importance in the living cell. Vitamin C is known to be one of these molecules. In this study we have analyzed the reactivity of vitamin C toward the $ \cdot OH $ and $ \cdot OOH $ ROS species, in all acidic, neutral and basic media. In order to do so, density functional theory (DFT) have been used. More concretely, the meta-GGA functional MPW1B95 have been used. Two reaction types have been studied in each case: addition to the ring atoms, and hydrogen/proton abstraction. Our results show that $ \cdot OH $ is the most reactive species, while $ \cdot OOH $ displays low reactivity. In all three media, vitamin C reactions with two hydroxyl radicals show a wide variety of possible products.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

Enhanced accuracy of kinetic information from CT-CPMG experiments by transverse rotating-frame spectroscopy

Micro-to-millisecond motions of proteins transmit pivotal signals for protein function. A powerful technique for the measurement of these motions is nuclear magnetic resonance spectroscopy. One of the most widely used methodologies for this purpose is the constant-time Carr–Purcell–Meiboom–Gill (CT-CPMG) relaxation dispersion experiment where kinetic and structural information can be obtained at atomic resolution. Extraction of accurate kinetics determined from CT-CPMG data requires refocusing frequencies that are much larger than the nuclei’s exchange rate between states. We investigated the effect when fast processes are probed by CT-CPMG experiments via simulation and show that if the intrinsic relaxation rate $ \left( {R_{2,0}^{CT - CPMG} } \right) $ ( R 2 , 0 CT ? CPMG ) is not known a priori the extraction of accurate kinetics is hindered. Errors on the order of 50 % in the exchange rate are attained when processes become fast, but are minimized to 5 % with a priori $ R_{2,0}^{CT - CPMG} $ R 2 , 0 CT ? CPMG information. To alleviate this shortcoming, we developed an experimental scheme probing $ R_{2,0}^{CT - CPMG} $ R 2 , 0 CT ? CPMG with large amplitude spin-lock fields, which specifically contains the intrinsic proton longitudinal Eigenrelaxation rate. Our approach was validated with ubiquitin and the Oscillatoria agardhii agglutinin (OAA). For OAA, an underestimation of 66 % in the kinetic rates was observed if $ R_{2,0}^{CT - CPMG}\, $ R 2 , 0 CT ? CPMG is not included during the analysis of CT-CPMG data and result in incorrect kinetics and imprecise amplitude information. This was overcome by combining CT-CPMG with $ R_{2,0}^{CT - CPMG} $ R 2 , 0 CT ? CPMG measured with a high power R_1ρ experiment. In addition, the measurement of $ R_{2,0}^{CT - CPMG} $ R 2 , 0 CT ? CPMG removes the ambiguities in choosing between different models that describe CT-CPMG data.     

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.     

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.     

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.     

The effects of temperature and exercise training on swimming performance in juvenile qingbo (Spinibarbus sinensis)

To investigate the effects of temperature and exercise training on swimming performance in juvenile qingbo (Spinibarbus sinensis), we measured the following: (1) the resting oxygen consumption rate $ \left( {{\dot{\text{M}}\text{O}}_{{ 2 {\text{rest}}}} } \right) $ , critical swimming speed (U _crit) and active oxygen consumption rate $ \left( {{\dot{\text{M}}\text{O}}_{{ 2 {\text{active}}}} } \right) $ of fish at acclimation temperatures of 10, 15, 20, 25 and 30 °C and (2) the $ \dot{M}{\text{O}}_{{ 2 {\text{rest}}}} $ , U _crit and $ \dot{M}{\text{O}}_{{ 2 {\text{active}}}} $ of both exercise-trained (exhaustive chasing training for 14 days) and control fish at both low and high acclimation temperatures (15 and 25 °C). The relationship between U _crit and temperature (T) approximately followed a bell-shaped curve as temperature increased: U _crit = 8.21/{1 + [(T ? 27.2)/17.0]²} (R ² = 0.915, P < 0.001, N = 40). The optimal temperature for maximal U _crit (8.21 BL s^?1) in juvenile qingbo was 27.2 °C. Both the $ \dot{M}{\text{O}}_{{ 2 {\text{active}}}} $ and the metabolic scope (MS, $ \dot{M}{\text{O}}_{{ 2 {\text{active}}}} - \dot{M}{\text{O}}_{{ 2 {\text{rest}}}} $ ) of qingbo increased with temperature from 10 to 25 °C (P < 0.05), but there were no significant differences between fish acclimated to 25 and 30 °C. The relationships between $ \dot{M}{\text{O}}_{{ 2 {\text{active}}}} $ or MS and temperature were described as $ {\dot{\text{M}}\text{O}}_{{ 2 {\text{active}}}} = 1,214.29/\left\{ {1 + \left[ {\left( {T - 28.8} \right)/10.6} \right]^{2} } \right\}\;\left( {R^{2} = 0.911,\;P < 0.001,\;N = 40} \right) $ and MS = 972.67/{1 + [(T ? 28.0)/9.34]²} (R ² = 0.878, P < 0.001, N = 40). The optimal temperatures for $ \dot{M}{\text{O}}_{{ 2 {\text{active}}}} $ and MS in juvenile qingbo were 28.8 and 28.0 °C, respectively. Exercise training resulted in significant increases in both U _crit and $ \dot{M}{\text{O}}_{{ 2 {\text{active}}}} $ at a low temperature (P < 0.05), but training exhibited no significant effect on either U _crit or $ \dot{M}{\text{O}}_{{ 2 {\text{active}}}} $ at a high temperature. These results suggest that exercise training had different effects on swimming performance at different temperatures. These differences may be related to changes in aerobic metabolic capability, arterial oxygen delivery, available dissolved oxygen, imbalances in ion fluxes and stimuli to remodel tissues with changes in temperature.     

On the page number of RNA secondary structures with pseudoknots

Let ${\mathcal {S}}$ denote the set of (possibly noncanonical) base pairs {i, j} of an RNA tertiary structure; i.e. ${\{i, j\} \in \mathcal {S}}$ if there is a hydrogen bond between the ith and jth nucleotide. The page number of ${\mathcal {S}}$ , denoted ${\pi(\mathcal {S})}$ , is the minimum number k such that ${\mathcal {S}}$ can be decomposed into a disjoint union of k secondary structures. Here, we show that computing the page number is NP-complete; we describe an exact computation of page number, using constraint programming, and determine the page number of a collection of RNA tertiary structures, for which the topological genus is known. We describe an approximation algorithm from which it follows that ${\omega(\mathcal {S}) \leq \pi(\mathcal {S}) \leq \omega(\mathcal {S}) \cdot \log n}$ , where the clique number of ${\mathcal {S}, \omega(\mathcal {S})}$ , denotes the maximum number of base pairs that pairwise cross each other.