A New Computational Method for Membrane Compressibility: Bilayer Mechanical Thickness Revisited

Because lipid bilayers can bend and stretch in ways similar to thin elastic sheets, physical models of bilayer deformation have utilized mechanical constants such as the moduli for bending rigidity $\left({\kappa }_C\right)$ and area compressibility (K_A). However, the use of these models to quantify the energetics of membrane deformation associated with protein-membrane interactions, and the membrane response to stress is often hampered by the shortage of experimental data suitable for the estimation of the mechanical constants of various lipid mixtures. Although computational tools such as molecular dynamics simulations can provide alternative means to estimate K_A values, current approaches suffer significant technical limitations. Here, we present a novel, to our knowledge, computational framework that allows for a direct estimation of K_A values for individual bilayer leaflets. The theory is based on the concept of elasticity and derives K_A from real-space analysis of local thickness fluctuations sampled in molecular dynamics simulations. We explore and validate the model on a large set of single and multicomponent bilayers of different lipid compositions and sizes, simulated at different temperatures. The calculated bilayer compressibility moduli agree with values estimated previously from experiments and those obtained from a standard computational method based on a series of constrained tension simulations. We further validate our framework in a comparison with an existing polymer brush model and confirm the polymer brush model’s predicted linear relationship with proportionality coefficient of 24, using elastic parameters calculated from the simulation trajectories. The robustness of the results that emerge from the method allows us to revisit the origins of the bilayer mechanical (compressible) thickness and in particular its dependence on acyl-chain unsaturation and the presence of cholesterol.     

The Mechanics of Metastasis: Insights from a Computational Model

Although it may seem obvious that mechanical forces are required to drive metastatic cell movements, understanding of the mechanical aspects of metastasis has lagged far behind genetic and biochemical knowledge. The goal of this study is to learn about the mechanics of metastasis using a cell-based finite element model that proved useful for advancing knowledge about the forces that drive embryonic cell and tissue movements. Metastasis, the predominant cause of cancer-related deaths, involves a series of mechanical events in which one or more cells dissociate from a primary tumour, migrate through normal tissue, traverse in and out of a multi-layer circulatory system vessel and resettle. The present work focuses on the dissemination steps, from dissociation to circulation. The model shows that certain surface tension relationships must be satisfied for cancerous cells to dissociate from a primary tumour and that these equations are analogous to those that govern dissociation of embryonic cells. For a dissociated cell to then migrate by invadopodium extension and contraction and exhibit the shapes seen in experiments, the invadopodium must generate a contraction equal to approximately twice that produced by the interfacial tension associated with surrounding cells. Intravasation through the wall of a vessel is governed by relationships akin to those in the previous two steps, while release from the vessel wall is governed by equations that involve surface and interfacial tensions. The model raises a number of potential research questions. It also identifies how specific mechanical properties and the sub-cellular structural components that give rise to them might be changed so as to thwart particular metastatic steps and thereby block the spread of cancer.     

Intuitive Logic Revisited: New Data and a Bayesian Mixed Model Meta-Analysis

Recent research on syllogistic reasoning suggests that the logical status (valid vs. invalid) of even difficult syllogisms can be intuitively detected via differences in conceptual fluency between logically valid and invalid syllogisms when participants are asked to rate how much they like a conclusion following from a syllogism (Morsanyi & Handley, 2012). These claims of an intuitive logic are at odds with most theories on syllogistic reasoning which posit that detecting the logical status of difficult syllogisms requires effortful and deliberate cognitive processes. We present new data replicating the effects reported by Morsanyi and Handley, but show that this effect is eliminated when controlling for a possible confound in terms of conclusion content. Additionally, we reanalyze three studies () without this confound with a Bayesian mixed model meta-analysis (i.e., controlling for participant and item effects) which provides evidence for the null-hypothesis and against Morsanyi and Handley''s claim.     

Chaste: An Open Source C++ Library for Computational Physiology and Biology

