Implications of different membrane compartmentalization models in particle-based in silico studies.

Studying membrane dynamics is important to understand the cellular response to environmental stimuli. A decisive spatial characteristic of the plasma membrane is its compartmental structure created by the actin-based membrane-skeleton (fences) and anchored transmembrane proteins (pickets). Particle-based reaction-diffusion simulation of the membrane offers a suitable temporal and spatial resolution to analyse its spatially heterogeneous and stochastic dynamics. Fences have been modelled via hop probabilities, potentials or explicit picket fences. Our study analyses the different approaches' constraints and their impact on simulation results and performance. Each of the methods comes with its own constraints; the picket fences require small timesteps, potential fences might induce a bias in diffusion in crowded systems, and probabilistic fences, in addition to carefully scaling the probability with the timesteps, induce higher computational costs for each propagation step.

Long-term stimulation with alternating electric fields modulates the differentiation and mineralization of human pre-osteoblasts.

Biophysical stimulation by electric fields can promote bone formation in bone defects of critical size. Even though, long-term effects of alternating electric fields on the differentiation of osteoblasts are not fully understood. Human pre-osteoblasts were stimulated over 31 days to gain more information about these cellular processes. An alternating electric field with 0.7 Vrms and 20 Hz at two distances was applied and viability, mineralization, gene expression, and protein release of differentiation factors were analyzed. The viability was enhanced during the first days of stimulation. A higher electric field resulted in upregulation of typical osteogenic markers like osteoprotegerin, osteopontin, and interleukin-6, but no significant changes in mineralization. Upregulation of the osteogenic markers could be detected with a lower electric field after the first days of stimulation. As a significant increase in the mineralized matrix was identified, an enhanced osteogenesis due to low alternating electric fields can be assumed.

Receptor/Raft Ratio Is a Determinant for LRP6 Phosphorylation and WNT/ß-Catenin Signaling.

Microdomains or lipid rafts greatly affect the distribution of proteins and peptides in the membrane and play a vital role in the formation and activation of receptor/protein complexes. A prominent example for the decisive impact of lipid rafts on signaling is LRP6, whose localization to the same lipid rafts domain as the kinase CK1γ is crucial for its successful phosphorylation and the subsequent activation of the signalosome, hence WNT/ß-catenin signaling. However, according to various experimental measurements, approximately 25 to 35 % of the cell plasma membrane is covered by nanoscopic raft domains with diameters ranging between 10 to 200 nm. Extrapolating/Translating these values to the membrane of a "normal sized" cell yields a raft abundance, that, by far, outnumbers the membrane-associated pathway components of most individual signaling pathway, such as receptor and kinases. To analyze whether and how the quantitative ratio between receptor and rafts affects LRP6 phosphorylation and WNT/ß-catenin pathway activation, we present a computational modeling study, that for the first time employs realistic raft numbers in a compartment-based pathway model. Our simulation experiments indicate, that for receptor/raft ratios smaller than 1, i.e., when the number of raft compartments clearly exceeds the number of pathway specific membrane proteins, we observe significant decrease in LRP6 phosphorylation and downstream pathway activity. Our results suggest that pathway specific targeting and sorting mechanism are required to significantly narrow down the receptor/raft ratio and to enable the formation of the LRP6 signalosome, hence signaling.

Multi-Transcript Level Profiling Revealed Distinct mRNA, miRNA, and tRNA-Derived Fragment Bio-Signatures for Coping Behavior Linked Haplotypes in HPA Axis and Limbic System.

