Development of a novel MATLAB-based framework for implementing mechanical joint stability constraints within OpenSim musculoskeletal models

The Static Optimization (SO) solver in OpenSim estimates muscle activations and forces that only equilibrate applied moments. In this study, SO was enhanced through an open-access MATLAB interface, where calculated muscle activations can additionally satisfy crucial mechanical stability requirements. This Stability-Constrained SO (SCSO) is applicable to many OpenSim models and can potentially produce more biofidelic results than SO alone, especially when antagonistic muscle co-contraction is required to stabilize body joints. This hypothesis was tested using existing models and experimental data in the literature. Muscle activations were calculated by SO and SCSO for a spine model during two series of static trials (i.e. simulation 1 and 2), and also for a lower limb model (supplementary material 2). In simulation 1, symmetric and asymmetric flexion postures were compared, while in simulation 2, various external load heights were compared, where increases in load height did not change the external lumbar flexion moment, but necessitated higher EMG activations. During the tasks in simulation 1, the predicted muscle activations by SCSO demonstrated less average deviation from the EMG data (6.8% −7.5%) compared to those from SO (10.2%). In simulation 2, SO predicts constant muscle activations and forces, while SCSO predicts increases in the average activations of back and abdominal muscles that better match experimental data. Although the SCSO results are sensitive to some parameters (e.g. musculotendon stiffness), when considering the strategy of the central nervous system in distributing muscle forces and in activating antagonistic muscles, the assigned activations by SCSO are more biofidelic than SO.     

Does a linear position transducer placed on a stick and belt provide sufficient validity and reliability of countermovement jump performance outcomes?

Manufacturers recommend that linear position transducers (LPTs) should be placed on the side of a barbell (or wooden dowel) to measure countermovement jump (CMJ) height, but the validity and reliability of this placement have not been compared to other attachment sites. Since this recommended attachment site is far from the centre of mass, a belt attachment where the LPT is placed between the feet may increase the validity and reliability of CMJ data. Thirty-six physical education students participated in the study (24.6 ± 4.3 years; 177.0 ± 7.7 cm; 77.2 ± 9.0 kg). Parameters from the two LPT attachments (barbell and belt) were simultaneously validated to force plate data, where the nature of bias was analysed (systematic vs random). The within-session and between-session reliability of both attachment sites were compared to force plate data using a test-retest protocol of two sets of 5 CMJs separated by 7 days. The LPT provided highly reliable and valid measures of peak force, mean force, mean power, and jump height, where the bias was mostly systematic (r² > 0.7; ICC > 0.9). Peak velocity, mean velocity, and peak power were in very good agreement with the force plate and were highly reliable (r² > 0.5; ICC > 0.7). Therefore, both attachment sites produced similar results with a systematic bias compared to force plate data. Thus, both attachment sites seem to be valid for assessing CMJs when the measuring tool and site remain consistent across measurements. However, if LPT data are to be compared to force plate data, recalculation equations should be used.     

Score Tests for Exploring Complex Models: Application to HIV Dynamics Models

In biostatistics, more and more complex models are being developed. This is particularly the case in system biology. Fitting complex models can be very time‐consuming, since many models often have to be explored. Among the possibilities are the introduction of explanatory variables and the determination of random effects. The particularity of this use of the score test is that the null hypothesis is not itself very simple; typically, some random effects may be present under the null hypothesis. Moreover, the information matrix cannot be computed, but only an approximation based on the score. This article examines this situation with the specific example of HIV dynamics models. We examine the score test statistics for testing the effect of explanatory variables and the variance of random effect in this complex situation. We study type I errors and the statistical powers of this score test statistics and we apply the score test approach to a real data set of HIV‐infected patients.     

An antitumor promoter from Moringa oleifera Lam.

In the course of studies on the isolation of bioactive compounds from Philippine plants, the seeds of Moringa oleifera Lam. were examined and from the ethanol extract were isolated the new O-ethyl-4-(α- -rhamnosyloxy)benzyl carbamate (1) together with seven known compounds, 4(α- -rhamnosyloxy)-benzyl isothiocyanate (2), niazimicin (3), niazirin (4), β-sitosterol (5), glycerol-1-(9-octadecanoate) (6), 3-O-(6′-O-oleoyl-β- -glucopyranosyl)-β-sitosterol (7), and β-sitosterol-3-O-β- -glucopyranoside (8). Four of the isolates (2, 3, 7, and 8), which were obtained in relatively good yields, were tested for their potential antitumor promoting activity using an in vitro assay which tested their inhibitory effects on Epstein–Barr virus-early antigen (EBV-EA) activation in Raji cells induced by the tumor promoter, 12-O-tetradecanoyl-phorbol-13-acetate (TPA). All the tested compounds showed inhibitory activity against EBV-EA activation, with compounds 2, 3 and 8 having shown very significant activities. Based on the in vitro results, niazimicin (3) was further subjected to in vivo test and found to have potent antitumor promoting activity in the two-stage carcinogenesis in mouse skin using 7,12-dimethylbenz(a)anthracene (DMBA) as initiator and TPA as tumor promoter. From these results, niazimicin (3) is proposed to be a potent chemo-preventive agent in chemical carcinogenesis.     

Enzymatic synthesis of lysophosphatidic acid and phosphatidic acid

