Simulation studies of correlation functions and relaxation in polymeric systems

We investigate the statics and dynamics of a glassy,non-entangled, short bead-spring polymer melt with moleculardynamics simulations. Temperature ranges from slightlyabove the mode-coupling critical temperature to the liquidregime where features of a glassy liquid are absent. Ouraim is to work out the polymer specific effects on therelaxation and particle correlation. We find the intra-chain static structure unaffected bytemperature, it depends only on the distance of monomersalong the backbone. In contrast, the distinct inter-chainstructure shows pronounced site-dependence effects at thelength-scales of the chain and the nearest neighbordistance. There, we also find the strongest temperaturedependence which drives the glass transition. Both the siteaveraged coupling of the monomer and center of mass (CM) andthe CM-CM coupling are weak and presumably not responsiblefor a peak in the coherent relaxation time at the chain'slength scale. Chains rather emerge as soft, easilyinterpenetrating objects. Three particle correlations arewell reproduced by the convolution approximation with theexception of model dependent deviations. In the spatially heterogeneous dynamics of our system weidentify highly mobile monomers which tend to follow eachother in one-dimensional paths forming ``strings''. Thesestrings have an exponential length distribution and aregenerally short compared to the chain length. Thus, arelaxation mechanism in which neighboring mobile monomersmove along the backbone of the chain seems unlikely.However, the correlation of bonded neighbors is enhanced. When liquids are confined between two surfaces in relativesliding motion kinetic friction is observed. We study ageneric model setup by molecular dynamics simulations for awide range of sliding speeds, temperatures, loads, andlubricant coverings for simple and molecular fluids. Instabilities in the particle trajectories are identified asthe origin of kinetic friction. They lead to high particlevelocities of fluid atoms which are gradually dissipatedresulting in a friction force. In commensurate systemsfluid atoms follow continuous trajectories for sub-monolayercoverings and consequently, friction vanishes at low slidingspeeds. For incommensurate systems the velocity probabilitydistribution exhibits approximately exponential tails. Weconnect this velocity distribution to the kinetic frictionforce which reaches a constant value at low sliding speeds. This approach agrees well with the friction obtaineddirectly from simulations and explains Amontons' law on themicroscopic level. Molecular bonds in commensurate systemslead to incommensurate behavior, but do not change thequalitative behavior of incommensurate systems. However,crossed chains form stable load bearing asperities whichstrongly increase friction.

Modeling and simulation of a synthetic genetic circuit that implements a multicelled behavior in a growing microcolony of E. coli

Synthetic Biology is a relatively new discipline, born at the beginning of the New Millennium, that brings the typical engineering approach (abstraction, modularity and standardization) to biotechnology. These principles aim to tame the extreme complexity of the various components and aid the construction of artificial biological systems with specific functions, usually by means of synthetic genetic circuits implemented in bacteria or simple eukaryotes like yeast. The cell becomes a programmable machine and its low-level programming language is made of strings of DNA. This work was performed in collaboration with researchers of the Department of Electrical Engineering of the University of Washington in Seattle and also with a student of the Corso di Laurea Magistrale in Ingegneria Biomedica at the University of Bologna: Marilisa Cortesi. During the collaboration I contributed to a Synthetic Biology project already started in the Klavins Laboratory. In particular, I modeled and subsequently simulated a synthetic genetic circuit that was ideated for the implementation of a multicelled behavior in a growing bacterial microcolony. In the first chapter the foundations of molecular biology are introduced: structure of the nucleic acids, transcription, translation and methods to regulate gene expression. An introduction to Synthetic Biology completes the section. In the second chapter is described the synthetic genetic circuit that was conceived to make spontaneously emerge, from an isogenic microcolony of bacteria, two different groups of cells, termed leaders and followers. The circuit exploits the intrinsic stochasticity of gene expression and intercellular communication via small molecules to break the symmetry in the phenotype of the microcolony. The four modules of the circuit (coin flipper, sender, receiver and follower) and their interactions are then illustrated. In the third chapter is derived the mathematical representation of the various components of the circuit and the several simplifying assumptions are made explicit. Transcription and translation are modeled as a single step and gene expression is function of the intracellular concentration of the various transcription factors that act on the different promoters of the circuit. A list of the various parameters and a justification for their value closes the chapter. In the fourth chapter are described the main characteristics of the gro simulation environment, developed by the Self Organizing Systems Laboratory of the University of Washington. Then, a sensitivity analysis performed to pinpoint the desirable characteristics of the various genetic components is detailed. The sensitivity analysis makes use of a cost function that is based on the fraction of cells in each one of the different possible states at the end of the simulation and the wanted outcome. Thanks to a particular kind of scatter plot, the parameters are ranked. Starting from an initial condition in which all the parameters assume their nominal value, the ranking suggest which parameter to tune in order to reach the goal. Obtaining a microcolony in which almost all the cells are in the follower state and only a few in the leader state seems to be the most difficult task. A small number of leader cells struggle to produce enough signal to turn the rest of the microcolony in the follower state. It is possible to obtain a microcolony in which the majority of cells are followers by increasing as much as possible the production of signal. Reaching the goal of a microcolony that is split in half between leaders and followers is comparatively easy. The best strategy seems to be increasing slightly the production of the enzyme. To end up with a majority of leaders, instead, it is advisable to increase the basal expression of the coin flipper module. At the end of the chapter, a possible future application of the leader election circuit, the spontaneous formation of spatial patterns in a microcolony, is modeled with the finite state machine formalism. The gro simulations provide insights into the genetic components that are needed to implement the behavior. In particular, since both the examples of pattern formation rely on a local version of Leader Election, a short-range communication system is essential. Moreover, new synthetic components that allow to reliably downregulate the growth rate in specific cells without side effects need to be developed. In the appendix are listed the gro code utilized to simulate the model of the circuit, a script in the Python programming language that was used to split the simulations on a Linux cluster and the Matlab code developed to analyze the data.

