Muscle moment-arms: a key element in muscle-force estimation.

A clear and rigorous definition of muscle moment-arms in the context of musculoskeletal systems modelling is presented, using classical mechanics and screw theory. The definition provides an alternative to the tendon excursion method, which can lead to incorrect moment-arms if used inappropriately due to its dependency on the choice of joint coordinates. The definition of moment-arms, and the presented construction method, apply to musculoskeletal models in which the bones are modelled as rigid bodies, the joints are modelled as ideal mechanical joints and the muscles are modelled as massless, frictionless cables wrapping over the bony protrusions, approximated using geometric surfaces. In this context, the definition is independent of any coordinate choice. It is then used to solve a muscle-force estimation problem for a simple 2D conceptual model and compared with an incorrect application of the tendon excursion method. The relative errors between the two solutions vary between 0% and 100%.

Prediction of Childhood Asthma Using Conditional Probability and Discrete Event Simulation

Abstract: Asthma prevalence in children and adolescents in Spain is 10-17%. It is the most common chronic illness during childhood. Prevalence has been increasing over the last 40 years and there is considerable evidence that, among other factors, continued exposure to cigarette smoke results in asthma in children. No statistical or simulation model exist to forecast the evolution of childhood asthma in Europe. Such a model needs to incorporate the main risk factors that can be managed by medical authorities, such as tobacco (OR = 1.44), to establish how they affect the present generation of children. A simulation model using conditional probability and discrete event simulation for childhood asthma was developed and validated by simulating realistic scenario. The parameters used for the model (input data) were those found in the bibliography, especially those related to the incidence of smoking in Spain. We also used data from a panel of experts from the Hospital del Mar (Barcelona) related to actual evolution and asthma phenotypes. The results obtained from the simulation established a threshold of a 15-20% smoking population for a reduction in the prevalence of asthma. This is still far from the current level in Spain, where 24% of people smoke. We conclude that more effort must be made to combat smoking and other childhood asthma risk factors, in order to significantly reduce the number of cases. Once completed, this simulation methodology can realistically be used to forecast the evolution of childhood asthma as a function of variation in different risk factors.

The clc element and related genomic islands in Proteobacteria

Genomic islands are large DNA segments, present in most bacterial genome, that are acquired via horizontal gene transfer and contribute to the rapid bacterial evolution and adaptation of the host cell. Here we focus on the clc element (or ICEclc), a 103‑kb genomic island first discovered in Pseudomonas knackmussii B13, as a model of this diverse group of mobile genetic elements. ICEclc is normally integrated in the host bacterial chromosome but can excise and transfer to a new host by conjugation. In this chapter we review the basic features of ICEclc, the mechanisms of its life‑style as well as evolutionary relationships with other known and unknown elements in a variety of Proteobacteria.

A qualitative continuous model of cellular auxin and brassinosteroid signaling and their crosstalk.

Motivation: Hormone pathway interactions are crucial in shaping plant development, such as synergism between the auxin and brassinosteroid pathways in cell elongation. Both hormone pathways have been characterized in detail, revealing several feedback loops. The complexity of this network, combined with a shortage of kinetic data, renders its quantitative analysis virtually impossible at present.Results: As a first step towards overcoming these obstacles, we analyzed the network using a Boolean logic approach to build models of auxin and brassinosteroid signaling, and their interaction. To compare these discrete dynamic models across conditions, we transformed them into qualitative continuous systems, which predict network component states more accurately and can accommodate kinetic data as they become available. To this end, we developed an extension for the SQUAD software, allowing semi-quantitative analysis of network states. Contrasting the developmental output depending on cell type-specific modulators enabled us to identify a most parsimonious model, which explains initially paradoxical mutant phenotypes and revealed a novel physiological feature.

Trace element chemistry and U-Pb dating of zircons from oceanic gabbros and their relationship with whole rock composition (Lanzo, Italian Alps)

The U-Pb ages and the trace element content of zircon U-Pb along with major and trace element whole rock data on gabbroic dikes from the Lanzo lherzolitic massif, N-Italy, have been determined to constrain crustal accretion in ocean-continent transition zones. Three Fe-Ti gabbros were dated from the central and the southern part of the massif providing middle Jurassic ages of 161 +/- 2, 158 +/- 2 and 163 +/- 1 Ma, which argue for magmatic activity over few millions of years. Zircon crystals are characterized by high but variable Th/U ratios, rare earth element patterns enriched in heavy rare earths, pronounced positive Ce and negative Eu-anomalies consistent with crystallization after substantial plagioclase fractionation. The zircon trace element composition coupled with whole rock chemistry was used to reconstruct the crystallization history of the gabbros. A number of gabbros crystallized in situ, and zircon precipitated from trapped, intercumulus liquid, while other gabbros represent residual liquids that were extracted from a cumulus pile and crystallized along syn-magmatic shear zones. We propose a model in which the emplacement mechanism of gabbroic rocks in ocean-continent transition zones evolves from in situ crystallization to stratified crystallization with efficient extraction of residual liquid along syn-magmatic shear zones. Such an evolution of the crystallization history is probably related to the thermal evolution of the underlying mantle lithosphere.