The molecular basis of porcine coping behavior (CB) relies on a sophisticated interplay of genetic and epigenetic features. Deep sequencing technologies allowed the identification of a plethora of new regulatory small non-coding RNA (sncRNA). We characterized mRNA and sncRNA profiles of central parts of the physiological stress response system including amygdala, hippocampus, hypothalamus and adrenal gland using systems biology for integration. Therefore, ten each of high- (HR) and low- (LR) reactive pigs (n = 20) carrying a CB associated haplotype in a prominent QTL-region on SSC12 were selected for mRNA and sncRNA expression profiling. The molecular markers related to the LR group included ATP1B2, MPDU1, miR-19b-5p, let-7g-5p, and 5'-tiRNA Leu in the adrenal gland, miR-194a-5p, miR-125a-5p, miR-7-1-5p, and miR-107-5p in the hippocampus and CBL and PVRL1 in the hypothalamus. Interestingly, amygdalae of the LR group showed 5'-tiRNA and 5'-tRF (5'-tRF Lys , 5'-tiRNA Lys , 5'-tiRNA Cys , and 5'-tiRNA Gln ) enrichment. Contrarily, molecular markers associated with the HR group encompassed miR-26b-5p, tRNA Arg , tRNA GlyiF in the adrenal gland, IGF1 and APOD in the amygdala and PBX1, TOB1, and C18orf1 in the hippocampus and miR-24 in the hypothalamus. In addition, hypothalami of the HR group were characterized by 3'-tiRNA enrichment (3'-tiRNAGln, 3'-tiRNA Asn , 3'-tiRNA Val , 3'-tRF Pro , 3'-tiRNA Cys , and 3'-tiRNA Ala ) and 3'-tRFs enrichment (3'-tRF Asn , 3'-tRF Glu , and 3'-tRF Val ). These evidence suggest that tRNA-derived fragments and their cleavage activity are a specific marker for coping behavior. Data integration revealed new bio-signatures of important molecular interactions on a multi-transcript level in HPA axis and limbic system of pigs carrying a CB-associated haplotype.

Relating simulation studies by provenance-Developing a family of Wnt signaling models.

For many biological systems, a variety of simulation models exist. A new simulation model is rarely developed from scratch, but rather revises and extends an existing one. A key challenge, however, is to decide which model might be an appropriate starting point for a particular problem and why. To answer this question, we need to identify entities and activities that contributed to the development of a simulation model. Therefore, we exploit the provenance data model, PROV-DM, of the World Wide Web Consortium and, building on previous work, continue developing a PROV ontology for simulation studies. Based on a case study of 19 Wnt/ß-catenin signaling models, we identify crucial entities and activities as well as useful metadata to both capture the provenance information from individual simulation studies and relate these forming a family of models. The approach is implemented in WebProv, a web application for inserting and querying provenance information. Our specialization of PROV-DM contains the entities Research Question, Assumption, Requirement, Qualitative Model, Simulation Model, Simulation Experiment, Simulation Data, and Wet-lab Data as well as activities referring to building, calibrating, validating, and analyzing a simulation model. We show that most Wnt simulation models are connected to other Wnt models by using (parts of) these models. However, the overlap, especially regarding the Wet-lab Data used for calibration or validation of the models is small. Making these aspects of developing a model explicit and queryable is an important step for assessing and reusing simulation models more effectively. Exposing this information helps to integrate a new simulation model within a family of existing ones and may lead to the development of more robust and valid simulation models. We hope that our approach becomes part of a standardization effort and that modelers adopt the benefits of provenance when considering or creating simulation models.

Wnt signaling related transcripts and their relationship to energy metabolism in C2C12 myoblasts under temperature stress.

Temperature stress is one of the main environmental stressors affecting the welfare, health and productivity of livestock. Temperature changes can modify cell membrane components, disrupting the crosstalk between the cell and its surroundings by affecting signaling pathways including Wnt signaling pathway, which subsequently disrupts cell energy metabolism. The present study aims to understand the effect of temperature stress on the expression of genes involved in Wnt signaling pathways, and their interaction with energy metabolism in C2C12 myoblasts cells. The C2C12 cells were exposed to cold stress (35 °C), mild heat stress (39 °C) and severe heat stress (41 °C), whereas 37 °C was used as control temperature. Transcript levels of important genes involved in Wnt signaling including Axin2, Tnks2, Sfrp1, Dkk1, Dact1, Cby1, Wnt5a, Wnt7a, Wnt11, Porcn, Ror2, Daam1, and Ppp3ca were significantly altered under severe heat stress (41 °C), whereas eight Wnt signaling-related transcripts (Daam1, Ppp3ca, Fzd7, Wnt5a, Porcn, Tnks2, Lrp6, and Aes) were significantly altered under cold stress (35 °C) compared to control. Under heat stress transcripts of the Wnt/ß-catenin inhibitors (Sfrp1, Dkk1, and Cby1) and negative regulators (Dact1 and Axin2) are activated. A positive correlation between oxidative phosphorylation and Wnt-related transcripts was found under high temperatures. Transcripts of the cell membrane receptors, including Lrp6 and Fzd7, and the members of Wnt/Ca+2 signaling pathway, including Ppp3ca and Porcn were downregulated under cold stress. Many Wnt signaling-related transcripts were positively correlated with glycolysis under cold stress. These findings indicate a cross-talk between Wnt signaling and energy metabolism under thermal stress.

