Coarse-graining methods and polymer dynamics

This thesis studies molecular dynamics simulations on two levels of resolution: the detailed level of atomistic simulations, where the motion of explicit atoms in a many-particle system is considered, and the coarse-grained level, where the motion of superatoms composed of up to 10 atoms is modeled. While atomistic models are capable of describing material specific effects on small scales, the time and length scales they can cover are limited due to their computational costs. Polymer systems are typically characterized by effects on a broad range of length and time scales. Therefore it is often impossible to atomistically simulate processes, which determine macroscopic properties in polymer systems. Coarse-grained (CG) simulations extend the range of accessible time and length scales by three to four orders of magnitude. However, no standardized coarse-graining procedure has been established yet. Following the ideas of structure-based coarse-graining, a coarse-grained model for polystyrene is presented. Structure-based methods parameterize CG models to reproduce static properties of atomistic melts such as radial distribution functions between superatoms or other probability distributions for coarse-grained degrees of freedom. Two enhancements of the coarse-graining methodology are suggested. Correlations between local degrees of freedom are implicitly taken into account by additional potentials acting between neighboring superatoms in the polymer chain. This improves the reproduction of local chain conformations and allows the study of different tacticities of polystyrene. It also gives better control of the chain stiffness, which agrees perfectly with the atomistic model, and leads to a reproduction of experimental results for overall chain dimensions, such as the characteristic ratio, for all different tacticities. The second new aspect is the computationally cheap development of nonbonded CG potentials based on the sampling of pairs of oligomers in vacuum. Static properties of polymer melts are obtained as predictions of the CG model in contrast to other structure-based CG models, which are iteratively refined to reproduce reference melt structures. The dynamics of simulations at the two levels of resolution are compared. The time scales of dynamical processes in atomistic and coarse-grained simulations can be connected by a time scaling factor, which depends on several specific system properties as molecular weight, density, temperature, and other components in mixtures. In this thesis the influence of molecular weight in systems of oligomers and the situation in two-component mixtures is studied. For a system of small additives in a melt of long polymer chains the temperature dependence of the additive diffusion is predicted and compared to experiments.

Effects of low-density lipoprotein receptor-related protein 4 and apolipoprotein E isoforms on bone metabolism in vivo

LRP4, member of the LDLR family, is a multifunctional membrane-bound receptor that is expressed in various tissues. The expression of LRP4 by osteoblasts, its novel interaction with Wnt-signaling inhibitors Dkk1 and SOST, and the lower levels of activated beta-catenin in different bone locations described here, adds another player to the long list of established factors that modulate canonical Wnt-signaling in bone. By demonstrating that in addition to Wise, LRP4 is able to interact with two additional important modulators of Wnt- and BMP-signaling, our perspective of the complexity of the integration of BMP and Wnt-signaling pathways on the osteoblast surface has expanded further. Nevertheless the recently described association of both the SOST and LRP4 genes with BMD in humans, together with our findings suggest that LRP4 plays a physiologically important role in the skeletal development and bone metabolism not only in rodents, but in humans as well. The efficiency with which LRP4 binds both SOST and Dkk1, presumably at the osteoblastic surface, LRP4 may act as a sink and competes with LRP5/6 for the binding of these Wnt antagonists, which then are no longer available for suppression of the signal through the LRP5/6 axis. rnApoE, a 299 amino acid glycoprotein, is a crucial regulator in the uptake of triglyceride, phospholipids, cholesteryl esters, and cholesterol into cells. ApoE has been linked to osteoporosis, and such a role is further strengthened by the present of a high bone mass phenotype in ApoE null mice. Until recently, the effects of respective ApoE isoforms E2, E3, and E4, and their impact on bone metabolism, have been unclear. Here we report that respective human ApoE knockin mice display diverse effects on bone metabolism. ApoE2 mice show decreased trabecular bone volume per total volume in femoral bone and lumbar spine in comparison to ApoE3 and E4 animals. In this context, urinary bone resorption marker DPD is increased in these animals, which is accompanied by a low ratio of osteoclastogenesis markers OPG/RANKL. Interestingly, serum bone formation markers ALP and OCN are diminished in ApoE4 mice. In contrast to this finding, ApoE2 mice show the lowest bone formation of all groups in vivo. These findings cannot be explained by the low receptor-affinity of ApoE2 and subsequent decreased uptake of triglyceride-rich lipoproteins by osteoblasts, resulting in elevated levels of undercarboxylated osteocalcin. Thus, other crucial pathways relevant for bone metabolism, e. g. Wnt/beta-catenin-signaling pathways, must be, compared to the ApoE3/4 isoforms, more affected by the ApoE2 isoform.