An operon of three transcriptional regulators controls horizontal gene transfer of the integrative and conjugative element ICEclc in Pseudomonas knackmussii B13.

The integrative and conjugative element ICEclc is a mobile genetic element in Pseudomonas knackmussii B13, and an experimental model for a widely distributed group of elements in Proteobacteria. ICEclc is transferred from specialized transfer competent cells, which arise at a frequency of 3-5% in a population at stationary phase. Very little is known about the different factors that control the transfer frequency of this ICE family. Here we report the discovery of a three-gene operon encoded by ICEclc, which exerts global control on transfer initiation. The operon consists of three consecutive regulatory genes, encoding a TetR-type repressor MfsR, a MarR-type regulator and a LysR-type activator TciR. We show that MfsR autoregulates expression of the operon, whereas TciR is a global activator of ICEclc gene expression, but no clear role was yet found for MarR. Deletion of mfsR increases expression of tciR and marR, causing the proportion of transfer competent cells to reach almost 100% and transfer frequencies to approach 1 per donor. mfsR deletion also caused a two orders of magnitude loss in population viability, individual cell growth arrest and loss of ICEclc. This indicates that autoregulation is an important feature maintaining ICE transfer but avoiding fitness loss. Bioinformatic analysis showed that the mfsR-marR-tciR operon is unique for ICEclc and a few highly related ICE, whereas tciR orthologues occur more widely in a large variety of suspected ICE among Proteobacteria.

Analysis of a finite volume method for a cross-diffusion model in population dynamics

The main goal of this paper is to propose a convergent finite volume method for a reactionâeuro"diffusion system with cross-diffusion. First, we sketch an existence proof for a class of cross-diffusion systems. Then the standard two-point finite volume fluxes are used in combination with a nonlinear positivity-preserving approximation of the cross-diffusion coefficients. Existence and uniqueness of the approximate solution are addressed, and it is also shown that the scheme converges to the corresponding weak solution for the studied model. Furthermore, we provide a stability analysis to study pattern-formation phenomena, and we perform two-dimensional numerical examples which exhibit formation of nonuniform spatial patterns. From the simulations it is also found that experimental rates of convergence are slightly below second order. The convergence proof uses two ingredients of interest for various applications, namely the discrete Sobolev embedding inequalities with general boundary conditions and a space-time $L^1$ compactness argument that mimics the compactness lemma due to Kruzhkov. The proofs of these results are given in the Appendix.

An active strain electromechanical model for cardiac tissue

We propose a finite element approximation of a system of partial differential equations describing the coupling between the propagation of electrical potential and large deformations of the cardiac tissue. The underlying mathematical model is based on the active strain assumption, in which it is assumed that a multiplicative decomposition of the deformation tensor into a passive and active part holds, the latter carrying the information of the electrical potential propagation and anisotropy of the cardiac tissue into the equations of either incompressible or compressible nonlinear elasticity, governing the mechanical response of the biological material. In addition, by changing from an Eulerian to a Lagrangian configuration, the bidomain or monodomain equations modeling the evolution of the electrical propagation exhibit a nonlinear diffusion term. Piecewise quadratic finite elements are employed to approximate the displacements field, whereas for pressure, electrical potentials and ionic variables are approximated by piecewise linear elements. Various numerical tests performed with a parallel finite element code illustrate that the proposed model can capture some important features of the electromechanical coupling, and show that our numerical scheme is efficient and accurate.

A semi-grand canonical Monte Carlo simulation model for ion binding to ionizable surfaces: Proton binding of carboxylated latex particles as a case study

In this paper, we present a computer simulation study of the ion binding process at an ionizable surface using a semi-grand canonical Monte Carlo method that models the surface as a discrete distribution of charged and neutral functional groups in equilibrium with explicit ions modelled in the context of the primitive model. The parameters of the simulation model were tuned and checked by comparison with experimental titrations of carboxylated latex particles in the presence of different ionic strengths of monovalent ions. The titration of these particles was analysed by calculating the degree of dissociation of the latex functional groups vs. pH curves at different background salt concentrations. As the charge of the titrated surface changes during the simulation, a procedure to keep the electroneutrality of the system is required. Here, two approaches are used with the choice depending on the ion selected to maintain electroneutrality: counterion or coion procedures. We compare and discuss the difference between the procedures. The simulations also provided a microscopic description of the electrostatic double layer (EDL) structure as a function of pH and ionic strength. The results allow us to quantify the effect of the size of the background salt ions and of the surface functional groups on the degree of dissociation. The non-homogeneous structure of the EDL was revealed by plotting the counterion density profiles around charged and neutral surface functional groups. © 2011 American Institute of Physics.

