Thermal activation and ATP dependence of the cytoskeleton remodeling dynamics

The cytoskeleton (CSK) is a nonequilibrium polymer network that uses hydrolyzable sources of free energy such as adenosine triphosphate (ATP) to remodel its internal structure. As in inert nonequilibrium soft materials, CSK remodeling has been associated with structural rearrangements driven by energy-activated processes. We carry out particle tracking and traction microscopy measurements of alveolar epithelial cells at various temperatures and ATP concentrations. We provide the first experimental evidence that the remodeling dynamics of the CSK is driven by structural rearrangements over free-energy barriers induced by thermally activated forces mediated by ATP. The measured activation energy of these forces is ~40kBTr (kB being the Boltzmann constant and Tr being the room temperature). Our experiments provide clues to understand the analogy between the dynamics of the living CSK and that of inert nonequilibrium soft materials.

Effect of expiratory flow limitation on respiratory mechanical impedance: a model study

Large phasic variations of respiratory mechanical impedance (Zrs) have been observed during induced expiratory flow limitation (EFL) (M. Vassiliou, R. Peslin, C. Saunier, and C. Duvivier. Eur. Respir. J. 9: 779-786, 1996). To clarify the meaning of Zrs during EFL, we have measured from 5 to 30 Hz the input impedance (Zin) of mechanical analogues of the respiratory system, including flow-limiting elements (FLE) made of easily collapsible rubber tubing. The pressures upstream (Pus) and downstream (Pds) from the FLE were controlled and systematically varied. Maximal flow (Vmax) increased linearly with Pus, was close to the value predicted from wave-speed theory, and was obtained for Pus-Pds of 4-6 hPa. The real part of Zin started increasing abruptly with flow (V) >85%Vmax and either further increased or suddenly decreased in the vicinity of V¿max. The imaginary part of Zin decreased markedly and suddenly above 95%Vmax. Similar variations of Zin during EFL were seen with an analogue that mimicked the changes of airway transmural pressure during breathing. After pressure andV measurements upstream and downstream from the FLE were combined, the latter was analyzed in terms of a serial (Zs) and a shunt (Zp) compartment. Zs was consistent with a large resistance and inertance, and Zp with a mainly elastic element having an elastance close to that of the tube walls. We conclude that Zrs data during EFL mainly reflect the properties of the FLE.

Genome sequence of the pea aphid Acyrthosiphon pisum

Aphids are important agricultural pests and also biological models for studies of insect-plant interactions, symbiosis, virus vectoring, and the developmental causes of extreme phenotypic plasticity. Here we present the 464 Mb draft genome assembly of the pea aphid Acyrthosiphon pisum. This first published whole genome sequence of a basal hemimetabolous insect provides an outgroup to the multiple published genomes of holometabolous insects. Pea aphids are host-plant specialists, they can reproduce both sexually and asexually, and they have coevolved with an obligate bacterial symbiont. Here we highlight findings from whole genome analysis that may be related to these unusual biological features. These findings include discovery of extensive gene duplication in more than 2000 gene families as well as loss of evolutionarily conserved genes. Gene family expansions relative to other published genomes include genes involved in chromatin modification, miRNA synthesis, and sugar transport. Gene losses include genes central to the IMD immune pathway, selenoprotein utilization, purine salvage, and the entire urea cycle. The pea aphid genome reveals that only a limited number of genes have been acquired from bacteria; thus the reduced gene count of Buchnera does not reflect gene transfer to the host genome. The inventory of metabolic genes in the pea aphid genome suggests that there is extensive metabolite exchange between the aphid and Buchnera, including sharing of amino acid biosynthesis between the aphid and Buchnera. The pea aphid genome provides a foundation for post-genomic studies of fundamental biological questions and applied agricultural problems.

Modeling the Emergence of Social Structure from a Phylogenetic Point of View

Based on provious (Hemelrijk 1998; Puga-González, Hildenbrant & Hemelrijk 2009), we have developed an agent-based model and software, called A-KinGDom, which allows us to simulate the emergence of the social structure in a group of non-human primates. The model includes dominance and affiliative interactions and incorporate s two main innovations (preliminary dominance interactions and a kinship factor), which allow us to define four different attack and affiliative strategies. In accordance with these strategies, we compared the data obtained under four simulation conditions with the results obtained in a provious study (Dolado & Beltran 2012) involving empirical observations of a captive group of mangabeys (Cercocebus torquatus)

Two different 3D biogeochemical models for the NW mediterranean sea: validation with Meris data and intercomparison

Report for the scientific sojourn carried out at the University of New South Wales from February to June the 2007. Two different biogeochemical models are coupled to a three dimensional configuration of the Princeton Ocean Model (POM) for the Northwestern Mediterranean Sea (Ahumada and Cruzado, 2007). The first biogeochemical model (BLANES) is the three-dimensional version of the model described by Bahamon and Cruzado (2003) and computes the nitrogen fluxes through six compartments using semi-empirical descriptions of biological processes. The second biogeochemical model (BIOMEC) is the biomechanical NPZD model described in Baird et al. (2004), which uses a combination of physiological and physical descriptions to quantify the rates of planktonic interactions. Physical descriptions include, for example, the diffusion of nutrients to phytoplankton cells and the encounter rate of predators and prey. The link between physical and biogeochemical processes in both models is expressed by the advection-diffusion of the non-conservative tracers. The similarities in the mathematical formulation of the biogeochemical processes in the two models are exploited to determine the parameter set for the biomechanical model that best fits the parameter set used in the first model. Three years of integration have been carried out for each model to reach the so called perpetual year run for biogeochemical conditions. Outputs from both models are averaged monthly and then compared to remote sensing images obtained from sensor MERIS for chlorophyll.