Funktionelle, strukturelle und phylogenetische Untersuchungen an Lipoproteinen des Flusskrebses Astacus leptodactylus

Diese Arbeit untersucht zwei Lipoproteine, das discoidale High density-Lipoprotein (dHDL) und das β-Glukan-Bindeprotein (BGBP) aus dem Flusskrebs Astacus leptodactylus in funktioneller, struktureller und phylogenetischer Hinsicht. Die Nukleotid-Sequenz des BGBP konnte nahezu vollständig entschlüsselt werden. Dabei errechnet sich aus der abgeleiteten Aminosäure-Sequenz ein Molekulargewicht von 153 kDa. Das reife BGBP hat nur eine molekulare Masse von 105 kDa. Vermutlich kommt es durch eine Furin-ähnliche Protease zu einer post-translationalen N- und C-terminalen Prozessierung: zwei bisher nicht beschriebene, aber auch in der BGBP-Sequenz von anderen höheren Krebsen vorhandene, typische Furin-Schnittstellen (RAKR, bzw. RARR) wurden anhand von Sequenzvergleichen identifiziert. BGBP hat zwei Funktionen: zum Einen ist es für den Transport und die Aktivierung des proPhenoloxidase-Systems zuständig, zum Anderen für die Versorgung der Organe mit Lipiden, welche vermutlich der Energiegewinnung dienen. Eine 100 kDa große, BGBP-bindende Rezeptor-Fraktion konnte in Hämocyten-Membranen identifiziert werden. Das Vorkommen von dHDL war aus eigenen Befunden bisher ausschließlich in Astacus leptodactylus bekannt, doch konnte in dieser Arbeit ein mit dem dHDL-Antikörper reagierendes Protein erstmalig auch in anderen Arthropoden-Spezies nachgewiesen werden. Die discoidale Form und das Untereinheiten-Muster (240 + 85 kDa) sind typisch für die bei Vertretern ursprünglicher Tiergruppen gefundenen Lipoproteine (z.B. beim Cheliceraten Limulus und beim Polychaeten Nereis). Eventuell handelt es sich bei dHDL also um einen ‚Prototypen’ in der Lipoprotein-Evolution. Obwohl die Sequenz des dHDL auf Nukleotid-Ebene unbekannt ist, wurden die Sequenzen einiger dHDL-Peptide aus massenspektroskopischen Analysen gewonnen. Überraschenderweise befinden sich diese Sequenzen in der Aminosäuresequenz des BGBP. Dabei liegen alle Peptide am N- und/oder am C-Terminus der abgeleiteten BGBP-Aminosäure-Sequenz, und zwar in den Bereichen, die vermutlich durch das erwähnte Furin vom BGBP abgeschnitten werden, im reifen BGBP also gar nicht mehr vorkommen. Deshalb ist zu vermuten, dass BGBP und dHDL ein gemeinsames Vorläuferprotein haben und durch Genduplikation entstanden sind, oder dass es sich beim dHDL- und beim BGBP-Gen um ein und dasselbe Gen handelt. Das Genprodukt wird dann auf unterschiedliche Weise prozessiert und es entstehen die beiden Proteine dHDL und BGBP. Die Funktion von dHDL ist noch nicht eindeutig geklärt, es ließen sich aber dHDL-bindende Rezeptor-Fraktionen mit einer molekularen Masse von 160 kDa in somatischen Geweben (Muskel, Darm, Hepatopankreas, Kiemen und Samenleiter) sowie in Oocyten und Hämocyten nachweisen. Deshalb wird vermutet, dass dHDL als Energielieferant in Stoffwechsel-aktiven Organen und als Speicherprotein in Oocyten dient. Eine endocytotische Aufnahme konnte gezeigt werden.

