Fractal-like Distributions over the Rational Numbers in High-throughput Biological and Clinical Data

Recent developments in extracting and processing biological and clinical data are allowing quantitative approaches to studying living systems. High-throughput sequencing, expression profiles, proteomics, and electronic health records are some examples of such technologies. Extracting meaningful information from those technologies requires careful analysis of the large volumes of data they produce. In this note, we present a set of distributions that commonly appear in the analysis of such data. These distributions present some interesting features: they are discontinuous in the rational numbers, but continuous in the irrational numbers, and possess a certain self-similar (fractal-like) structure. The first set of examples which we present here are drawn from a high-throughput sequencing experiment. Here, the self-similar distributions appear as part of the evaluation of the error rate of the sequencing technology and the identification of tumorogenic genomic alterations. The other examples are obtained from risk factor evaluation and analysis of relative disease prevalence and co-mordbidity as these appear in electronic clinical data. The distributions are also relevant to identification of subclonal populations in tumors and the study of the evolution of infectious diseases, and more precisely the study of quasi-species and intrahost diversity of viral populations

ADD force field for sugars and polyols: predicting the additivity of protein-osmolyte interaction

The protein-osmolyte interaction has been shown experimentally to follow an additive construct, where the individual osmolyte-backbone and osmolyte-side-chain interactions contribute to the overall conformational stability of proteins. Here, we computationally reconstruct this additive relation using molecular dynamics simulations, focusing on sugars and polyols, including sucrose and sorbitol, as model osmolytes. A new set of parameters (ADD) is developed for this purpose, using the individual Kirkwood-Buff integrals for sugar-backbone and sugar-side-chain interactions as target experimental data. We show that the ADD parameters can reproduce the additivity of protein-sugar interactions and correctly predict sucrose and sorbitol self-association and their interaction with water. The accurate description of the separate osmolyte-backbone and osmolyte-side-chain contributions also automatically translates into a good prediction of preferential exclusion from the surface of ribonuclease A and α-chymotrypsinogen A. The description of sugar polarity is improved compared to previous force fields, resulting in closer agreement with the experimental data and better compatibility with charged groups, such as the guanidinium moiety. The ADD parameters are developed in combination with the CHARMM36m force field for proteins, but good compatibility is also observed with the AMBER 99SB-ILDN and the OPLS-AA force fields. Overall, exploiting the additivity of protein-osmolyte interactions is a promising approach for the development of new force fields

Relativistic Proton Production During the 14 July 2000 Solar Event: The Case for Multiple Source Mechanisms

Protons accelerated to relativistic energies by transient solar and interplanetary phenomena caused a ground-level cosmic ray enhancement on 14 July 2000, Bastille Day. Near-Earth spacecraft measured the proton flux directly and ground-based observatories measured the secondary responses to higher energy protons. We have modelled the arrival of these relativistic protons at Earth using a technique which deduces the spectrum, arrival direction and anisotropy of the high-energy protons that produce increased responses in neutron monitors. To investigate the acceleration processes involved we have employed theoretical shock and stochastic acceleration spectral forms in our fits to spacecraft and neutron monitor data. During the rising phase of the event (10:45 UT and 10:50 UT) we find that the spectrum between 140 MeV and 4 GeV is best fitted by a shock acceleration spectrum. In contrast, the spectrum at the peak (10:55 UT and 11:00 UT) and in the declining phase (11:40 UT) is best fitted with a stochastic acceleration spectrum. We propose that at least two acceleration processes were responsible for the production of relativistic protons during the Bastille Day solar event: (1) protons were accelerated to relativistic energies by a shock, presumably a coronal mass ejection (CME). (2) protons were also accelerated to relativistic energies by stochastic processes initiated by magnetohydrodynamic (MHD) turbulence.Comment: 38 pages, 9 figures, accepted for publication in the Astrophysical Journal, January, 200

Motif Statistics and Spike Correlations in Neuronal Networks

Motifs are patterns of subgraphs of complex networks. We studied the impact of such patterns of connectivity on the level of correlated, or synchronized, spiking activity among pairs of cells in a recurrent network model of integrate and fire neurons. For a range of network architectures, we find that the pairwise correlation coefficients, averaged across the network, can be closely approximated using only three statistics of network connectivity. These are the overall network connection probability and the frequencies of two second-order motifs: diverging motifs, in which one cell provides input to two others, and chain motifs, in which two cells are connected via a third intermediary cell. Specifically, the prevalence of diverging and chain motifs tends to increase correlation. Our method is based on linear response theory, which enables us to express spiking statistics using linear algebra, and a resumming technique, which extrapolates from second order motifs to predict the overall effect of coupling on network correlation. Our motif-based results seek to isolate the effect of network architecture perturbatively from a known network state

Modeling the Mutualistic Interactions between Tubeworms and Microbial Consortia

The deep-sea vestimentiferan tubeworm Lamellibrachia luymesi forms large aggregations at hydrocarbon seeps in the Gulf of Mexico that may persist for over 250 y. Here, we present the results of a diagenetic model in which tubeworm aggregation persistence is achieved through augmentation of the supply of sulfate to hydrocarbon seep sediments. In the model, L. luymesi releases the sulfate generated by its internal, chemoautotrophic, sulfide-oxidizing symbionts through posterior root-like extensions of its body. The sulfate fuels sulfate reduction, commonly coupled to anaerobic methane oxidation and hydrocarbon degradation by bacterial–archaeal consortia. If sulfate is released by the tubeworms, sulfide generation mainly by hydrocarbon degradation is sufficient to support moderate-sized aggregations of L. luymesi for hundreds of years. The results of this model expand our concept of the potential benefits derived from complex interspecific relationships, in this case involving members of all three domains of life

Effect of Sting Geometry on Axial Force Calculation for the Space Launch System

The primary purpose of this study is to determine the extent to which the size and shape of the wind tunnel sting affect the accuracy of the base pressure corrections applied to measured axial force. The study also includes an assessment of the overall accuracy of the corrections. To accomplish these goals, Computational Fluid Dynamics is used to simulate a simplified version of the geometry of the Space Launch System Block 1B Cargo configuration, paired with a range of wind tunnel sting sizes, over a variety of ascent flight conditions. The base pressure correction method used in the wind tunnel is emulated on the base pressures from the simulated flows and results are compared to direct integration of the base pressures. Differences in results between the two methods provides an assessment of the accuracy of the base force correction method and how that accuracy is affected by sting size

Dispelling the myths of online education: learning via the information superhighway

There continues to be a perception that online education is inferior to traditional education. In the U.S. online learning is more developed than in the U.K. This paper provides insights into a U.S. provision and takes a close look at what are perceived as weaknesses of on line learning and argues that these are not necessarily inherent weaknesses of this form of educational delivery. Then, results of two major studies, undertaken in the U.S. are provided comparing the effectiveness of online education to traditional education as perceived by current MBA students and past graduates. Results of these studies suggest that students of MBA modules and MBA graduates perceive the quality and effectiveness of online education to be similar to, if not higher than, the quality and effectiveness of traditional modules and programmes