Responses to Central Oxotremorine and Scopolamine Support the Cholinergic Control of Male Mating Behavior in Hamsters

The responses of hamsters to intracranial injections of the cholinergic agonist oxotremorine (OXO) implicate cholinergic mechanisms in the medial preoptic area (MPOA) in the control of male mating behavior. To extend these observations, we ran three studies of responses to cholinergic drugs delivered singly or in combination to the vicinity of the MPOA. The first tested responses to OXO, confirming its ability to reduce the postejaculatory interval. The second complemented the first by examining responses to MPOA microinjections of the cholinergic antagonist scopolamine (SCO). These caused several changes revolving around intromission. These included increases in intromission frequency and ejaculation latency. They also included a change in the patterning of intromissions, marked by continuous strings without the usual separation by dismounts. The final study resembled the others in examining the effects of MPOA injections of OXO and SCO but focused on the ability of each drug to antagonize responses to the other. Most of the responses to OXO and SCO individually replicated earlier findings, though the measures examined here also permitted the description of effects on some noncopulatory sexual behaviors, specifically the male's inspection of the female. However, the most interesting results may be those suggesting asymmetry in the responses to the addition of the second drug: Whereas responses to OXO tended to be antagonized by SCO, OXO was less effective at counteracting responses to SCO. Though the explanation of this asymmetry is not completely clear, it is consistent with previous suggestions of differences in the affinities of these drugs for subtypes of muscarinic receptors. Therefore, it suggests that the cholinergic synapses and circuits controlling distinct elements of male behavior could differ in their dependence on these receptors. Copyright 2013 Elsevier Inc. All rights reserved.

Analysis of gene order data supports vertical inheritance of the leukotoxin operon and genome rearrangements in the 5' flanking region in genus Mannheimia

BACKGROUND: The Mannheimia subclades belong to the same bacterial genus, but have taken divergent paths toward their distinct lifestyles. For example, M. haemolytica + M. glucosida are potential pathogens of the respiratory tract in the mammalian suborder Ruminantia, whereas M. ruminalis, the supposed sister group, lives as a commensal in the ovine rumen. We have tested the hypothesis that vertical inheritance of the leukotoxin (lktCABD) operon has occurred from the last common ancestor of genus Mannheimia to any ancestor of the diverging subclades by exploring gene order data. RESULTS: We examined the gene order in the 5' flanking region of the leukotoxin operon and found that the 5' flanking gene strings, hslVU-lapB-artJ-lktC and xylAB-lktC, are peculiar to M. haemolytica + M. glucosida and M. granulomatis, respectively, whereas the gene string hslVU-lapB-lktC is present in M. ruminalis, the supposed sister group of M. haemolytica + M. glucosida, and in the most ancient subclade M. varigena. In M. granulomatis, we found remnants of the gene string hslVU-lapB-lktC in the xylB-lktC intergenic region. CONCLUSION: These observations indicate that the gene string hslVU-lapB-lktC is more ancient than the hslVU-lapB-artJ-lktC and xylAB-lktC gene strings. The presence of (remnants of) the ancient gene string hslVU-lapB-lktC among any subclades within genus Mannheimia supports that it has been vertically inherited from the last common ancestor of genus Mannheimia to any ancestor of the diverging subclades, thus reaffirming the hypothesis of vertical inheritance of the leukotoxin operon. The presence of individual 5' flanking regions in M. haemolytica + M. glucosida and M. granulomatis reflects later genome rearrangements within each subclade. The evolution of the novel 5' flanking region in M. haemolytica + M. glucosida resulted in transcriptional coupling between the divergently arranged artJ and lkt promoters. We propose that the chimeric promoter have led to high level expression of the leukotoxin operon which could explain the increased potential of certain M. haemolytica + M. glucosida strains to cause a particular type of infection.