The influence of individual and colony level variation in behavior on colony performance in ants

In my dissertation I investigated the influence of behavioral variation between and within ant colonies on group performance. In particular, I analyzed how evolution shapes behavior in response to ecological conditions, and whether within-group diversity improves productivity as suggested by theory. Our field and laboratory experiments showed that behavioral diverse groups are more productive. Different aggression levels within colonies were beneficial under competitive field situations, whereas diversity in brood care and exploratory behavior were favored in non-competitive laboratory situations. We then examined whether population density and social parasite presence shape aggression through phenotypic plasticity and/or natural selection. The importance of selection was indicated by the absence of density or parasite effects on aggression in a field manipulation. Indeed, more aggressive colonies fared better under high density and during parasite attack. When analyzing the proximate causes of individual behavioral variation, ovarian development was shown to be linked to division of labor and aggressiveness. Finally, our studies show that differences in the collective behavior can be linked to immune defense and productivity. My dissertation demonstrates that behavioral variation should be studied on multiple scales and when possible combined with physiological analyses to better understand the evolution of animal personalities in social groups.rn

Natural compounds Camptothecin and Triptolide: highly specific enzyme inhibitors and tools to dissect transcriptional functions

With this work I elucidated new and unexpected mechanisms of two strong and highly specific transcription inhibitors: Triptolide and Campthotecin. Triptolide (TPL) is a diterpene epoxide derived from the Chinese plant Trypterigium Wilfoordii Hook F. TPL inhibits the ATPase activity of XPB, a subunit of the general transcription factor TFIIH. In this thesis I found that degradation of Rbp1 (the largest subunit of RNA Polymerase II) caused by TPL treatments, is preceded by an hyperphosphorylation event at serine 5 of the carboxy-terminal domain (CTD) of Rbp1. This event is concomitant with a block of RNA Polymerase II at promoters of active genes. The enzyme responsible for Ser5 hyperphosphorylation event is CDK7. Notably, CDK7 downregulation rescued both Ser5 hyperphosphorylation and Rbp1 degradation triggered by TPL. Camptothecin (CPT), derived from the plant Camptotheca acuminata, specifically inhibits topoisomerase 1 (Top1). We first found that CPT induced antisense transcription at divergent CpG islands promoter. Interestingly, by immunofluorescence experiments, CPT was found to induce a burst of R loop structures (DNA/RNA hybrids) at nucleoli and mitochondria. We then decided to investigate the role of Top1 in R loop homeostasis through a short interfering RNA approach (RNAi). Using DNA/RNA immunoprecipitation techniques coupled to NGS I found that Top1 depletion induces an increase of R loops at a genome-wide level. We found that such increase occurs on the entire gene body. At a subset of loci R loops resulted particularly stressed after Top1 depletion: some of these genes showed the formation of new R loops structures, whereas other loci showed a reduction of R loops. Interestingly we found that new peaks usually appear at tandem or divergent genes in the entire gene body, while losses of R loop peaks seems to be a feature specific of 3’ end regions of convergent genes.

The CFTR frameshift mutation 3905insT and its effect at transcript and protein level