Chaste — Cancer, Heart And Soft Tissue Environment — is an open source C++ library for the computational simulation of mathematical models developed for physiology and biology. Code development has been driven by two initial applications: cardiac electrophysiology and cancer development. A large number of cardiac electrophysiology studies have been enabled and performed, including high-performance computational investigations of defibrillation on realistic human cardiac geometries. New models for the initiation and growth of tumours have been developed. In particular, cell-based simulations have provided novel insight into the role of stem cells in the colorectal crypt. Chaste is constantly evolving and is now being applied to a far wider range of problems. The code provides modules for handling common scientific computing components, such as meshes and solvers for ordinary and partial differential equations (ODEs/PDEs). Re-use of these components avoids the need for researchers to ‘re-invent the wheel’ with each new project, accelerating the rate of progress in new applications. Chaste is developed using industrially-derived techniques, in particular test-driven development, to ensure code quality, re-use and reliability. In this article we provide examples that illustrate the types of problems Chaste can be used to solve, which can be run on a desktop computer. We highlight some scientific studies that have used or are using Chaste, and the insights they have provided. The source code, both for specific releases and the development version, is available to download under an open source Berkeley Software Distribution (BSD) licence at http://www.cs.ox.ac.uk/chaste, together with details of a mailing list and links to documentation and tutorials.

This is a PLOS Computational Biology Software Article

     

Development of Current Concepts of Physiology of Cerebral Circulation: A Comparative Analysis

The paper presents historical review on evolution of the concepts about the structure and functioning of the cerebral circulation system, which are formed both on the basis of experimental physiological data and as a result of observations on human. The advance in this field was, in many aspects, determined by development of the methodical base that provided, at first, only performance of experiments on animals. However, at present, a possibility has appeared to carry out parallel studies on animals and human with use of computer methods of recording and processing of the information, which increases their scientific and applied significance. At present, it can be believed firmly established that the specific intensity of the total cerebral circulation rises in the phylogenetic line of vertebrates and this regularity is repeated in ontogeny of both higher mammals and human. Several factors take part in the cerebral circulation regulation. The neurogenic regulatory mechanism is realized in full measure in adult individuals, whereas the metabolic control is already developed in early ontogeny. The high degree of the local circulation heterogeneity, its rhythmic pulsation, and the latent period liability of functional responses are characteristic features of the cerebral circulation system in higher animals. In man, of particular development is close relationship of the vascular and cerebrospinal fluid systems in the single craniospinal cavity that provides additional mechanisms of maintenance of the optimal blood flow through the brain.     

A Patient-Specific Computational Fluid Dynamic Model of Middle Cerebral Artery Aneurysm Before and One Year After Surgery

Computational fluid dynamics (CFD) has been widely used for studying intracranial aneurysm hemodynamics, while its use for guiding clinical strategy is still in development. In this study, CFD simulations helped informtreatment decision for a middle cerebral artery (MCA) aneurysm case was investigated. A patient with a 10.4 × 9.8 mm aneurysm attached with a small aneurysmat the edge of the trifurcation in the left MCA was included in this study. Forremoving the MCA aneurysm, two scenarios were considered: Plan-A involvedclipping the small aneurysm and Plan-B involved clipping the whole aneurysm.A suitable treatment plan was decided by comparing the clinical measurementsand CFD analysis between these two plans. One-year after the surgery, theCFD analysis was conducted again on the post-operative aneurysm model to verify the selected surgical plan in terms of morphometric and hemodynamic properties changes in the aneurysm. Based on the CFD simulation and clinicalexperience, surgical Plan-A was adopted. One-year after the surgery, both thehemodynamic and morphological properties improved in the post-operativeaneurysm model, indicating the recovery of the patient. The patient-specificaneurysm CFD analysis can help to determine a better surgical plan for patientswith special cerebral aneurysms. This study showed how CFD analysis can beused to aid clinical diagnosis and treatment.     

A New Computational Model for Neuro-Glio-Vascular Coupling: Astrocyte Activation Can Explain Cerebral Blood Flow Nonlinear Response to Interictal Events