The proportional distribution in a cooperative model with external opportunities

[cat] En aquest article es considera un problema de cooperació entre agents on cada agent realitza una contribució (diners, capital, treball, esforç) per tal d'obtenir un benefici comú a repartir. El repartiment proporcional respecte a les contribucions és una distribució que pertany al nucli del joc cooperatiu associat. A partir d'aquest model bàsic s'introdueix un agent extern que pot realitzar una determinada aportació que serveix per avaluar el potencial benefici de cada subcoalició d'agents si aquest nou agent finalment entrés. Aquesta anàlisi pot produir que el poder relatiu dels agents hagi variat. en concret s'avalua si la distribució proporcional és encara robusta des del punt de vista de la seva pertinença al conjunt de negociació. Amb aquest objectiu, analitzem el problema utilitzant el model de joc cooperatius amb estructura de coalició. Donat que, en general, la distribució proporcional, no pertany al conjunt de negociació, s'estudia una condició suficient per a que així sigui. També enunciem una condició necessària, i finalment es proposa una condició suficient que garanteix que el repartiment proporcional és la única distribució existent dins del conjunt de negociació.

The proportional distribution in a cooperative model with external opportunities

[cat] En aquest article es considera un problema de cooperació entre agents on cada agent realitza una contribució (diners, capital, treball, esforç) per tal d'obtenir un benefici comú a repartir. El repartiment proporcional respecte a les contribucions és una distribució que pertany al nucli del joc cooperatiu associat. A partir d'aquest model bàsic s'introdueix un agent extern que pot realitzar una determinada aportació que serveix per avaluar el potencial benefici de cada subcoalició d'agents si aquest nou agent finalment entrés. Aquesta anàlisi pot produir que el poder relatiu dels agents hagi variat. en concret s'avalua si la distribució proporcional és encara robusta des del punt de vista de la seva pertinença al conjunt de negociació. Amb aquest objectiu, analitzem el problema utilitzant el model de joc cooperatius amb estructura de coalició. Donat que, en general, la distribució proporcional, no pertany al conjunt de negociació, s'estudia una condició suficient per a que així sigui. També enunciem una condició necessària, i finalment es proposa una condició suficient que garanteix que el repartiment proporcional és la única distribució existent dins del conjunt de negociació.

Robust Estimates of the Negative Binomial Model

We consider robust parametric procedures for univariate discrete distributions, focusing on the negative binomial model. The procedures are based on three steps: ?First, a very robust, but possibly inefficient, estimate of the model parameters is computed. ?Second, this initial model is used to identify outliers, which are then removed from the sample. ?Third, a corrected maximum likelihood estimator is computed with the remaining observations. The final estimate inherits the breakdown point (bdp) of the initial one and its efficiency can be significantly higher. Analogous procedures were proposed in [1], [2], [5] for the continuous case. A comparison of the asymptotic bias of various estimates under point contamination points out the minimum Neyman's chi-squared disparity estimate as a good choice for the initial step. Various minimum disparity estimators were explored by Lindsay [4], who showed that the minimum Neyman's chi-squared estimate has a 50% bdp under point contamination; in addition, it is asymptotically fully efficient at the model. However, the finite sample efficiency of this estimate under the uncontaminated negative binomial model is usually much lower than 100% and the bias can be strong. We show that its performance can then be greatly improved using the three step procedure outlined above. In addition, we compare the final estimate with the procedure described in

A cognitive hierarchy model of behavior in the action commitment game

We apply the cognitive hierarchy model of Camerer et al. (Q J Econ 119(3):861-898, 2004)-where players have different levels of reasoning-to Huck et al. (Games Econ Behav 38:240-264, 2002) discrete version of Hamilton and Slutsky (Games Econ Behav 2:29-46, 1990) action commitment game-a duopoly with endogenous timing of entry. We show that, for an empirically reasonable average number of thinking steps, the model rules out Stackelberg equilibria, generates Cournot outcomes including delay, and outcomes where the first mover commits to a quantity higher than Cournot but lower than Stackelberg leader. We show that a cognitive hierarchy model with quantal responses can explain the most important features of the experimental data on the action commitment game in (2002). In order to gauge the success of the model in fitting the data, we compare it to a noisy Nash model. We find that the cognitive hierarchy model with quantal responses fits the data better than the noisy Nash model.

El Turisme de memòria a Catalunya: realitat i perspectives d'un model incipient