Cystic fibrosis (CF) is one of the most common genetic diseases in the Caucasian population and is characterized by chronic obstructive pulmonary disease, exocrine pancreatic insufficiency, and elevation of sodium and chloride concentrations in the sweat and infertility in men. The disease is caused by mutations in the CF transmembrane conductance regulator (CFTR) gene, which encodes a protein that functions as chloride channel at the apical membrane of different epithelia. Owing to the high genotypic and phenotypic disease heterogeneity, effects and consequences of the majority of the CFTR mutations have not yet been studied. Recently, the frameshift mutation 3905insT was identified as the second most frequent mutation in the Swiss population and found to be associated with a severe phenotype. The frameshift mutation produces a premature termination codon (PTC) in exon 20, and transcripts bearing this PTC are potential targets for degradation through nonsense-mediated mRNA decay (NMD) and/or for exon skipping through nonsense-associated alternative splicing (NAS). Using RT-PCR analysis in lymphocytes and different tissue types from patients carrying the mutation, we showed that the PTC introduced by the mutation does neither elicit a degradation of the mRNA through NMD nor an alternative splicing through NAS. Moreover, immunocytochemical analysis in nasal epithelial cells revealed a significantly reduced amount of CFTR at the apical membrane providing a possible molecular explanation for the more severe phenotype observed in F508del/3905insT compound heterozygotes compared with F508del homozygotes. However, further experiments are needed to elucidate the fate of the 3905insT CFTR in the cell after its biosynthesis.

Hydrodynamic correlation functions of a driven granular fluid in steady state

We study a homogeneously driven granular fluid of hard spheres at intermediate volume fractions and focus on time-delayed correlation functions in the stationary state. Inelastic collisions are modeled by incomplete normal restitution, allowing for efficient simulations with an event-driven algorithm. The incoherent scattering function Fincoh(q,t ) is seen to follow time-density superposition with a relaxation time that increases significantly as the volume fraction increases. The statistics of particle displacements is approximately Gaussian. For the coherent scattering function S(q,ω), we compare our results to the predictions of generalized fluctuating hydrodynamics, which takes into account that temperature fluctuations decay either diffusively or with a finite relaxation rate, depending on wave number and inelasticity. For sufficiently small wave number q we observe sound waves in the coherent scattering function S(q,ω) and the longitudinal current correlation function Cl(q,ω). We determine the speed of sound and the transport coefficients and compare them to the results of kinetic theory.

Prime and Semiprime Inner Functions

This paper studies the structure of inner functions under the operation of composition, and in particular the notions or primeness and semiprimeness. Results proved include the density of prime finite Blaschke products in the set of finite Blaschke products, the semiprimeness of finite products of thin Blaschke products and their approximability by prime Blaschke products. An example of a nonsemiprime Blaschke product that is a Frostman Blaschke product is also provided.

Prime and semiprime inner functions

Abstract This paper studies the structure of inner functions under the operation of composition, and in particular the notions or primeness and semiprimeness. Results proved include the density of prime finite Blaschke products in the set of finite Blaschke products, the semiprimeness of finite products of thin Blaschke products and their approximability by prime Blaschke products. An example of a nonsemiprime Blaschke product that is a Frostman Blaschke product is also provided.

Stocking of captive-bred fish can cause long-term population decline and gene pool replacement: predictions from a population dynamics model incorporating density-dependent mortality

Releasing captive-bred fish into natural environments (stocking) is common in fisheries worldwide. Although stocking is believed to have a positive effect on fish abundance over the short term, little is known about the long-term consequences of recurrent stocking and its influence on natural populations. In fact, there are growing concerns that genetically maladapted captive-bred fish can eventually reduce the abundance of natural population. In this study, we develop a simple model to quantitatively investigate the condition under which recurrent stocking has long-term effects on the natural population. Using a population dynamics model that takes into account a density-dependent recruitment, a gene responsible for the fitness difference between wild and captive-bred fish, and hybridization between them, we show that there is little or no contribution of recurrent stocking to the stock enhancement without a replacement of the wild gene pool by the captive-bred gene pool. The model further predicted that stocking of an intermediate level causes a reduction, rather than enhancement, of population size over the long term. The population decline due to stocking was attributed to the fitness disadvantage of captive-bred fish and strong overcompensation at recruitment stage. These results suggest that it would be difficult to simultaneously attain population size recovery and conservation of the local gene pool when captive-bred fish have fitness disadvantage in the wild, although caution is needed when applying the predictions from the simplified model to a specific species or population.

