Local stress-strain distribution and load transfer across cartilage matrix at micro-scale using combined microscopy-based finite element method

Articular cartilage is the load-bearing tissue that consists of proteoglycan macromolecules entrapped between collagen fibrils in a three-dimensional architecture. To date, the drudgery of searching for mathematical models to represent the biomechanics of such a system continues without providing a fitting description of its functional response to load at micro-scale level. We believe that the major complication arose when cartilage was first envisaged as a multiphasic model with distinguishable components and that quantifying those and searching for the laws that govern their interaction is inadequate. To the thesis of this paper, cartilage as a bulk is as much continuum as is the response of its components to the external stimuli. For this reason, we framed the fundamental question as to what would be the mechano-structural functionality of such a system in the total absence of one of its key constituents-proteoglycans. To answer this, hydrated normal and proteoglycan depleted samples were tested under confined compression while finite element models were reproduced, for the first time, based on the structural microarchitecture of the cross-sectional profile of the matrices. These micro-porous in silico models served as virtual transducers to produce an internal noninvasive probing mechanism beyond experimental capabilities to render the matrices micromechanics and several others properties like permeability, orientation etc. The results demonstrated that load transfer was closely related to the microarchitecture of the hyperelastic models that represent solid skeleton stress and fluid response based on the state of the collagen network with and without the swollen proteoglycans. In other words, the stress gradient during deformation was a function of the structural pattern of the network and acted in concert with the position-dependent compositional state of the matrix. This reveals that the interaction between indistinguishable components in real cartilage is superimposed by its microarchitectural state which directly influences macromechanical behavior.

Selecting volatility forecasting models for portfolio allocation purposes

Techniques for evaluating and selecting multivariate volatility forecasts are not yet understood as well as their univariate counterparts. This paper considers the ability of different loss functions to discriminate between a set of competing forecasting models which are subsequently applied in a portfolio allocation context. It is found that a likelihood-based loss function outperforms its competitors, including those based on the given portfolio application. This result indicates that considering the particular application of forecasts is not necessarily the most effective basis on which to select models.

Kinematic models for manoeuvring and seakeeping of marine vessels

The motion of marine vessels has traditionally been studied using two different approaches: manoeuvring and seakeeping. These two approaches use different reference frames and coordinate systems to describe the motion. This paper derives the kinematic models that characterize the transformation of motion variables (position, velocity, accelerations) and forces between the different coordinate systems used in these theories. The derivations hereby presented are done in terms of the formalism adopted in robotics. The advantage of this formulation is the use of matrix notation and operations. As an application, the transformation of linear equations of motion used in seakeeping into body-fixed coordinates is considered for both zero and forward speed.

Myogel, a novel, basement membrane-rich, extracellular matrix derived from skeletal muscle, is highly adipogenic in vivo and in vitro

Background/Aims Biological and synthetic scaffolds play important roles in tissue engineering and are being developed towards human clinical applications. Based on previous work from our laboratory, we propose that extracellular matrices from skeletal muscle could be developed for adipose tissue engineering. Methods Extracellular matrices (Myogels) extracted from skeletal muscle of various species were assessed using biochemical assays including ELISA and Western blotting. Biofunctionality was assessed using an in vitro differentiation assay and a tissue engineering construct model in the rat. Results Myogels were successfully extracted from mice, rats, pigs and humans. Myogels contained significant levels of laminin α4- and α2-subunits and collagen I compared to Matrigel™, which contains laminin 1 (α1β1γ1) and collagen IV. Levels of growth factors such as fibroblast growth factor 2 were significantly higher than Matrigel, vascular endothelial growth factor-A levels were significantly lower and all other growth factors were comparable. Myogels reproducibly stimulated adipogenic differentiation of preadipocytes in vitro and the growth of adipose tissue in the rat. Conclusions We found Myogel induces adipocyte differentiation in vitroand shows strong adipogenic potential in vivo, inducing the growth of well-vascularised adipose tissue. Myogel offers an alternative for current support scaffolds in adipose tissue engineering, allowing the scaling up of animal models towards clinical adipose tissue engineering applications.

An MMP13-selective inhibitor delays primary tumor growth and the onset of tumor-associated osteolytic lesions in experimental models of breast cancer