Immobilised 1,3-specific lipase from Rhizopus arrhizus was used as catalyst for the esterification of -glycero-3-phosphate and fatty acid or fatty acid vinyl ester in a solvent-free system. With lauric acid vinyl ester as acyl donor, a_w<0.53 favored the synthesis of lysophosphatidic acid (1-acyl-rac-glycero-3-phosphate, LPA1) and the spontaneous acyl migration of the fatty acid on the molecule. Subsequent acylation by the enzyme resulted in high phosphatidic acid (1,2-diacyl-rac-glycero-3-phosphate, PA) formation and high total conversions (>95%). With oleic acid, maximum conversions of 55% were obtained at low water activities. Temperatures below melting point of the product favored precipitation and resulted in high final conversion and high product ratio [LPA/(PA+LPA)]. Thus, LPA was the only product with lauric acid vinyl ester as acyl donor at 25°C. Increased substrate ratio ( -glycero-3-phosphate/fatty acid) from 0.05 to 1 resulted in a higher ratio of LPA to PA formed, but a lower total conversion of -glycero-3-phosphate. Increased amounts of enzyme preparation did not result in higher esterification rates, probably due to high mass-transfer limitations.     

Glucose-induced period-doubling cascade in the electrical activity of pancreatic β-cells

 In the presence of stimulatory concentrations of glucose, the membrane potential of pancreatic β-cells may experience a transition from periods of rapid spike-like oscillations alternating with a pseudo-steady state to spike-only oscillations. Insulin secretion from β-cells closely correlates the periods of spike-like oscillations. The purpose of this paper is to study the mathematical structure which underlines this transitional stage in a pancreatic β-cell model. It is demonstrated that the transition can be chaotic but becomes more and more regular with increase in glucose. In particular, the system undergoes a reversed period-doubling cascade leading to the spike-only oscillations as the glucose concentration crosses a threshold. The transition interval in glucose concentration is estimated to be extremely small in terms of the rate of change for the calcium dynamics in the β-cells. The methods are based on the theory of unimodal maps and the geometric and asymptotic theories of singular perturbations. Received: 25 October 1996/Revised version: 18 August 1997     

Stature of sub-arctic birch in relation to growth rate, lifespan and tree form

BACKGROUND AND AIMS: Sub-arctic mountain birch Betula pubescens var. pumila communities in the North Atlantic region are of variable stature, ranging from prostrate scrubs to forests with trees up to 12 m high. Four hypotheses were tested, relating growth and population characteristics of sub-arctic birch woodland and scrub to tree stature; i.e. the variable stature of birch woods is due to differences in (1) the mean growth rate; (2) the age-related patterns of growth rate; (3) the life expectancy of stems; or (4) the tree form. Methods: A stratified random sample of 300 birch trees was drawn from the total population of indigenous birch woodlands and scrub in Iceland, yielding 286 valid sample genets. The population was divided into three sub-populations with dominant trees 0-2, 2-4 and 4-12 m tall, referred to as birch scrub, birch scrub-woodland and birch forest, respectively. KEY RESULTS: Trees in the scrub population were of more contorted growth form than birch in the scrub-woodland and forest populations. Mean growth rates, mean age and median life expectancies increased significantly with sub-population of greater tree stature. At the population level, annual increment and longevity of birch stems was apparently interrelated as the stems in vigorously growing birch sub-populations had a longer life expectancy than those of slower growth. However, no difference was observed between sub-populations in age-related patterns of extension growth rate. CONCLUSIONS: The results were consistent with hypotheses (1), (3) and (4), but hypothesis (2) was rejected. Hence, mountain birch of more vigorous growth attains a greater stature than birch of lesser increment due to faster extension growth rate and a longer lifespan. In addition, the more contorted stem form of scrub populations contributes to their low stature.     

Hydraulic responses to height growth in maritime pine trees

As trees grow taller, decreased xylem path conductance imposes a major constraint on plant water and carbon balance, and is thus a key factor underlying forest productivity decline with age. The responses of stomatal conductance, leaf area: sapwood area ratio (A_L : A_S) and soil–leaf water potential gradient (ΔΨ_S–L) to height growth were investigated in maritime pine trees. Extensive measurements of in situ sap flow, stomatal conductance and (non‐gravitational) needle water potential (_L = Ψ_L ? ρ_wgh) were made during 2 years in a chronosequence of four even‐aged stands, under both wet and dry soil conditions. Under wet soil conditions, _L was systematically lower in taller trees on account of differences in gravitational potential. In contrast, under dry soil conditions, our measurements clearly showed that _L was maintained above a minimum threshold value of ?2.0 MPa independently of tree height, thus limiting the range of compensatory change in ΔΨ_S–L. Although a decrease in the A_L : A_S ratio occurred with tree height, this compensation was not sufficient to prevent a decline in leaf‐specific hydraulic conductance, K_L (50% lower in 30 m trees than in 10 m trees). An associated decline in stomatal conductance with tree height thus occurred to maintain a balance between water supply and demand. Both the increased investment in non‐productive versus productive tissues (A_S : A_L) and stomatal closure may have contributed to the observed decrease in tree growth efficiency with increasing tree height (by a factor of three from smallest to tallest trees), although other growth‐limiting responses (e.g. soil nutrient sequestration, increased respiratory costs) cannot be excluded.