ROS Dependent Wnt/ß-Catenin Pathway and Its Regulation on Defined Micro-Pillars-A Combined In Vitro and In Silico Study.

The physico-chemical surface design of implants influences the surrounding cells. Osteoblasts on sharp-edged micro-topographies revealed an impaired cell phenotype, function and Ca2+ mobilization. The influence of edges and ridges on the Wnt/ß-catenin pathway in combination with the cells' stress response has not been clear. Therefore, MG-63 osteoblasts were studied on defined titanium-coated micro-pillars (5 × 5 × 5 µm) in vitro and in silico. MG-63s on micro-pillars indicated an activated state of the Wnt/ß-catenin pathway. The ß-catenin protein accumulated in the cytosol and translocated into the nucleus. Gene profiling indicated an antagonism mechanism of the transcriptional activity of ß-catenin due to an increased expression of inhibitors like ICAT (inhibitor of ß-catenin and transcription factor-4). Cells on pillars produced a significant reactive oxygen species (ROS) amount after 1 and 24 h. In silico analyses provided a detailed view on how transcriptional activity of Wnt signaling is coordinated in response to the oxidative stress induced by the micro-topography. Based on a coordinated expression of regulatory elements of the Wnt/ß-catenin pathway, MG-63s are able to cope with an increased accumulation of ß-catenin on micro-pillars and suppress an unintended target gene expression. Further, ß-catenin may be diverted into other signaling pathways to support defense mechanisms against ROS.

Exploring the mechanistic and temporal regulation of LRP6 endocytosis in canonical WNT signaling.

Endocytosis plays a pivotal regulatory role in canonical WNT signaling. Internalization of the low-density lipoprotein receptor-related protein 6 (LRP6) receptor complex can either promote or attenuate canonical WNT signaling, depending on the employed internalization pathway. Detailed analysis of the mechanism of LRP6 internalization and its temporal regulation is crucial for understanding the different cellular responses to WNT stimulation under varying conditions and in various cell types. Here, we elucidate the mechanisms involved in the internalization of LRP6 and re-evaluate existing, partly contradicting, theories on the regulation of LRP6 receptor internalization. We utilize a computational approach that aims at finding a set of mechanisms that accounts for the temporal dynamics of LRP6 receptor internalization upon WNT stimulation. Starting with a simple simulation model, we successively extend and probe the model's behavior based on quantitative measurements. The final model confirms that LRP6 internalization is clathrin independent in vertebrates, is not restricted to microdomains, and that signalosome formation delays LRP6 internalization within the microdomains. These findings partly revise the current understanding of LRP6 internalization in vertebrates.

Deep sequencing of small non-coding RNA highlights brain-specific expression patterns and RNA cleavage.

With the advance of high-throughput sequencing technology numerous new regulatory small RNAs have been identified, that broaden the variety of processing mechanisms and functions of non-coding RNA. Here we explore small non-coding RNA (sncRNA) expression in central parts of the physiological stress and anxiety response system. Therefore, we characterize the sncRNA profile of tissue samples from Amygdala, Hippocampus, Hypothalamus and Adrenal Gland, obtained from 20 pigs. Our analysis reveals that all tissues but Amygdala and Hippocampus possess distinct, tissue-specific expression pattern of miRNA that are associated with Hypoxia, stress responses as well as memory and fear conditioning. In particular, we observe marked differences in the expression profile of limbic tissues compared to those associated to the HPA/stress axis, with a surprisingly high aggregation of 3´-tRNA halves in Amygdala and Hippocampus. Since regulation of sncRNA and RNA cleavage plays a pivotal role in the central nervous system, our work provides seminal insights in the role/involvement of sncRNA in the transcriptional and post-transcriptional regulation of negative emotion, stress and coping behaviour in pigs, and mammals in general.

Haplotypes of coping behavior associated QTL regions reveal distinct transcript profiles in amygdala and hippocampus.