ADVENTURES IN FRIEDMANN COSMOLOGIES---INTERACTION OF POSITIVE ENERGY DENSITIES WITH NEGATIVE ENERGY DENSITIES AND CURVATURE OF THE UNIVERSE

In this report we will investigate the effect of negative energy density in a classic Friedmann cosmology. Although never measured and possibly unphysical, the evolution of a Universe containing a significant cosmological abundance of any of a number of hypothetical stable negative energy components is explored. These negative energy (Ω < 0) forms include negative phantom energy (w<-1), negative cosmological constant (w=-1), negative domain walls (w=-2/3), negative cosmic strings (w= -1/3), negative mass (w=0), negative radiation (w=1/3), and negative ultra-light (w > 1/3). Assuming that such universe components generate pressures as perfect fluids, the attractive or repulsive nature of each negative energy component is reviewed. The Friedmann equations can only be balanced when negative energies are coupled to a greater magnitude of positive energy or positive curvature, and minimal cases of both of these are reviewed. The future and fate of such universes in terms of curvature, temperature, acceleration, and energy density are reviewed including endings categorized as a Big Crunch, Big Void, or Big Rip and further qualified as "Warped", "Curved", or "Flat", "Hot" versus "Cold", "Accelerating" versus" Decelerating" versus "Coasting". A universe that ends by contracting to zero energy density is termed a Big Poof. Which contracting universes ``bounce" in expansion and which expanding universes ``turnover" into contraction are also reviewed. The name by which the ending of the Universe is mentioned is our own nomenclature.

Monitoring the Time Course of Perception without Awareness: A Comparison of Mirror Masked Words and Nonwords

Mirror masked words are embedded into a context that makes them appear as senseless patterns or as strings of unfamiliar letters. Thus, mirror masked words can be shown for several hundreds of milliseconds without being recognised as words. We sought to further investigate effects of nonsconscious reading by monitoring event-related brain potentials (ERPs) while participants observed mirror masked letter strings. ERPs were recorded while participants observed mirror masked words and nonwords. Data of 15 participants was segmented into periods of quasi-stable field topography (microstates). Microstates for masked words and nonwords were compared using randomization tests, statistical parametric scalp maps and Low Resolution Electromagnetic Tomography (LORETA). ERPs to masked words and nonwords showed significant topographic differences between 136 and 256 ms, indicating that stimuli were nonconsciously discriminated. A LORETA model localised sources of activation discriminating between masked words and nonwords in left operculum, the right superior parietal lobe and right superior temporal gyrus indicating higher current density for nonwords than for words in these areas. ERPs of mirror masked stimuli can indicate unconscious discrimination even in cases where behavioural priming is unreliable. This approach might be useful for investigating differences in early, nonconscious stages of word perception.

Dirichlet-Neumann inverse spectral problem for a star graph of Stieltjes strings

We solve two inverse spectral problems for star graphs of Stieltjes strings with Dirichlet and Neumann boundary conditions, respectively, at a selected vertex called root. The root is either the central vertex or, in the more challenging problem, a pendant vertex of the star graph. At all other pendant vertices Dirichlet conditions are imposed; at the central vertex, at which a mass may be placed, continuity and Kirchhoff conditions are assumed. We derive conditions on two sets of real numbers to be the spectra of the above Dirichlet and Neumann problems. Our solution for the inverse problems is constructive: we establish algorithms to recover the mass distribution on the star graph (i.e. the point masses and lengths of subintervals between them) from these two spectra and from the lengths of the separate strings. If the root is a pendant vertex, the two spectra uniquely determine the parameters on the main string (i.e. the string incident to the root) if the length of the main string is known. The mass distribution on the other edges need not be unique; the reason for this is the non-uniqueness caused by the non-strict interlacing of the given data in the case when the root is the central vertex. Finally, we relate of our results to tree-patterned matrix inverse problems.