We investigated the effects of the matrix metalloproteinase 13 (MMP13)-selective inhibitor, 5-(4-{4-[4-(4-fluorophenyl)-1,3-oxazol-2-yl]phenoxy}phenoxy)-5-(2-methoxyethyl) pyrimidine-2,4,6(1H,3H,5H)-trione (Cmpd-1), on the primary tumor growth and breast cancer-associated bone remodeling using xenograft and syngeneic mouse models. We used human breast cancer MDA-MB-231 cells inoculated into the mammary fat pad and left ventricle of BALB/c Nu/Nu mice, respectively, and spontaneously metastasizing 4T1.2-Luc mouse mammary cells inoculated into mammary fat pad of BALB/c mice. In a prevention setting, treatment with Cmpd-1 markedly delayed the growth of primary tumors in both models, and reduced the onset and severity of osteolytic lesions in the MDA-MB-231 intracardiac model. Intervention treatment with Cmpd-1 on established MDA-MB-231 primary tumors also significantly inhibited subsequent growth. In contrast, no effects of Cmpd-1 were observed on soft organ metastatic burden following intracardiac or mammary fat pad inoculations of MDA-MB-231 and 4T1.2-Luc cells respectively. MMP13 immunostaining of clinical primary breast tumors and experimental mice tumors revealed intra-tumoral and stromal expression in most tumors, and vasculature expression in all. MMP13 was also detected in osteoblasts in clinical samples of breast-to-bone metastases. The data suggest that MMP13-selective inhibitors, which lack musculoskeletal side effects, may have therapeutic potential both in primary breast cancer and cancer-induced bone osteolysis.

Models for studying cellular invasion of basement membranes

Invasion of extracellular matrices is crucial to a number of physiological and pathophysiological states, including tumor cell metastasis, arthritis, embryo implantation, wound healing, and early development. To isolate invasion from the additional complexities of these scenarios a number of in vitro invasion assays have been developed over the years. Early studies employed intact tissues, like denuded amniotic membrane (1) or embryonic chick heart fragments (2), however recently, purified matrix components or complex matrix extracts have been used to provide more uniform and often more rapid analyses (for examples, see the following integrin studies). Of course, the more holistic view of invasion offered in the earlier assays is valuable and cannot be fully reproduced in these more rapid assays, but advantages of reproducibility among replicates, ease of preparation and analysis, and overall high throughput favor the newer assays. In this chapter, we will focus on providing detailed protocols for Matrigel-based assays (Matrigel=reconstituted basement membrane; reviewed in ref. (3)). Matrigel is an extract from the transplantable Engelbreth-Holm-Swarm murine sarcoma that deposits a multilammelar basement membrane. Matrigel is available commercially (Becton Dickinson, Bedford, MA), and can be manipulated as a liquid at 4°C into a variety of different formats. Alternatively, cell culture inserts precoated with Matrigel can be purchased for even greater simplicity.

Three dimensional in-vitro models for studying cancer angiogenesis

Introduction Hydrogels prepared from star-shaped poly(ethylene glycol) (PEG) and maleimide-functionalized heparin provide a potential matrix for use in developing three dimensional (3D) models. We have previously demonstrated that these hydrogels support the cultivation of human umbilical vein endothelial cells (HUVECs). We extend this body of work to study the ability to create an extracellular matrix (ECM)-like model to study breast and prostate cancer cell growth in 3D. Also, we investigate the ability to produce a tri-culture mimicking tumour angiogenesis with cancer spheroids, HUVECs and mesenchymal stem cells (MSCs). Materials and Methods The breast cancer cell lines, MCF-7 and MDA-MB-231, and prostate cancer cell lines, LNCaP and PC3, were seeded into starPEG-heparin hydrogels and grown for 14 Days to analyze the effects of varying hydrogel stiffness on spheroid development. Resulting hydrogel constructs were analyzed via proliferation assays, light microscopy, and immunostaining. Cancer cell lines were then seeded into starPEG-heparin hydrogels functionalized with growth factors as spheroids with HUVECs and MSCs and grown as a tri-culture. Cultures were analyzed via immunostaining and observed using confocal microscopy. Results Cultures prepared in MMP-cleavable starPEG-heparin hydrogels display spheroid formation in contrast to adherent growth on tissue culture plastic. Small differences were visualized in cancer spheroid growth between different gel stiffness across the range of cell lines. Cancer cell lines were able to be co-cultivated with HUVECs and MSC. Interaction was visualized between tumours and HUVECs via confocal microscopy. Further studies intend to further optimize and mimic the ECM environment of in-situ tumour angiogenesis. Discussion Our results confirm the suitability of hydrogels constructed from starPEG-heparin for HUVEC and MSC co-cultivation with cancer cell lines to study cell-cell and cell-matrix interactions in a 3D environment. This represents a step forward in the development of 3D culture models to study the pathomechanisms of breast and prostate cancer.