Stress response and coping behavior in pigs are largely shaped by hypothalamic-pituitary-adrenal axis and sympatho-adrenomedullary system action. However, the dynamic interaction between amygdala and hippocampus crucially modulates the behavioral response towards significant emotional events. While this functional relationship is well documented, the molecular underpinnings still remain insufficiently understood. Our study used transcriptome profiling of porcine amygdala and hippocampus to identify molecular pathways that are differentially activated depending on the haplotype of a significantly coping behavior-associated region on pig chromosome 12 (SSC12). The pigs were classified into two groups based on the haplotype information of this QTL-region discovered in our previous genome-wide association study. Ten each of high- (HR) and low- (LR) reactive pigs (n = 20) were selected for differential gene expression analysis and weighted gene co-expression analysis with subsequent pathway analysis. Differentially expressed genes identified in the amygdala include SELL, CXCR7 and NTS, while TRAF3, PTGS2 and CFI were detected in the hippocampus indicating a role of neuroinflammation and immunological processes. Pathway analysis revealed IL-8 signaling, NF-κB signaling, glutamate and GABA metabolism, glucocorticoid receptor signaling and chemokine signaling in the amygdala and ephrin receptor signaling, as well as NF-κB signaling in the hippocampus. We discovered candidate genes in regions detected by genome-wide association study including ARRB2, ADRBK2, THRB, NEK7 and ACVR2B, which relate to dopaminergic and other monoaminergic neurotransmitter systems, neuroimmunomodulation, neuroinflammation and GABA-ergic neurotransmission. These findings provide insights into the molecular underpinning of divergent coping behavior and associated haplotypes in limbic forebrain system in pig.

Genetic architecture and regulatory impact on hepatic microRNA expression linked to immune and metabolic traits.

Regulation of microRNA (miRNA) expression contributes to a wide range of target gene expression and phenotypes. The miRNA expression in the liver, the central metabolic organ, was examined in 209 pigs, and integrated with haematological and clinical biomarkers of metabolic and overall health, mRNA-target expression levels and single-nucleotide polymorphism (SNP) genotypes. The expression levels of 426 miRNA species correlated with plasma haematological or biochemical traits (r² = |0.19-0.45|, false discovery rate < 5%). Pairs of these miRNAs and their predicted target mRNAs showing expressing levels associated with the identical traits were examined to understand how immune and metabolic traits are affected by miRNA-mediated regulatory networks derived by mapping miRNA abundance as an expression quantitative trait. In total, 221 miRNA-expression-QTL correspond to 164 SNPs and 108 miRNAs, including miR-34a, miR-30e, miR-148-3p, miR-204, miR-181-5p, miR-143-5p and let-7 g that also correlate with the biomarkers. Sixty-one SNPs were simultaneously associated with 29 miRNA and 41 mRNA species. The expression levels of 13 out of 29 miRNA were correlated with one of the biochemical or haematological traits. For example, the expression levels of miR-34a were correlated with serum phosphorus and cholesterin levels; miR-204, miR-15a and miR-16b were correlated with triglyceride. For haematological traits, the expression levels of miR-652 and miR-204 were correlated with the mean corpuscular haemoglobin concentration, and the expression of miR-143 was correlated with plateletcrit. Pleiotropic association analyses revealed genetic links between mRNA and miRNA on SSC6 for miR-34a, SSC9 for miR-708 and SSC14 for miR-652. Our analysis of miRNA and mRNA transcript profiles, their correlation with clinically important plasma parameters of hepatic functions as well as information on their genetic regulation provide novel regulatory networks and potential new biomarkers for immune and metabolic traits.

Genetically regulated hepatic transcripts and pathways orchestrate haematological, biochemical and body composition traits.