Reduced Cerebral Blood Flow Within the Default-Mode Network and Within Total Gray Matter in Major Depression

The default-mode network (DMN) was shown to have aberrant blood oxygenation-level-dependent (BOLD) activity in major depressive disorder (MDD). While BOLD is a relative measure of neural activity, cerebral blood flow (CBF) is an absolute measure. Resting-state CBF alterations have been reported in MDD. However, the association of baseline CBF and CBF fluctuations is unclear in MDD. Therefore, the aim was to investigate the CBF within the DMN in MDD, applying a strictly data-driven approach. In 22 MDD patients and 22 matched healthy controls, CBF was acquired using arterial spin labeling (ASL) at rest. A concatenated independent component analysis was performed to identify the DMN within the ASL data. The perfusion of the DMN and its nodes was quantified and compared between groups. The DMN was identified in both groups with high spatial similarity. Absolute CBF values within the DMN were reduced in MDD patients (p<0.001). However, after controlling for whole-brain gray matter CBF and age, the group difference vanished. In patients, depression severity was correlated with reduced perfusion in the DMN in the posterior cingulate cortex and the right inferior parietal lobe. Hypoperfusion within the DMN in MDD is not specific to the DMN. Still, depression severity was linked to DMN node perfusion, supporting a role of the DMN in depression pathobiology. The finding has implications for the interpretation of BOLD functional magnetic resonance imaging data in MDD.

Expression of antagonists of WNT and BMP signaling after non-rigid fixation of osteotomies

Delayed fracture healing and non-unions represent rare but severe complications in orthopedic surgery. Further knowledge on the mechanisms of the bone repair process and of the development of a pseudoarthrosis is essential to predict and prevent impaired healing of fractures. The present study aimed at elucidating differences in gene expression during the repair of rigidly and non-rigidly fixed osteotomies. For this purpose, the MouseFix™ and the FlexiPlate™ systems (AO Development Institute, Davos, CH), allowing the creation of well defined osteotomies in mouse femora, were employed. A time course following the healing process of the osteotomy was performed and bones and periimplant tissues were analyzed by high-resolution X-ray, MicroCT and by histology. For the assessment of gene expression, Low Density Arrays (LDA) were done. In animals with rigid fixation, X-ray and MicroCT revealed healing of the osteotomy within 3 weeks. Using the FlexiPlate™ system, the osteotomy was still visible by X-ray after 3 weeks and a stabilizing cartilaginous callus was formed. After 4.5 weeks, the callus was remodeled and the osteotomy was, on a histological level, healed. Gene expression studies revealed levels of transcripts encoding proteins associated with inflammatory processes not to be altered in tissues from bones with rigid and non-rigid fixation, respectively. Levels of transcripts encoding proteins of the extracellular matrix and essential for bone cell functions were not increased in the rigidly fixed group when compared to controls without osteotomy. In the FlexiPlate™ group, levels of transcripts encoding the same set of genes were significantly increased 3 weeks after surgery. Expression of transcripts encoding BMPs and BMP antagonists was increased after 3 weeks in repair tissues from bones fixed with FlexiPlate™, as were inhibitors of the WNT signaling pathways. Little changes only were detected in transcript levels of tissues from rigidly fixed bones. The data of the present study suggest that rigid fixation enables accelerated healing of an experimental osteotomy as compared to non-rigid fixation. The changes in the healing process after non-rigid fixation are accompanied by an increase in the levels of transcripts encoding inhibitors of osteogenic pathways and, probably as a consequence, by temporal changes in bone matrix synthesis.

The future of osteoporosis treatment - a research update