Developing a clear understanding of the relationship between cerebral blood flow (CBF) response and neuronal activity is of significant importance because CBF increase is essential to the health of neurons, for instance through oxygen supply. This relationship can be investigated by analyzing multimodal (fMRI, PET, laser Doppler…) recordings. However, the important number of intermediate (non-observable) variables involved in the underlying neurovascular coupling makes the discovery of mechanisms all the more difficult from the sole multimodal data. We present a new computational model developed at the population scale (voxel) with physiologically relevant but simple equations to facilitate the interpretation of regional multimodal recordings. This model links neuronal activity to regional CBF dynamics through neuro-glio-vascular coupling. This coupling involves a population of glial cells called astrocytes via their role in neurotransmitter (glutamate and GABA) recycling and their impact on neighboring vessels. In epilepsy, neuronal networks generate epileptiform discharges, leading to variations in astrocytic and CBF dynamics. In this study, we took advantage of these large variations in neuronal activity magnitude to test the capacity of our model to reproduce experimental data. We compared simulations from our model with isolated epileptiform events, which were obtained in vivo by simultaneous local field potential and laser Doppler recordings in rats after local bicuculline injection. We showed a predominant neuronal contribution for low level discharges and a significant astrocytic contribution for higher level discharges. Besides, neuronal contribution to CBF was linear while astrocytic contribution was nonlinear. Results thus indicate that the relationship between neuronal activity and CBF magnitudes can be nonlinear for isolated events and that this nonlinearity is due to astrocytic activity, highlighting the importance of astrocytes in the interpretation of regional recordings.     

Allele Frequencies at Microsatellite Loci: The Stepwise Mutation Model Revisited

We summarize available data on the frequencies of alleles at microsatellite loci in human populations and compare observed distributions of allele frequencies to those generated by a simulation of the stepwise mutation model. We show that observed frequency distributions at 108 loci are consistent with the results of the model under the assumption that mutations cause an increase or decrease in repeat number by one and under the condition that the product Nu, where N is the effective population size and u is the mutation rate, is larger than one. We show that the variance of the distribution of allele sizes is a useful estimator of Nu and performs much better than previously suggested estimators for the stepwise mutation model. In the data, there is no correlation between the mean and variance in allele size at a locus or between the number of alleles and mean allele size, which suggests that the mutation rate at these loci is independent of allele size.     

Socio-Ecological Webs and Sites of Sociality:Levins' Strategy of Model Building Revisited

This essay extends Levins' 1966 analysis of modelbuilding in ecology and evolutionary biology. Amodel, as the product of modeling, might bevalued according to its correspondence to reality. Yet Levins' emphasis on provisionality and changeredirects attention to the processes ofmodeling, through which scientists select and generatetheir problems, define their categories, collect theirdata, compare competing models, and present theirfindings. I identify several points where decisionsare required that are not determined by nature. Thisinvites examination of the social considerationsmodelers are reacting to at the sites of sociality.Modelers must weave socio-ecological webs so thatthe models can be seen to represent their subjectmatter at the same time as the modelers secure thesupport of colleagues, collaborators and institutions,and enjoin others to act upon their conclusions. Notonly do theory justification and theory generationmerge, but the joint project becomes simultaneouslyphilosophical and sociological.     

Dynamic Mechanisms of Cell Rigidity Sensing: Insights from a Computational Model of Actomyosin Networks

Cells modulate themselves in response to the surrounding environment like substrate elasticity, exhibiting structural reorganization driven by the contractility of cytoskeleton. The cytoskeleton is the scaffolding structure of eukaryotic cells, playing a central role in many mechanical and biological functions. It is composed of a network of actins, actin cross-linking proteins (ACPs), and molecular motors. The motors generate contractile forces by sliding couples of actin filaments in a polar fashion, and the contractile response of the cytoskeleton network is known to be modulated also by external stimuli, such as substrate stiffness. This implies an important role of actomyosin contractility in the cell mechano-sensing. However, how cells sense matrix stiffness via the contractility remains an open question. Here, we present a 3-D Brownian dynamics computational model of a cross-linked actin network including the dynamics of molecular motors and ACPs. The mechano-sensing properties of this active network are investigated by evaluating contraction and stress in response to different substrate stiffness. Results demonstrate two mechanisms that act to limit internal stress: (i) In stiff substrates, motors walk until they exert their maximum force, leading to a plateau stress that is independent of substrate stiffness, whereas (ii) in soft substrates, motors walk until they become blocked by other motors or ACPs, leading to submaximal stress levels. Therefore, this study provides new insights into the role of molecular motors in the contraction and rigidity sensing of cells.     