Gene expression in sea urchin development: Study of {\it LpS1\/} gene regulation and cloning and characterization of the {\it SpKrox}-{\it 1\/} gene

Cell differentiation are associated with activation of cell lineage-specific genes. The $LpS{\it 1}\beta$ gene of Lytechinus pictus is activated at the late cleavage stage. $LpS{\it 1}\beta$ transcripts accumulate exclusively in aboral ectoderm lineages. Previous studies demonstrated two G-string DNA-elements, proximal and distal G-strings, which bind to an ectoderm-enriched nuclear factor. In order to define the cis-elements which control positive expression of the $LpS{\it 1}\beta$ gene, the regulatory region from $-$108 to +17 bp of the $LpS{\it 1}\beta$ gene promoter was characterized. The ectoderm G-string factor binds to a G/C-rich region larger than the G-string itself and the binding of the G-string factor requires sequences immediately downstream from the G-string. These downstream sequences are essential for full promoter activity. In addition, only 108 bp of $LpS{\it 1}\beta\ 5\sp\prime$ flanking DNA drives $LpS{\it 1}\beta$ gene expression in aboral ectoderm/mesenchyme cells. Therefore, for positive control of $LpS{\it 1}\beta$ gene expression, two regions of 5$\sp\prime$ flanking DNA are required: region I from base pairs $-$762 to $-$511, and region II, which includes the G/C-rich element, from base pairs $-$108 to $-$61. A mesenchyme cell repressor element is located within region I.^ DNA-binding proteins play key roles in determination of cell differentiation. The zinc finger domain is a DNA-binding domain present in many transcription factors. Based on homologies in zinc fingers, a zinc finger-encoding gene, SpKrox-1, was cloned from S. purpuratus. The putative SpKrox-1 protein has all structural characteristics of a transcription factor: four zinc fingers for DNA binding; acidic domain for transactivation; basic domain for nuclear targeting; and leucine zipper for dimerization. SpKrox-1 RNA transcripts showed a transient expression pattern which correlates largely with early embryonic development. The spatial expression of SpKrox-1 mRNA was distributed throughout the gastrula and larva ectodermal wall. However, SpKrox-1 was not expressed in pigment cells. The SpKrox-1 gene is thus a marker of a subset of SMCs or ectoderm cells. The structural features, and the transient temporal and restricted spatial expression patterns suggest that SpKrox-1 plays a role in a specific developmental event. ^

Volume of natural gas within gas-hydrate and free-gas occurrences in ODP Leg 164 sites

Leg 164 of the Ocean Drilling Program was designed to investigate the occurrence of gas hydrate in the sedimentary section beneath the Blake Ridge on the southeastern continental margin of North America. Sites 994, 995, and 997 were drilled on the Blake Ridge to refine our understanding of the in situ characteristics of natural gas hydrate. Because gas hydrate is unstable at surface pressure and temperature conditions, a major emphasis was placed on the downhole logging program to determine the in situ physical properties of the gas hydrate-bearing sediments. Downhole logging tool strings deployed on Leg 164 included the Schlumberger quad-combination tool (NGT, LSS/SDT, DIT, CNT-G, HLDT), the Formation MicroScanner (FMS), and the Geochemical Combination Tool (GST). Electrical resistivity (DIT) and acoustic transit-time (LSS/SDT) downhole logs from Sites 994, 995, and 997 indicate the presence of gas hydrate in the depth interval between 185 and 450 mbsf on the Blake Ridge. Electrical resistivity log calculations suggest that the gas hydrate-bearing sedimentary section on the Blake Ridge may contain between 2 and 11 percent bulk volume (vol%) gas hydrate. We have determined that the log-inferred gas hydrates and underlying free-gas accumulations on the Blake Ridge may contain as much as 57 trillion m**3 of gas.

Benthic and substrate cover data derived from a time series of photo-transect surveys for the Eastern Banks, Moreton Bay Australia, 2004-2015