Reaction-diffusion wave fronts: Multigeneration biological species under climate change

A generalization of reaction-diffusion models to multigeneration biological species is presented. It is based on more complex random walks than those in previous approaches. The new model is developed analytically up to infinite order. Our predictions for the speed agree to experimental data for several butterfly species better than existing models. The predicted dependence for the speed on the number of generations per year allows us to explain the change in speed observed for a specific invasion

From SNPs to pathways: integration of functional effect of sequence variations on models of cell signalling pathways

Background: Single nucleotide polymorphisms (SNPs) are the most frequent type of sequence variation between individuals, and represent a promising tool for finding genetic determinants of complex diseases and understanding the differences in drug response. In this regard, it is of particular interest to study the effect of non-synonymous SNPs in the context of biological networks such as cell signalling pathways. UniProt provides curated information about the functional and phenotypic effects of sequence variation, including SNPs, as well as on mutations of protein sequences. However, no strategy has been developed to integrate this information with biological networks, with the ultimate goal of studying the impact of the functional effect of SNPs in the structure and dynamics of biological networks. Results: First, we identified the different challenges posed by the integration of the phenotypic effect of sequence variants and mutations with biological networks. Second, we developed a strategy for the combination of data extracted from public resources, such as UniProt, NCBI dbSNP, Reactome and BioModels. We generated attribute files containing phenotypic and genotypic annotations to the nodes of biological networks, which can be imported into network visualization tools such as Cytoscape. These resources allow the mapping and visualization of mutations and natural variations of human proteins and their phenotypic effect on biological networks (e.g. signalling pathways, protein-protein interaction networks, dynamic models). Finally, an example on the use of the sequence variation data in the dynamics of a network model is presented. Conclusion: In this paper we present a general strategy for the integration of pathway and sequence variation data for visualization, analysis and modelling purposes, including the study of the functional impact of protein sequence variations on the dynamics of signalling pathways. This is of particular interest when the SNP or mutation is known to be associated to disease. We expect that this approach will help in the study of the functional impact of disease-associated SNPs on the behaviour of cell signalling pathways, which ultimately will lead to a better understanding of the mechanisms underlying complex diseases.

Identifying the preferred subset of enzymatic profiles in nonlinear kinetic metabolic models via multiobjective global optimization and pareto filters

Optimization models in metabolic engineering and systems biology focus typically on optimizing a unique criterion, usually the synthesis rate of a metabolite of interest or the rate of growth. Connectivity and non-linear regulatory effects, however, make it necessary to consider multiple objectives in order to identify useful strategies that balance out different metabolic issues. This is a fundamental aspect, as optimization of maximum yield in a given condition may involve unrealistic values in other key processes. Due to the difficulties associated with detailed non-linear models, analysis using stoichiometric descriptions and linear optimization methods have become rather popular in systems biology. However, despite being useful, these approaches fail in capturing the intrinsic nonlinear nature of the underlying metabolic systems and the regulatory signals involved. Targeting more complex biological systems requires the application of global optimization methods to non-linear representations. In this work we address the multi-objective global optimization of metabolic networks that are described by a special class of models based on the power-law formalism: the generalized mass action (GMA) representation. Our goal is to develop global optimization methods capable of efficiently dealing with several biological criteria simultaneously. In order to overcome the numerical difficulties of dealing with multiple criteria in the optimization, we propose a heuristic approach based on the epsilon constraint method that reduces the computational burden of generating a set of Pareto optimal alternatives, each achieving a unique combination of objectives values. To facilitate the post-optimal analysis of these solutions and narrow down their number prior to being tested in the laboratory, we explore the use of Pareto filters that identify the preferred subset of enzymatic profiles. We demonstrate the usefulness of our approach by means of a case study that optimizes the ethanol production in the fermentation of Saccharomyces cerevisiae.

On the Rule of Mixtures for Predicting Stress-Softening and Residual Strain Effects in Biological Tissues and Biocompatible Materials

In this work, we use the rule of mixtures to develop an equivalent material model in which the total strain energy density is split into the isotropic part related to the matrix component and the anisotropic energy contribution related to the fiber effects. For the isotropic energy part, we select the amended non-Gaussian strain energy density model, while the energy fiber effects are added by considering the equivalent anisotropic volumetric fraction contribution, as well as the isotropized representation form of the eight-chain energy model that accounts for the material anisotropic effects. Furthermore, our proposed material model uses a phenomenological non-monotonous softening function that predicts stress softening effects and has an energy term, derived from the pseudo-elasticity theory, that accounts for residual strain deformations. The model’s theoretical predictions are compared with experimental data collected from human vaginal tissues, mice skin, poly(glycolide-co-caprolactone) (PGC25 3-0) and polypropylene suture materials and tracheal and brain human tissues. In all cases examined here, our equivalent material model closely follows stress-softening and residual strain effects exhibited by experimental data

The Role of the Delay Time in the Modeling of Biological Range Expansions

The time interval between successive migrations of biological species causes a delay time in the reaction-diffusion equations describing their space-time dynamics. This lowers the predicted speed of the waves of advance, as compared to classical models. It has been shown that this delay-time effect improves the modeling of human range expansions. Here, we demonstrate that it can also be important for other species. We present two new examples where the predictions of the time-delayed and the classical (Fisher) approaches are compared to experimental data. No free or adjustable parameters are used. We show that the importance of the delay effect depends on the dimensionless product of the initial growth rate and the delay time. We argue that the delay effect should be taken into account in the modeling of range expansions for biological species