Simplifying a life cycle assessment of a mobile phone

Purpose

The possibilities for full life cycle assessment (LCA) of new Information and Communication Technology (ICT) products are often limited, so simplification approaches are needed. The aim of this paper is to investigate possible simplifications in LCA of a mobile phone and to use the results to discuss the possibilities of LCA simplifications for ICT products in a broader sense. Another aim is to identify processes and data that are sensitive to different methodological choices and assumptions related to the environmental impacts of a mobile phone.

Methods

Different approaches to a reference LCA of a mobile phone was tested: (1) excluding environmental impact categories, (2) excluding life cycle stages/processes, (3) using secondary process data from generic databases, (4) using input-output data and (5) using a simple linear relationship between mass and embodied emissions.

Results and discussion

It was not possible to identify one or a few impact categories representative of all others. If several impact categories would be excluded, information would be lost. A precautionary approach of not excluding impact categories is therefore recommended since impacts from the different life cycle stages vary between impact categories. Regarding use of secondary data for an ICT product similar to that studied here, we recommend prioritising collection of primary (specific) data on energy use during production and use, key component data (primarily integrated circuits) and process-specific data regarding raw material acquisition of specific metals (e.g. gold) and air transport. If secondary data are used for important processes, the scaling is crucial. The use of input-output data can be a considerable simplification and is probably best used to avoid data gaps when more specific data are lacking.

Conclusions

Further studies are needed to provide for simplified LCAs for ICT products. In particular, the end-of-life treatment stage need to be further addressed, as it could not be investigated here for all simplifications due to data gaps.     

Cerebral Processing of Prosodic Emotional Signals: Evaluation of a Network Model Using rTMS

A great number of functional imaging studies contributed to developing a cerebral network model illustrating the processing of prosody in the brain. According to this model, the processing of prosodic emotional signals is divided into three main steps, each related to different brain areas. The present study sought to evaluate parts of the aforementioned model by using low-frequency repetitive transcranial magnetic stimulation (rTMS) over two important brain regions identified by the model: the superior temporal cortex (Experiment 1) and the inferior frontal cortex (Experiment 2). The aim of both experiments was to reduce cortical activity in the respective brain areas and evaluate whether these reductions lead to measurable behavioral effects during prosody processing. However, results obtained in this study revealed no rTMS effects on the acquired behavioral data. Possible explanations for these findings are discussed in the paper.     

Computational Optogenetics: Empirically-Derived Voltage- and Light-Sensitive Channelrhodopsin-2 Model

Channelrhodospin-2 (ChR2), a light-sensitive ion channel, and its variants have emerged as new excitatory optogenetic tools not only in neuroscience, but also in other areas, including cardiac electrophysiology. An accurate quantitative model of ChR2 is necessary for in silico prediction of the response to optical stimulation in realistic tissue/organ settings. Such a model can guide the rational design of new ion channel functionality tailored to different cell types/tissues. Focusing on one of the most widely used ChR2 mutants (H134R) with enhanced current, we collected a comprehensive experimental data set of the response of this ion channel to different irradiances and voltages, and used these data to develop a model of ChR2 with empirically-derived voltage- and irradiance- dependence, where parameters were fine-tuned via simulated annealing optimization. This ChR2 model offers: 1) accurate inward rectification in the current-voltage response across irradiances; 2) empirically-derived voltage- and light-dependent kinetics (activation, deactivation and recovery from inactivation); and 3) accurate amplitude and morphology of the response across voltage and irradiance settings. Temperature-scaling factors (Q₁₀) were derived and model kinetics was adjusted to physiological temperatures. Using optical action potential clamp, we experimentally validated model-predicted ChR2 behavior in guinea pig ventricular myocytes. The model was then incorporated in a variety of cardiac myocytes, including human ventricular, atrial and Purkinje cell models. We demonstrate the ability of ChR2 to trigger action potentials in human cardiomyocytes at relatively low light levels, as well as the differential response of these cells to light, with the Purkinje cells being most easily excitable and ventricular cells requiring the highest irradiance at all pulse durations. This new experimentally-validated ChR2 model will facilitate virtual experimentation in neural and cardiac optogenetics at the cell and organ level and provide guidance for the development of in vivo tools.