This paper describes seagrass species and percentage cover point-based field data sets derived from georeferenced photo transects. Annually or biannually over a ten year period (2004-2015) data sets were collected using 30-50 transects, 500-800 m in length distributed across a 142 km**2 shallow, clear water seagrass habitat, the Eastern Banks, Moreton Bay, Australia. Each of the eight data sets include seagrass property information derived from approximately 3000 georeferenced, downward looking photographs captured at 2-4 m intervals along the transects. Photographs were manually interpreted to estimate seagrass species composition and percentage cover (Coral Point Count excel; CPCe). Understanding seagrass biology, ecology and dynamics for scientific and management purposes requires point-based data on species composition and cover. This data set, and the methods used to derive it are a globally unique example for seagrass ecological applications. It provides the basis for multiple further studies at this site, regional to global comparative studies, and, for the design of similar monitoring programs elsewhere.

Synchronization between prefrontal and posterior association cortex during human working memory

We measured coherence between the electroencephalogram at different scalp sites while human subjects performed delayed response tasks. The tasks required the retention of either verbalizable strings of characters or abstract line drawings. In both types of tasks, a significant enhancement in coherence in the θ range (4–7 Hz) was found between prefrontal and posterior electrodes during 4-s retention intervals. During 6-s perception intervals, far fewer increases in θ coherence were found. Also in other frequency bands, coherence increased; however, the patterns of enhancement made a relevance for working memory processes seem unlikely. Our results suggest that working memory involves synchronization between prefrontal and posterior association cortex by phase-locked, low frequency (4–7 Hz) brain activity.

Natural genetic exchange between Haemophilus and Neisseria: Intergeneric transfer of chromosomal genes between major human pathogens

Members of the bacterial families Haemophilus and Neisseria, important human pathogens that commonly colonize the nasopharynx, are naturally competent for DNA uptake from their environment. In each genus this process is discriminant in favor of its own and against foreign DNA through sequence specificity of DNA receptors. The Haemophilus DNA uptake apparatus binds a 29-bp oligonucleotide domain containing a highly conserved 9-bp core sequence, whereas the neisserial apparatus binds a 10-bp motif. Each motif (“uptake sequence”, US) is highly over-represented in the chromosome of the corresponding genus, particularly concentrated with core sequences in inverted pairs forming gene terminators. Two Haemophilus core USs were unexpectedly found forming the terminator of sodC in Neisseria meningitidis (meningococcus), and sequence analysis strongly suggests that this virulence gene, located next to IS1106, arose through horizontal transfer from Haemophilus. By using USs as search strings in a computer-based analysis of genome sequence, it was established that while USs of the “wrong” genus do not occur commonly in Neisseria or Haemophilus, where they do they are highly likely to flag domains of chromosomal DNA that have been transferred from Haemophilus. Three independent domains of Haemophilus-like DNA were found in the meningococcal chromosome, associated respectively with the virulence gene sodC, the bio gene cluster, and an unidentified orf. This report identifies intergenerically transferred DNA and its source in bacteria, and further identifies transformation with heterologous chromosomal DNA as a way of establishing potentially important chromosomal mosaicism in these pathogenic bacteria.

Mapping species distributions : A comparison of skilled naturalist and lay citizen science recording

Acknowledgements We are grateful to Elaine O’Mahony, Imogen Pearce, Richard Comont, Anthony McCluskey and other BBCT staff for the many hours of BeeWatch species identification and for all people who submitted sightings to BeeWatch, OPAL, BWARS and the various local recording schemes and societies. We thank the NBN for allowing us to download the bumblebee records without strings attached, and the Essex, Greater London, Cumbria and Sussex based recording centres for providing records upon request. Finally, we are indebted to Tom August and two anonymous reviewers for their valuable critique on an earlier version of this work.

Improved DNA binding specificity from polyzinc finger peptides by using strings of two-finger units

Multizinc finger peptides are likely to reach increased prominence in the search for the “ideal” designer transcription factor for in vivo applications such as gene therapy. However, for these treatments to be effective and safe, the peptides must bind with high affinity and, more importantly, with great specificity. Our previous research has shown that zinc finger arrays can be made to bind 18 bp of DNA with picomolar affinity, but also has suggested that arrays of fingers also may bind tightly to related sequences. This work addresses the question of zinc finger DNA binding specificity. We show that by changing the way in which zinc finger arrays are constructed—by linking three two-finger domains rather than two three-finger units—far greater target specificity can be achieved through increased discrimination against mutated or closely related sequences. These new peptides have the added capability of being able to span two short gaps of unbound DNA, although still binding with picomolar affinity to their target sites. We believe that this new method of constructing zinc finger arrays will offer greater efficacy in the fields of gene therapy and in the production of transgenic organisms than previously reported zinc finger arrays.