Osteoporosis is characterised by a progressive loss of bone mass and microarchitecture which leads to increased fracture risk. Some of the drugs available to date have shown reductions in vertebral and non-vertebral fracture risk. However, in the ageing population of industrialised countries, still more fractures happen today than are avoided, which highlights the large medical need for new treatment options, models, and strategies. Recent insights into bone biology, have led to a better understanding of bone cell functions and crosstalk between osteoblasts, osteoclasts, and osteocytes at the molecular level. In the future, the armamentarium against osteoporotic fractures will likely be enriched by (1.) new bone anabolic substances such as antibodies directed against the endogenous inhibitors of bone formation sclerostin and dickkopf-1, PTH and PTHrp analogues, and possibly calcilytics; (2.) new inhibitors of bone resorption such as cathepsin K inhibitors which may suppress osteoclast function without impairing osteoclast viability and thus maintain bone formation by preserving the osteoclast-osteoblast crosstalk, and denosumab, an already widely available antibody against RANKL which inhibits osteoclast formation, function, and survival; and (3.) new therapeutic strategies based on an extended understanding of the pathophysiology of osteoporosis which may include sequential therapies with two or more bone active substances aimed at optimising the management of bone capital acquired during adolescence and maintained during adulthood in terms of both quantity and quality. Finally, one of the future challenges will be to identify those patients and patient populations expected to benefit the most from a given drug therapy or regimen. The WHO fracture risk assessment tool FRAX® and improved access to bone mineral density measurements by DXA will play a key role in this regard.

Sequential Quantitative Trait Locus Mapping in Experimental Crosses

The etiology of complex diseases is heterogeneous. The presence of risk alleles in one or more genetic loci affects the function of a variety of intermediate biological pathways, resulting in the overt expression of disease. Hence, there is an increasing focus on identifying the genetic basis of disease by sytematically studying phenotypic traits pertaining to the underlying biological functions. In this paper we focus on identifying genetic loci linked to quantitative phenotypic traits in experimental crosses. Such genetic mapping methods often use a one stage design by genotyping all the markers of interest on the available subjects. A genome scan based on single locus or multi-locus models is used to identify the putative loci. Since the number of quantitative trait loci (QTLs) is very likely to be small relative to the number of markers genotyped, a one-stage selective genotyping approach is commonly used to reduce the genotyping burden, whereby markers are genotyped solely on individuals with extreme trait values. This approach is powerful in the presence of a single quantitative trait locus (QTL) but may result in substantial loss of information in the presence of multiple QTLs. Here we investigate the efficiency of sequential two stage designs to identify QTLs in experimental populations. Our investigations for backcross and F2 crosses suggest that genotyping all the markers on 60% of the subjects in Stage 1 and genotyping the chromosomes significant at 20% level using additional subjects in Stage 2 and testing using all the subjects provides an efficient approach to identify the QTLs and utilizes only 70% of the genotyping burden relative to a one stage design, regardless of the heritability and genotyping density. Complex traits are a consequence of multiple QTLs conferring main effects as well as epistatic interactions. We propose a two-stage analytic approach where a single-locus genome scan is conducted in Stage 1 to identify promising chromosomes, and interactions are examined using the loci on these chromosomes in Stage 2. We examine settings under which the two-stage analytic approach provides sufficient power to detect the putative QTLs.

Imaging emotional brain functions: conceptual and methodological issues

This article reviews the psychophysiological and brain imaging literature on emotional brain function from a methodological point of view. The difficulties in defining, operationalising and measuring emotional activation and, in particular, aversive learning will be considered. Emotion is a response of the organism during an episode of major significance and involves physiological activation, motivational, perceptual, evaluative and learning processes, motor expression, action tendencies and monitoring/subjective feelings. Despite the advances in assessing the physiological correlates of emotional perception and learning processes, a critical appraisal shows that functional neuroimaging approaches encounter methodological difficulties regarding measurement precision (e.g., response scaling and reproducibility) and validity (e.g., response specificity, generalisation to other paradigms, subjects or settings). Since emotional processes are not only the result of localised but also of widely distributed activation, a more representative model of assessment is needed that systematically relates the hierarchy of high- and low-level emotion constructs with the corresponding patterns of activity and functional connectivity of the brain.