The liver is the central metabolic organ and exhibits fundamental functions in haematological traits. Hepatic expression, haematological, plasma biochemical, and body composition traits were assessed in a porcine model (n = 297) to establish tissue-specific genetic variations that influence the function of immune-metabolism-correlated expression networks. At FDR (false discovery rate) <1%, more than 3,600 transcripts were jointly correlated (r = |0.22-0.48|) with the traits. Functional enrichment analysis demonstrated common links of metabolic and immune traits. To understand how immune and metabolic traits are affected via genetic regulation of gene expression, eQTLs were assessed. 20517 significant (FDR < 5%) eQTLs for 1401 transcripts were identified, among which 443 transcripts were associated with at least one of the examined traits and had cis-eQTL (such as ACO1 (6.52 × 10-7) and SOD1 (6.41 × 10-30). The present study establishes a comprehensive view of hepatic gene activity which links together metabolic and immune traits in a porcine model for medical research.

Spatio-temporal model of endogenous ROS and raft-dependent WNT/beta-catenin signaling driving cell fate commitment in human neural progenitor cells.

Canonical WNT/ß-catenin signaling is a central pathway in embryonic development, but it is also connected to a number of cancers and developmental disorders. Here we apply a combined in-vitro and in-silico approach to investigate the spatio-temporal regulation of WNT/ß-catenin signaling during the early neural differentiation process of human neural progenitors cells (hNPCs), which form a new prospect for replacement therapies in the context of neurodegenerative diseases. Experimental measurements indicate a second signal mechanism, in addition to canonical WNT signaling, being involved in the regulation of nuclear ß-catenin levels during the cell fate commitment phase of neural differentiation. We find that the biphasic activation of ß-catenin signaling observed experimentally can only be explained through a model that combines Reactive Oxygen Species (ROS) and raft dependent WNT/ß-catenin signaling. Accordingly after initiation of differentiation endogenous ROS activates DVL in a redox-dependent manner leading to a transient activation of down-stream ß-catenin signaling, followed by continuous auto/paracrine WNT signaling, which crucially depends on lipid rafts. Our simulation studies further illustrate the elaborate spatio-temporal regulation of DVL, which, depending on its concentration and localization, may either act as direct inducer of the transient ROS/ß-catenin signal or as amplifier during continuous auto-/parcrine WNT/ß-catenin signaling. In addition we provide the first stochastic computational model of WNT/ß-catenin signaling that combines membrane-related and intracellular processes, including lipid rafts/receptor dynamics as well as WNT- and ROS-dependent ß-catenin activation. The model's predictive ability is demonstrated under a wide range of varying conditions for in-vitro and in-silico reference data sets. Our in-silico approach is realized in a multi-level rule-based language, that facilitates the extension and modification of the model. Thus, our results provide both new insights and means to further our understanding of canonical WNT/ß-catenin signaling and the role of ROS as intracellular signaling mediator.

Adaptive spectral clustering with application to tripeptide conformation analysis.

A decomposition of a molecular conformational space into sets or functions (states) allows for a reduced description of the dynamical behavior in terms of transition probabilities between these states. Spectral clustering of the corresponding transition probability matrix can then reveal metastabilities. The more states are used for the decomposition, the smaller the risk to cover multiple conformations with one state, which would make these conformations indistinguishable. However, since the computational complexity of the clustering algorithm increases quadratically with the number of states, it is desirable to have as few states as possible. To balance these two contradictory goals, we present an algorithm for an adaptive decomposition of the position space starting from a very coarse decomposition. The algorithm is applied to small data classification problems where it was shown to be superior to commonly used algorithms, e.g., k-means. We also applied this algorithm to the conformation analysis of a tripeptide molecule where six-dimensional time series are successfully analyzed.

Studying the role of lipid rafts on protein receptor bindings with cellular automata.

It is widely accepted that lipid rafts promote receptor clustering and thereby facilitate signaling transduction. The role of lipid rafts in inducing and promoting receptor accumulation within the cell membrane has been explored by several computational and experimental studies. However, it remains unclear whether lipid rafts influence the recruitment and binding of proteins from the cytosol as well. To provide an answer to this question a spatial membrane model has been developed based on cellular automata. Our results indicate that lipid rafts indeed influence protein receptor bindings. In particular processes with slow dissociation and binding kinetics are promoted by lipid rafts, whereas fast binding processes are slightly hampered. However, the impact depends on a variety of parameters, such as the size and mobility of the lipid rafts, the induced slow down of receptors within rafts, and also the dissociation and binding kinetics of the cytosolic proteins. Thus, for any individual signaling pathway the influence of lipid rafts on protein binding might be different. To facilitate analyzing this influence given a specific pathway, our approach has been generalized into LiRaM, a modeling and simulation tool for lipid rafts models.