Catalunya disposa a l’actualitat de nombrosos equipaments culturals destinats a la difusió memorial mitjançant la pràctica turística creant-se, tanmateix, les primeres estructures en xarxa per tal de fer que aquests espais disposin d’un relat conjunt, element indispensable en el qual recolzar-se a l’hora d’emprendre qualsevol estratègia de difusió. El model català, però, pot dir-se que està en fase embrionària, i des d’aquest moment en endavant ens correspon als professionals i acadèmics treballar de cara a la seva consolidació. Aquest estudi posa les bases per a la investigació en clau turística sobre els Espais de Memòria

Modeling the transport phenomena within arterial wall: porous media approach

The transport of macromolecules, such as low-density lipoprotein (LDL), and their accumulation in the layers of the arterial wall play a critical role in the creation and development of atherosclerosis. Atherosclerosis is a disease of large arteries e.g., the aorta, coronary, carotid, and other proximal arteries that involves a distinctive accumulation of LDL and other lipid-bearing materials in the arterial wall. Over time, plaque hardens and narrows the arteries. The flow of oxygen-rich blood to organs and other parts of the body is reduced. This can lead to serious problems, including heart attack, stroke, or even death. It has been proven that the accumulation of macromolecules in the arterial wall depends not only on the ease with which materials enter the wall, but also on the hindrance to the passage of materials out of the wall posed by underlying layers. Therefore, attention was drawn to the fact that the wall structure of large arteries is different than other vessels which are disease-resistant. Atherosclerosis tends to be localized in regions of curvature and branching in arteries where fluid shear stress (shear rate) and other fluid mechanical characteristics deviate from their normal spatial and temporal distribution patterns in straight vessels. On the other hand, the smooth muscle cells (SMCs) residing in the media layer of the arterial wall respond to mechanical stimuli, such as shear stress. Shear stress may affect SMC proliferation and migration from the media layer to intima. This occurs in atherosclerosis and intimal hyperplasia. The study of blood flow and other body fluids and of heat transport through the arterial wall is one of the advanced applications of porous media in recent years. The arterial wall may be modeled in both macroscopic (as a continuous porous medium) and microscopic scales (as a heterogeneous porous medium). In the present study, the governing equations of mass, heat and momentum transport have been solved for different species and interstitial fluid within the arterial wall by means of computational fluid dynamics (CFD). Simulation models are based on the finite element (FE) and finite volume (FV) methods. The wall structure has been modeled by assuming the wall layers as porous media with different properties. In order to study the heat transport through human tissues, the simulations have been carried out for a non-homogeneous model of porous media. The tissue is composed of blood vessels, cells, and an interstitium. The interstitium consists of interstitial fluid and extracellular fibers. Numerical simulations are performed in a two-dimensional (2D) model to realize the effect of the shape and configuration of the discrete phase on the convective and conductive features of heat transfer, e.g. the interstitium of biological tissues. On the other hand, the governing equations of momentum and mass transport have been solved in the heterogeneous porous media model of the media layer, which has a major role in the transport and accumulation of solutes across the arterial wall. The transport of Adenosine 5´-triphosphate (ATP) is simulated across the media layer as a benchmark to observe how SMCs affect on the species mass transport. In addition, the transport of interstitial fluid has been simulated while the deformation of the media layer (due to high blood pressure) and its constituents such as SMCs are also involved in the model. In this context, the effect of pressure variation on shear stress is investigated over SMCs induced by the interstitial flow both in 2D and three-dimensional (3D) geometries for the media layer. The influence of hypertension (high pressure) on the transport of lowdensity lipoprotein (LDL) through deformable arterial wall layers is also studied. This is due to the pressure-driven convective flow across the arterial wall. The intima and media layers are assumed as homogeneous porous media. The results of the present study reveal that ATP concentration over the surface of SMCs and within the bulk of the media layer is significantly dependent on the distribution of cells. Moreover, the shear stress magnitude and distribution over the SMC surface are affected by transmural pressure and the deformation of the media layer of the aorta wall. This work reflects the fact that the second or even subsequent layers of SMCs may bear shear stresses of the same order of magnitude as the first layer does if cells are arranged in an arbitrary manner. This study has brought new insights into the simulation of the arterial wall, as the previous simplifications have been ignored. The configurations of SMCs used here with elliptic cross sections of SMCs closely resemble the physiological conditions of cells. Moreover, the deformation of SMCs with high transmural pressure which follows the media layer compaction has been studied for the first time. On the other hand, results demonstrate that LDL concentration through the intima and media layers changes significantly as wall layers compress with transmural pressure. It was also noticed that the fraction of leaky junctions across the endothelial cells and the area fraction of fenestral pores over the internal elastic lamina affect the LDL distribution dramatically through the thoracic aorta wall. The simulation techniques introduced in this work can also trigger new ideas for simulating porous media involved in any biomedical, biomechanical, chemical, and environmental engineering applications.