Mixture models of nucleotide sequence evolution that account for heterogeneity in the substitution process across sites and across lineages

Molecular phylogenetic studies of homologous sequences of nucleotides often assume that the underlying evolutionary process was globally stationary, reversible, and homogeneous (SRH), and that a model of evolution with one or more site-specific and time-reversible rate matrices (e.g., the GTR rate matrix) is enough to accurately model the evolution of data over the whole tree. However, an increasing body of data suggests that evolution under these conditions is an exception, rather than the norm. To address this issue, several non-SRH models of molecular evolution have been proposed, but they either ignore heterogeneity in the substitution process across sites (HAS) or assume it can be modeled accurately using the distribution. As an alternative to these models of evolution, we introduce a family of mixture models that approximate HAS without the assumption of an underlying predefined statistical distribution. This family of mixture models is combined with non-SRH models of evolution that account for heterogeneity in the substitution process across lineages (HAL). We also present two algorithms for searching model space and identifying an optimal model of evolution that is less likely to over- or underparameterize the data. The performance of the two new algorithms was evaluated using alignments of nucleotides with 10 000 sites simulated under complex non-SRH conditions on a 25-tipped tree. The algorithms were found to be very successful, identifying the correct HAL model with a 75% success rate (the average success rate for assigning rate matrices to the tree's 48 edges was 99.25%) and, for the correct HAL model, identifying the correct HAS model with a 98% success rate. Finally, parameter estimates obtained under the correct HAL-HAS model were found to be accurate and precise. The merits of our new algorithms were illustrated with an analysis of 42 337 second codon sites extracted from a concatenation of 106 alignments of orthologous genes encoded by the nuclear genomes of Saccharomyces cerevisiae, S. paradoxus, S. mikatae, S. kudriavzevii, S. castellii, S. kluyveri, S. bayanus, and Candida albicans. Our results show that second codon sites in the ancestral genome of these species contained 49.1% invariable sites, 39.6% variable sites belonging to one rate category (V1), and 11.3% variable sites belonging to a second rate category (V2). The ancestral nucleotide content was found to differ markedly across these three sets of sites, and the evolutionary processes operating at the variable sites were found to be non-SRH and best modeled by a combination of eight edge-specific rate matrices (four for V1 and four for V2). The number of substitutions per site at the variable sites also differed markedly, with sites belonging to V1 evolving slower than those belonging to V2 along the lineages separating the seven species of Saccharomyces. Finally, sites belonging to V1 appeared to have ceased evolving along the lineages separating S. cerevisiae, S. paradoxus, S. mikatae, S. kudriavzevii, and S. bayanus, implying that they might have become so selectively constrained that they could be considered invariable sites in these species.

Predictive models for water sources with high susceptibility for bromine-containing disinfection by-product formation: Implications for water treatment

This study examines a matrix of synthetic water samples designed to include conditions that favour brominated disinfection by-product (Br-DBP) formation, in order to provide predictive models suitable for high Br-DBP forming waters such as salinity-impacted waters. Br-DBPs are known to be more toxic than their chlorinated analogues, in general, and their formation may be favoured by routine water treatment practices such as coagulation/flocculation under specific conditions; therefore, circumstances surrounding their formation must be understood. The chosen factors were bromide concentration, mineral alkalinity, bromide to dissolved organic carbon (Br/DOC) ratio and Suwannee River natural organic matter concentration. The relationships between these parameters and DBP formation were evaluated by response surface modelling of data generated using a face-centred central composite experimental design. Predictive models for ten brominated and/or chlorinated DBPs are presented, as well as models for total trihalomethanes (tTHMs) and total dihaloacetonitriles (tDHANs), and bromide substitution factors for the THMs and DHANs classes. The relationships described revealed that increasing alkalinity and increasing Br/DOC ratio were associated with increasing bromination of THMs and DHANs, suggesting that DOC lowering treatment methods that do not also remove bromide such as enhanced coagulation may create optimal conditions for Br-DBP formation in waters in which bromide is present.

Designing innovative business models with a framework that promotes experimentation

Purpose – Business models to date have remained the creation of management, however, it is the belief of the authors that designers should be critically approaching, challenging and creating new business models as part of their practice. This belief portrays a new era where business model constructs become the new design brief of the future and fuel design and innovation to work together at the strategic level of an organisation. Design/methodology/approach – The purpose of this paper is to explore and investigate business model design. The research followed a deductive structured qualitative content analysis approach utilizing a predetermined categorization matrix. The analysis of forty business cases uncovered commonalities of key strategic drivers behind these innovative business models. Findings – Five business model typologies were derived from this content analysis, from which quick prototypes of new business models can be created. Research limitations/implications – Implications from this research suggest there is no “one right” model, but rather through experimentation, the generation of many unique and diverse concepts can result in greater possibilities for future innovation and sustained competitive advantage. Originality/value – This paper builds upon the emerging research and exploration into the importance and relevance of dynamic, design-driven approaches to the creation of innovative business models. These models aim to synthesize knowledge gained from real world examples into a tangible, accessible and provoking framework that provide new prototyping templates to aid the process of business model experimentation.

Discussion of “Generalized estimating equations: Notes on the choice of the working correlation matrix”

Objective To discuss generalized estimating equations as an extension of generalized linear models by commenting on the paper of Ziegler and Vens "Generalized Estimating Equations. Notes on the Choice of the Working Correlation Matrix". Methods Inviting an international group of experts to comment on this paper. Results Several perspectives have been taken by the discussants. Econometricians have established parallels to the generalized method of moments (GMM). Statisticians discussed model assumptions and the aspect of missing data Applied statisticians; commented on practical aspects in data analysis. Conclusions In general, careful modeling correlation is encouraged when considering estimation efficiency and other implications, and a comparison of choosing instruments in GMM and generalized estimating equations, (GEE) would be worthwhile. Some theoretical drawbacks of GEE need to be further addressed and require careful analysis of data This particularly applies to the situation when data are missing at random.

The accuracy of varietal selection using factor analytic models for multi-environment plant breeding trials

Modeling of cultivar x trial effects for multienvironment trials (METs) within a mixed model framework is now common practice in many plant breeding programs. The factor analytic (FA) model is a parsimonious form used to approximate the fully unstructured form of the genetic variance-covariance matrix in the model for MET data. In this study, we demonstrate that the FA model is generally the model of best fit across a range of data sets taken from early generation trials in a breeding program. In addition, we demonstrate the superiority of the FA model in achieving the most common aim of METs, namely the selection of superior genotypes. Selection is achieved using best linear unbiased predictions (BLUPs) of cultivar effects at each environment, considered either individually or as a weighted average across environments. In practice, empirical BLUPs (E-BLUPs) of cultivar effects must be used instead of BLUPs since variance parameters in the model must be estimated rather than assumed known. While the optimal properties of minimum mean squared error of prediction (MSEP) and maximum correlation between true and predicted effects possessed by BLUPs do not hold for E-BLUPs, a simulation study shows that E-BLUPs perform well in terms of MSEP.

Estimation in a multiplicative mixed model involving a genetic relationship matrix.

Genetic models partitioning additive and non-additive genetic effects for populations tested in replicated multi-environment trials (METs) in a plant breeding program have recently been presented in the literature. For these data, the variance model involves the direct product of a large numerator relationship matrix A, and a complex structure for the genotype by environment interaction effects, generally of a factor analytic (FA) form. With MET data, we expect a high correlation in genotype rankings between environments, leading to non-positive definite covariance matrices. Estimation methods for reduced rank models have been derived for the FA formulation with independent genotypes, and we employ these estimation methods for the more complex case involving the numerator relationship matrix. We examine the performance of differing genetic models for MET data with an embedded pedigree structure, and consider the magnitude of the non-additive variance. The capacity of existing software packages to fit these complex models is largely due to the use of the sparse matrix methodology and the average information algorithm. Here, we present an extension to the standard formulation necessary for estimation with a factor analytic structure across multiple environments.

Oxidative stability of phytosterols in food models and foods

An important safety aspect to be considered when foods are enriched with phytosterols and phytostanols is the oxidative stability of these lipid compounds, i.e. their resistance to oxidation and thus to the formation of oxidation products. This study concentrated on producing scientific data to support this safety evaluation process. In the absence of an official method for analyzing of phytosterol/stanol oxidation products, we first developed a new gas chromatographic - mass spectrometric (GC-MS) method. We then investigated factors affecting these compounds' oxidative stability in lipid-based food models in order to identify critical conditions under which significant oxidation reactions may occur. Finally, the oxidative stability of phytosterols and stanols in enriched foods during processing and storage was evaluated. Enriched foods covered a range of commercially available phytosterol/stanol ingredients, different heat treatments during food processing, and different multiphase food structures. The GC-MS method was a powerful tool for measuring the oxidative stability. Data obtained in food model studies revealed that the critical factors for the formation and distribution of the main secondary oxidation products were sterol structure, reaction temperature, reaction time, and lipid matrix composition. Under all conditions studied, phytostanols as saturated compounds were more stable than unsaturated phytosterols. In addition, esterification made phytosterols more reactive than free sterols at low temperatures, while at high temperatures the situation was the reverse. Generally, oxidation reactions were more significant at temperatures above 100°C. At lower temperatures, the significance of these reactions increased with increasing reaction time. The effect of lipid matrix composition was dependent on temperature; at temperatures above 140°C, phytosterols were more stable in an unsaturated lipid matrix, whereas below 140°C they were more stable in a saturated lipid matrix. At 140°C, phytosterols oxidized at the same rate in both matrices. Regardless of temperature, phytostanols oxidized more in an unsaturated lipid matrix. Generally, the distribution of oxidation products seemed to be associated with the phase of overall oxidation. 7-ketophytosterols accumulated when oxidation had not yet reached the dynamic state. Once this state was attained, the major products were 5,6-epoxyphytosterols and 7-hydroxyphytosterols. The changes observed in phytostanol oxidation products were not as informative since all stanol oxides quantified represented hydroxyl compounds. The formation of these secondary oxidation products did not account for all of the phytosterol/stanol losses observed during the heating experiments, indicating the presence of dimeric, oligomeric or other oxidation products, especially when free phytosterols and stanols were heated at high temperatures. Commercially available phytosterol/stanol ingredients were stable during such food processes as spray-drying and ultra high temperature (UHT)-type heating and subsequent long-term storage. Pan-frying, however, induced phytosterol oxidation and was classified as a rather deteriorative process. Overall, the findings indicated that although phytosterols and stanols are stable in normal food processing conditions, attention should be paid to their use in frying. Complex interactions between other food constituents also suggested that when new phytosterol-enriched foods are developed their oxidative stability must first be established. The results presented here will assist in choosing safe conditions for phytosterol/stanol enrichment.

Top quark mass measurement in the tt̅ all hadronic channel using a matrix element technique in pp̅ collisions at √s=1.96 TeV

We present a measurement of the top quark mass in the all-hadronic channel (\tt $\to$ \bb$q_{1}\bar{q_{2}}q_{3}\bar{q_{4}}$) using 943 pb$^{-1}$ of \ppbar collisions at $\sqrt {s} = 1.96$ TeV collected at the CDF II detector at Fermilab (CDF). We apply the standard model production and decay matrix-element (ME) to $\ttbar$ candidate events. We calculate per-event probability densities according to the ME calculation and construct template models of signal and background. The scale of the jet energy is calibrated using additional templates formed with the invariant mass of pairs of jets. These templates form an overall likelihood function that depends on the top quark mass and on the jet energy scale (JES). We estimate both by maximizing this function. Given 72 observed events, we measure a top quark mass of 171.1 $\pm$ 3.7 (stat.+JES) $\pm$ 2.1 (syst.) GeV/$c^{2}$. The combined uncertainty on the top quark mass is 4.3 GeV/$c^{2}$.